function compar(s0,r0)
% compare les valeurs exactes avec les deux approximations : Euler et
% Runge-Kutta

clf

pas = 0.01 ;
a = 2 ;
b = 0.4 ;
c = 1 ;
d = 0.1 ;
T=5;


[X,Y] = meshgrid(0.1:pas:30, 0.1:pas:15) ;
Z = b*Y - a*log(Y) + d*X - c*log(X) ;

v = b*8 - a*log(8) + d*22 - c*log(22) 


    A = LVEuler(s0,r0,a,b,c,d,T,pas) ;
    B = LVRK(s0,r0,a,b,c,d,T,pas) ;
   
 hold on   
   
  
    plot(A(:,2),A(:,3),'r');
 % pause

    plot(B(:,2),B(:,3), 'b') ;

    pause
contour(X,Y,Z,[v v],'g');
