Question: How to avoid a singularity

Hi all,

I deal with the following ODE system which present a singularity when v(x)=0. I want to inject the condition "if v(x)=0 then u(x)=0" as may be a solution to avoid this singularity. If someone could help me to do that. Thanks.


ICS := { u(-2) = 4,   v(-2) = 4, (D(v))(-2) = 0 };
p:= dsolve( SYS union ICS,method=rkf45_dae,numeric,output=listprocedure);
solp := plots[odeplot](p, [[x, u(x)],[x, v(x)]], x = -5 .. 0);

Please Wait...