Advanced Search

Question: Errors in code

Question:Errors in code

iman 35 Maple 2015
dsolve stringtools
June 19 2016
0

hi.please help me for solve this equation

i encounter with error''

Error, (in StringTools:-IsPrefix) second argument must be a string''

equations which be solved attached as pdf file

thanks

Kernel4.mw

root.pdf
 
> 

restart
> 

with(LinearAlgebra):
> 

Typesetting:-Settings(functionassign=false):
> 

NULL

Constants
 
> 

landa := 0.404e11; -1; mu := 0.27e11; -1; alpha := 0.23e-4; -1; rho := 2707; -1; k := 204; -1; c := 903; -1; nu := .3; -1; E := 0.70e11; -1; T0 := 293; -1; omega := 0.1e-1

0.1e-1

 (1.1.1)
> 

beta := alpha*(3*landa+2*mu):
> 

NULL
> 

varpi := 0.1e-1; -1; No := 15

15

 (1.1.2)

 

 
> 

Eq[1] := besselj(0, xi*b)*(eval(diff(bessely(0, xi*r), r), r = a))-(eval(diff(besselj(0, xi*r), r), r = a))*bessely(0, xi*b):
> 

 

> 

wf1 := unapply(Eq[1], xi):

1

  

1.794010904

  

1

  

2

  

1.794010904

  

1

  

3

  

4.802060761

  

2

  

4

  

4.802060761

  

2

  

5

  

4.802060761

  

2

  

6

  

7.908961712

  

3

  

7

  

7.908961712

  

3

  

8

  

7.908961712

  

3

  

9

  

11.03509457

  

4

  

10

  

11.03509457

  

4

  

11

  

11.03509457

  

4

  

12

  

11.03509457

  

4

  

13

  

14.16798650

  

5

  

14

  

14.16798650

  

5

  

15

  

14.16798650

  

5

  

16

  

17.30400975

  

6

 (1.2.1)
> 

Eq[2] := MTM:-besselj(1, eta*b)*(eval(diff(MTM:-bessely(1, eta*r), r), r = a))-(eval(diff(MTM:-besselj(1, eta*r), r), r = a))*MTM:-bessely(1, eta*b):
> 

wf2 := unapply(Eq[2], eta):

1

  

1.958510605

  

1

  

2

  

1.958510605

  

1

  

3

  

4.857021628

  

2

  

4

  

4.857021628

  

2

  

5

  

4.857021628

  

2

  

6

  

7.941288451

  

3

  

7

  

7.941288451

  

3

  

8

  

7.941288451

  

3

  

9

  

11.05802155

  

4

  

10

  

11.05802155

  

4

  

11

  

11.05802155

  

4

  

12

  

11.05802155

  

4

  

13

  

14.18576207

  

5

  

14

  

14.18576207

  

5

  

15

  

14.18576207

  

5

  

16

  

17.31852918

  

6

 (1.2.2)
> 

 

> 

for m to MM do K0[m] := proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc; KK0[m] := proc (r, m) options operator, arrow; diff(K0[m](r, m), r) end proc; K1[n] := proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc; KK1[n] := proc (r, n) options operator, arrow; diff(K1[n](r, n), r) end proc end do

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

  

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

  

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

  

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

  

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

  

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

  

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

  

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

  

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

  

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

  

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

  

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

  

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

  

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

  

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

  

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

  

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

  

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

  

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

  

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

  

proc (r, m) options operator, arrow; BesselJ(0, xi[m]*r)*BesselY(0, xi[m]*b)-BesselJ(0, xi[m]*b)*BesselY(0, xi[m]*r) end proc

  

proc (r, m) options operator, arrow; MTM:-diff(K0[m](r, m), r) end proc

  

proc (r, n) options operator, arrow; BesselJ(1, eta__n*r)*BesselY(1, eta__n*b)-BesselJ(1, eta__n*b)*BesselY(1, eta__n*r) end proc

  

proc (r, n) options operator, arrow; MTM:-diff(K1[n](r, n), r) end proc

 (1.2.3)
> 

> 

U1 := -(int(r*K0[m]*(diff(K1[n], r)+K1[n]/r), r = a .. b))/(int(r*K0[m]^2, r = a .. b)); -1; U2 := -(int(r*K1[n]*(diff(K0[m], r)), r = a .. b))/(int(r*K1[n]^2, r = a .. b)); -1; U3 := (int(r^2*omega^2*K1[n], r = a .. b))/(int(r*K1[n]^2, r = a .. b))

0.1555555555e-3/K1[0]

 (1.2.4)
> 

m := 0; -1; for m to MM do M__m := int(r*K1[m](r, m)^2, r = a .. b); bb__m := 1/M__m end do

int(r*K1[1](r, 1)^2, r = 1 .. 2)

  

1/(int(r*K1[1](r, 1)^2, r = 1 .. 2))

  

int(r*K1[2](r, 2)^2, r = 1 .. 2)

  

1/(int(r*K1[2](r, 2)^2, r = 1 .. 2))

  

int(r*K1[3](r, 3)^2, r = 1 .. 2)

  

1/(int(r*K1[3](r, 3)^2, r = 1 .. 2))

  

int(r*K1[4](r, 4)^2, r = 1 .. 2)

  

1/(int(r*K1[4](r, 4)^2, r = 1 .. 2))

  

int(r*K1[5](r, 5)^2, r = 1 .. 2)

  

1/(int(r*K1[5](r, 5)^2, r = 1 .. 2))

  

int(r*K1[6](r, 6)^2, r = 1 .. 2)

  

1/(int(r*K1[6](r, 6)^2, r = 1 .. 2))

 (1.2.5)
> 

MM; 1; n; 1; m; 1; U1; 1; U2; 1; U3; 1; xi

6

  

0

  

7

  

-(2/3)*K1[0]/K0[7]

  

0

  

0.1555555555e-3/K1[0]

  

xi

 (1.2.6)
> 

for m to MM do for n to MM do dsys := {diff(S[m][n](t), t, t, t)+xi^2*[m]*(diff(S[m][n](t), t, t))+(-U1*U2+`η__η__n__`^2)*(diff(S[m][n](t), t))+xi[m]^2*`η__η__n__`^2*S[m][n](t) = -(2*U2*bb[m]/(Pi*xi[m])*(-BesselJ(0, xi[m]*b)/BesselJ(1, xi[m]*a)))*q+xi^2*[m]*U3} end do end do; sol := dsolve(dsys)

Error, (in StringTools:-IsPrefix) second argument must be a string

  
> 

 

> 

NULL
> 

for m to MM do for n to MM do dsys2 := {diff(Q__mn(t), t, t, t)+xi[m]^2*(diff(Q__mn(t), t, t))+(-U1*U2+eta__n^2)*(diff(Q__mn(t), t))+xi[m]^2*eta__n^2*Q__mn(t) = -2*BesselJ(0, xi[m]*b)*U1*U2*b__m*(1-exp(-xi[m]^2*t))/(BesselJ(1, xi[m]*a)*Pi*xi[m]^3)} end do end do;
> 

sol2 := dsolve(dsys2)

Error, (in dsolve) invalid input: `PDEtools/sdsolve` expects its 1st argument, SYS, to be of type Or(set({`<>`, `=`, algebraic}), list({`<>`, `=`, algebraic}), `casesplit/ans`(list, list)), but received [{Q__mn(t)*pochhammer(1-n, n)+(1497143767/5000000)*(diff(Q__mn(t), [`$`(t, t)]))+eta__n^2*(diff(Q__mn(t), t))+(1497143767/5000000)*eta__n^2*Q__mn(t) = 0}]

  
> 

``

NULL

NULL

 

Download Kernel4.mw
Please Wait...