iman

Error, (in dsolve/numeric/bvp/convertsys...

Asked: iman 35 Product: Maple 2016

August 07 2016

1 0

hi...

pleaase help me for dsolve equations...

$sys := {-436.394689888941*(diff(f1(x), x, x))-217.044912608998*(diff(f2(x), x))-272.746681180588*(diff(f3(x), x))+103.631312144012*f1(x) = 0, -61.3519912916472*(diff(f2(x), x, x))+103.604272941402*f2(x)+81.3706056883046*(diff(f1(x), x))+325.482422753220*f3(x) = 0, -.368424397601154/((1.-.921060994002885*f3(x))*ln(20.-18.4212198800577*f3(x))^2)-.392026631136497/((1.-.980066577841242*f3(x))*ln(20.-19.6013315568248*f3(x))^2)-0.731318801575311e-1/(1.-.877582561890373*f3(x))^4-.150074517058779/(1.-.900447102352677*f3(x))^4-.166458376732494/(1.-.998750260394966*f3(x))^4-0.767550828335738e-1/(1.-.921060994002885*f3(x))^4-.248643987308266/((1.-.621609968270664*f3(x))*ln(20.-12.4321993654133*f3(x))^2)-.108060461173628/((1.-.540302305868140*f3(x))*ln(20.-10.8060461173628*f3(x))^2)-.398001666111210/((1.-.995004165278026*f3(x))*ln(20.-19.9000833055605*f3(x))^2)-.799000208315973/((1.-.998750260394966*f3(x))*ln(20.-19.9750052078993*f3(x))^2)-0.637368489403740e-1/(1.-.764842187284488*f3(x))^4-0.829170137731688e-1/(1.-.995004165278026*f3(x))^4-.164795179656007/(1.-.988771077936042*f3(x))^4-.751498170277903/((1.-.939372712847379*f3(x))*ln(20.-18.7874542569476*f3(x))^2)-.382134595650242/((1.-.955336489125606*f3(x))*ln(20.-19.1067297825121*f3(x))^2)-.682019617647605/((1.-.852524522059506*f3(x))*ln(20.-17.0504904411901*f3(x))^2)-.351033024756149/((1.-.877582561890373*f3(x))*ln(20.-17.5516512378075*f3(x))^2)-.720357681882142/((1.-.900447102352677*f3(x))*ln(20.-18.0089420470535*f3(x))^2)-.305936874913795/((1.-.764842187284488*f3(x))*ln(20.-15.2968437456898*f3(x))^2)-.132680633091509/(1.-.796083798549056*f3(x))^4-0.687779679091399e-1/(1.-.825335614909678*f3(x))^4-.142087420343251/(1.-.852524522059506*f3(x))^4-0.416666666666667e-1/(1.-1.*f3(x))^4-.330134245963871/((1.-.825335614909678*f3(x))*ln(20.-16.5067122981936*f3(x))^2)-.527986516707986/((1.-.659983145884982*f3(x))*ln(20.-13.1996629176996*f3(x))^2)-.278682683738866/((1.-.696706709347165*f3(x))*ln(20.-13.9341341869433*f3(x))^2)-.585351095099057/((1.-.731688868873821*f3(x))*ln(20.-14.6337773774764*f3(x))^2)-0.816722148201035e-1/(1.-.980066577841242*f3(x))^4-.109997190980830/(1.-.659983145884982*f3(x))^4-0.225125960778392e-1/(1.-.540302305868140*f3(x))^4-0.518008306892220e-1/(1.-.621609968270664*f3(x))^4-0.580588924455971e-1/(1.-.696706709347165*f3(x))^4-.121948144812303/(1.-.731688868873821*f3(x))^4-.921167219057882*(diff(f3(x), x, x))+.727324483148235*(diff(f3(x), x, x, x, x))+545.493362361176*(diff(f1(x), x))+5455.22529158978*f3(x)+1736.35930087199*f2(x)-.465346471571107/((1.-.581683089463883*f3(x))*ln(20.-11.6336617892777*f3(x))^2)-.156562118807896/(1.-.939372712847379*f3(x))^4-0.796113740938005e-1/(1.-.955336489125606*f3(x))^4-.161485403618441/(1.-.968912421710645*f3(x))^4-.636867038839245/((1.-.796083798549056*f3(x))*ln(20.-15.9216759709811*f3(x))^2)-.200000000000000/((1.-1.*f3(x))*ln(20.-20.*f3(x))^2)-.775129937368516/((1.-.968912421710645*f3(x))*ln(20.-19.3782484342129*f3(x))^2)-0.969471815773139e-1/(1.-.581683089463883*f3(x))^4-.791016862348834/((1.-.988771077936042*f3(x))*ln(20.-19.7754215587208*f3(x))^2) = 0}; ode := diff(sys[2], x); sys2 := {ode, sys[1], sys[3]}; bcs := {D@@1*f1(0) = 0, D@@3*f3(1) = 0, f1(0) = 0, f2(0) = 0, f3(0) = 0, (D(f2))(0) = 0, (D(f3))(0) = 0, ((D@@2)(f3))(1) = 0}; bcs11 := eval[recurse](convert(sys[2], D), `union`({x = 1}, bcs)); res := dsolve(`union`(`union`(sys2, {bcs11}), bcs), 'maxmesh' = 2024, numeric, range = 0 .. 1, abserr = 0.1e-3, output = listprocedure)$

Error, (in dsolve/numeric/bvp/convertsys) too few boundary conditions: expected 10, got 9

Download expected.mw

error ''"initial Newton iteration is not...

Asked: iman 35 Product: Maple 2015

July 03 2016

0 0

hi .i have a problem with dsolve this equation.please help me for removing error

"initial Newton iteration is not converging"

thanks..

SGT-FDM.mw

Download SGT-FDM.mw

Errors in code ...

Asked: iman 35 Product: Maple 2015

June 19 2016

0 11

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

>

>	Typesetting:-Settings(functionassign=false):

>

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

(1.1.1)

>

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

>

(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.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.2.2)

>

(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))

(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

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

(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))+(-U1U2+`η__η__n__`^2)(diff(S[m][n](t), t))+xi[m]^2`η__η__n__`^2S[m][n](t) = -(2U2bb[m]/(Pixi[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

>

>	$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))+(-U1U2+eta__n^2)(diff(Q__mn(t), t))+xi[m]^2eta__n^2Q__mn(t) = -2BesselJ(0, xi[m]b)U1U2b__m(1-exp(-xi[m]^2t))/(BesselJ(1, xi[m]a)Pi*xi[m]^3)} end do end do;$

>

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}]

>

Download Kernel4.mw

From RootOf to explicit form...

Asked: iman 35 Product: Maple 2016

June 14 2016

2 4

hi.how convert root of to explicit form.

w is a imaginary..

thanks

123.mw

Download 123.mw

Roots of the characteristic polynomial ...

Asked: iman 35 Product: Maple 2015

June 02 2016

1 0

how i can calculate roots of the characteristic polynomial equations {dsys and dsys2}
and dsolve them with arbitrary initial condition for differennt amont of m and n?
thanks
Kr.mw