Advanced Search

Question: Animation of the plots

Question:Animation of the plots

SHIVAS 35 Maple 2018
ode animation homework dsolve numeric odeplot differential-equations plotting
May 31 2024
1

How to get the animation graphs for eta =0..10

NULL

> 

restart; with(plots)
> 

``
> 

ga := .2; Gc := .2; n := 2; Sc := .5; Kp := .2; Q := 0.5e-1; Gr := .2
> 

A := .1Pr := 6.2; Nt := .2; alpha := .1; Rd := 1.5; M := .5; E1 := .3; Ec := .3; Thetap := .2; Nb := .2
> 

NULL
> 

a1 := 1.301348831
NULL
> 

a2 := 1.298194584
a3 := .9728927630; a4 := .9161173998
> 

a5 := 1.316893419
> 

a6 := 1.333333333
 

 
> 

OdeSys := a1*((diff(f(eta), eta, eta, eta))*(2*eta*ga+1)+2*(diff(f(eta), eta, eta))*ga)/a2+A^2-(diff(f(eta), eta))^2+f(eta)*(diff(f(eta), eta, eta))-a1*Kp*(diff(f(eta), eta))/a2-a6*M*(diff(f(eta), eta))/a2+a4*(Theta(eta)*Gr+Phi(eta)*Gc)/a2, f(eta)*(diff(Theta(eta), eta))+a5*(1+4*Rd*(1/3))*((diff(Theta(eta), eta, eta))*(2*eta*ga+1)+2*(diff(Theta(eta), eta))*ga)/(a3*Pr)+(diff(f(eta), eta, eta))^2*(2*eta*ga+1)*Ec*a1/a3+Theta(eta)*Q/a3, f(eta)*(diff(Phi(eta), eta))+((diff(Phi(eta), eta, eta))*(2*eta*ga+1)+2*(diff(Phi(eta), eta))*ga)/Sc+Nt*((diff(Theta(eta), eta, eta))*(2*eta*ga+1)+2*(diff(Theta(eta), eta))*ga)/(Nb*Sc)-Kr*(1+Theta(eta)*(Thetap-1))^n*exp(-E1/(1+Theta(eta)*(Thetap-1))); Cond := f(0) = 0, (D(f))(0) = 1, a5*(D(Theta))(0) = -alpha*(Theta(0)-1), Phi(0) = 1, (D(f))(10) = A, Theta(10) = 0, Phi(10) = 0
> 

KrVals := [0.1e-1, .1, .2, .3]
> 

for j to numelems(KrVals) do Ans[j] := dsolve(eval([OdeSys, Cond], Kr = KrVals[j]), numeric, output = listprocedure) end do
> 

``
> 

``
> 

with(plots):
  cols := [red, blue,green, black]:

 plotA:= display
  ( [ seq
      ( odeplot
        ( Ans[k],[eta,D(f)(eta)],
          eta=0..10,
          color=cols[k]
        ),
        k=1..numelems(KrVals)
      )
    ],
    'axes'= 'boxed',labels=[eta,'f(eta)']
  );

 
> 

with(plots):
  cols := [red, blue, green,black]:

plotC:= display( [ seq( odeplot
        ( Ans[k],[eta,Theta(eta)],
          eta=0..10,
          color=cols[k]
        ),
        k=1..numelems(KrVals)
      )
    ],
    'axes'= 'boxed',labels=[eta,'Theta(eta)']
  );

 
> 

 

``

> 

plotA1 := display(seq(plot3d(r*(eval(f(:-eta), Ans[k]))(eta), eta = 0 .. 10, r = -5 .. 5, color = cols[k]), k = 1 .. nops(KrVals)), linestyle = "solid", style = contour, thickness = 1)

 
 

 

Download ode_plots_animation_graphs.mw

i need like this demo plot  
Please Wait...