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

N = 

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


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)

lam =

   0.100997929685749

format bank
vuoden2lopussa=round(N(2,lam)) % Pyoristetaan kokonaiseen ihmiseen.

vuoden2lopussa =

    2187939.00

format   % Palautetaan oletusformaatti.

publish('H2T9R','pdf') % Komentoikkunassa, ei skriptissa (muuten ikuinen silmukka)