Advanced Search

Question: Error, (in fproc) unable to store '[0.183848448391796e-1]+0.389981557194716e-2' when datatype=float[8]

Question:Error, (in fproc) unable to store '[0.183848448391796e-1]+0.389981557194716e-2' when datatype=float[8]

DimaBrt 105 Maple
dsolve numeric bvp
November 05 2017
1

Hello. Please help me to correct the error. Thank you

lambda := proc (r) options operator, arrow; 3*r end proc;

proc (r) options operator, arrow; 3*r end proc

 (1)

mu := proc (r) options operator, arrow; 4*r end proc;

proc (r) options operator, arrow; 4*r end proc

 (2)

r1 := 5;

5

 (3)

r2 := 3;

3

 (4)

omega := 7;

7

 (5)

alpha := 1;

1

 (6)

beta := 1;

1

 (7)

rho := 1;

1

 (8)

n := 1;

1

 (9)

rho1 := 4;

4

 (10)

A := 8;

8

 (11)

f11 := proc (r) options operator, arrow; (lambda(r)+2*mu(r))*r^2 end proc;

proc (r) options operator, arrow; (lambda(r)+2*mu(r))*r^2 end proc

 (12)

f12 := proc (r) options operator, arrow; [(diff(lambda(r), r)+2*(diff(mu(r), r)))*r^2+(lambda(r)+2*mu(r))*r] end proc

proc (r) options operator, arrow; [(diff(lambda(r), r)+2*(diff(mu(r), r)))*r^2+(lambda(r)+2*mu(r))*r] end proc

 (13)

f13 := proc (r) options operator, arrow; -n*(lambda(r)+mu(r))*r end proc;

proc (r) options operator, arrow; -n*(lambda(r)+mu(r))*r end proc

 (14)

NULL

f14 := proc (r) options operator, arrow; -alpha*(lambda(r)+mu(r))*r^2 end proc;

proc (r) options operator, arrow; -alpha*(lambda(r)+mu(r))*r^2 end proc

 (15)

f15 := proc (r) options operator, arrow; [(diff(lambda(r), r))*r-lambda(r)-(n^2+2+alpha^2*r^2)*mu(r)+omega^2*rho*r^2] end proc;

proc (r) options operator, arrow; [(diff(lambda(r), r))*r-lambda(r)-(n^2+2+alpha^2*r^2)*mu(r)+omega^2*rho*r^2] end proc

 (16)

f16 := proc (r) options operator, arrow; -n*[(diff(lambda(r), r))*r-lambda(r)-3*mu(r)] end proc;

proc (r) options operator, arrow; -n*[(diff(lambda(r), r))*r-lambda(r)-3*mu(r)] end proc

 (17)

f17 := proc (r) options operator, arrow; -alpha*(diff(lambda(r), r))*r^2 end proc;

proc (r) options operator, arrow; -alpha*(diff(lambda(r), r))*r^2 end proc

 (18)

f21 := proc (r) options operator, arrow; (lambda(r)+2*mu(r))*r^2 end proc;

proc (r) options operator, arrow; (lambda(r)+2*mu(r))*r^2 end proc

 (19)

f22 := proc (r) options operator, arrow; [(diff(lambda(r), r)+2*(diff(mu(r), r)))*r^2+(lambda(r)+2*mu(r))*r] end proc

proc (r) options operator, arrow; [(diff(lambda(r), r)+2*(diff(mu(r), r)))*r^2+(lambda(r)+2*mu(r))*r] end proc

 (20)

f23 := proc (r) options operator, arrow; n*(lambda(r)+mu(r))*r end proc;

proc (r) options operator, arrow; n*(lambda(r)+mu(r))*r end proc

 (21)

f24 := proc (r) options operator, arrow; alpha*(lambda(r)+mu(r))*r^2 end proc;

proc (r) options operator, arrow; alpha*(lambda(r)+mu(r))*r^2 end proc

 (22)

f25 := proc (r) options operator, arrow; [(diff(lambda(r), r))*r-lambda(r)-(n^2+2+alpha^2*r^2)*mu(r)+omega^2*rho*r^2] end proc;

proc (r) options operator, arrow; [(diff(lambda(r), r))*r-lambda(r)-(n^2+2+alpha^2*r^2)*mu(r)+omega^2*rho*r^2] end proc

 (23)

f26 := proc (r) options operator, arrow; n*[(diff(lambda(r), r))*r-lambda(r)-3*mu(r)] end proc;

proc (r) options operator, arrow; n*[(diff(lambda(r), r))*r-lambda(r)-3*mu(r)] end proc

 (24)

f27 := proc (r) options operator, arrow; alpha*(diff(lambda(r), r))*r^2 end proc;

proc (r) options operator, arrow; alpha*(diff(lambda(r), r))*r^2 end proc

 (25)

f31 := proc (r) options operator, arrow; mu(r)*r^2 end proc;

proc (r) options operator, arrow; mu(r)*r^2 end proc

 (26)

f32 := proc (r) options operator, arrow; -n*(lambda(r)+mu(r))*r end proc;

proc (r) options operator, arrow; -n*(lambda(r)+mu(r))*r end proc

 (27)

f33 := proc (r) options operator, arrow; (diff(mu(r), r))*r^2+mu(r)*r end proc;

proc (r) options operator, arrow; (diff(mu(r), r))*r^2+mu(r)*r end proc

 (28)

f34 := proc (r) options operator, arrow; -n*((diff(mu(r), r))*r+lambda(r)+3*mu(r)) end proc;

proc (r) options operator, arrow; -n*((diff(mu(r), r))*r+lambda(r)+3*mu(r)) end proc

 (29)

f35 := proc (r) options operator, arrow; [-(diff(mu(r), r))*r-n^2*lambda(r)-(2*n^2+alpha^2*r^2+1)*mu(r)+omega^2*rho*r^2] end proc;

proc (r) options operator, arrow; [-(diff(mu(r), r))*r-n^2*lambda(r)-(2*n^2+alpha^2*r^2+1)*mu(r)+omega^2*rho*r^2] end proc

 (30)

f36 := proc (r) options operator, arrow; -n*alpha*(lambda(r)+mu(r))*r end proc;

proc (r) options operator, arrow; -n*alpha*(lambda(r)+mu(r))*r end proc

 (31)

f41 := proc (r) options operator, arrow; mu(r)*r^2 end proc;

proc (r) options operator, arrow; mu(r)*r^2 end proc

 (32)

f42 := proc (r) options operator, arrow; n*(lambda(r)+mu(r))*r end proc;

proc (r) options operator, arrow; n*(lambda(r)+mu(r))*r end proc

 (33)

f43 := proc (r) options operator, arrow; (diff(mu(r), r))*r^2+mu(r)*r end proc;

proc (r) options operator, arrow; (diff(mu(r), r))*r^2+mu(r)*r end proc

 (34)

f44 := proc (r) options operator, arrow; n*((diff(mu(r), r))*r+lambda(r)+3*mu(r)) end proc;

proc (r) options operator, arrow; n*((diff(mu(r), r))*r+lambda(r)+3*mu(r)) end proc

 (35)

f45 := proc (r) options operator, arrow; [-(diff(mu(r), r))*r-n^2*lambda(r)-(2*n^2+alpha^2*r^2+1)*mu(r)+omega^2*rho*r^2] end proc;

proc (r) options operator, arrow; [-(diff(mu(r), r))*r-n^2*lambda(r)-(2*n^2+alpha^2*r^2+1)*mu(r)+omega^2*rho*r^2] end proc

 (36)

f46 := proc (r) options operator, arrow; -n*alpha*(lambda(r)+mu(r))*r end proc;

proc (r) options operator, arrow; -n*alpha*(lambda(r)+mu(r))*r end proc

 (37)

f51 := proc (r) options operator, arrow; mu(r)*r^2 end proc;

proc (r) options operator, arrow; mu(r)*r^2 end proc

 (38)

f52 := proc (r) options operator, arrow; -alpha*(lambda(r)+mu(r))*r^2 end proc;

proc (r) options operator, arrow; -alpha*(lambda(r)+mu(r))*r^2 end proc

 (39)

f53 := proc (r) options operator, arrow; (diff(mu(r), r))*r^2+mu(r)*r end proc;

proc (r) options operator, arrow; (diff(mu(r), r))*r^2+mu(r)*r end proc

 (40)

f54 := proc (r) options operator, arrow; -alpha*[(diff(mu(r), r))*r^2+(lambda(r)+mu(r))*r] end proc;

proc (r) options operator, arrow; -alpha*[(diff(mu(r), r))*r^2+(lambda(r)+mu(r))*r] end proc

 (41)

f55 := proc (r) options operator, arrow; -n*alpha*(lambda(r)+mu(r))*r end proc;

proc (r) options operator, arrow; -n*alpha*(lambda(r)+mu(r))*r end proc

 (42)

f56 := proc (r) options operator, arrow; [-(n^2+2*alpha^2*r^2)*mu(r)-alpha^2*lambda(r)*r^2+omega^2*rho*r^2] end proc;

proc (r) options operator, arrow; [-(n^2+2*alpha^2*r^2)*mu(r)-alpha^2*lambda(r)*r^2+omega^2*rho*r^2] end proc

 (43)

f61 := proc (r) options operator, arrow; mu(r)*r^2 end proc;

proc (r) options operator, arrow; mu(r)*r^2 end proc

 (44)

f62 := proc (r) options operator, arrow; alpha*(lambda(r)+mu(r))*r^2 end proc;

proc (r) options operator, arrow; alpha*(lambda(r)+mu(r))*r^2 end proc

 (45)

f63 := proc (r) options operator, arrow; (diff(mu(r), r))*r^2+mu(r)*r end proc;

proc (r) options operator, arrow; (diff(mu(r), r))*r^2+mu(r)*r end proc

 (46)

f64 := proc (r) options operator, arrow; alpha*[(diff(mu(r), r))*r^2+(lambda(r)+mu(r))*r] end proc;

proc (r) options operator, arrow; alpha*[(diff(mu(r), r))*r^2+(lambda(r)+mu(r))*r] end proc

 (47)

f65 := proc (r) options operator, arrow; -n*alpha*(lambda(r)+mu(r))*r end proc;

proc (r) options operator, arrow; -n*alpha*(lambda(r)+mu(r))*r end proc

 (48)

f66 := proc (r) options operator, arrow; [-(n^2+2*alpha^2*r^2)*mu(r)-alpha^2*lambda(r)*r^2+omega^2*rho*r^2] end proc;

proc (r) options operator, arrow; [-(n^2+2*alpha^2*r^2)*mu(r)-alpha^2*lambda(r)*r^2+omega^2*rho*r^2] end proc

 (49)

``

DJ := proc (x, n) options operator, arrow; diff(BesselJ(n, x), x) end proc;

proc (x, n) options operator, arrow; diff(BesselJ(n, x), x) end proc

 (50)

DY := proc (x, n) options operator, arrow; diff(BesselY(n, x), x) end proc;

proc (x, n) options operator, arrow; diff(BesselY(n, x), x) end proc

 (51)

``

g11 := proc (r) options operator, arrow; (lambda(r)+2*mu(r))*subs(x = beta*r1, DJ(x, n)) end proc;

proc (r) options operator, arrow; (lambda(r)+2*mu(r))*subs(x = beta*r1, DJ(x, n)) end proc

 (52)

g12 := proc (r) options operator, arrow; -(lambda(r)+2*mu(r))*subs(x = beta*r1, DY(x, n)) end proc;

proc (r) options operator, arrow; -(lambda(r)+2*mu(r))*subs(x = beta*r1, DY(x, n)) end proc

 (53)

````

g13 := proc (r) options operator, arrow; lambda(r)*subs(x = beta*r1, DJ(x, n))/r+omega^2*rho1*subs(x = beta*r1, BesselJ(n, x))/beta end proc;

proc (r) options operator, arrow; lambda(r)*subs(x = beta*r1, DJ(x, n))/r+omega^2*rho1*subs(x = beta*r1, BesselJ(n, x))/beta end proc

 (54)

g14 := proc (r) options operator, arrow; -lambda(r)*subs(x = beta*r1, DY(x, n))/r-omega^2*rho1*subs(x = beta*r1, BesselY(n, x))/beta end proc;

proc (r) options operator, arrow; -lambda(r)*subs(x = beta*r1, DY(x, n))/r-omega^2*rho1*subs(x = beta*r1, BesselY(n, x))/beta end proc

 (55)

g15 := proc (r) options operator, arrow; -n*lambda(r)*subs(x = beta*r1, DY(x, n))/r end proc;

proc (r) options operator, arrow; -n*lambda(r)*subs(x = beta*r1, DY(x, n))/r end proc

 (56)

BesselY(0, 5)-(1/5)*BesselY(1, 5)

 (57)

g16 := proc (r) options operator, arrow; -n*lambda(r)*subs(x = beta*r1, DJ(x, n))/r end proc;

proc (r) options operator, arrow; -n*lambda(r)*subs(x = beta*r1, DJ(x, n))/r end proc

 (58)

 

g17 := proc (r) options operator, arrow; -alpha*lambda(r)*subs(x = beta*r1, DY(x, n)) end proc;

proc (r) options operator, arrow; -alpha*lambda(r)*subs(x = beta*r1, DY(x, n)) end proc

 (59)

g18 := proc (r) options operator, arrow; -alpha*lambda(r)*subs(x = beta*r1, DJ(x, n)) end proc

proc (r) options operator, arrow; -alpha*lambda(r)*subs(x = beta*r1, DJ(x, n)) end proc

 (60)

q1 := proc (r) options operator, arrow; piecewise(`mod`(n, 2) = 0, 2*i^n*A*omega*rho1/(Pi*beta*r1), `mod`(n, 2) <> 0, 0) end proc

proc (r) options operator, arrow; piecewise(`mod`(n, 2) = 0, 2*i^n*A*omega*rho1/(Pi*beta*r1), `mod`(n, 2) <> 0, 0) end proc

 (61)

g21 := proc (r) options operator, arrow; (lambda(r)+2*mu(r))*subs(x = beta*r1, DY(x, n)) end proc;

proc (r) options operator, arrow; (lambda(r)+2*mu(r))*subs(x = beta*r1, DY(x, n)) end proc

 (62)

g22 := proc (r) options operator, arrow; (lambda(r)+2*mu(r))*subs(x = beta*r1, DJ(x, n)) end proc;

proc (r) options operator, arrow; (lambda(r)+2*mu(r))*subs(x = beta*r1, DJ(x, n)) end proc

 (63)

g23 := proc (r) options operator, arrow; lambda(r)*subs(x = beta*r1, DY(x, n))/r+omega^2*rho1*subs(x = beta*r1, BesselY(n, x))/beta end proc;

proc (r) options operator, arrow; lambda(r)*subs(x = beta*r1, DY(x, n))/r+omega^2*rho1*subs(x = beta*r1, BesselY(n, x))/beta end proc

 (64)

g24 := proc (r) options operator, arrow; -lambda(r)*subs(x = beta*r1, DJ(x, n))/r-omega^2*rho1*subs(x = beta*r1, BesselJ(n, x))/beta end proc;

proc (r) options operator, arrow; -lambda(r)*subs(x = beta*r1, DJ(x, n))/r-omega^2*rho1*subs(x = beta*r1, BesselJ(n, x))/beta end proc

 (65)

g25 := proc (r) options operator, arrow; n*lambda(r)*subs(x = beta*r1, DJ(x, n))/r end proc;

proc (r) options operator, arrow; n*lambda(r)*subs(x = beta*r1, DJ(x, n))/r end proc

 (66)

g26 := proc (r) options operator, arrow; -n*lambda(r)*subs(x = beta*r1, DY(x, n))/r end proc;

proc (r) options operator, arrow; -n*lambda(r)*subs(x = beta*r1, DY(x, n))/r end proc

 (67)

g27 := proc (r) options operator, arrow; alpha*lambda(r)*subs(x = beta*r1, DJ(x, n)) end proc;

proc (r) options operator, arrow; alpha*lambda(r)*subs(x = beta*r1, DJ(x, n)) end proc

 (68)

g28 := proc (r) options operator, arrow; -alpha*lambda(r)*subs(x = beta*r1, DY(x, n)) end proc

proc (r) options operator, arrow; -alpha*lambda(r)*subs(x = beta*r1, DY(x, n)) end proc

 (69)

q2 := proc (r) options operator, arrow; piecewise(`mod`(n, 2) <> 0, 2*i^n*A*omega*rho1/(i*Pi*beta*r1), `mod`(n, 2) = 0, 0) end proc

proc (r) options operator, arrow; piecewise(`mod`(n, 2) <> 0, 2*i^n*A*omega*rho1/(i*Pi*beta*r1), `mod`(n, 2) = 0, 0) end proc

 (70)

g31 := proc (r) options operator, arrow; mu(r) end proc;

proc (r) options operator, arrow; mu(r) end proc

 (71)

g32 := proc (r) options operator, arrow; -n*mu(r)/r end proc;

proc (r) options operator, arrow; -n*mu(r)/r end proc

 (72)

g33 := proc (r) options operator, arrow; -mu(r)/r end proc;

proc (r) options operator, arrow; -mu(r)/r end proc

 (73)

g41 := proc (r) options operator, arrow; mu(r) end proc;

proc (r) options operator, arrow; mu(r) end proc

 (74)

g42 := proc (r) options operator, arrow; n*mu(r)/r end proc;

proc (r) options operator, arrow; n*mu(r)/r end proc

 (75)

g43 := proc (r) options operator, arrow; -mu(r)/r end proc;

proc (r) options operator, arrow; -mu(r)/r end proc

 (76)

g51 := proc (r) options operator, arrow; mu(r) end proc;

proc (r) options operator, arrow; mu(r) end proc

 (77)

g52 := proc (r) options operator, arrow; -alpha*mu(r) end proc;

proc (r) options operator, arrow; -alpha*mu(r) end proc

 (78)

g61 := proc (r) options operator, arrow; mu(r) end proc;

proc (r) options operator, arrow; mu(r) end proc

 (79)

g62 := proc (r) options operator, arrow; alpha*mu(r) end proc;

proc (r) options operator, arrow; alpha*mu(r) end proc

 (80)

``

sys := (diff(x1(r), `$`(r, 2)))*f11(r)+(diff(x1(r), r))*f12(r)+(diff(y2(r), r))*f13(r)+(diff(y3(r), y3))*f14(r)+x1(r)*f15(r)+y2(r)*f16(r)+y3(r)*f17(r) = 0, (diff(y1(r), `$`(r, 2)))*f21(r)+(diff(y1(r), r))*f22(r)+(diff(x2(r), `$`(r, 1)))*f23(r)+(diff(x3(r), r))*f24(r)+y1(r)*f25(r)+x2(r)*f26(r)+x3(r)*f27(r) = 0, (diff(x2(r), `$`(r, 2)))*f31(r)+(diff(y1(r), r))*f32(r)+(diff(x2(r), `$`(r, 1)))*f33(r)+y1(r)*f34(r)+x2(r)*f35(r)+x3(r)*f36(r) = 0, (diff(y2(r), `$`(r, 2)))*f41(r)+(diff(x1(r), r))*f42(r)+(diff(y2(r), r))*f43(r)+x1(r)*f44(r)+y2(r)*f45(r)+y3(r)*f46(r) = 0, (diff(x3(r), `$`(r, 2)))*f51(r)+(diff(y1(r), r))*f52(r)+(diff(x3(r), r))*f53(r)+y1(r)*f54(r)+x2(r)*f55(r)+x3(r)*f56(r) = 0, (diff(y3(r), `$`(r, 2)))*f61(r)+(diff(x1(r), r))*f62(r)+(diff(y3(r), r))*f63(r)+x1(r)*f64(r)+y2(r)*f65(r)+y3(r)*f66(r) = 0;

11*(diff(diff(x1(r), r), r))*r^3+(diff(x1(r), r))*[22*r^2]-7*(diff(y2(r), r))*r^2+x1(r)*[-4*(3+r^2)*r+49*r^2]+y2(r)*[12*r]-3*y3(r)*r^2 = 0, 11*(diff(diff(y1(r), r), r))*r^3+(diff(y1(r), r))*[22*r^2]+7*(diff(x2(r), r))*r^2+7*(diff(x3(r), r))*r^3+y1(r)*[-4*(3+r^2)*r+49*r^2]+x2(r)*[-12*r]+3*x3(r)*r^2 = 0, 4*(diff(diff(x2(r), r), r))*r^3-7*(diff(y1(r), r))*r^2+8*(diff(x2(r), r))*r^2-19*y1(r)*r+x2(r)*[-7*r-4*(3+r^2)*r+49*r^2]-7*x3(r)*r^2 = 0, 4*(diff(diff(y2(r), r), r))*r^3+7*(diff(x1(r), r))*r^2+8*(diff(y2(r), r))*r^2+19*x1(r)*r+y2(r)*[-7*r-4*(3+r^2)*r+49*r^2]-7*y3(r)*r^2 = 0, 4*(diff(diff(x3(r), r), r))*r^3-7*(diff(y1(r), r))*r^3+8*(diff(x3(r), r))*r^2+y1(r)*[-11*r^2]-7*x2(r)*r^2+x3(r)*[-4*(1+2*r^2)*r-3*r^3+49*r^2] = 0, 4*(diff(diff(y3(r), r), r))*r^3+7*(diff(x1(r), r))*r^3+8*(diff(y3(r), r))*r^2+x1(r)*[11*r^2]-7*y2(r)*r^2+y3(r)*[-4*(1+2*r^2)*r-3*r^3+49*r^2] = 0

 (81)

a11 := (D(x1))(r1);

(D(x1))(5)

 (82)

a12 := (D(y1))(r1);

(D(y1))(5)

 (83)

a21 := (D(x2))(r1);

(D(x2))(5)

 (84)

a22 := (D(y2))(r1);

(D(y2))(5)

 (85)

a31 := (D(x3))(r1);

(D(x3))(5)

 (86)

a32 := (D(y3))(r1);

(D(y3))(5)

 (87)

Inits := a11*evalf(g11(r1))+a12*evalf(g12(r1))+x1(r1)*evalf(g13(r1))+y1(r1)*evalf(g14(r1))+x2(r1)*evalf(g15(r1))+y2(r1)*evalf(g16(r1))+x3(r1)*evalf(g17(r1))+y3(r1)*evalf(g18(r1)) = evalf(q1(r1)), a11*evalf(g21(r1))+a12*evalf(g22(r1))+x1(r1)*evalf(g23(r1))+y1(r1)*evalf(g24(r1))+x2(r1)*evalf(g25(r1))+y2(r1)*evalf(g26(r1))+x3(r1)*evalf(g27(r1))+y3(r1)*evalf(g28(r1)) = evalf(q2(r1)), a21*g31(r1)+y1(r1)*g32(r1)+x2(r1)*g33(r1) = 0, a22*g41(r1)+x1(r1)*g42(r1)+y2(r1)*g43(r1) = 0, a31*g51(r1)+y1(r1)*g52(r1) = 0, a32*g61(r1)+x1(r1)*g62(r1) = 0, x1(r2) = 0, x2(r2) = 0, x3(r2) = 0, y1(r2) = 0, y2(r2) = 0, y3(r2) = 0;

-6.164451908*(D(x1))(5)+18.59496397*(D(y1))(5)-64.54175380*x1(5)-27.96690534*y1(5)+1.014270762*x2(5)+.3362428313*y2(5)+5.071353808*x3(5)+1.681214157*y3(5) = 0., -18.59496397*(D(x1))(5)-6.164451908*(D(y1))(5)+27.96690534*x1(5)+64.54175380*y1(5)-.3362428313*x2(5)+1.014270762*y2(5)-1.681214157*x3(5)+5.071353808*y3(5) = 28.52056579, 20*(D(x2))(5)-4*y1(5)-4*x2(5) = 0, 20*(D(y2))(5)+4*x1(5)-4*y2(5) = 0, 20*(D(x3))(5)-20*y1(5) = 0, 20*(D(y3))(5)+20*x1(5) = 0, x1(3) = 0, x2(3) = 0, x3(3) = 0, y1(3) = 0, y2(3) = 0, y3(3) = 0

 (88)

dde := dsolve({Inits, sys}, numeric)

Error, (in fproc) unable to store '[0.183848448391796e-1]+0.389981557194716e-2' when datatype=float[8]

  

NULL


 

Download ODU.mw

 
Please Wait...