clear all
global k2 k4 K M e0 f0 

t0=0
t1=10
p0=1
s0=2

ei=[s0;p0];

k2=1
k4=1
f0=1
e0=1
K=1
M=1


[t,et]=ode45('inclass',[t0,t1], ei);
plot(t,et (:,1), 'r')
hold on
plot(t,et(:,2)), legend ('s','p');