Question: implicitplot with empty output

November 02 2012 J4James 175
Maple 16
0

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,

output=listprocedure,parameters=[M,phi,beta]);

F,F1,F2:=op(subs(respar,[f(eta),diff(f(eta),eta),diff(f(eta),eta,eta)]));

p:=proc(N,M,phi,beta) if not type([_passed],list(numeric)) then return 'procname(_passed)' 
end if;
respar(parameters=[M,phi,beta]);
F1(N)-1
end proc;
implicitplot(p(5,M,phi,1)=0,phi=0..0.2,M=-2..2,grid=[35,35],gridrefine=2,resolution=1000,
crossingrefine=3);
animate(implicitplot,[p(N, M,phi,1) =0, phi=0..0.2, M= 0..2, color = black,grid = [35, 35], 
gridrefine = 2, resolution = 1000, crossingrefine = 3],N=1.2..10,frames=10);
 

 

Please Wait...