Skip to main content
added 104 characters in body
Source Link
Nasser
  • 150.6k
  • 12
  • 162
  • 376

This is just to confirm @user21 result using Maple. Too smallI used bvp[middefer] (For BVP problem like this one) in Maple to postget solution. I am not sure what this as commentmaps to in NDSolve now if any.

restart;
ode := diff(X(x),x$2)=X(x)^(3/2)/sqrt(x);
bc  := X(0)=1,X(10)=0;
sol := dsolve([ode,bc],X(x),`numeric`,method=bvp[middefer],abserr=0.001,maxmesh=8192*4);
plots:-odeplot(sol, x=0..10);

Similar (but not exact) looking solution is generated.

Mathematica graphics

Note that Maple does not have FEM solver, only FDM.

This is just to confirm @user21 result using Maple. Too small to post this as comment.

restart;
ode := diff(X(x),x$2)=X(x)^(3/2)/sqrt(x);
bc  := X(0)=1,X(10)=0;
sol := dsolve([ode,bc],X(x),`numeric`,method=bvp[middefer],abserr=0.001,maxmesh=8192*4);
plots:-odeplot(sol, x=0..10);

Similar (but not exact) looking solution is generated.

Mathematica graphics

Note that Maple does not have FEM solver, only FDM.

This is just to confirm @user21 result using Maple. I used bvp[middefer] (For BVP problem like this one) in Maple to get solution. I am not sure what this maps to in NDSolve now if any.

restart;
ode := diff(X(x),x$2)=X(x)^(3/2)/sqrt(x);
bc  := X(0)=1,X(10)=0;
sol := dsolve([ode,bc],X(x),`numeric`,method=bvp[middefer],abserr=0.001,maxmesh=8192*4);
plots:-odeplot(sol, x=0..10);

Similar (but not exact) looking solution is generated.

Mathematica graphics

Note that Maple does not have FEM solver, only FDM.

Source Link
Nasser
  • 150.6k
  • 12
  • 162
  • 376

This is just to confirm @user21 result using Maple. Too small to post this as comment.

restart;
ode := diff(X(x),x$2)=X(x)^(3/2)/sqrt(x);
bc  := X(0)=1,X(10)=0;
sol := dsolve([ode,bc],X(x),`numeric`,method=bvp[middefer],abserr=0.001,maxmesh=8192*4);
plots:-odeplot(sol, x=0..10);

Similar (but not exact) looking solution is generated.

Mathematica graphics

Note that Maple does not have FEM solver, only FDM.