Advanced Search

Question: multiple plots in single layout

Question:multiple plots in single layout

Gabriel samaila 35 Maple
differential-equation homework multiple plotting
November 22 2018
0

Hi guys;

can one please show me how to plot graphs in on layout. i want like soundalgekar2004-deka.pdf using this code

variable.mw
 

> 

restart
> 

with*plots; -1; ode1 := diff(f(eta), eta, eta, eta)-(diff(theta(eta), eta))*(diff(f(eta), eta, eta))/(theta(eta)-R)-(theta(eta)-R)*f(eta)*(diff(f(eta), eta, eta))/(2*R) = 0

diff(diff(diff(f(eta), eta), eta), eta)-(diff(theta(eta), eta))*(diff(diff(f(eta), eta), eta))/(theta(eta)-R)-(1/2)*(theta(eta)-R)*f(eta)*(diff(diff(f(eta), eta), eta))/R = 0

 (1)
> 

ode2 := diff(theta(eta), eta, eta)-N*pr*(diff(f(eta), eta))*theta(eta)+(1/2)*pr*f(eta)*(diff(theta(eta), eta)) = 0

diff(diff(theta(eta), eta), eta)-N*pr*(diff(f(eta), eta))*theta(eta)+(1/2)*pr*f(eta)*(diff(theta(eta), eta)) = 0

 (2)
> 

bcs1 := f(0) = 0, (D(f))(0) = 1, (D(f))(38) = 0;

f(0) = 0, (D(f))(0) = 1, (D(f))(38) = 0

 (3)
> 

fixedparameter := [pr = 0.2e-1, N = 1];

[pr = 0.2e-1, N = 1]

 (4)
> 

ode3 := eval(ode1, fixedparameter);

diff(diff(diff(f(eta), eta), eta), eta)-(diff(theta(eta), eta))*(diff(diff(f(eta), eta), eta))/(theta(eta)-R)-(1/2)*(theta(eta)-R)*f(eta)*(diff(diff(f(eta), eta), eta))/R = 0

 (5)
> 

bcs2 := theta(38) = 0, theta(0) = 1;

theta(38) = 0, theta(0) = 1

 (6)
> 

ode4 := eval(ode2, fixedparameter);

diff(diff(theta(eta), eta), eta)-0.2e-1*(diff(f(eta), eta))*theta(eta)+0.1000000000e-1*f(eta)*(diff(theta(eta), eta)) = 0

 (7)
> 

L := [2, 10, 15];

[2, 10, 15]

 (8)
> 

``
> 

for k to 3 do sol_All := dsolve(eval({bcs1, bcs2, ode3, ode4}, R = L[k]), [f(eta), theta(eta)], numeric, output = listprocedure); Y_sol || k := rhs(sol_All[5]); YP_sol || k := -rhs(sol_All[6]); ftheta || k := rhs(sol_All[3]); pftheta || k := -rhs(sol_All[4]) end do:
> 

for k to 3 do L[k], [(Y_sol || k)(0), (YP_sol || k)(0)] end do;

2, [HFloat(1.0000000000000002), HFloat(0.07429586344532203)]

  

10, [HFloat(1.0000000000000007), HFloat(0.06509471738660902)]

  

15, [HFloat(0.9999999999999997), HFloat(0.06453272061469233)]

 (9)
> 

NULL
> 

NULL
> 

  plot( [ seq((Y_sol||j)(eta), j = 1..3)],
         eta = 0 .. 6,
         labels = [eta, theta(eta)],
         axes = boxed
      );
plot( [ seq((YP_sol||j)(eta), j = 1..3)],
         eta = 0 .. 6,
         labels = [eta, thetaprime(eta)],
         axes = boxed
      );

 plot( [ seq((ftheta||j)(eta), j = 1..3)],
         eta = 0 .. 6,
         labels = [eta, f(eta)],
         axes = boxed
      );
  plot( [ seq((pftheta||j)(eta), j = 1..3)],
         eta = 0 .. 6,
         labels = [eta, fprime(eta)],
         axes = boxed
      );

 

 

 

 
> 

NULL
> 

``


 

Download variable.mw

 
Please Wait...