Question: Error in dsolve/numeric/bvp

i have run higer order nonlinear ode bvp

but cant solve the error

please help me

NULL

restart

with(plots)

Pr := .71; beta := .5; alpha := .1; S := .1; Du := .1; Nb := .1; Nt := .1; Sc := .67; Sr := .1; omega := 1.0; Lb := 1.0; Pe := 1.0; delta := 1.0; Nc := 1.0; p := .1; q := 5; r := 5; s := 5; a := 1; b := 2; epsilon := .1

Eq1 := (101-100*lambda)*(diff(f(eta), `$`(eta, 3)))-theta(eta)*beta*(diff(f(eta), `$`(eta, 2)))/(1+theta(eta)*beta)+(1+theta(eta)*beta)*f(eta)*(diff(f(eta), `$`(eta, 2)))-(1+theta(eta)*beta)*(diff(f(eta), eta))^2-(M-alpha)*(1+theta(eta)*beta)*(diff(f(eta), eta))

(101-100*lambda)*(diff(diff(diff(f(eta), eta), eta), eta))-.5*theta(eta)*(diff(diff(f(eta), eta), eta))/(1+.5*theta(eta))+(1+.5*theta(eta))*f(eta)*(diff(diff(f(eta), eta), eta))-(1+.5*theta(eta))*(diff(f(eta), eta))^2-(M-.1)*(1+.5*theta(eta))*(diff(f(eta), eta))

(1)

Eq2 := (1+epsilon*theta(eta))*(diff(theta(eta), `$`(eta, 2)))+f(eta)*(diff(theta(eta), eta))+epsilon*(diff(theta(eta), eta))^2+Pr*S*theta(eta)+Pr*Du*(diff(phi(eta), `$`(eta, 2)))+Pr*Nb*(diff(theta(eta), eta))*(diff(phi(eta), eta))+Pr*Nt*(diff(theta(eta), eta))^2

(1+.1*theta(eta))*(diff(diff(theta(eta), eta), eta))+f(eta)*(diff(theta(eta), eta))+.171*(diff(theta(eta), eta))^2+0.71e-1*theta(eta)+0.71e-1*(diff(diff(phi(eta), eta), eta))+0.71e-1*(diff(theta(eta), eta))*(diff(phi(eta), eta))

(2)

Eq3 := (diff(phi(eta), `$`(eta, 2)))/Sc+f(eta)*(diff(phi(eta), eta))-omega*(diff(theta(eta), eta))*(diff(phi(eta), eta))-omega*(diff(theta(eta), `$`(eta, 2)))*phi(eta)-delta*phi(eta)+omega*Nc*(diff(theta(eta), `$`(eta, 2)))+Sr*(diff(theta(eta), `$`(eta, 2)))

1.492537313*(diff(diff(phi(eta), eta), eta))+f(eta)*(diff(phi(eta), eta))-1.0*(diff(theta(eta), eta))*(diff(phi(eta), eta))-1.0*(diff(diff(theta(eta), eta), eta))*phi(eta)-1.0*phi(eta)+1.10*(diff(diff(theta(eta), eta), eta))

(3)

Eq4 := diff(chi(eta), `$`(eta, 2))+Lb*f(eta)*(diff(chi(eta), eta))-Pe*(diff(phi(eta), eta))*(diff(chi(eta), eta))-Pe*(diff(phi(eta), `$`(eta, 2)))*chi(eta)

diff(diff(chi(eta), eta), eta)+1.0*f(eta)*(diff(chi(eta), eta))-1.0*(diff(phi(eta), eta))*(diff(chi(eta), eta))-1.0*(diff(diff(phi(eta), eta), eta))*chi(eta)

(4)

VM := [0., .5, 1.0]

etainf := 9

bcs := f(0) = 0, (D(f))(0) = p*(D@@2)(f)*0, theta(0) = 1+q*(D(theta))(0), phi(0) = 1+r*(D(phi))(0), chi(0) = 1+s*(D(chi))(0), D(f)*etainf = b/a, theta(etainf) = 0, phi(etainf) = 0, chi(etainf) = 0

f(0) = 0, (D(f))(0) = 0., theta(0) = 1+5*(D(theta))(0), phi(0) = 1+5*(D(phi))(0), chi(0) = 1+5*(D(chi))(0), 9*D(f) = 2, theta(9) = 0, phi(9) = 0, chi(9) = 0

(5)

dsys := {Eq1, Eq2, Eq3, Eq4, bcs}

for i to 3 do M := VM[i]; dsol[i] := dsolve(dsys, numeric, continuation = lambda); print(M); print(dsol[i](0)) end do

 

 

 

 
Please Wait...