Advanced Search

Question: How to DEplot solutions that go through singularities ? (problem : "cannot evaluate the solution further")

Question:How to DEplot solutions that go through singularities ? (problem : "cannot evaluate the solution further")

Kaluza 10 Maple 14
February 05 2011
1

Hello everybody !

I am trying to display some solutions for the differential equation :

eqdiffyx75:=

diff(y(x), x) = y(x)*(4.384477959*10^15/(1/(y(x)*x^2))^.25-1.330000000*10^16/x)/(x*(1.102240000*10^15*y(x)^2-2.192238980*10^15/(1/(y(x)*x^2))^.25))

When I make a DEplot, with the command : DEplot(eqdiffyx75, y(x), x = 0 .. 5, y = 0 .. 5)

the result is a field plot, which is perfect.

The problem occurs when I try to plot some of the solutions with the command :

DEplot(eqdiffyx75, y(x), x = 0 .. 5, y = 0 .. 5, {[y(2) = 2], [y(2) = 1.6], [y(3) = 2], [y(3) = 2.1], [y(.5) = 1.2], [y(.5) = 1.3]}, stepsize = 0.1e-1, linecolor = [green, blue, yellow, yellow, orange, orange], arrows = medium, size = 1, dirfield = [30, 30], resolution = 1000, scaling = constrained)

And the result is :


Issue

My question is : what is the correct syntax to plot the solutions that are in 2 parts and that go through a singularity ?

Thank you very much !


 
Please Wait...