# Question:Improving Phase Plot Smoothness Without Exceeding Output Limits

## Question:Improving Phase Plot Smoothness Without Exceeding Output Limits

Maple 2024

I attempted to create a phase plot using the `Odeplot` command in Maple, but the resulting plot isn't as smooth as I'd like. To improve the smoothness, I tried increasing the `numpoints` parameter to 20,000, but this resulted in the error message: `[Length of output exceeds limit of 1000000]`. What adjustments should I make to obtain a smoother phase plot without encountering this error?

```restart;
r := 1;
b := 1;
c := 0.01;
a := 0.36;
mu := 0.4;
beta := 0.75;

sys := {diff(x(t), t) = r*x(t) - b*x(t)^2 - c*x(t)*y(t) - beta*x(t)*y(t)/(a + x(t)), diff(y(t), t) = -mu*y(t) + beta*x(t)*y(t)/(a + x(t))};

ics := {x(0) = 0.2, y(0) = 0.05};

sol := dsolve(sys union ics, {x(t), y(t)}, numeric, range = 0 .. 20000, maxfun = 0, output = listprocedure, abserr = 0.1e-7, relerr = 0.1e-7)

with(DEtools);
DEplot3d(sys, {x(t), y(t)}, t = 0 .. 20000, [[x(0) = 0.2, y(0) = 0.05]], maxfun = 0, numpoints = 35000, thickness = 1, linestyle = solid);
```

```with(plots);
phaseplot := odeplot(sol, [x(t), y(t)], 0 .. 20000, numpoints = 10000, color = red, thickness = 2, axes = boxed, gridlines, title = "Phase-space Diagram");```