main.m

%Valores Tabelados
x = [0.5 1 2 3 4 5 6 7 9]';
y = [6 4.4 3.2 2.7 2 1.9 1.7 1.4 1.1]';
%Variavel para facilitar no uso do numero de euler
e = exp(1);

%Equacao eh:
%P(t) = Ae^-1.5t + Be^-0.3t + Ce^-0.05t
%Valor dos coeficientes em T = 0???
%Criacao da Matriz Z
z = [e.^(-1.5*x) e.^(-0.3*x) e.^(-0.05*x)];

fprintf('Matriz Z\n');
disp(z);

%Calculo dos Coeficientes
coef = z'*z;
fprintf('Matriz de coeficientes\n');
disp(coef);

%Calculo das Variaveis independendes
inde = z'*y;
fprintf('Variaveis independentes\n');
disp(inde);

%Calculo das variaveis desconhecidas
a = coef\inde;
fprintf('Valores encontrados\n');
disp(a);

%Calculo do Erro Padrao
Sr = sum((y - z*a).^2);
syx = sqrt(Sr/(length(x) - length(a)));
fprintf('Erro Padrao:\n');
disp(syx);

%Calculo do coeficiente de determinacao
St = sum((y - mean(y)).^2);
Determinacao = (St - Sr)/St;
fprintf('Valor do coeficiente de determinacao em porcentagem: %f%%\n', Determinacao*100);

%como o tempo eh 0, e quando euler for elevado a 0 o resultado sera 1.
%neste caso, cada coeficiente será multiplicado por 1, entao a equacao fica
%reduzida a soma dos coeficientes e do erro.
resultado = a(1) + a(2) + a(3);
fprintf('Valor estimado do numero de organismos no tempo 0: %f\n', resultado);