Items tagged with ode

Feed

Hi

 

I have an ODE which is based on a seperate function, and I would like to make a plot with the information

dsolve([diff(X(W), W) = (0.536000000000000e-3*(1-X(W)))*(1+X(W)), X(0) = 0], numeric)

and

C_A:= C_A0*(1-X(W))*(1+X(W))

which has been used as part of the ODE.

I would really like to plot C_A as a function of W. I have no problem plotting X as a function to W using odeplot. Ideally I would like to plot C_A and X vs W in the same plot.

Regards

the question is 

ODE5:= y(x)*diff(y(x),x,x) + (diff(y(x),x)^2=0

dsolve({ODE5,y(0)=4,D(y)(0)=7},y(x))

and my answer appears to be an integration! which is wrong

the correct answer : 2*(4+14*x)^(1/2)

Could someone tell me what did I do wrong? And how could I get to this result?

Thanks a lot!

 

Please i need help to plot the graph of f'' against episoln using the below BVP

 

HELP.mw

Dears

Hope everything fine with you. I want to solve the attached problem by numarically and want to plot it but failed. Please see the attachement and correct it. I am waiting your positive respone. 

System_of_ODEs.mw

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

i don't know much about maple, i need to solve the following odes system... I study a little on the help page of maple about numeric[midrich] that takes bvp and deal singularity as well but dint know how to used in the following system

odes.mw

Respected member!
Please help me to find the solution of attached problem, I am a new user so pleaes forgive any mistakes.
 

``


``


``

NULL

NULL

restart

R := 2.0

2.0

(1)

ODEforNum := r^3*((D@@4)(F))(r)+r^2*R*((D@@3)(F))(r)*(F(r)-2/R)+R*((D(F))(r)-r*((D@@2)(F))(r))*(r*(D(F))(r)+3*(F(r)-1/R)) = 0:

numsol := dsolve({BCSforNum, ODEforNum}, numeric, output = listprocedure)

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

 

``


 

Download mplprimes.mw

G := 6.6743*10^(-8);

a := 1.9501*10^24;

b := .3306;

c := 2.99792458*10^10;

d := 2.035;

pi := 3.143;

eps := 3.8220*10^35;

g(r) = 1-s(r)/0.06123;

j(r) = e^(-(1/2)*w(r))*(1-2*G*v(r)/(r*c^2))^.5

sys := diff(v(r), r) = 4*pi*r^2*eps/c^2, ics=v(0)=0

diff(u(r), r) = -G*(eps+u(r))*(v(r)+4*Pi*r^3*u(r)/c^2)/(c^2*(r^2-2*G*r*v(r)/c^2)),u(0)=1.3668*10^34

diff(w(r), r) = 1.485232054*10^(-28)*(v(r)+4.450600224*10^(-21)*pi*r^3*u(r))/(r^2-2*G*r*v(r)/c^2), conditions: w(0)=0,iterate it to find w(688240)=-2.05684, it solve must satistfy the both conditions.

diff(r^4*j(r)*(diff(g(r), r)), r)+4*r^3*g(r)*(diff(j(r), r)) = 0, conditions dg(r)/dr =0  at r=0, g(688240) =0.87214

diff(J(r), r) = (8*pi*(1/3))*(eps/c^2+u(r)/c^2)*(g(r)*j(r).(r^4))/(1-2*G*v(r)/(r*c^2)) condition J(0)=0.

$$\textbf{x}' = \begin{bmatrix} -4 & -2 \\ 3 & 1 \end{bmatrix} \begin{bmatrix} x \\ y \end{bmatrix}+\begin{bmatrix} -t \\ -2t-1 \end{bmatrix},\textbf{x}(0)=\begin{bmatrix} 3 \\ -5 \end{bmatrix}$$

As I know firstly, when the matrix is denoted by $A$, we must compute $e^{At}$ by diagonalizing $A$: if $A=PDP^{-1}$ for a diagonal $D$ then $e^{At} = P e^{Dt} P^{-1}$ where $e^{Dt}$ is a diagonal matrix with $(e^{Dt})_{ii} = e^{D_{ii} t}$...
 
How can I write The Maple code? maple.stackexchange)

restart: with(LinearAlgebra):

A := Matrix(2,2,[-4,-2,3,1]);

....

Dear sir in this problem should accept five boundaryconditions but it is not working for five boundary conditions and showing the following error please can you tell why it is like this ??

Error, (in dsolve/numeric/bvp/convertsys) too many boundary conditions: expected 4, got 5
Error, (in plots:-display) expecting plot structures but received: [fplt[1], fplt[2], fplt[3], fplt[4], fplt[5], fplt[6], fplt[7]]
Error, (in plots:-display) expecting plot structures but received: [tplt[1], tplt[2], tplt[3], tplt[4], tplt[5], tplt[6], tplt[7]]
 

and for the progam please check the following link

stretching_cylinder_new1.mw

Dear all

I am trying to solve system of ODEs by TWS command for traveling wave solution, but an error is showing. When I enter sinlge ODE or PDE the command does not show any error. Why it is showing error for system of ODEs ?


 

with(PDEtools, TWSolutions, declare):

with(DEtools):

Sys := {125*xi^3*(diff(f(xi), xi, xi, xi))-90*f(xi)*xi*(diff(h(xi), xi))-180*h(xi)*xi*(diff(f(xi), xi))+750*xi^2*(diff(f(xi), xi, xi))-180*f(xi)*h(xi)+830*xi*(diff(f(xi), xi))+80*f(xi)-108*(diff(h(xi), xi)), 15*f(xi)*xi*(diff(f(xi), xi))+6*f(xi)^2+10*xi*(diff(g(xi), xi))+8*g(xi)+6*(diff(f(xi), xi)) = 0, 5*f(xi)*xi*(diff(g(xi), xi))+10*g(xi)*xi*(diff(f(xi), xi))+8*f(xi)*g(xi)+10*xi*(diff(h(xi), xi))+12*h(xi)+6*(diff(g(xi), xi)) = 0}

TWSolutions(Sys)

Error, (in pdsolve/sys/info) required an indication of the solving variables for the given system

 

``


 

Download ODEs.mw

I try to solve numerically a boundary VP for ODE with different order of discontinuity of right part.

Say, the following BVP is given:

y''(x)+y'(x)+y(x)=F(x)

y(0)=1, y(2)=1

Let's use piecewise right part

F  := piecewise(x<=1, -x, x>1, 2*x+(x-1)^2)

plot(piecewise(x<=1, -x, x>1, 2*x+(x-1)^2), x=0..2,thickness=5)

The function

piecewise(x<=1, 1-x, x>1, (x-1)^2)

plot(piecewise(x<=1, 1-x, x>1, (x-1)^2), x=0..2, color=blue,thickness=5)

as obviuos, satisfies the BVP exclung the point x=1, where its 1st and 2nd derivatives are discontinuos.

Numerical solution

N0:=6:
As:=dsolve([diff(y(x), x$2)+diff(y(x), x)+y(x)=F,  y(0)=1, y(2)=1], y(x), type=numeric, output = Array([seq(2.0*k/N0, k=0..N0)]), 'maxmesh'=500, 'abserr'=1e-3):

provides the solution essentially different to exact one described above:

But if to use the right part

F := piecewise(x<=1, x^2+x+2, x>1, -x^2+x)

plot(piecewise(x<=1, x^2+x+2, x>1, -x^2+x), x=0..2, color=blue,thickness=5)

for which the function

piecewise(x<=1, 1-x+x^2, x>1, -1+3*x-x^2)

plot(piecewise(x<=1, 1-x+x^2, x>1, -1+3*x-x^2), x=0..2, thickness=5)

satisfies the BVP excluding x=1, where this function has discontinuity of 2nd derivative only, the corresponding numerical solution is very similar to this exact solution:

This reason of the difference between these two cases is clear. In the first case both 1st and 2nd derivatives are discontiuos, while in the second one -- 1st derivative is contiuos.

I wonder, if there are numerical methods, implemeted in Maple, for numerical solution of the first type BVP with non-smooth right part?

what does it means and what it will do. Can some one help me for solving this

Shootlib := "C:/Shoot9"; libname := Shootlib, libname; with(Shoot);

while i m receiving the following message:

"Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received Shoot "

Full program is :

restart; Shootlib := "C:/Shoot9"; libname := Shootlib, libname; with(Shoot);
with(plots):
N1 := 1.0; N2 := 2.0; N3 := .5; Bt := 6; Re_m := N1*Bt; gamma1 := 1;
FNS := {f(eta), fp(eta), fpp(eta), g(eta), gp(eta), m(eta), mp(eta), n(eta), np(eta), fppp(eta)};
ODE := {diff(f(eta), eta) = fp(eta),
        diff(fp(eta), eta) = fpp(eta),
        diff(fpp(eta), eta) = fppp(eta),
        diff(g(eta), eta) = gp(eta),
        diff(gp(eta), eta) = N1*(2.*g(eta)+(eta-2.*f(eta))*gp(eta)+2.*g(eta)*fp(eta)+2.*N2*N3*(m(eta)*np(eta)-n(eta)*mp(eta))),
        diff(m(eta), eta) = mp(eta),
        diff(mp(eta), eta) = Re_m*(m(eta)+(eta-2.*f(eta))*mp(eta)+2.*m(eta)*fp(eta)),
        diff(n(eta), eta) = np(eta),
        diff(np(eta), eta) = Re_m*(2.*n(eta)+(eta-2.*f(eta))*np(eta)+2.*N2/N3*m(eta)*gp(eta)),
        diff(fppp(eta), eta) = N1*(3.*fpp(eta)+(eta-2.*f(eta))*fppp(eta)-2.*N2*N2*m(eta)*(diff(mp(eta), eta)))
       }:
   
blt := 1.0;
IC := { f(0) = 0,
        fp(0) = 0,
        fpp(0) = alpha1,
        g(0) = 1,
        gp(0) = beta1,
        m(0) = 0,
        mp(0) = beta2,
        n(0) = 0,
        np(0) = beta3,
        fppp(0) = alpha2
      };
BC := { f(blt) = .5,
        fp(blt) = 0,
        g(blt) = 0,
        m(blt) = 1,
        n(blt) = 1};
infolevel[shoot] := 1;
 

Hi, 

I'm being unable to plot the solution to the equation

schro := {-(diff(psi(x), x, x))+(2*a*b*x^4+a^2*x^6+(b^2-a*(2*p+3))*x^2-(2*p+1)*b)*psi(x) = 0};

for the special cases of a=b=1 and p=0 and p=1

I've used dsolve and am getting Heun functions :-(

The claim is that the solutions come out to be exponentials of the form:

psi(x)=(x^p)*exp(-(a*(x^4))/4 - (b*(x^2))/2)

thanks in advance

Am trying to valid a research work done by kuiken(1968)

Kuiken_(1968).pdf

where we have this two eauations:

restart;
Digits := 35;
with(ODETools);
with(student);
with(plots);
inf := 4;
equ1 := diff(f[0](eta), `$`(eta, 3))+theta[0](eta);
equ2 := diff(theta[0](eta), `$`(eta, 2))+3*f[0](eta)*(diff(theta[0](eta), eta));
Bcs1 := f[0](0) = 0, (D(f[0]))(0) = 0, theta[0](0) = 1, theta[0](inf) = 0, (D(D(f[0])))(inf) = 0;
S1 := dsolve({Bcs1, equ1, equ2}, {f[0](eta), theta[0](eta)}, type = numeric, method = bvp[midrich]);
proc(x_bvp)  ...  end;
S1(0);
[                            d                   
[eta = 0., f[0](eta) = 0., ----- f[0](eta) = 0., 
[                           deta                 

    d   /  d            \                                          
  ----- |----- f[0](eta)| = 0.82449782146165697398999365896678734, 
   deta \ deta          /                                          

  theta[0](eta) = 1.0000000000000000000000000000000000, 

    d                                                         ]
  ----- theta[0](eta) = -0.71098574970825563256340736114251047]
   deta                                                       ]
S1(inf);
[                                                            
[eta = 4., f[0](eta) = 1.7815670728545914261072119522795076, 
[                                                            

    d                                                      
  ----- f[0](eta) = 0.51061876174095320088291844433043562, 
   deta                                                    

    d   /  d            \                           
  ----- |----- f[0](eta)| = 0., theta[0](eta) = 0., 
   deta \ deta          /                           

    d                                                             
  ----- theta[0](eta) = -0.000054818176138173095945902421930470836
   deta                                                           

  ]
  ]
  ]
 

 

Pls, I need to find the function of the limit of f[0](eta) at eta tend to infinity. checked equation 45 of the attached document and for the two equation pls checked equation 36 and 37 for the ODE equation solved above.

Kuiken_solution for equation 36 and 37.pdf

staganation_point11.mw
 

``

restart

l := 1:

1

 

1.5

 

.5

 

[blue, green, red, yellow]

(1)

``

for j to nops(A) do R1 := 2*n/(n+1); R2 := 2*p/(n+1); R3 := 2/(n+1); sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A[j]*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(20) = 1, theta(0) = 1+s*(D(theta))(0), theta(20) = 0], numeric, method = bvp); plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol1, [eta, diff(f(eta), eta)], color = K[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal) end do:

 

 

``

l := 1:

1

 

[0, 1, 1.5]

 

.5

 

[blue, green, red, yellow]

(2)

for j to nops(n1) do R4 := 2*n1[j]/(n1[j]+1); R5 := 2*p1/(n1[j]+1); R6 := 2/(n1[j]+1); sol2 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R4*(1-(diff(f(eta), eta))^2)-M1*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr1*f(eta)*(diff(theta(eta), eta))-R5*pr1*(diff(f(eta), eta))*theta(eta)+R6*(A1*(diff(f(eta), eta))+B1*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L1+b1*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s1*(D(theta))(0), theta(7) = 0], numeric, method = bvp); plots[odeplot](sol2, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol2, sol1, [[eta, diff(f(eta), eta)], [eta, diff(f(eta), eta)]], color = K1[j], axes = boxed); tplt[j] := plots[odeplot](sol2, [[eta, theta(eta)]], color = K1[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n1))]); plots:-display([seq(tplt[j], j = 1 .. nops(n1))])

Error, (in plots/odeplot) invalid argument: sol1

 

 

 

``

``

sol1(0);

[eta = 0., f(eta) = 1.00000000000000044, diff(f(eta), eta) = .899599635987511914, diff(diff(f(eta), eta), eta) = -.200800728024976643, theta(eta) = 1.18657688243172332, diff(theta(eta), eta) = .373153764863445037]

(3)

sol1(.1)

[eta = .1, f(eta) = 1.08902313464617162, diff(f(eta), eta) = .881503890693141945, diff(diff(f(eta), eta), eta) = -.162309073227910134, theta(eta) = 1.22040489003745489, diff(theta(eta), eta) = .304232440930656767]

(4)

sol1(.2)

[eta = .2, f(eta) = 1.17641763749368966, diff(f(eta), eta) = .866916683092940898, diff(diff(f(eta), eta), eta) = -.130454301210374102, theta(eta) = 1.24759607023709362, diff(theta(eta), eta) = .240488988787701030]

(5)

sol1(.3)

[eta = .3, f(eta) = 1.19045803452309284, diff(f(eta), eta) = .579367537136023514, diff(diff(f(eta), eta), eta) = -.285675511621370782, theta(eta) = 1.18389221591022696, diff(theta(eta), eta) = .146013974567769960]

(6)

sol1(.4)

[eta = .4, f(eta) = 1.40000000000000034, diff(f(eta), eta) = 1.00000000000000022, diff(diff(f(eta), eta), eta) = -0.243774513041384287e-17, theta(eta) = .625958972186505536, diff(theta(eta), eta) = -.314549395236395634]

(7)

sol1(.5)

[eta = .5, f(eta) = 1.37026161183094430, diff(f(eta), eta) = .874752886901313142, diff(diff(f(eta), eta), eta) = .345911467377074400, theta(eta) = .432494259338694842, diff(theta(eta), eta) = -.382764248064397461]

(8)

sol1(.6)

[eta = .6, f(eta) = 1.36678221814533528, diff(f(eta), eta) = .771028661281065508, diff(diff(f(eta), eta), eta) = .407805382194403932, theta(eta) = .876413930517023876, diff(theta(eta), eta) = -.197648778495384870]

(9)

sol1(2)

[eta = 2., f(eta) = 2.66120522956795602, diff(f(eta), eta) = .991532161353848585, diff(diff(f(eta), eta), eta) = 0.251405465681268682e-1, theta(eta) = .635967939441598018, diff(theta(eta), eta) = -.144641270049362308]

(10)

``

``

``

 

``

``

NULL


 

Download staganation_point11.mw
 

``

restart

l := 1:

1

 

1.5

 

.5

 

[blue, green, red, yellow]

(1)

``

for j to nops(A) do R1 := 2*n/(n+1); R2 := 2*p/(n+1); R3 := 2/(n+1); sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A[j]*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(20) = 1, theta(0) = 1+s*(D(theta))(0), theta(20) = 0], numeric, method = bvp); plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol1, [eta, diff(f(eta), eta)], color = K[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal) end do:

 

 

``

l := 1:

1

 

[0, 1, 1.5]

 

.5

 

[blue, green, red, yellow]

(2)

for j to nops(n1) do R4 := 2*n1[j]/(n1[j]+1); R5 := 2*p1/(n1[j]+1); R6 := 2/(n1[j]+1); sol2 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R4*(1-(diff(f(eta), eta))^2)-M1*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr1*f(eta)*(diff(theta(eta), eta))-R5*pr1*(diff(f(eta), eta))*theta(eta)+R6*(A1*(diff(f(eta), eta))+B1*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L1+b1*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s1*(D(theta))(0), theta(7) = 0], numeric, method = bvp); plots[odeplot](sol2, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol2, sol1, [[eta, diff(f(eta), eta)], [eta, diff(f(eta), eta)]], color = K1[j], axes = boxed); tplt[j] := plots[odeplot](sol2, [[eta, theta(eta)]], color = K1[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n1))]); plots:-display([seq(tplt[j], j = 1 .. nops(n1))])

Error, (in plots/odeplot) invalid argument: sol1

 

 

 

``

``

sol1(0);

[eta = 0., f(eta) = 1.00000000000000044, diff(f(eta), eta) = .899599635987511914, diff(diff(f(eta), eta), eta) = -.200800728024976643, theta(eta) = 1.18657688243172332, diff(theta(eta), eta) = .373153764863445037]

(3)

sol1(.1)

[eta = .1, f(eta) = 1.08902313464617162, diff(f(eta), eta) = .881503890693141945, diff(diff(f(eta), eta), eta) = -.162309073227910134, theta(eta) = 1.22040489003745489, diff(theta(eta), eta) = .304232440930656767]

(4)

sol1(.2)

[eta = .2, f(eta) = 1.17641763749368966, diff(f(eta), eta) = .866916683092940898, diff(diff(f(eta), eta), eta) = -.130454301210374102, theta(eta) = 1.24759607023709362, diff(theta(eta), eta) = .240488988787701030]

(5)

sol1(.3)

[eta = .3, f(eta) = 1.19045803452309284, diff(f(eta), eta) = .579367537136023514, diff(diff(f(eta), eta), eta) = -.285675511621370782, theta(eta) = 1.18389221591022696, diff(theta(eta), eta) = .146013974567769960]

(6)

sol1(.4)

[eta = .4, f(eta) = 1.40000000000000034, diff(f(eta), eta) = 1.00000000000000022, diff(diff(f(eta), eta), eta) = -0.243774513041384287e-17, theta(eta) = .625958972186505536, diff(theta(eta), eta) = -.314549395236395634]

(7)

sol1(.5)

[eta = .5, f(eta) = 1.37026161183094430, diff(f(eta), eta) = .874752886901313142, diff(diff(f(eta), eta), eta) = .345911467377074400, theta(eta) = .432494259338694842, diff(theta(eta), eta) = -.382764248064397461]

(8)

sol1(.6)

[eta = .6, f(eta) = 1.36678221814533528, diff(f(eta), eta) = .771028661281065508, diff(diff(f(eta), eta), eta) = .407805382194403932, theta(eta) = .876413930517023876, diff(theta(eta), eta) = -.197648778495384870]

(9)

sol1(2)

[eta = 2., f(eta) = 2.66120522956795602, diff(f(eta), eta) = .991532161353848585, diff(diff(f(eta), eta), eta) = 0.251405465681268682e-1, theta(eta) = .635967939441598018, diff(theta(eta), eta) = -.144641270049362308]

(10)

``

``

``

 

``

``

NULL


 

Download staganation_point11.mw
 

``

restart

l := 1:

1

 

1.5

 

.5

 

[blue, green, red, yellow]

(1)

``

for j to nops(A) do R1 := 2*n/(n+1); R2 := 2*p/(n+1); R3 := 2/(n+1); sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A[j]*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(20) = 1, theta(0) = 1+s*(D(theta))(0), theta(20) = 0], numeric, method = bvp); plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol1, [eta, diff(f(eta), eta)], color = K[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal) end do:

 

 

``

l := 1:

1

 

[0, 1, 1.5]

 

.5

 

[blue, green, red, yellow]

(2)

for j to nops(n1) do R4 := 2*n1[j]/(n1[j]+1); R5 := 2*p1/(n1[j]+1); R6 := 2/(n1[j]+1); sol2 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R4*(1-(diff(f(eta), eta))^2)-M1*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr1*f(eta)*(diff(theta(eta), eta))-R5*pr1*(diff(f(eta), eta))*theta(eta)+R6*(A1*(diff(f(eta), eta))+B1*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L1+b1*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s1*(D(theta))(0), theta(7) = 0], numeric, method = bvp); plots[odeplot](sol2, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol2, sol1, [[eta, diff(f(eta), eta)], [eta, diff(f(eta), eta)]], color = K1[j], axes = boxed); tplt[j] := plots[odeplot](sol2, [[eta, theta(eta)]], color = K1[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n1))]); plots:-display([seq(tplt[j], j = 1 .. nops(n1))])

Error, (in plots/odeplot) invalid argument: sol1

 

 

 

``

``

sol1(0);

[eta = 0., f(eta) = 1.00000000000000044, diff(f(eta), eta) = .899599635987511914, diff(diff(f(eta), eta), eta) = -.200800728024976643, theta(eta) = 1.18657688243172332, diff(theta(eta), eta) = .373153764863445037]

(3)

sol1(.1)

[eta = .1, f(eta) = 1.08902313464617162, diff(f(eta), eta) = .881503890693141945, diff(diff(f(eta), eta), eta) = -.162309073227910134, theta(eta) = 1.22040489003745489, diff(theta(eta), eta) = .304232440930656767]

(4)

sol1(.2)

[eta = .2, f(eta) = 1.17641763749368966, diff(f(eta), eta) = .866916683092940898, diff(diff(f(eta), eta), eta) = -.130454301210374102, theta(eta) = 1.24759607023709362, diff(theta(eta), eta) = .240488988787701030]

(5)

sol1(.3)

[eta = .3, f(eta) = 1.19045803452309284, diff(f(eta), eta) = .579367537136023514, diff(diff(f(eta), eta), eta) = -.285675511621370782, theta(eta) = 1.18389221591022696, diff(theta(eta), eta) = .146013974567769960]

(6)

sol1(.4)

[eta = .4, f(eta) = 1.40000000000000034, diff(f(eta), eta) = 1.00000000000000022, diff(diff(f(eta), eta), eta) = -0.243774513041384287e-17, theta(eta) = .625958972186505536, diff(theta(eta), eta) = -.314549395236395634]

(7)

sol1(.5)

[eta = .5, f(eta) = 1.37026161183094430, diff(f(eta), eta) = .874752886901313142, diff(diff(f(eta), eta), eta) = .345911467377074400, theta(eta) = .432494259338694842, diff(theta(eta), eta) = -.382764248064397461]

(8)

sol1(.6)

[eta = .6, f(eta) = 1.36678221814533528, diff(f(eta), eta) = .771028661281065508, diff(diff(f(eta), eta), eta) = .407805382194403932, theta(eta) = .876413930517023876, diff(theta(eta), eta) = -.197648778495384870]

(9)

sol1(2)

[eta = 2., f(eta) = 2.66120522956795602, diff(f(eta), eta) = .991532161353848585, diff(diff(f(eta), eta), eta) = 0.251405465681268682e-1, theta(eta) = .635967939441598018, diff(theta(eta), eta) = -.144641270049362308]

(10)

``

``

``

 

``

``

NULL


 

Download staganation_point11.mw

 

in this program im trying to combine the result, but it showing some error can help me please

 

 

 

1 2 3 4 5 6 7 Last Page 3 of 29