%% mlNl09R [2014 Harj 2 teht. 9, vaestonkasvumalli]
%  H2T9R
%        (Skandit ovat ikuinen ongelma, valtellaan niita.)
clear;close all
N0=10^6;
v=435000;
N1=1564000;

% Huom! 
% Kaikilla funktiomaaritelyssa olevilla argumenttilistan ulkopuolisilla
% parametreilla pitaa olla arvo. 
%%
N=@(t,lam)N0*exp(lam*t)+(v/lam)*(exp(lam*t)-1)

% Maarattava lam niin, etta f(lam)=0, kun f on:
f=@(lam)N(1,lam)-N1
%%
fplot(f,[0,1])    % fplot ja ezplot lienevat identtisia.
grid on;shg
format long
lam=fzero(f,0.2)
%%
format bank
vuoden2lopussa=round(N(2,lam)) % Pyoristetaan kokonaiseen ihmiseen.
%%
format   % Palautetaan oletusformaatti.
%%
% publish('H2T9R','pdf')  % Komentoikkunassa, ei skriptissa 
%                           (muuten ikuinen silmukka)