% Delfiinit1bratk.m
%  HA 5.1.2021 
% Tehtävä 1 ratkaisu (Vakiorajoitteinen kasvu)
% Annetut vakioparametrit:
r=0.06;
t0=2017.5; y0=100;tloppu=2025.5;
% Rajoitusvakion B eri arvoja.
% Suorita skriptin ajo ensin muutoksitta, sitten
% poistamalla kommentteja vuorollaan.
% Riittää poistaa yksi kommentti, sillä viimeksi sijoitettu
% B:n arvo jää voimaan.
B=0;
% B=3;
% B=6;
% B=8
% Ratkaisufunktio:
% t:n on sallittava olla vektori.
% Käytä varmuuden vuoksi .* ja ./, vaikkei tässä tarttis.
% (Jakolakussa skalaarit, kertolaskussa toinen on skalaari.)    
 y=@(t) (y0-B./r).*exp(r.*(t-t0)) + B./r ;
yloppu=round(y(tloppu));
% Piirretään ratkaisufunktio
t=linspace(t0,tloppu);  % 100:aan osaan jako on linspace:n oletus
plot(t,y(t),'LineWidth',2)
hold on; grid on
plot(tloppu,yloppu,'.','MarkerSize',10)
% Seuraavassa title-komennossa pieni kikkailu: 
% Numeeriset B ja  yloppu muunnetaan "stringiksi", jolloin 
% ne voidaan liittää tekstiin.
title(['B=',num2str(B),', Vuonna 2025.5 Delfiinejä ',num2str(yloppu),' kpl'])
display(B)
display('Loppuarvo vuonna 2025.5:');
display(yloppu)
hold off


