Helo every one,

I'm trying to find the positive real root from:

and substitute it in:

my prog is:

eq5:=u->-3*b*k*u^2-3*b*k^2*u^4-k^3*u^6*b+k*u+5*k^2*u^3-b:  

                   (solve(eq5(u),u)): S:=array([...

Dear all,

need to take data(a,u) from eq1 and plug it in eq2

eq1 is:

eq2 is:

Then need to plot data (a,v)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

I did that:

points:=seq([a,evalf(solve((alpha+(l+alpha)*u+alpha*k*u^2)*a=u*(alpha+(l+alpha)*u+alpha*k*u^2...

Hi there i have function as:

f(B,F)=

where: F,B>0 and 

,

such that: 

h=B^2;

I need to plot the surface (contour plot) f(B,F) at 0 and -ve points i.e f(B,F)<0.

Thanks

Hi experts,

I would like to ask what is the defult method's name for the numerical integration?

Thanks,

Hello experts in Maple,

I'm going to solve IVP numerically by this command bellow:

var:={u(t),v(t),w(t)}: 
dsys1 :=diff(u(t),t)=-(Gamma+AM*cos(q*t+phi))*u(t)+(Delta-AM*sin(q*t+phi))*v(t)+Omega*w(t), diff(v(t),t)=-(Delta+AM*sin(q*t+phi))*u(t)-(Gamma-AM*cos(q*t+phi))*v(t), diff(w(t),t)=-0.5-2*Gamma*w(t)-Omega*u(t):
dsol1 :=dsolve({dsys1,ini1},numeric,method=lsode,var, abserr=1e-9, relerr=1e-8,output=listprocedure):
dsolu:=subs(dsol1,u(t)):dsolv:=subs(dsol1,v(t)):dsolw:=subs(dsol1,w(t)):...