Hello,

I need to plot M vs phi from the following ode. For this I am using implicitplot but unfortunately

I am getting empty plot. Please guide me through this issue.

Thanks 

eq1 := 1/(1-phi)^2.5/(1+4117/783*phi)*diff(f(eta),`$`(eta,3))+1/2*f(eta)*diff(f(eta),`$`(eta,2))

+M*(1-diff(f(eta),eta)) = 0;

bc:=f(0)=0,D(f)(0)=beta/(1-phi)^(2.5)*(D@D)(f)(0),D(f)(N)=1;

respar:=dsolve({eq1,f(0)=0,D(f)(0)=beta/(1-phi)^(2.5)*(D@D)(f)(0)},numeric,

How to include/code/type such type of boundary conditions in maple?

Thanks

Hello,

I need help on plotting in for loop.

for M from 1.1 by .2 to 1.7 do 
for F from -3 to 3 by 0.1 do F2[F]:=dsolve(d1, numeric,output=listprocedure):  end do:
for F from -3 to 3 by 0.1 do 
P2[F]:=subs(F2[F],P1);
P3[F]:=evalf(Int(P2[F],0..1)); 
end do:
IQ:=Vector([seq(F,F=-3..3,0.1)]):
DT:=Vector([seq(P3[F],F=-3..3,0.1)]):
plot(IQ,DT);
end do;

The above code gives me 4 separate plots, but I want them all in just one plot.

I have a differential equation having a parameter say 'a'.

If I increase 'a' by 10%, got change in the solution as well.

How to calculate change in percentage in the solution?.

Thanks

restart:with(plots):

Nb:=0.5:Nt:=0.5:G:= 2:B:= 2:beta:= .2:Q:=2: alpha:= 2:

de := diff(alpha*f(y), y, y, y, y)+G*(diff(theta(y), y, y))+B*(diff(phi(y), y, y))

+6*beta*(diff(f(y), y, y))*(diff(f(y), y, y, y))^2+3*beta*(diff(f(y), y, y, y, y))*

(diff(f(y), y, y))^2 = 0, diff(theta(y), y, y)+Nb*(diff(theta(y), y))*(diff(phi(y), y))

+Nt*(diff(theta(y), y))^2 = 0, diff(phi(y), y, y)+Nb*(diff(theta(y), y, y))/Nt = 0: