Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

HI.please help me for solve differenrtial equation with finite difference method not dsolve numeric solver in maple

thanks

FDM2.mw

L := 1/50000000; -1; eta := 1; -1; PDE[111] := 7.65692307692309*10^(-8)*(diff(f1(x), x, x, x, x))-1.56784615384616*10^12*(diff(f1(x), x, x))+220.592307692308*(diff(f2(x), x, x, x))-3.52947692307693*10^21*(diff(f2(x), x))+43.7538461538462*(diff(f3(x), x, x, x))+4.81292307692309*10^20*(diff(f3(x), x))+6.50473846153848*10^30*f1(x)-7.90000000000000*10^(-8)*eta*f1(x)

0.7656923077e-7*(diff(diff(diff(diff(f1(x), x), x), x), x))-0.1567846154e13*(diff(diff(f1(x), x), x))+220.592307692308*(diff(diff(diff(f2(x), x), x), x))-0.3529476923e22*(diff(f2(x), x))+43.7538461538462*(diff(diff(diff(f3(x), x), x), x))+0.4812923077e21*(diff(f3(x), x))+0.6504738462e31*f1(x)

(1)

PDE[222] := 2.14211538461539*10^(-8)*(diff(f2(x), x, x, x, x))-1.64988461538462*10^12*(diff(f2(x), x, x))+7.90486153846156*10^30*f2(x)-220.592307692308*(diff(f1(x), x, x, x))+3.52947692307693*10^21*(diff(f1(x), x))-5.94323076923080*10^11*(diff(f3(x), x, x))+5.13378461538463*10^30*f3(x)-7.90000000000000*10^(-8)*eta*f2(x)

0.2142115385e-7*(diff(diff(diff(diff(f2(x), x), x), x), x))-0.1649884615e13*(diff(diff(f2(x), x), x))+0.7904861538e31*f2(x)-220.592307692308*(diff(diff(diff(f1(x), x), x), x))+0.3529476923e22*(diff(f1(x), x))-0.5943230769e12*(diff(diff(f3(x), x), x))+0.5133784615e31*f3(x)

(2)

PDE[333] := -6.38076923076924*10^(-31)*(diff(f3(x), x, x, x, x, x, x))+9.66537046153848*10^(-8)*(diff(f3(x), x, x, x, x))-3.10154753538461*10^12*(diff(f3(x), x, x))-43.7538461538462*(diff(f1(x), x, x, x))-4.81292307692309*10^20*(diff(f1(x), x))-5.94323076923080*10^11*(diff(f2(x), x, x))+5.13378461538463*10^30*f2(x)+2.29989058707693*10^31*f3(x)-7.90105333333333*10^(-8)*omega^2*f3(x)+6.58333333333333*10^(-31)*eta*(diff(f3(x), x, x))

-0.6380769231e-30*(diff(diff(diff(diff(diff(diff(f3(x), x), x), x), x), x), x))+0.9665370462e-7*(diff(diff(diff(diff(f3(x), x), x), x), x))-0.3101547535e13*(diff(diff(f3(x), x), x))-43.7538461538462*(diff(diff(diff(f1(x), x), x), x))-0.4812923077e21*(diff(f1(x), x))-0.5943230769e12*(diff(diff(f2(x), x), x))+0.5133784615e31*f2(x)+0.2299890587e32*f3(x)-0.7901053333e-7*omega^2*f3(x)

(3)

bcs := {f1(0) = 0, f1(L) = 0, f2(0) = 0, f2(L) = 0, f3(0) = 0, f3(L) = 0, ((D@@1)(f1))(0) = 0, ((D@@1)(f1))(L) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(L) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(L) = 0, ((D@@2)(f3))(0) = 0, ((D@@2)(f3))(L) = 0}

{f1(0) = 0, f1(1/50000000) = 0, f2(0) = 0, f2(1/50000000) = 0, f3(0) = 0, f3(1/50000000) = 0, (D(f1))(0) = 0, (D(f1))(1/50000000) = 0, (D(f2))(0) = 0, (D(f2))(1/50000000) = 0, (D(f3))(0) = 0, (D(f3))(1/50000000) = 0, ((D@@2)(f3))(0) = 0, ((D@@2)(f3))(1/50000000) = 0}

(4)

``


Download FDM2.mw

hi

i have a problem with convert mathematica code in to maple one

i attached maple file .some line have error.please help me for remove error

thanks

fdm.mw

fdm.pdf

hello.

i have two problem in maple file, that is attached..

one of them is RootOf...note that i suppose that [varepsilon := -2.3650203724313] for i can going on following calculation

and second is  Float(undefined) in calculation integral

please help me

thanks

(m=1_n=6)2.mw

from mathematica,

 

n = 5;
CalabiYau[z_, k1_, k2_] := Module[{z1 = Exp[2Pi I k1/n]Cosh[z]^(2/n), z2 = Exp[2Pi I k2/n]Sinh[z]^(2/n)}, {Re[z1], Re[z2], Cos[alpha]Im[z1] + Sin[alpha]Im[z2]}];
Do[alpha = (0.25 + t)Pi; Show[Graphics3D[Table[ParametricPlot3D[CalabiYau[x + I y, k1, k2], {x, -1, 1}, {y, 0, Pi/2}, DisplayFunction -> Identity, Compiled ->False][[1]], {k1, 0, n - 1}, {k2, 0, n - 1}], PlotRange -> 1.5{{-1, 1}, {-1, 1}, {-1, 1}}, ViewPoint -> {1, 1, 0}]], {t, 0, 1, 0.1}];

 

n := 5;

z1 := exp(2*3.14*I*k1/n)*cosh(z)^(2/n);
z2 := exp(2*3.14*I*k2/n)*sinh(z)^(2/n);

alpha = (0.25 + t)Pi;

xx := Re(z1);
yy := Re(z2);
uu := cos(alpha)*Im(z1) + sin(alpha)*Im(z2);

 

where k1, k2, alpha are variables

print([xx,yy,uu]);

i find algcurve has implicitize

how to use this implicitize to find 3d surface?

is there any other method to find?

 

i searched groebner basis can do this, but in mathematica is different from maple example

i got this error in window 8 in surface 2  then follow this post and install again still error

https://www.maplesoft.com/support/faqs/detail.aspx?sid=139020

then follow

http://www.maplesoft.com/support/faqs/detail.aspx?sid=32607

then follow and install again same error

http://www.maplesoft.com/support/faqs/detail.aspx?sid=32631

and install again same eror

then i add option -f c:\Program File (x86)\MapleXX in cmd and then no error any more 

but no install succeed 

where it go, it still not install

then i try again, there is no room enough to install,  hard disk do not have enough space, then i go to c:\Windows\Temp, after deleted file in it, still not enough space

 

i find 

https://www.maplesoft.com/support/install/maple15_install.html

but template do not state how to activate later

how to write this template and how to clear the temp file created by previous failed cmd install method

hello friend

please run this code and resend result to me.

thank you.

 

restart;
with(LinearAlgebra); with(plots);
c[1] := 0; c[2] := 1/2;
a[2, 1] := 1/8;
b[-1] := -1/3; b[1] := 2/3; b[2] := 5/6;
bb[2] := 5/12;
s := 2;
e[1] := Matrix(s+1, 1); e[2] := Matrix(s+1, 1); e[3] := Matrix(s+1, 1); e[4] := Matrix(s+1, 1);
c[0] := 0;
for ii to s+1 do e[1][ii, 1] := 1; e[2][ii, 1] := -1; e[3][ii, 1] := c[ii-1] end do;
e[1][1, 1] := 0; e[2][2, 1] := 0; e[4] := -e[3];
A := Matrix(s+1, s+1);
for i from 3 to s+1 do for j from 2 to i-1 do A[i, j] := a[i-1, j-1] end do; A[i, 1] := A[i, 2] end do;
Id := Matrix(s+1, s+1, shape = identity);
N := Matrix(s+1, s+1);
for i to s+1 do for j to s+1 do N[i, j] := -H*A[i, j]+Id[i, j] end do end do;
Bb := Matrix(s+1, 1);
for i from 3 to s+1 do Bb[i, 1] := bb[i-1] end do;
B := Matrix(s+1, 1);
for i from 2 to s+1 do B[i, 1] := b[i-1] end do; B[1, 1] := b[-1];
Z := Multiply(Transpose(Bb), 1/N);
for i to s+1 do Z[1, i] := H*Z[1, i] end do;
Q := Multiply(Transpose(B), 1/N);
for i to s+1 do Q[1, i] := H*Q[1, i] end do;
Dh := proc (H) options operator, arrow; Matrix([[1+Multiply(Z, e[1]), Multiply(Z, e[2]), 3/2+Multiply(Z, e[3]), -1/2+Multiply(Z, e[4])], [1, 0, 0, 0], [Multiply(Q, e[1]), Multiply(Q, e[2]), 1+Multiply(Q, e[3]), Multiply(Q, e[4])], [0, 0, 1, 0]]) end proc;
k := Matrix(4, 4);
for i to 4 do for j to 4 do if i = j then k[i, j] := ep-Dh(H)[i, j] else k[i, j] := -Dh(H)[i, j] end if end do end do;
sp := unapply(collect(Determinant(k), ep), H, ep);
sol := solve(sp(H, ep), ep, explicit);
sol := simplify(map(allvalues, {sol}));
l := numelems(sol);
inequal({evalc(abs(subs(H = x+I*y, sol[1]))) <= 1, evalc(abs(subs(H = x+I*y, sol[2]))) <= 1, evalc(abs(subs(H = x+I*y, sol[3]))) <= 1, evalc(abs(subs(H = x+I*y, sol[4]))) <= 1}, x = -3 .. 3, y = -3 .. 3, color = "Nautical 1");

 

hi every one
please help me for this problem
f

restart; Digits := 10; Ha := 2; R := 2; `&theta;r` := .5; Rt := 1; B := 1.5; Xi := 0; U[0] := 0; U[1] := alpha; U[2] := -6+(1/16)*Rt*Xi; U[3] := (1/6)*beta; Theta[0] := -(1/2)*Rt Theta[1]:=phi:   delta:=(k)->`if`(k=0, 1,0);   for k  from 0 to 20do    U[k+4]:=((Ha^(2)+R)*U[k+2]+Xi*Theta[k+2])/(4(k+3)*(k+4)) 

Theta[k+2] := (-4*R*`&theta;r`*(`&theta;r`*Theta[k]+1)^2*(sum((i+1)*Theta[i+1]*(k-i+1)*Theta[k-i+1], i = 0 .. k))-B*(sum((i+1)*U[i+1]*(k-i+1)*U[k-i+1], i = 0 .. k))+(1/4)*Ha^2*B*(sum(U[i]*(k-i)*U[k-i], i = 0 .. k)))/((1+(4/3)*R*(`&theta;r`*Theta[k]+1)^3)*(k+1)*(k+2))

 od:  u:=0:  theta:=0:  for k from 0 to 20 do  u:=u+U[k]*y^k:  theta:=theta+Theta[k]*y^k:  od:  print(expand(u)):

with(numapprox); pade(u, y, [3, 3])

pade(diff(u, `$`(y, 2)), y, [3, 3])

pade(theta, y, [3, 3])

solve({limit(pade(theta, y, [3, 3]), y = 1) = (1/2)*Rt, limit(pade(u, y, [3, 3]), y = 1) = 0, limit(pade(diff(u, `$`(y, 2)), y, [3, 3]), y = 1) = -12-(1/8)*Rt*Xi}, [alpha, beta, phi])

This is a specific Maple 2016 question: I find that plots no longer remember their attributes (like size on the screen) upon re-execution of a worksheet, even though I have set the appropriate flag in the Interface pane of the preferences.

Can someone confirm this; or am I just being dense? Is this an OS X issue? An El Capitan issue?? Maple 2015 on Yosemite behaves as expected.

Maple 2016.1 on Mac OS X 10.11.5.

Thanks,

M.D.

hi .please help me for this error

Error, (in fprintf) number expected for floating point format...

float.mw

T(r, t) := -(1.082744105*(.4426072302*BesselJ(1, 1.511248492*r)+.5201541352*BesselY(1, 1.511248492*r)))*Besselj(0, 3.588021808)*(1.-1.*exp(-3.218475123670897216*t))/Besselj(1, 1.794010904)-(0.7489242453e-1*(-.4693819312*BesselJ(1, 2.197141326*r)+.6841159909*BesselY(1, 2.197141326*r)))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)-(0.1290812036e-1*(-.6766364675*BesselJ(1, 3.111788255*r)-.7242359729*BesselY(1, 3.111788255*r)))*Besselj(0, 15.817923424)*(1.-1.*exp(-62.551675361881970944*t))/Besselj(1, 7.908961712)-(0.5633715672e-1*(.8435152108*BesselJ(1, 4.692507174*r)+.8812474286*BesselY(1, 4.692507174*r)))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)-(0.6254981294e-1*(-.8481299777*BesselJ(1, 5.429681040*r)+1.002082035*BesselY(1, 5.429681040*r)))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)

-1.082744105*(.4426072302*BesselJ(1, 1.511248492*r)+.5201541352*BesselY(1, 1.511248492*r))*Besselj(0, 3.588021808)*(1.-1.*exp(-3.218475123670897216*t))/Besselj(1, 1.794010904)-0.7489242453e-1*(-.4693819312*BesselJ(1, 2.197141326*r)+.6841159909*BesselY(1, 2.197141326*r))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)-0.1290812036e-1*(-.6766364675*BesselJ(1, 3.111788255*r)-.7242359729*BesselY(1, 3.111788255*r))*Besselj(0, 15.817923424)*(1.-1.*exp(-62.551675361881970944*t))/Besselj(1, 7.908961712)-0.5633715672e-1*(.8435152108*BesselJ(1, 4.692507174*r)+.8812474286*BesselY(1, 4.692507174*r))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)-0.6254981294e-1*(-.8481299777*BesselJ(1, 5.429681040*r)+1.002082035*BesselY(1, 5.429681040*r))*Besselj(0, 9.604121522)*(1.-1.*exp(-23.059787552335899121*t))/Besselj(1, 4.802060761)

(1)

``

F53 := proc (t) options operator, arrow; eval(T(r, t), r = 1.5) end proc; G53 := proc (t) options operator, arrow; Re(evalf[20](F53(t))) end proc; Temp3 := cat("T-t (N=", No, ",r=1.5).dat"); fd53 := fopen(Temp3, WRITE, TEXT); i := 0; for i from 0 by 0.5e-1 to 500 do tt := i; y := eval(G53(t), t = tt); fprintf(fd53, "%10.5e %10.5e  \n", tt, y) end do; fclose(fd53)

proc (t) options operator, arrow; eval(T(r, t), r = 1.5) end proc

 

Error, (in fprintf) number expected for floating point format

 

``

``


Download float.mw

i work on Maple 18 and recently maple standard worksheet interface not run properly.

i already have  uninstall and again i have same problems

if i install maple 2015 in amazon cloud, activate with my license

and i do not use amazon after several trial, 

will the license be used and can not activate in future if i use amazon cloud again in future?

any command or configure needed when using large amount memory such as 100GB memory?

hi everyone..i would like to ask,why do we need to type in (101-100λ) in our maple program..why do we need the continuation..thanks in advance for answering..

Hi 

Sometimes this happen to me.

why ?  and how can I fix it?

Please recover the content of the attached file if it is possible

Thanks

Check.mw

hi.i have a problem with gain Eigenvectors

please help me

thanks

lorenz_1.mw

restart; with(LinearAlgebra[Modular]); A := Matrix(3, 3, [[-a, a, 0], [0, 0, -sqrt(l.b.c.(h+k))/(h+k)], [(2.*h.sqrt(l.b.c.(h+k)))/(l.(h+k)), (2.*k.sqrt(l.b.c.(h+k)))/(l.(h+k)), -c]])

A := Matrix(3, 3, {(1, 1) = -a, (1, 2) = a, (1, 3) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = -(`.`(l, b, c, h+k))^(1/2)/(h+k), (3, 1) = 2.*(h.((`.`(l, b, c, h+k))^(1/2)))/(l.(h+k)), (3, 2) = 2.*(k.((`.`(l, b, c, h+k))^(1/2)))/(l.(h+k)), (3, 3) = -c})

(1)

Sol := LinearAlgebra[Eigenvectors](A)

Error, (in evala/preproc3) floats not handled yet

 

``


Download lorenz_1.mw

First 139 140 141 142 143 144 145 Last Page 141 of 2097