clear all
close all
global a b d
a=0.1; b=10; d=0.1;
u=linspace(0.1,4,100);
v=u;
[f0, g0]=izokliny(u,v,'fun_f', 'fun_g' );

plot(u,f0,u,g0)
hold on
ylim([0 3]);
x=u(1:10:end);
y=v(1:10:end);
[X,Y]=meshgrid(x,y);
ff=fun_f(Y,X);
gg=fun_g(Y,X);
quiver(X,Y,ff./sqrt(ff.^2+gg.^2),gg./sqrt(ff.^2+gg.^2),0.5)

for i=1:5
disp('wybierz mysza punkt startowy')
%[T,Y] = solver('F',tspan,y0,options,p1,p2...)
options = odeset('AbsTol',1e-9);
p=ginput(1)


[t,p] = ode45('rhs',[0 20],p(1,:),options); 
plot(p(:,1),p(:,2),'r')
p(1,:)=p(end,:);
pause(0.01)
end