Question: How to extract and plot data from a loop?

August 31 2012 J4James 175
0

Hi,

 

eq1:=diff(f(eta),eta,eta,eta)+f(eta)*diff(f(eta),eta,eta)+beta*(1-diff(f(eta),eta)^2)=0;

eq2:=diff(theta(eta),eta,eta)+pr*f(eta)*diff(theta(eta),eta)+pr*Ec*(diff(f(eta),eta,eta))^2+(2-beta)*pr*Q*theta(eta)=0;

bc:=f(0)=gama,D(f)(0)=lamda,D(f)(N)=1,theta(N)=0,D(theta)(0)=-Nc*(1-theta(0)):

gama:=-1.5:Ec:=0.5:
lamda:=1:
Nc:=3:
pr:=5:
beta:=1:
N:=5:for Q from -1.5 to 1.5 by 0.1 do A[Q]:=dsolve({eq1,eq2,bc},numeric, 
output=array([seq(i, i=0..5,0.1 )]))end do:

for Q from -1.5 to 1.5 by 0.1 do Dtheta[Q]:=-simplify((sqrt(1/(2-beta))*A[Q][2,1][1,6]))
 end do;
 
I need help to plot and extract data from above loop.
The plot and data should be of Dtheta[Q] vs Q
Thanks
 
 
 
Please Wait...