clear

fplot(@(t)75000*exp(-0.045*t)+100000, [0,40]);
hold on
grid on
fplot(@(t)(1.2*300000)/(1+29*exp(-0.08*t)), [0,40]);
legend('Linna keskuse elanikud','Äärlinna elanikud')
xlabel('aastat')
ylabel('elanike arv')

f=@(x)(1.2*300000)/(1+29*exp(-0.08*x))-75000*exp(-0.045*x)-100000;
x=33; %jooniselt
ytuletis=75000*(-0.045)*exp(-0.045*x)-360000*(-1)*(1+29*exp(-0.08*x))^-2*29*exp(-0.08*x)*(-0.08);

%newtoni meetod
xvana=33-f(x)/ytuletis
i=0;

%lõikajate meetod
while abs(x-xvana)>10^-5
    xuus=x-(x-xvana)/(f(x)-f(xvana))*f(x);
    xvana=x;
    x=xuus;
    i=i+1
end
xvana
x
i

Pu=75000*exp(-0.045*x)+100000;
Ps=300000/(1+29*exp(-0.08*x));

Pu
Ps