% Delfiinit1b.m
%   HA 5.1.2021 
% Solmu: Diffyhtälöitä, ohjelmointa, delfiinejä
% Tehtävä 1 (Vakiorajoitteinen kasvu)

% close all  % Poista alkukommentti,
%              jos haluat sulkea edelliset grafiikkaikkunat.
%    
% Annetut vakioparametrit:
r=0.06;
t0=2017.5; y0=100;tloppu=2025.5;

% Rajoitusvakion B arvoja voit muutella.
%
% Suorita skriptin ajo ensin muutoksitta, sitten
% poistamalla kommentteja vuorollaan.
% Riittää poistaa yksi kommentti, sillä viimeksi sijoitettu
% B:n arvo jää voimaan.

B=0;   % Rajoittamaton kasvu (vanha vitsi)
% B=3;
% B=6;
% B=8
%    
% Seuraavaksi ratkaisufunktion määrittely.
% Ainoa syöteargumentti on t.
% t:n on sallittava olla vektori.
% Käytä  .* ja ./ 
%
% Tässä ainoa oma koodirivisi:
% y=@(t) .... Täydennä (ja poista kommentti edestä)
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:');
display(yloppu)
hold off  % Jos haluat seuraavat samaan kuvaan, lisää kommentti alkuun.


