Advanced Search

Question: how to solve BVP

Question:how to solve BVP

adel-00 135 Maple
ode dsolve numeric differential-equations bvp singularity
February 23 2021
0

Here the code for BVP, but there is a singularity how to avoid it

a:=-0.1:alpha:=0.1:
dsys:={diff(u(r),r,r)+(1-alpha)/r+r^(alpha-2)*diff(u(r),r)=a};

res1:=dsolve(dsys union {u(-1)=0,u(1)=0},numeric,output=listprocedure);

plots[odeplot](res1,[[r,(u(r))]],-5..5,axes=boxed,tickmarks = [3, 2], color = black, thickness = 1, linestyle = solid, titlefont = [Helvetica, roman, 18], labeldirections = [horizontal, vertical], labelfont = [Helvetica, roman, 24]);

 
Please Wait...