I need the python code for this system of differential equations. I have the octave code.

In 1995, 21 wolves were initially released in the Yellowstone park and their numbers have risen. In 2007, biologists estimated their number to get to 171. Study the interaction with elks.

E(t) , W(t) elk/ wolf population at time t, we set 0 to be 1995

The estimated growth rate for elks will be r= 0.0325.

This is the Octave code, I need it converted in python.

>> function xdot = f(x,t)

xdot(1) = 0.0325 * x(1) - 0.8 * x(1) * x(2);

xdot(2) = -0.6 * x(2) + 0.05 * x(1) * x(2);

endfunction

>> x = 1sode (f, [18; 0.021], (t = 1inspace (0, 50, 200)));

>> plot (t,x)

dE 0.0325E 0.8EW dt