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