Hi every one,

I am trying to find a way to plot a dependent variable against a parameter. 

The data I want to use for this ploting should be from the numerical results

obtained by using shooting method.

restart:ShootLib := "C:/MP/Shoot9/":
> libname := ShootLib, libname:
> with( Shoot ):
> with(plots):with(LinearAlgebra):
> A:=0:
> FNS:={ f(z), G(z), H(z)}:
> ODE:={diff(f(z),z)  = G(z),

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, 

Hi every one,

How can I get bouble Curl of this equation

In the above equation, the first term has the Laplacian, the second one has the gradient and in the term "k" is the unit vector and also u=(u1,u2,u3). sqrt(R) is a constant.