Advanced Search

SHIVAS

35 Reputation

4 Badges

1 years, 92 days
Contact SHIVAS

MaplePrimes Activity


These are questions asked by SHIVAS

Animation of the plots...

SHIVAS 35 Maple 2018
ode animation homework
May 31 2024
2 3

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  

Ode exact solutions...

SHIVAS 35 Maple 2022
homework dsolve numeric
May 27 2024
1 5
 
> 

restart
> 

with(PDEtools); with(plots)
> 

OdeSys := D1*(diff(U0(y), y, y))+A*S*(diff(U0(y), y))-(M1+D1/K)*U0(y)+B*G1*Theta0(y) = 0, D1*(diff(U1(y), y, y))+A*S*(diff(U1(y), y))-(M1+D1/K+A*i*omega)*U1(y)+B*G1*Theta0(y) = 0, E1*(diff(Theta0(y), y, y))+Pr*S*C*(diff(Theta0(y), y))-F*Theta0(y)+Du*(diff(Phi0(y), y, y))+Pr*C*QL*Phi0(y) = 0, E1*(diff(Theta1(y), y, y))+Pr*S*C*(diff(Theta1(y), y))-F*Pr*C*i*omega*Theta1(y)+Du*(diff(Phi1(y), y, y))+Pr*C*QL*Phi1(y) = 0, diff(Phi0(y), y, y)+S*Sc*(diff(Phi0(y), y))-Kr*Sc*Phi0(y) = 0, diff(Phi1(y), y, y)+S*Sc*(diff(Phi1(y), y))-(i*omega+Kr)*Sc*Phi1(y) = 0
> 

Cond := U0(0) = 1, U1(0) = 0, Theta0(0) = 1, Theta1(0) = 1, Phi0(0) = 1, Phi1(0) = 1, U0(inf) = 0, U1(inf) = 0, Theta0(inf) = 0, Theta1(inf) = 0, Phi0(inf) = 0, Phi1(inf) = 0
> 

solutions := dsolve({Cond, OdeSys})
> 

U0_exact := eval(U0(y), solutions)

Error, invalid input: eval received NULL, which is not valid for its 2nd argument, eqns

  
> 

U1_exact := eval(U1(y), solutions)

Error, invalid input: eval received NULL, which is not valid for its 2nd argument, eqns

  
> 

Theta0_exact := eval(Yheta0(y), solutions)

Error, invalid input: eval received NULL, which is not valid for its 2nd argument, eqns

  
> 

Theta1_exact := eval(Theta1(y), solutions)

Error, invalid input: eval received NULL, which is not valid for its 2nd argument, eqns

  
> 

Phi0_exact := eval(Phi0(y), solutions)

Error, invalid input: eval received NULL, which is not valid for its 2nd argument, eqns

  
> 

Phi1_exact := eval(Phi1(y), solutions)

Error, invalid input: eval received NULL, which is not valid for its 2nd argument, eqns

  
 

 

Download ode_exact.mw

How to find the exact solution of ODE's

Ode equations solution ...

SHIVAS 35 Maple 2018
ode homework numeric
December 15 2023
0 4

How to solve this type of ode in maple

need the value of S, Q, E

how to solve Pde using LaPlace Transfor...

SHIVAS 35 Maple 2018
syntax homework pde
October 13 2023
1 7

Please help me how to impliment to solve this problem 

Download dust_practice_code.mw

dust_paper.pdf

time dependent 2d plot...

SHIVAS 35 Maple 2018
homework pde numeric
September 12 2023
1 4
> 

restart:
> 

  ra:=2: b1:=1.41: na:=0.7: we:=0.5: eta[1]:=4*0.1: d:=0.5:
  xi:=0.1: m:=na: ea:=0.5: pr:=21: gr:=0.1: R:=0.9323556933:
> 

  PDE1:=ra*(diff(f(x,t),t))=+b1*(1+ea*cos(t))+(1/(R^2))*((diff(f(x,t),x,x))+(1/x)*diff(f(x,t),x));
  IBC:= {D[1](f)(0,t)=0,f(1,t)=0,f(x,0)=0};

2*(diff(f(x, t), t)) = 1.41+.705*cos(t)+1.150367877*(diff(diff(f(x, t), x), x))+1.150367877*(diff(f(x, t), x))/x

  

{f(1, t) = 0, f(x, 0) = 0, (D[1](f))(0, t) = 0}

 (1)
> 

sol := pdsolve({PDE1}, IBC, numeric); sol:-plot(f(x, t), t = 1.2, linestyle = "solid", title = "Velocity Profile", labels = ["r", "f"])

 

``

Download pde.mw

for different time plot of f(x,t) in single plot with different color 

1 2 Page 1 of 2