x0=[2;2];
h=0.1;
ti=0;
tf=100;
[t,x]=traprule(@lotka,x0,ti,tf,h);
plot(t,x)