%%  y'-y=cos x

%    -> y' = y+cos x
%   dy funktiona:
dy = @(x,y) y+cos(x);

%   alkuarvo y(0)=1
y0 = [1];

T= 0:0.25:5;
[t,y]=ode45(dy,T,y0); % huom: inline funktio ei tarvitse @-merkkiä eteensä, toisin kuin m-tiedostoon viittaava funktio.
          
plot(t,y);

figure, hold on;
for c=-0.9:0.1:0.0
    [t,y]=ode45(dy,T,[c]);
    plot(t,y);
end
hold off;