## Error in Display the plots...

Dear Maple user I am facing error while running the codes  to plot the graph for two data sets .

I am attaching the files.

Error_in_Display_1.mw

## Error in procedure...

Dear maple user i have a partial differential equation whose first approximation solution is w(r,z)=(r^2-h^2)/4 , i want to calculate the two approximations. I am attaching the file please see

# I dont want to solve directly.

HP_Help.mw

## System of PDE's error in solve and plot the grap...

dear maple user help me to correct this error .

I want to solve the coupled pdes  plot the graphs but unable to get it. please help me. Here are the maple codes

restart:
L:=0.4:
d1:=0.3: d2:=0.3:

h:=z->piecewise( z<=d1,   1,
z<=d1+L,   1-(delta/(2))*(1 + cos((2*(Pi)/L)*(z - d1 - L/2))),
z<=d1+L+d2,  1 ):

a:=x4*S*Gr/(4*x1*x5)*sin(alpha):
b:=1/Da+(x3*M/(x1*(1+m^2))):
c:=(1/x1)*DP:
x1:=1/((1-phi1)^2.5*(1-phi2)^2.5):
x5:=khnf/kf:
x3:=shnf/sf:
x4:=(1-phi2)*((1-phi1)+phi1*(RBs1)/(RBf))+phi2*RBs2/RBf:
shnf:=sbf*((ss2+2*sbf-2*phi2*(sbf-ss2))/(ss2+2*sbf+phi2*(sbf-ss2))):
sbf:=sf*((ss1+2*sf-2*phi1*(sf-ss1))/(ss1+2*sf+phi1*(sf-ss1))):
ss2:=2.7*10^(-8):
ss1:=59.6*10^(6):
sf:=6.67*10^(-1):
khnf:=kbf*((ks2+2*kbf-2*phi2*(kbf-ks2))/(ks2+2*kbf+phi2*(kbf-ks2))):
kbf:=kf*((ks1+2*kf-2*phi1*(kf-ks1))/(ks1+2*kf+phi1*(kf-ks1))):
ks1:=401:
ks2:=76.5:
kf:=0.4972:
RBs1:=8933*16.7*10^(6):
RBs2:=6320*18*10^(6):
RBf:=1063*1.8*10^(6):
DP:=(1/192)*(18432*F-77*b^2*h(z)^8+308*b*h(z)^8*a+768*h(z)^6*a+12*b*h(z)^8)/(h(z)^4*(b*h(z)^2-6)):
conds:= w(h(z),0)=0, theta(h(z),0)=0, D[1](w)(r,0)=0, D[1](theta)(r,0)=0:
pdes:= diff(w(r, z), r, r)+(1/r)*(diff(w(r, z), r))-(1/Da+(x3/x1)*M)*w(r, z)-(1/x1)*(DP+theta*x4*sin(alpha))=0,
diff(theta(r, z), r,r)+(1/r)*(diff(theta(r, z), r))+(S/x5)=0:
pars:= { F=1.5,Da=0.1,phi1=0.01,phi2=0.02,alpha=Pi/4,S=0.1,Gr=2,m=0.5}:

PrVals:=[2, 5, 7,9]:
colors:=[red, green, blue, black]:
for j from 1 to numelems(PrVals) do
pars1:=`union`( pars, {M=PrVals[j]}):
pdSol:= pdsolve( eval([pdes], pars1),
eval([conds], pars1),
numeric
);
plt[j]:=pdSol:-plot( w(r,z), z=1, r=0..1, numpoints=200, color=colors[j]);
od:
plots:-display( [seq(plt[j], j=1..numelems(PrVals))]);
Error, final value in for loop must be numeric or character
Error, unable to execute seq
PrVals:=[2, 5, 7,9]:
colors:=[red, green, blue, black]:
for j from 1 to numelems(PrVals) do
pars1:=`union`( pars, {M=PrVals[j]}):
pdSol:= pdsolve( eval([pdes], pars1),
eval([conds], pars1),
numeric
);
plt[j]:=pdSol:-plot( theta(r,z), z=1, r=0..1, numpoints=200, color=colors[j]);
od:
plots:-display( [seq(plt[j], j=1..numelems(PrVals))]);
Error, final value in for loop must be numeric or character
Error, unable to execute seq

PrVals:=[2, 5, 7,9]:
colors:=[red, green, blue, black]:
for j from 1 to numelems(PrVals) do
pars1:=`union`( pars, {M=PrVals[j]}):
pdSol:= pdsolve( eval([pdes], pars1),
eval([conds], pars1),
numeric
);
plt[j]:=pdSol:-plot( DP,  z=0..1, numpoints=200, color=colors[j]);
od:
plots:-display( [seq(plt[j], j=1..numelems(PrVals))]);

## unable to get the desire plot...

Dear Maple users Help me to  get the desire graph for this codes.

restart:
with(plots):
with(IntegrationTools):
h:=z->piecewise( z<=d+1,   1,
z<=d+4,   1-(delta/(2))(1 + cos(2(Pi)*(z - 1 - 1/2))),                                                           z<=d+6,   1 ):
w0:=(-c*h(z)^2/4)+(3/64)(b*c-4*a)*h(z)^4+(19/2304)*b(b-4*a)*h(z)^6:
w1:=(c/4)+(1/16)*(4*a-b*c)*h(z)^2:
w2:=(1/256)(4(b*c-4*a)-b*h(z)^2):
w3:=(1/2304)b(b-4*a):
a:=(x4*S*Gr)*sin(alpha)/(4*x1*x5):
b:=(1/Da)+(x3*M/(x1*(1+m^2))):
c:=(1/x1)*Dp:
Dp:=96*x1/((6-b*h(z)^2)h(z)^4)(F+(a*h(z)/24)-((11/6144)b(b-4*a)*h(z)^8)):
x1:=1/((1-phi1)^2.5*(1-phi2)^2.5):
x2:=(1-phi2)((1-phi1)+phi1*Rs1/Rf)+phi2(Rs2/Rf):
x3:=(shnf)/(sf):
x4:=(1-phi2)((1-phi1)+phi1(RBs1)/(RBf))+phi2*((RBs2)/(RBf)):
x5:=khnf/kf:
shnf:=sbf*((ss2+2*sbf-2*phi2*(sbf-ss2))/(ss2+2*sbf+phi2*(sbf-ss2))):
sbf:=sf*((ss1+2*sf-2*phi1*(sf-ss1))/(ss1+2*sf+phi1*(sf-ss1))):
khnf:=kbf*((ks2+2*kbf-2*phi2*(kbf-ks2))/(ks2+2*kbf+phi2*(kbf-ks2))):
kbf:=kf*((ks1+2*kf-2*phi1*(kf-ks1))/(ks1+2*kf+phi1*(kf-ks1))):
RBs1:=(8933*16.7*10^6):
RBf:=(1063*1.8*10^6):
RBs2:=6320*18*10^6:
kf:=0.492:
sbf:=6.67*10^(-1): ss2:=2.7*10^(-8):
sf:=6.67*10^(-1):ss1:=59.6*10^(6):
ks2:=76.5:kf:=0.492: ks1:=401:
phi1:=0.01: phi2:=0.02:alpha:=Pi/4:m:=0.5:Da:=0.1:Gr:=5:delta:=1:S:=0.5:  d:=1:

W1:=w0+w1*r^2+w2*r^4+w3*r^6:

by varing M =2,5,7 and r varies from 0 to 1 i want this type of graphs.  please see the sample graphs

## How to find abundant/deficient numbers...

Hi MaplePrimes,

I have decided to look for an odd perfect number.
I want Maple code to determine if there is a perfect number beween bounds on a factor.
For definitions of deficient and perfect and abundant numbers, see
https://mathworld.wolfram.com/DeficientNumber.html
https://mathworld.wolfram.com/PerfectNumber.html
https://mathworld.wolfram.com/AbundantNumber.html

I need some sort of loop that closes in on the edge of abundant/deficient numbers.

Specifically, an IF() statement is needed about wheather the function

sigma(a)-2*a is positive or negative.

Regards,

Matt

abundant_edge_30.mw

## Sum of proper divisors procedure...

Hi all,

I wrote a little simple minded procedure for finding the sum of divisors for a positive integer.  I thought I would share.

sum_of_proper_divisors.mw

sum_of_proper_divisors.pdf

Has someone written better code for this task?

Regards,
Matt

## idk what this line should do...

H:= Fit(a*B+b*sinh(c*B), Bdata, Hdata, B)

can someone help me ? i have never used maple13 and i have to describe a code but i cant understand what this line does. can someone help me, please?!

## Maple unter Window 11...

My Maple Program run perfectly  with Window 10.  My new Laptop has Window 11 and I get system  crash wenn I applied

plot3d(...) .

So I think It is somethink with the Grapics Redering ?

Can somebody help ?

Betriebsystemname    Microsoft Windows 11 Home
Version    10.0.22621 Build 22621
Weitere Betriebsystembeschreibung     Nicht verfügbar
Betriebsystemhersteller    Microsoft Corporation
Systemname    ELIPAN
Systemhersteller    SAMSUNG ELECTRONICS CO., LTD.
Systemmodell    750XED
Systemtyp    x64-basierter PC
Prozessor    12th Gen Intel(R) Core(TM) i5-1235U, 1300 MHz, 10 Kern(e), 12 logische(r) Prozessor(en)
BIOS-Version/-Datum    American Megatrends International, LLC. P08CFG.033.220913.HQ, 13.09.2022
SMBIOS-Version    3.4
Version des eingebetteten Controllers    255.255
BIOS-Modus    UEFI
BaseBoard-Hersteller    SAMSUNG ELECTRONICS CO., LTD.
BaseBoard-Produkt    NP750XED-KC5DE
BaseBoard-Version    SAMSUNG_SW_REVISION_12345+10.0.22000
Plattformrolle    Mobil
Sicherer Startzustand    Ein
PCR7-Konfiguration    Erweiterung zum Anzeigen erforderlich
Windows-Verzeichnis    C:\WINDOWS
Systemverzeichnis    C:\WINDOWS\system32
Startgerät    \Device\HarddiskVolume1
Gebietsschema    Deutschland
Hardwareabstraktionsebene    Version = "10.0.22621.819"
Benutzername    ELIPAN\Pan
Zeitzone    Mitteleuropäische Zeit
Installierter physischer Speicher (RAM)    16,0 GB
Gesamter physischer Speicher    15,7 GB
Verfügbarer physischer Speicher    8,28 GB
Gesamter virtueller Speicher    16,7 GB
Verfügbarer virtueller Speicher    6,98 GB
Größe der Auslagerungsdatei    1,00 GB
Auslagerungsdatei    C:\pagefile.sys
Kernel-DMA-Schutz    Ein
Virtualisierungsbasierte Sicherheit    Wird ausgeführt...
Virtualisierungsbasierte Sicherheit – erforderliche Sicherheitseigenschaften
Virtualisierungsbasierte Sicherheit – verfügbare Sicherheitseigenschaften    Allgemeine Virtualisierungsunterstützung, Sicherer Start, DMA-Schutz, UEFI-Code Readonly, SMM Security Mitigations 1.0, Modusbasierte Ausführungssteuerung, APIC-Virtualisierung
Virtualisierungsbasierte Sicherheit – konfigurierte Dienste    Durch Hypervisor erzwungene Codeintegrität
Virtualisierungsbasierte Sicherheit – ausgeführte Dienste    Durch Hypervisor erzwungene Codeintegrität
Windows Defender-Anwendungssteuerungsrichtlinie    Erzwungen
Windows Defender-Anwendungssteuerungs-Richtlinie für den Benutzermodus    Aus
Unterstützung der Geräteverschlüsselung    Erweiterung zum Anzeigen erforderlich
Es wurde ein Hypervisor erkannt. Features, die für Hyper-V erforderlich sind, werden nicht angezeigt.

## create a lib . Problem with Windows 11...

I get this error unter WINDOW 11, with WINDOW10 it was o.k

> lib := cat(kernelopts(mapledir), "/lib/mylibn.mla"); savelib('MYEIGW', lib)

Error, cannot lock repository, C:\Program Files\Maple 13/lib/mylibn.mla

 >
 >
 (1)
 >

## How do I solve a differential equation in Maple?...

By printing each processing step for a parametric surface entered in 3-dimensional space,

i- Finding whether the surface is a ruled(regle) surface and printing it in the ruled surface form,

ii- Calculating the dispersion parameter of the surface and the throat line,

iii- Draws the surface and the throat line on the surface together,

## How can one put 2 print() statements on the same l...

Simple question

Got stuck

tup_10_g_trying.mw

Want output to look like -

tuple10(50,12,14,0,0,0,0,0,0,0,0)

1 5 17 19
2 17 29 31
3 29 41 43
4 47 59 61

Regards,
Matt A

## Matt's procedure for pairs of prime numbers...

Hi again everybody,
I hope this can help in some classrooms.

This procedure has some warnings.  Someone else can clean it up.
Also, the ithprime() command can make it a bit slicker.
I did work with prime constelations project, and you can see my web page at mattanderson.fun  .

Refference oeis.org/A40 .
Some prime constelations, and k-tuples, are not of 'general intrest' yet, so they are not in the OEIS encyclopedia YET!

But prime numbers are exciting.  I like doing calculations.

Hopefully someone in the next generation will take up an intrest and do what I did.

pairs_procedure.mw

pairs_procedure.pdf

I wish the best on all of you.

Regards,

Matt C. Anderson

## second question about a simple two column matrix...

Thank you very much, Maple Community !

I have another (hopefully) easy question about a data set

matrix_start_2_trim.pdf

I want to imitate Online Encyclopedia of Integer Sequences .b file

reference oeis.ort/A40

Cheers,

Matt