% Delfiinit2.m
% Parannetun delfiinimallin ratkaisutiedosto
%
% Tehtävän data:
t0=2017.5; tl=2025.5;r=0.06; % tl: "t-loppu" 
y0=100;                         % Alkuehto
b=@(t) -0.032*(2 - exp(-5*(sin(pi*t)).^2))
f=@(t,y) (r + b(t)).*y;        % y' = f(t,y)

% Lasketaan ratkaisuapproksimaatio Eulerin menetelmällä:
N=90;   % Euler-askelten lkm., muuttele huviksesi!
[T,Y]=eulerS(f,[t0,tl],y0,N); 
YL=round(Y(end)) %  Viimeinen Y-vektorin alkio.  
% Piirretään:
plot(T,Y,'.-','LineWidth',2,'Markersize',8); grid on;
xlabel('t'); ylabel('y(t)')
title(['N=',num2str(N),', Delfiinit 2025 heinäkuussa = ',num2str(YL)])
display(['N= ',num2str(N)])
display(['Delfiinit heinäkuussa 2025: ',num2str(YL)])

% Puolitetaan Euler-askel, ei piirretä uudestaan.

N=2*N;   
[T,Y]=eulerS(f,[t0,tl],y0,N); 
YL=round(Y(end));
display(['N= ',num2str(N)])
display(['Delfiinit heinäkuussa 2025: ',num2str(YL)])