%EX413
clear all; 
clc;

t = [0 1 2 3 4]';
h= [0.70 0.61 0.55 0.45 0.4]';

A = [ones(size(t)) t];
b = log(h);
x = A\b;

x0 = exp(x(1));
k= -x(2);

step = 4.0/10;
tplot = 4:-step:0;
plot(tplot, x0*exp(-k*tplot),'-', t, h, '*r');