% väestönkasvu

clear;close all
N0=10^6;
v=435000;
N1=1564000;

% Huom! Kaikilla funktiomääritelyssä olevilla parametreillä pitää olla
% arvo, elleivät ne esiinny argumenttilistassa. (Kokeile %-merkillä.)

N=@(t,lam)N0*exp(lam*t)+(v/lam)*(exp(lam*t)-1)

% Määrättävä lam niin, että f(lam)=0, kun f on:
f=@(lam)N(1,lam)-N1

fplot(f,[0,1])
grid on;shg
format long
lam=fzero(f,0.2)

figure
fplot(@(t)N(t,lam),[0 2]);grid on;shg

vuoden2lopussa=N(2,lam)

format   % Palautetaan oletusformaatti.
%%
