error in dsolve.. - MaplePrimes

Ask a Question

Create a Post

Question:error in dsolve..

Posted: torabi 90 Product: Maple 2017

September 16 2017

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)