# how to plot solution of differential system

Hi, everyone! I got a question. I have a system of three Newton's equations

```
x''=-x/(x*x+y*y+z*z)^3/2-F(t)
y''=-x/(x*x+y*y+z*z)^3/2
z''=-x/(x*x+y*y+z*z)^3/2
```

I use `ode_solver()`

to solve numerically the system by using `y[0]`

for x, `y[1]`

for y, `y[2]`

for z like in example in the manual without using jacobian:

```
j_1= lambda t,y: [ y[3], y[4], y[5], -y[0]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5)-F*cos(w*t), -y[1]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5), -y[2]/(y[0]*y[0]+y[1]*y[1]+y[2]*y[2])^(1.5)]
```

How may I plot the curve (x(t), y(t)) or (x(t), z(t))?
`T.plot_solution(i=0, filename=outfile)`

allows to plot only (x(t), t), (y(t), t) or (z(t), t)
Thank you