torabi

error in dsolve '''too few boundary cond...

Asked: torabi 90 Product: Maple 2017

October 19 2017

1 0

I have a problem for dsolve for equations.How can I remove it?

"too few boundary conditions: expected %1, got %2\"",20,19

thanks

FEW_BCS.mw

restart; Digits := 15; c := 30*10^(-3); b := 50*10^(-3); Et := 210*10^9; Eb := 210*10^9; Ec := 56*10^6; Gc := 22*10^6; dt := .2*10^(-2); db := .2*10^(-2); mc := 60*b*c; mb := 7900*b*db; mt := 7900*b*dt; LL := 300*10^(-3); ks := (1/12)*Pi^2; v := .272727272727272; Gt := Et/(2.6); Gb := Eb/(2.6)

$newsys2 := {1.46666666666667*10^6*(diff(g4(y), y))-1.46666666666667*10^6*(diff(g5(y), y))-6.00000000000000*10^12*g2(y)*omega3-6.00000000000000*10^12*g1(y)*omega3+1466.66666666667*(diff(g7(y), y))-24444.4444444444*(diff(g2(y), y, y))-29333.3333333333*(diff(g9(y), y))+1466.66666666667*(diff(g6(y), y))-1.95555555555556*10^5*(diff(g3(y), y, y))-24444.4444444444*(diff(g1(y), y, y))+2.48888888888889*10^8*g1(y)+2.48888888888889*10^8*g2(y)-4.97777777777778*10^8*g3(y)-4.80000000000000*10^13*g3(y)*omega3, -1.22222222222222*10^5*(diff(g4(y), y))-2.32222222222222*10^6*(diff(g5(y), y))+3.11111111111111*10^7*g1(y)+2.17777777777778*10^8*g2(y)-6.00000000000000*10^12*g3(y)*omega3+3.00000000000000*10^12*g1(y)*omega3-2.48888888888889*10^8*g3(y)-2.21434654297945*10^7*(diff(g7(y), y))+8.14814814814813*10^7*(diff(g8(y), y))-1.20000000000000*10^13*g2(y)*omega3-7.38600329141297*10^7*(diff(g2(y), y, y))+12222.2222222222*(diff(g1(y), y, y))+14666.6666666667*(diff(g9(y), y))-122.222222222222*(diff(g6(y), y))-24444.4444444444*(diff(g3(y), y, y)), -1.20000000000000*10^13*g1(y)*omega3-7.38783662474630*10^7*(diff(g1(y), y, y))+2.32222222222222*10^6*(diff(g4(y), y))+1.22222222222222*10^5*(diff(g5(y), y))-2.21410209853500*10^7*(diff(g6(y), y))+2322.22222222222*(diff(g7(y), y))+2.17777777777778*10^8*g1(y)+3.11111111111111*10^7*g2(y)-2.48888888888889*10^8*g3(y)-8.14814814814813*10^7*(diff(g8(y), y))-6.00000000000000*10^12*g3(y)*omega3-24444.4444444444*(diff(g3(y), y, y))+30555.5555555556*(diff(g2(y), y, y))+3.00000000000000*10^12*g2(y)*omega3+14666.6666666667*(diff(g9(y), y)), -9.77777777777778*10^7*g4(y)-9.77777777777778*10^7*g5(y)+1.95555555555556*10^8*g8(y)-97777.7777777777*g6(y)-97777.7777777777*g7(y)-4.27777777777777*10^6*(diff(g2(y), y))+4.27777777777777*10^6*(diff(g1(y), y))-6.00000000000000*10^12*g4(y)*omega3-6.00000000000000*10^12*g5(y)*omega3+6.00000000000000*10^9*g7(y)*omega3-6.00000000000000*10^9*g6(y)*omega3-4.80000000000000*10^13*g8(y)*omega3-5.37777777777778*10^5*(diff(g8(y), y, y))-67222.2222222222*(diff(g4(y), y, y))-67222.2222222222*(diff(g5(y), y, y))-67.2222222222222*(diff(g6(y), y, y))+67.2222222222222*(diff(g7(y), y, y)), 880.000000000000*g6(y)+880.000000000000*g7(y)+26400.*g9(y)-25666.6666666666*(diff(g1(y), y))+3.85714285714286*10^7*g7(y)*omega3+8.80000*10^5*g4(y)-8.80000*10^5*g5(y)-1.54285714285714*10^9*g9(y)*omega3+3.85714285714286*10^10*g4(y)*omega3-3.85714285714286*10^10*g5(y)*omega3+51333.3333333333*(diff(g3(y), y))+3.85714285714286*10^7*g6(y)*omega3-25666.6666666666*(diff(g2(y), y))-17.2857142857143*(diff(g9(y), y, y))-702.232142857142*(diff(g4(y), y, y))+.432142857142857*(diff(g6(y), y, y))-432.142857142858*(diff(g5(y), y, y))+.432142857142857*(diff(g7(y), y, y)), -8.24095238095238*10^6*g7(y)*omega3+1.28571428571429*10^6*g6(y)*omega3+3.85714285714286*10^7*g9(y)*omega3+1.28571428571429*10^9*g4(y)*omega3+1.27500000000000*10^10*g8(y)*omega3+7.71428571428571*10^9*g5(y)*omega3+17.1111111111112*g6(y)+0.144047619047619e-1*(diff(g6(y), y, y))+6.64302007338279*10^6*g7(y)-77.8642063492063*(diff(g7(y), y, y))-336.111111111107*(diff(g1(y), y))-2566.66666666666*(diff(g3(y), y))+2.21420293186834*10^7*(diff(g2(y), y))-1.14888888888889*10^5*g5(y)+86.4285714285714*(diff(g5(y), y, y))+97777.7777777777*g8(y)+67.2222222222222*(diff(g8(y), y, y))+17111.1111111111*g4(y)+14.4047619047619*(diff(g4(y), y, y))+880.000000000000*g9(y)+.432142857142857*(diff(g9(y), y, y)), 1.28571428571429*10^6*g7(y)*omega3-8.24095238095238*10^6*g6(y)*omega3+3.85714285714286*10^7*g9(y)*omega3-7.71428571428571*10^9*g4(y)*omega3-6.00000000000000*10^9*g8(y)*omega3-1.28571428571429*10^9*g5(y)*omega3+6.64311785116057*10^6*g6(y)-77.8642063492063*(diff(g6(y), y, y))+114.888888888889*g7(y)+0.144047619047619e-1*(diff(g7(y), y, y))+2.21420293186834*10^7*(diff(g1(y), y))-2566.66666666666*(diff(g3(y), y))+213.888888888888*(diff(g2(y), y))-17111.1111111111*g5(y)-14.4047619047619*(diff(g5(y), y, y))-97777.7777777777*g8(y)-67.2222222222222*(diff(g8(y), y, y))+1.14888888888889*10^5*g4(y)-86.4285714285714*(diff(g4(y), y, y))+880.000000000000*g9(y)+.432142857142857*(diff(g9(y), y, y)), 1.28571428571429*10^9*g7(y)*omega3-7.71428571428571*10^9*g6(y)*omega3+3.85714285714286*10^10*g9(y)*omega3-7.97714285714286*10^14*g4(y)*omega3-6.00000000000000*10^12*g8(y)*omega3-1.28571428571429*10^12*g5(y)*omega3+1.14888888888889*10^5*g6(y)-86.4285714285714*(diff(g6(y), y, y))+1.14888888888889*10^5*g7(y)+14.4047619047619*(diff(g7(y), y, y))-1.31388888888889*10^6*(diff(g1(y), y))-2.56666666666666*10^6*(diff(g3(y), y))+2.13888888888889*10^5*(diff(g2(y), y))-1.71111111111111*10^7*g5(y)+1.51250000000000*10^5*(diff(g5(y), y, y))-9.77777777777778*10^7*g8(y)+1.00833333333333*10^5*(diff(g8(y), y, y))+1.14888888888889*10^8*g4(y)-2.33587817460318*10^8*(diff(g4(y), y, y))+8.80000*10^5*g9(y)+432.142857142857*(diff(g9(y), y, y)), 7.71428571428571*10^9*g7(y)*omega3-1.28571428571429*10^9*g6(y)*omega3-3.85714285714286*10^10*g9(y)*omega3-1.28571428571429*10^12*g4(y)*omega3-6.00000000000000*10^12*g8(y)*omega3-7.97714285714286*10^14*g5(y)*omega3-17111.1111111111*g6(y)-14.4047619047619*(diff(g6(y), y, y))-17111.1111111111*g7(y)+86.4285714285714*(diff(g7(y), y, y))+3.36111111111110*10^5*(diff(g1(y), y))+1.74166666666667*10^6*(diff(g3(y), y))+1.31388888888889*10^6*(diff(g2(y), y))+1.14888888888889*10^8*g5(y)-2.33419761904762*10^8*(diff(g5(y), y, y))-9.77777777777778*10^7*g8(y)+2.68888888888889*10^5*(diff(g8(y), y, y))-1.71111111111111*10^7*g4(y)-14404.7619047619*(diff(g4(y), y, y))-8.80000*10^5*g9(y)-432.142857142858*(diff(g9(y), y, y))}$

$bcs3 := {(8/15)*Ec*b*c*((D@@1)(g8))(0)/(-v^2+1)+(1/15)*Ec*b*c*((D@@1)(g4))(0)/(-v^2+1)+(1/15)*Ec*b*c*((D@@1)(g5))(0)/(-v^2+1)+(1/30)*Ec*b*c*dt*((D@@1)(g6))(0)/(-v^2+1)-(1/30)*Ec*b*c*db*((D@@1)(g7))(0)/(-v^2+1)+(2/3)*Ec*b*v*((D@@0)(g2))(0)/(-v^2+1)-(2/3)*Ec*b*v*((D@@0)(g1))(0)/(-v^2+1) = 0, (8/15)*Ec*b*c*((D@@1)(g8))(1)/(-v^2+1)+(1/15)*Ec*b*c*((D@@1)(g4))(1)/(-v^2+1)+(1/15)*Ec*b*c*((D@@1)(g5))(1)/(-v^2+1)+(1/30)*Ec*b*c*dt*((D@@1)(g6))(1)/(-v^2+1)-(1/30)*Ec*b*c*db*((D@@1)(g7))(1)/(-v^2+1)+(2/3)*Ec*b*v*((D@@0)(g2))(1)/(-v^2+1)-(2/3)*Ec*b*v*((D@@0)(g1))(1)/(-v^2+1) = 0, (2/105)*Ec*b*c^3*((D@@1)(g9))(0)/(-v^2+1)+(13/560)*Ec*b*c^2*((D@@1)(g4))(0)/(-v^2+1)-(1/140)*Ec*b*c^2*dt*((D@@1)(g6))(0)/(-v^2+1)+(1/70)*Ec*b*c^2*((D@@1)(g5))(0)/(-v^2+1)-(1/140)*Ec*b*c^2*db*((D@@1)(g7))(0)/(-v^2+1)+(2/15)*Ec*b*c*v*((D@@0)(g1))(0)/(-v^2+1)+(2/15)*Ec*b*c*v*((D@@0)(g2))(0)/(-v^2+1)-(4/15)*Ec*b*c*`#mfrac(mi("v"),mrow(mn("1"),mo("−"),msup(mi("v"),mn("2"))),linethickness = "1")`*((D@@0)(g3))(0) = 0, (2/105)*Ec*b*c^3*((D@@1)(g9))(1)/(-v^2+1)+(13/560)*Ec*b*c^2*((D@@1)(g4))(1)/(-v^2+1)-(1/140)*Ec*b*c^2*dt*((D@@1)(g6))(1)/(-v^2+1)+(1/70)*Ec*b*c^2*((D@@1)(g5))(1)/(-v^2+1)-(1/140)*Ec*b*c^2*db*((D@@1)(g7))(1)/(-v^2+1)+(2/15)*Ec*b*c*v*((D@@0)(g1))(1)/(-v^2+1)+(2/15)*Ec*b*c*v*((D@@0)(g2))(1)/(-v^2+1)-(4/15)*Ec*b*c*v*((D@@0)(g3))(1)/(-v^2+1) = 0, Eb*b*c*((D@@1)(g5))(0)+(1/15)*Ec*b*c*((D@@1)(g8))(0)/(-v^2+1)+(1/70)*Ec*b*c^2*((D@@1)(g9))(0)/(-v^2+1)+(1/70)*Ec*b*c*((D@@1)(g4))(0)/(-v^2+1)+(3/35)*Ec*b*c*((D@@1)(g5))(0)/(-v^2+1)+(1/140)*Ec*b*c*dt*((D@@1)(g6))(0)/(-v^2+1)-(3/70)*Ec*b*c*db*((D@@1)(g7))(0)/(-v^2+1)+(11/30)*Ec*b*v*((D@@0)(g2))(0)/(-v^2+1)+(1/30)*Ec*b*v*((D@@0)(g1))(0)/(-v^2+1)-(2/5)*Ec*b*v*((D@@0)(g3))(0)/(-v^2+1) = 0, Eb*b*c*((D@@1)(g5))(1)+(1/15)*Ec*b*c*((D@@1)(g8))(1)/(-v^2+1)+(1/70)*Ec*b*c^2*((D@@1)(g9))(1)/(-v^2+1)+(1/70)*Ec*b*c*((D@@1)(g4))(1)/(-v^2+1)+(3/35)*Ec*b*c*((D@@1)(g5))(1)/(-v^2+1)+(1/140)*Ec*b*c*dt*((D@@1)(g6))(1)/(-v^2+1)-(3/70)*Ec*b*c*db*((D@@1)(g7))(1)/(-v^2+1)+(11/30)*Ec*b*v*((D@@0)(g2))(1)/(-v^2+1)+(1/30)*Ec*b*v*((D@@0)(g1))(1)/(-v^2+1)-(2/5)*Ec*b*v*((D@@0)(g3))(1)/(-v^2+1) = 0, (1/12)*Eb*b*db^3*((D@@1)(g7))(0)-(1/30)*Ec*b*c*db*((D@@1)(g8))(0)/(-v^2+1)-(1/140)*Ec*b*c^2*db*((D@@1)(g9))(0)/(-v^2+1)-(1/140)*Ec*b*c*db*((D@@1)(g4))(0)/(-v^2+1)-(3/70)*Ec*b*c*db*((D@@1)(g5))(0)/(-v^2+1)+(3/140)*Ec*b*c*db^2*((D@@1)(g7))(0)/(-v^2+1)-(1/280)*Ec*b*c*dt*db*((D@@1)(g6))(0)/(-v^2+1)-(11/60)*Ec*b*dt*v*((D@@0)(g2))(0)/(-v^2+1)+(1/12)*Ec*b*dt*v*((D@@0)(g1))(0)/(-v^2+1)+(1/5)*Ec*b*db*v*((D@@0)(g3))(0)/(-v^2+1) = 0, (1/12)*Eb*b*db^3*((D@@1)(g7))(1)-(1/30)*Ec*b*c*db*((D@@1)(g8))(1)/(-v^2+1)-(1/140)*Ec*b*c^2*db*((D@@1)(g9))(1)/(-v^2+1)-(1/140)*Ec*b*c*db*((D@@1)(g4))(1)/(-v^2+1)-(3/70)*Ec*b*c*db*((D@@1)(g5))(1)/(-v^2+1)+(3/140)*Ec*b*c*db^2*((D@@1)(g7))(1)/(-v^2+1)-(1/280)*Ec*b*c*dt*db*((D@@1)(g6))(1)/(-v^2+1)-(11/60)*Ec*b*dt*v*((D@@0)(g2))(1)/(-v^2+1)+(1/12)*Ec*b*dt*v*((D@@0)(g1))(1)/(-v^2+1)+(1/5)*Ec*b*db*v*((D@@0)(g3))(1)/(-v^2+1) = 0, Et*b*c*((D@@1)(g4))(0)+(1/15)*Ec*b*c*((D@@1)(g8))(0)/(-v^2+1)-(1/70)*Ec*b*c^2*((D@@1)(g9))(0)/(-v^2+1)+(3/35)*Ec*b*c*((D@@1)(g4))(0)/(-v^2+1)-(3/20)*Ec*b*c*((D@@1)(g5))(0)/(-v^2+1)+(3/70)*Ec*b*c*dt*((D@@1)(g6))(0)/(-v^2+1)-(1/140)*Ec*b*c*db*((D@@1)(g7))(0)/(-v^2+1)-(1/30)*Ec*b*v*((D@@0)(g2))(0)/(-v^2+1)-(11/30)*Ec*b*v*((D@@0)(g1))(0)/(-v^2+1)+(2/5)*Ec*b*v*((D@@0)(g3))(0)/(-v^2+1) = 0, Et*b*c*((D@@1)(g4))(1)+(1/15)*Ec*b*c*((D@@1)(g8))(1)/(-v^2+1)-(1/70)*Ec*b*c^2*((D@@1)(g9))(1)/(-v^2+1)+(3/35)*Ec*b*c*((D@@1)(g4))(1)/(-v^2+1)-(3/20)*Ec*b*c*((D@@1)(g5))(1)/(-v^2+1)+(3/70)*Ec*b*c*dt*((D@@1)(g6))(1)/(-v^2+1)-(1/140)*Ec*b*c*db*((D@@1)(g7))(1)/(-v^2+1)-(1/30)*Ec*b*v*((D@@0)(g2))(1)/(-v^2+1)-(11/30)*Ec*b*v*((D@@0)(g1))(1)/(-v^2+1)+(2/5)*Ec*b*v*((D@@0)(g3))(1)/(-v^2+1) = 0, (1/12)*Et*b*dt^3*((D@@1)(g6))(0)+(1/30)*Ec*b*c*dt*((D@@1)(g8))(0)/(-v^2+1)-(1/140)*Ec*b*c^2*dt*((D@@1)(g9))(0)/(-v^2+1)+(3/70)*Ec*b*c*dt*((D@@1)(g4))(0)/(-v^2+1)+(1/140)*Ec*b*c*dt*((D@@1)(g5))(0)/(-v^2+1)+(3/140)*Ec*b*c*dt^2*((D@@1)(g6))(0)/(-v^2+1)-(1/280)*Ec*b*c*dt*db*((D@@1)(g7))(0)/(-v^2+1)-(1/60)*Ec*b*dt*v*((D@@0)(g2))(0)/(-v^2+1)-(11/60)*Ec*b*dt*v*((D@@0)(g1))(0)/(-v^2+1)+(1/5)*Ec*b*dt*v*((D@@0)(g3))(0)/(-v^2+1) = 0, (1/12)*Et*b*dt^3*((D@@1)(g6))(1)+(1/30)*Ec*b*c*dt*((D@@1)(g8))(1)/(-v^2+1)-(1/140)*Ec*b*c^2*dt*((D@@1)(g9))(1)/(-v^2+1)+(3/70)*Ec*b*c*dt*((D@@1)(g4))(1)/(-v^2+1)+(1/140)*Ec*b*c*dt*((D@@1)(g5))(1)/(-v^2+1)+(3/140)*Ec*b*c*dt^2*((D@@1)(g6))(1)/(-v^2+1)-(1/280)*Ec*b*c*dt*db*((D@@1)(g7))(1)/(-v^2+1)-(1/60)*Ec*b*dt*v*((D@@0)(g2))(1)/(-v^2+1)-(11/60)*Ec*b*dt*v*((D@@0)(g1))(1)/(-v^2+1)+(1/5)*Ec*b*dt*v*((D@@0)(g3))(1)/(-v^2+1) = 0, g1(0) = 0, g1(1) = 0, g2(0) = 0, g2(1) = 0, g3(0) = 0, g3(1) = 0}:$

$extra_bcs := `minus`({seq(seq(seq(((D@@i)(cat(g, j)))(a), j = 1 .. 2), i = 1 .. 1), a = 0 .. 1), seq(seq(seq(((D@@i)(cat(g, j)))(a), j = 3 .. 5), i = 1 .. 1), a = 0 .. 1), seq(seq(seq(((D@@i)(cat(g, j)))(a), j = 4 .. 6), i = 0 .. 1), a = 0 .. 1), seq(seq(seq(((D@@i)(cat(g, j)))(a), j = 7 .. 9), i = 1 .. 1), a = 0 .. 1)}, `~`[lhs](bcs3))$

(1)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

Download FEW_BCS.mw

Error, (in fsolve) {Eq1, Eq2} are in the...

Asked: torabi 90 Product: Maple 2017

October 03 2017

1 4

Hello.

I have a problem for determining a and b?

please help me.

thanks

PADE.mw

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/PADE.mw .

Download PADE.mw

error in dsolve.....

Asked: torabi 90 Product: Maple 2017

September 16 2017

1 0

I have a problem for converging dsolve?

"matrix is singular" and '"initial Newton iteration is not converging"'

thanks

dsolve.mw

>

restart; Digits := 15; c := 5*10^(-12); b := 5*10^(-12); Eb := 210*10^9; Ec := 56*10^6; Gc := 22*10^6; d := 2*10^(-12); L[0] := 17.6*10^(-6); L[1] := 17.6*10^(-6); L[2] := 17.6*10^(-6); `ρc` := 60; `ρs` := 7900; L := 20*10^(-12); `ζc` := 1; `ζs` := 1

>

wt := proc (x, t) options operator, arrow; f1(x)*exp(I*omega*t) end proc:

>

sys3 := {f11, f22, f33, f44, f55};

${-(125/393216)*(diff(diff(diff(diff(diff(diff(f3(x), x), x), x), x), x), x))+(15625/135168)*(diff(diff(diff(diff(f3(x), x), x), x), x))-(4375/44)*(diff(diff(diff(f1(x), x), x), x))-(30625/1936)*(diff(diff(f3(x), x), x)), -12000*(diff(diff(f4(x), x), x))-7500*f3(x)-(9540/79)*f4(x)*omega^2-(30/79)*f5(x)*omega^2+(3/79)*(diff(f1(x), x))*omega^2-(3/158)*(diff(f2(x), x))*omega^2, -12000*(diff(diff(f5(x), x), x))+7500*f3(x)-(9540/79)*f5(x)*omega^2-(30/79)*f4(x)*omega^2+(3/158)*(diff(f1(x), x))*omega^2-(3/79)*(diff(f2(x), x))*omega^2, (2/3)*(diff(diff(f2(x), x), x))+(30/79)*f2(x)*omega^2+(9540/79)*f1(x)*omega^2-(23/33792)*(diff(diff(diff(f3(x), x), x), x))+(1/147456)*(diff(diff(diff(diff(diff(f3(x), x), x), x), x), x))+(1079/72)*(diff(diff(diff(diff(f1(x), x), x), x), x))-(3/1580)*(diff(diff(f1(x), x), x))*omega^2+(3/79)*(diff(f4(x), x))*omega^2+(3/158)*(diff(f5(x), x))*omega^2-128*f1(x)+128*f2(x)+(3/3160)*(diff(diff(f2(x), x), x))*omega^2+(7/20)*(diff(f3(x), x))-(2/3)*(diff(diff(f1(x), x), x)), (9540/79)*f2(x)*omega^2+(30/79)*f1(x)*omega^2-(3/1580)*(diff(diff(f2(x), x), x))*omega^2-(3/158)*(diff(f4(x), x))*omega^2-(3/79)*(diff(f5(x), x))*omega^2+128*f1(x)-128*f2(x)+(3/3160)*(diff(diff(f1(x), x), x))*omega^2+15*(diff(diff(diff(diff(f2(x), x), x), x), x))+(7/20)*(diff(f3(x), x))-(2/3)*(diff(diff(f2(x), x), x))+(2/3)*(diff(diff(f1(x), x), x))-(7/4224)*(diff(diff(diff(f3(x), x), x), x))+(1/147456)*(diff(diff(diff(diff(diff(f3(x), x), x), x), x), x))+(1/64)*(diff(diff(diff(diff(f1(x), x), x), x), x))}$

(1)

>

(2)

>	$extra_bcs := `minus`({seq(seq(((D@@i)(f1))(a), i = 0 .. 3), a = 0 .. 1), seq(seq(((D@@i)(f2))(a), i = 0 .. 3), a = 0 .. 1), seq(seq(((D@@i)(f3))(a), i = 0 .. 5), a = 0 .. 1), seq(seq(((D@@i)(f4))(a), i = 0 .. 1), a = 0 .. 1), seq(seq(((D@@i)(f5))(a), i = 0 .. 1), a = 0 .. 1)}, `~`[lhs](bcs3))$

(3)

>

(2/3)*((D@@2)(f2))(0)-(23/33792)*((D@@3)(f3))(0)+(1/147456)*((D@@5)(f3))(0)+(1079/72)*((D@@4)(f1))(0)-(3/158)*((D@@2)(f1))(0)*omega3+(30/79)*(D(f4))(0)*omega3+(15/79)*(D(f5))(0)*omega3+(3/316)*((D@@2)(f2))(0)*omega3+(7/20)*(D(f3))(0)-(2/3)*((D@@2)(f1))(0)

(4)

>	$extra_bcs := {seq(seq(seq(((D@@i)(cat(f, j)))(a), i = 0 .. N[j]-1), j = 1 .. nops(N)), a = 0 .. 1)}$

(5)

>

Error, invalid input: indices received res, which is not valid for its 1st argument, t

Error, invalid input: subs received res[indx[1]](0), which is not valid for its 1st argument

(1/2)*10^(1/2)*P[1]^(1/2)/Pi

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

>

Download dsolve.mw

problem in dsolve...

Asked: torabi 90 Product: Maple 2017

September 13 2017

1 3

hi..i have a problem with dsolve.

please help me thanks

1.mw

>

restart; Digits := 15; c := 5*10^(-12); b := 5*10^(-12); E[s] := 210*10^9; E[c] := 56*10^6; G[c] := 22*10^6; d := 2*10^(-12); L[0] := 17.6*10^(-6); L[1] := 17.6*10^(-6); L[2] := 17.6*10^(-6); `ρ__c` := 60; `ρ__s` := 7900; L := 20*10^(-12); `ζ__c` := 1; `ζ__s` := 1

>

w[t] := proc (x, t) options operator, arrow; f1(x)*exp(I*omega*t) end proc:

>

$dsys3 := {f11, f22, f33, f44, f55, 12*(L/d)^3*((D@@1)(f5))(0) = 0, 12*(L/d)^3*((D@@1)(f5))(1) = 0, 12*(L/d)^3*((D@@1)(f4))(0)+(1/5)*Zeta*(diff(w__t(x, t), x, x, x, x))*((D@@4)(f1))(0) = 0, 12*(L/d)^3*((D@@1)(f4))(1)+(1/5)*Zeta*(diff(w__t(x, t), x, x, x, x))*((D@@4)(f1))(1) = 0, (1/120)*zeta__c*(c/d)^4*(c/L)^5*((D@@3)(f3))(0)-((1/12)*zeta__c*(c/d)^4*(c/L)^3*((1/2)*E__c/G__c))*((D@@1)(f3))(0)+(1/24)*zeta__c*(c/d)^4*c*((D@@2)(f1))(0)/L = 0, (1/120)*zeta__c*(c/d)^4*(c/L)^5*((D@@3)(f3))(1)-((1/12)*zeta__c*(c/d)^4*(c/L)^3*((1/2)*E__c/G__c))*((D@@1)(f3))(1)+(1/24)*zeta__c*(c/d)^4*c*((D@@2)(f1))(1)/L = 0, -(1/120)*zeta__c*(c/d)^4*(c/L)^5*((D@@4)(f3))(0)+(1/12)*zeta__c*(c/d)^4*(c/L)^3*(1+(1/2)*E__c/G__c)*((D@@2)(f3))(0)-(1/24)*zeta__c*(c/d)^4*c*(diff(w__t(x, t), x, x, x))*((D@@3)(f1))(0)/L = 0, -(1/120)*zeta__c*(c/d)^4*(c/L)^5*((D@@4)(f3))(1)+(1/12)*zeta__c*(c/d)^4*(c/L)^3*(1+(1/2)*E__c/G__c)*((D@@2)(f3))(1)-(1/24)*zeta__c*(c/d)^4*c*(diff(w__t(x, t), x, x, x))*((D@@3)(f1))(1)/L = 0, -((1/4)*zeta__c*E__c*(c/d)^4*c/(E__s*L)*(E__c/G__c))*((D@@1)(f3))(0)+(1/24)*zeta__c*E__c*(c/d)^4*(c/L)^3*((D@@3)(f3))(0)/E__s+(Zeta*c/d*(L/d)-(1/3)*zeta__c*E__c*(c/d)^3*L/(E__s*d))*((D@@2)(f1))(0)-L*((D@@3)(f1))(0)/d = 0, -((1/4)*zeta__c*E__c*(c/d)^4*c/(E__s*L)*(E__c/G__c))*((D@@1)(f3))(1)+(1/24)*zeta__c*E__c*(c/d)^4*(c/L)^3*((D@@3)(f3))(1)/E__s+(Zeta*c/d*(L/d)-(1/3)*zeta__c*E__c*(c/d)^3*L/(E__s*d))*((D@@2)(f1))(1)-L*((D@@3)(f1))(1)/d = 0, -L*((D@@2)(f2))(0)/d-((1/4)*zeta__c*E__c*(c/d)^4*c/(E__s*L)*(E__c/G__c))*((D@@1)(f3))(0)+(1/24)*zeta__c*E__c*(c/d)^4*(c/L)^3*((D@@3)(f3))(0)/E__s+(3/8)*zeta__c*E__c*(c/d)^3*L*((D@@2)(f1))(0)/(E__s*d)+(Zeta*c/d*(L/d))*((D@@2)(f2))(0) = 0, -L*((D@@2)(f2))(1)/d-((1/4)*zeta__c*E__c*(c/d)^4*c/(E__s*L)*(E__c/G__c))*((D@@1)(f3))(1)+(1/24)*zeta__c*E__c*(c/d)^4*(c/L)^3*((D@@3)(f3))(1)/E__s+(3/8)*zeta__c*E__c*(c/d)^3*L*((D@@2)(f1))(1)/(E__s*d)+(Zeta*c/d*(L/d))*((D@@2)(f2))(1) = 0, (1/120)*zeta__c*(c/d)^4*(c/L)^5*((D@@5)(f3))(0)-(1/12)*zeta__c*(c/d)^4*(c/L)^3*(1+(1/2)*E__c/G__c)*((D@@3)(f3))(0)+(1/24)*zeta__c*(c/d)^4*c*((D@@4)(f1))(0)/L+(1/4)*zeta__c*(c/d)^3*L*E__c*((D@@2)(f1))(0)/(d*G__c)+(1/4)*zeta__c*(c/d)^4*c*(E__c/G__c)^2*((D@@1)(f3))(0)/L = 0, (1/120)*zeta__c*(c/d)^4*(c/L)^5*((D@@5)(f3))(1)-(1/12)*zeta__c*(c/d)^4*(c/L)^3*(1+(1/2)*E__c/G__c)*((D@@3)(f3))(1)+(1/24)*zeta__c*(c/d)^4*c*((D@@4)(f1))(1)/L+(1/4)*zeta__c*(c/d)^3*L*E__c*((D@@2)(f1))(1)/(d*G__c)+(1/4)*zeta__c*(c/d)^4*c*(E__c/G__c)^2*((D@@1)(f3))(1)/L = 0, f1(0) = 0, f1(1) = 0, f2(0) = 0, f2(1) = 0}:$

>

solve(sys, {diff(f1(x), `$`(x, 4)), diff(f2(x), `$`(x, 4)), diff(f3(x), `$`(x, 6)), diff(f4(x), `$`(x, 2)), diff(f5(x), `$`(x, 2))}):

>

$"input system must be an ODE system, found %1", {f3(x), u__ob(x, t), u__ot(x, t), w__b(x, t), w__t(x, t)}$

Error, invalid input: indices received res, which is not valid for its 1st argument, t

Error, invalid input: subs received res[indx[1]](0), which is not valid for its 1st argument

>

Download 1.mw

how i can reduce the runtime of the prog...

Asked: torabi 90 Product: Maple 2017

August 14 2017

4 2

hi..how i can rewrite section of this code as another form i,e ''for section''

I have a lot of line as this and runnig cise is time consuming.

is there another way to write this section in order to the runtime of the program is reduced??

thanks

for.mw

>

restart;

>	with(LinearAlgebra):

>	with(VectorCalculus):

>	#Digits:=5: k:=6:

>

l:=0:

>

h:=1:

>

m:=4:

>

n:=4:

>

l1:=2*h:

>

l2:=2*h:

>

N:=0.5:

>

nu:=.3:

>	E_m:=70e9:

>	E_c:=380e9:

>	rho_m:=2702:

>	rho_c:=3800:

>	lambda_m:=nuE_m/((1+nu)(1-2*nu)):

>	lambda_c:=nuE_c/((1+nu)(1-2*nu)):

>	mu_m:=E_m/(2*(1+nu)):

>	mu_c:=E_c/(2*(1+nu)):

>	with(orthopoly):

>	for i from 0 to 5 do: L(i):=sqrt((2i+1)/2)P(i,z): end do:

>	Z:=rho_m+(rho_c-rho_m)*((1/2)+(z/h))^N;

(1)

>	U:=lambda_m+(lambda_c-lambda_m)*((1/2)+(z/h))^N;

(2)

>	S:=mu_m+(mu_c-mu_m)*((1/2)+(z/h))^N;

(3)

>	d:=Matrix([[0,0,0,0,0,0,0,0],[sqrt(3),0,0,0,0,0,0,0],[0,sqrt(15),0,0,0,0,0,0],[sqrt(7),0,sqrt(35),0,0,0,0,0],[0,sqrt(27),0,sqrt(63),0,0,0,0],[sqrt(11),0,sqrt(55),0,sqrt(99),0,0,0],[0,sqrt(39),0,sqrt(91),0,sqrt(143),0,0],[sqrt(15),0,sqrt(75),0,sqrt(135),0,sqrt(195),0]]);

$d := Matrix(8, 8, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (1, 7) = 0, (1, 8) = 0, (2, 1) = 3^(1/2), (2, 2) = 0, (2, 3) = 0, (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (2, 7) = 0, (2, 8) = 0, (3, 1) = 0, (3, 2) = 15^(1/2), (3, 3) = 0, (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (3, 7) = 0, (3, 8) = 0, (4, 1) = 7^(1/2), (4, 2) = 0, (4, 3) = 35^(1/2), (4, 4) = 0, (4, 5) = 0, (4, 6) = 0, (4, 7) = 0, (4, 8) = 0, (5, 1) = 0, (5, 2) = 3*3^(1/2), (5, 3) = 0, (5, 4) = 3*7^(1/2), (5, 5) = 0, (5, 6) = 0, (5, 7) = 0, (5, 8) = 0, (6, 1) = 11^(1/2), (6, 2) = 0, (6, 3) = 55^(1/2), (6, 4) = 0, (6, 5) = 3*11^(1/2), (6, 6) = 0, (6, 7) = 0, (6, 8) = 0, (7, 1) = 0, (7, 2) = 39^(1/2), (7, 3) = 0, (7, 4) = 91^(1/2), (7, 5) = 0, (7, 6) = 143^(1/2), (7, 7) = 0, (7, 8) = 0, (8, 1) = 15^(1/2), (8, 2) = 0, (8, 3) = 5*3^(1/2), (8, 4) = 0, (8, 5) = 3*15^(1/2), (8, 6) = 0, (8, 7) = 195^(1/2), (8, 8) = 0})$