Saha

40 Reputation

3 Badges

0 years, 125 days

MaplePrimes Activity


These are replies submitted by Saha

I can't plot because sv value is not available,Please help,What to do

restart; with(plots)

eq1 := diff(f(x), x, x, x)+(1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[fl])*(eta*`cosω`+f(x)*`sinω`)*(diff(f(x), x, x))+(1-phi)^2.5*M*sin^2*alpha*(1-(diff(f(x), x)))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[fl])*lambda[T]*theta;

diff(diff(diff(f(x), x), x), x)+(1/2)*(1-phi)^2.5*(1-phi+phi*rho[s]/rho[fl])*(eta*`cosω`+f(x)*`sinω`)*(diff(diff(f(x), x), x))+(1-phi)^2.5*M*sin^2*alpha*(1-(diff(f(x), x)))+(1-phi)^2.5*(1-phi+phi*`ρβ`[s]/`ρβ`[fl])*lambda[T]*theta

(1)

eq2 := K[nf]*(diff(theta(x), x, x))/K[f]+(1/2)*Pr*(eta*`cosω`+f*`sinω`)*(diff(theta(x), x));

K[nf]*(diff(diff(theta(x), x), x))/K[f]+(1/2)*Pr*(eta*`cosω`+f*`sinω`)*(diff(theta(x), x))

(2)

params := [M = 1, alpha = 0, eta = 0, phi = 0.5e-1, `cosω` = 1, `sinω` = 1, sin(alpha) = 0, lambda[T] = 0, Pr = 6.2, rho[s] = 5200, rho[fl] = 997.1, c[s] = 670, c[f] = 4179, `ρβ`[s] = 6500, `ρβ`[fl] = 20939.1, K[s] = 6, K[f] = .613, K[nf] = .6842];

[M = 1, alpha = 0, eta = 0, phi = 0.5e-1, `cosω` = 1, `sinω` = 1, sin(alpha) = 0, lambda[T] = 0, Pr = 6.2, rho[s] = 5200, rho[fl] = 997.1, c[s] = 670, c[f] = 4179, `ρβ`[s] = 6500, `ρβ`[fl] = 20939.1, K[s] = 6, K[f] = .613, K[nf] = .6842]

 

F(0) = 0, F(1) = val, F(10) = 1, Theta(0) = 1, Theta(10) = 0

(3)

odesys := eval([eq1, eq2], params);

[diff(diff(diff(f(x), x), x), x)+.5325197465*f(x)*(diff(diff(f(x), x), x)), 1.116150082*(diff(diff(theta(x), x), x))+3.100000000*f*(diff(theta(x), x))]

(4)

SM := proc (shoot) local ans; if type(shoot, numeric) then ans := dsolve*(eval*([odesys[], bcs], val = shoot), numeric); return rhs(ans(1)[2]) else return 'procname(args)' end if end proc:

fsolve(SM)

(5)

sol := dsolve*(eval*([odesys[], bcs], val = sv), numeric); sol(1)[1 .. 2]; odeplot*(sol, [[x, f(x)], [x, theta(x)]], x = 0 .. 10, color = [red, blue], legend = [typeset(f(x)), typeset(theta(x))])

Error, (in dsolve) expecting an ODE or a set or list of ODEs. Received 1

 

plots:-odeplot*(dsolve*(eval*([diff(diff(diff(f(x), x), x), x)+.5325197465*f(x)*(diff(diff(f(x), x), x)), 1.116150082*(diff(diff(theta(x), x), x))+3.100000000*f*(diff(theta(x), x)), F(0) = 0, F(1) = val, F(10) = 1, Theta(0) = 1, Theta(10) = 0], val = fsolve(SM)), numeric), [[x, f(x)], [x, theta(x)]], x = 0 .. 10, color = [red, blue], legend = [typeset(f(x)), typeset(theta(x))])

(6)

NULL


 

Download Numeric-Solution.mw

@tomleslie 

@tomleslie Sorry Sir,I typed this in Document mode.So it is like this. 

Page 1 of 1