How can I obtain a field plot instead of getting this particular solution plot only in Maple?

Viewed 7

I have a fourth order autonomous ordinary differential equation from which I need to obtain a field plot but I only get a line plot of the particular solution;

EQ1 := (-X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)+X(mu)^2*(diff(X(mu), $(mu, 2)))^2-X(mu)^2*(diff(X(mu), $(mu, 3)))(diff(X(mu), mu)-1)-2(diff(X(mu), mu))(diff(X(mu), mu)-1)^3+(8/3)(diff(X(mu), mu)-1)^3+2X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)^2)((3/2)(diff(X(mu), mu))(diff(X(mu), mu)-1)^2-(diff(X(mu), mu)-1)^2-3X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)+(3/2)X(mu)^2(diff(X(mu), $(mu, 2)))^2-(3/2)X(mu)^2(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1))-((3/2)(diff(X(mu), mu))(diff(X(mu), mu)-1)-(diff(X(mu), mu))+1-(3/2)X(mu)(diff(X(mu), $(mu, 2))))((2(diff(X(mu), mu))(diff(X(mu), mu)-1)-8/3(diff(X(mu), mu)-1)-2X(mu)(diff(X(mu), $(mu, 2))))(2(diff(X(mu), mu))(diff(X(mu), mu)-1)^3-(8/3)(diff(X(mu), mu)-1)^3-2X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)^2-X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)+X(mu)^2(diff(X(mu), $(mu, 2)))^2+X(mu)^2*(diff(X(mu), $(mu, 3)))(diff(X(mu), mu)-1))-2X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)^4+2X(mu)(diff(X(mu), $(mu, 2)))((diff(X(mu), mu)-1)^4+(diff(X(mu), mu)-1)^3-X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)^2)+2X(mu)^2*(diff(X(mu), $(mu, 3)))(diff(X(mu), mu)-1)^3-X(mu)^2(diff(X(mu), $(mu, 3)))(diff(X(mu), mu)-1)^2-2X(mu)(diff(X(mu), $(mu, 2)))((diff(X(mu), mu)-1)^3+(diff(X(mu), mu)-1)^2-X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1))+X(mu)(diff(X(mu), mu))(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)^2+2X(mu)^3*(diff(X(mu), $(mu, 2)))(diff(X(mu), $(mu, 3)))(diff(X(mu), mu)-1)+3X(mu)^2(diff(X(mu), $(mu, 2)))^2*((diff(X(mu), mu)-1)^2+diff(X(mu), mu)-1-X(mu)(diff(X(mu), $(mu, 2))))-X(mu)^2(diff(X(mu), mu))(diff(X(mu), $(mu, 2)))^2(diff(X(mu), mu)-1)-2X(mu)(diff(X(mu), $(mu, 3)))((diff(X(mu), mu)-1)^4+(diff(X(mu), mu)-1)^3-X(mu)(diff(X(mu), $(mu, 2)))(diff(X(mu), mu)-1)^2)-X(mu)^3(diff(X(mu), $(mu, 4)))*(diff(X(mu), mu)-1)^2) = 0

The code I run is;

> with(DEtools);
> DEplot(EQ1, X(mu), mu = 0 .. 1, X = 0 .. 1, {[X(0) = .5, (D(X))(0) = .5, ((D@@2)(X))(0) = -.5, ((D@@3)(X))(0) = -.5]})

The plot I get which turns out to be for one particular solution is this; DEplot output

0 Answers
Related