%% skrypt rozwiazuje rownianie rozniczkowe y'=ysin(x*y)

x=[0:0.01:pi];
y=sin(x);
poch=diff(y)./diff(x);
plot(x,y,'-r',[0.005:0.01:pi-0.005],poch,'-b');
