clear all
clf
global ep

xo=[0.01 0];
T=50;
t0=0;
N=1000;
ep=10;

ho=ruku(xo,t0,T,N);

plot(ho(:,1),ho(:,2),'r')