## How to execute a set of values in a program...

Hello, Sir

I tried to execute the program for a set values for more than one parameter but it is not existing, please can you do a favor for me in this case, that is how to write a program to execute set of values for more than one parameter at a time and how to plot the graph?

stretching_cylinder_new.mw

## Error, (in dsolve/numeric/BVPSolve) singularity en...

hi.

how i can dsolve this differential equations?

thanks

## ODE with boundary conditions involving limits...

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

 >

 >
 >
 >
 >
 >
 (1)
 >
 >

 >
 >
 >

## how do i solve the boundary conditions...

 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >

I found the solution of P(r) at P(0)=0, but could obtain the result of v(r) at v(1014030)=-0.4283, v(r) may have a graph such that i can goes from -0.4283 to 0.

## how to solve the challenging non linear coupled 1s...

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.

## Error in dsolve/numeric/bvp...

Hello every one ;

i have these flowing error :

 >
 >
 (1)
 >
 (2)
 >
 >
 (3)
 >
 >
 (4)
 >
 >
 >
 >
 >

any solution ?

thank you

Best regards !

## boundary conditions...

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]]

stretching_cylinder_new1.mw

## Error, (in fproc) unable to store 'HFloat(1.043007...

Please anyone, I have been battling with this problem for a while yet the error message keeps coming. Would be happy if responded to.

Thanks

 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 >
 (4)
 >
 >
 >
 (5)
 >
 (6)
 >
 (7)
 >
 (8)
 >
 (9)
 >
 >
 >
 >
 >
 >
 >
 (10)
 >
 (11)
 >
 (12)
 >
 >
 >
 (13)
 >
 (14)
 >
 >

 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 >
 (4)
 >
 >
 >
 (5)
 >
 (6)
 >
 (7)
 >
 (8)
 >
 (9)
 >
 >
 >
 >
 >
 >
 >
 (10)
 >
 (11)
 >
 (12)
 >
 >
 >
 (13)
 >
 (14)
 >
 >

## Error, (in dsolve/numeric/bvp/convertsys) too few ...

If i have boundary conditions with D(psi), i have no problem. But if i have condition with psi(infinity) (which i need), Maple says "too few boundary conditions". Maybe i make stupid mistakes, but i don't see.

restart;
assume(r, nonnegative);
ic_Re := `&psi;Re`(0) = 0, (D(`&psi;Re`))(0) = 0;
ic_Im := `&psi;Im`(0) = 0, (D(`&psi;Im`))(0) = 0;
V0 := 2.5; ERe := 1.5; EIm := 1.2; `&hbar;` := 6.582; mu := 938.27*(1/2); Q0 := 1.5; Rq := 4.5; Rv := 2.5;
Q := proc (r) options operator, arrow; -Q0*exp(-r/Rq) end proc;
V := proc (r) options operator, arrow; -V0*exp(-r/Rv) end proc;
`Eqn_&psi;Re` := -`&hbar;`^2*(diff(`&psi;Re`(r), r, r)+2*(diff(`&psi;Re`(r), r))/r)/(2*mu)-ERe*`&psi;Re`(r)+V(r)+EIm*`&psi;Re`(r) = Q(r);
`Eqn_&psi;Im` := -`&hbar;`^2*(diff(`&psi;Im`(r), r, r)+2*(diff(`&psi;Im`(r), r))/r)/(2*mu)-EIm*`&psi;Re`(r)-ERe*`&psi;Im`(r) = 0;
F := dsolve({ic_Im, ic_Re, `Eqn_&psi;Im`, `Eqn_&psi;Re`}, numeric);
plots[odeplot](F, [r, `&psi;Re`(r)], r = 0 .. 20, numpoints = 500);

plots[odeplot](F, [r, `&psi;Re`(r)], r = 0 .. 20, numpoints = 500);

restart;
assume(r, nonnegative);
ic_Re := `&psi;Re`(0) = 0, `&psi;Re`(infinity) = 0;
ic_Im := `&psi;Im`(0) = 0, `&psi;Im`(infinity) = 0;
V0 := 2.5; ERe := 1.5; EIm := 1.2; `&hbar;` := 6.582; mu := 938.27*(1/2); Q0 := 1.5; Rq := 4.5; Rv := 2.5;
Q := proc (r) options operator, arrow; -Q0*exp(-r/Rq) end proc;
V := proc (r) options operator, arrow; -V0*exp(-r/Rv) end proc;
`Eqn_&psi;Re` := -`&hbar;`^2*(diff(`&psi;Re`(r), r, r)+2*(diff(`&psi;Re`(r), r))/r)/(2*mu)-ERe*`&psi;Re`(r)+V(r)+EIm*`&psi;Re`(r) = Q(r);
`Eqn_&psi;Im` := -`&hbar;`^2*(diff(`&psi;Im`(r), r, r)+2*(diff(`&psi;Im`(r), r))/r)/(2*mu)-EIm*`&psi;Re`(r)-ERe*`&psi;Im`(r) = 0;
F := dsolve({ic_Im, ic_Re, `Eqn_&psi;Im`, `Eqn_&psi;Re`}, numeric);
Error, (in dsolve/numeric/bvp/convertsys) too few boundary conditions: expected 5, got 4

## Error, (in dsolve/numeric/bvp) singularity encount...

Hi

may every one help to me for dsolve this differentia1l equation?

error:

Error, (in dsolve/numeric/bvp) singularity encountered

Turbulent2-kw.mw

## Default method??...

Hello all..

Im sharena and i am solving ODE BVP by using maple. i used this command to solved the equation..

However, i dont know which method this programm solved my ODE. Is it RK45 method??

## Error, (in dsolve/numeric/bvp/convertsys) unable t...

restart;
Digits := 5;
with(ODETools);
with(student);
with(plots);
inf := 5;
with(LinearAlgebra);
equ1 := (1+2*n)*f(eta)*(diff(theta(eta), eta))/(1+3*n)-(diff(theta(eta), `\$`(eta, 2))) = 0;
equ2 := ((1+n)*(diff(f(eta), eta))^2/(1+3*n)-(1+2*n)*f(eta)*(diff(f(eta), eta, eta))/(1+3*n))/Bo+(diff(f(eta), `\$`(eta, 3)))^n-theta(eta) = 0;
Bcs := f(0) = 0, (D(f))(0) = 0, (D(f))(inf) = 0, theta(0) = 1, theta(inf) = 0;
Bo := 1; n := 2;
SolP1 := dsolve({Bcs, equ1, equ2}, numeric);
Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system
SolP1(0);
SolP1(0)

## fourth order equation...

> restart:with(plots):blt:=7:

> lambda:=2:m:=3:s:=1:

> Eq1:=(diff(f(eta),eta,eta,eta))+(f(eta)*(diff(f(eta),eta,eta)))-((diff(f(eta),eta))^2)+lambda*(((f(eta)*(diff(f(eta),eta,eta,eta))))-2*(diff(f(eta),eta))*(diff(f(eta),eta,eta,eta))^2)-(M/(1+m^2))*((diff(f(eta),eta)+ms))=0;

>

> Eq2:=(diff(h(eta),eta,eta))+(f(eta)*(diff(h(eta),eta)))-((diff(f(eta),eta))*(h(eta)))+lambda*(((f(eta)*(diff(h(eta),eta,eta,eta))))+(h(eta)*(diff(f(eta),eta,eta,eta)))+(diff(f(eta),eta,eta))*(diff(h(eta),eta))-2*(diff(f(eta),eta))*(diff(h(eta),eta,eta))+(M/(1+m^2)))*(m*(diff(f(eta),eta)-h))=0;

>

> Eq3:=((f(eta))*(diff(theta(eta),eta)))+Pr*((diff(theta(eta),eta,eta)))=0;

>

> bcs1 := f(0) = 0, (D(f))(0) = 1, (h)(0) = 0, (theta)(0) = 1, (D(f))(blt) = 0,  h(blt) = 0, theta(blt) = 0, D(D(f))(blt) = 0, (D(h))(blt)=0;

>

> L := [1,2,3];

> for k to 3 do R := dsolve(eval({Eq1, Eq2, Eq3, bcs1}, M = L[k]), [f(eta),h(eta),theta(eta)], numeric, maxmesh=10000, output = listprocedure);;Y || k := rhs(R[3]); YA || k := rhs(R[6]);YB || k := rhs(R[5]);YC || k := rhs(R[4]);YD || k := rhs(R[7]);end do:

>

>

>

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

hi