clf
clear all
x=linspace(-1,1,6);
y=linspace(1,3,6);
[X,Y]=meshgrid(x,y)
pause
plot(X,Y,'ok')
axis([-1.5,1.5,0.5,3.5])
shg
pause
hold on
%y'(x)=2+2x-y
XKAT=ones(size(X))
YKAT=2+2*X-Y
quiver(X,Y,XKAT,YKAT)
shg
pause
xfin=-1.1:0.01:1.4;
for c=[1,1.3,1.6,5/exp(1)]
    yfin=2*xfin+exp(-xfin)*c;
    plot(xfin,yfin,'r','linewidth',2)
    shg
    pause
end
