Este é o projecto:
http://www.igidl.ul.pt/pmiranda/LabNum/LN09_Proj01.pdf
O que fiz até agora foi:
function M01
a0=4
b0=3
Tol=1*10^-6
[a,b,k,L]=recorrencia(a0,b0,Tol)
fprintf('k =%2.0f a = %4.5f b = %4.5f\n', k, a, b);
end
function [a,b,k,L]=recorrencia(a0,b0,Tol)
n=1;
L=sqrt((a0^2)-(b0^2))/(2*log(a0/b0));
a(n)=sqrt(a0*(a0+b0)/2);
b(n)=sqrt(b0*(a0+b0)/2);
while( abs(a(n)-L)>Tol | abs(b(n)-L)>Tol );
n=n+1;
a(n)=sqrt(a(n-1)*((a(n-1)+b(n-1))/2));
b(n)=sqrt(b(n-1)*((a(n-1)+b(n-1))/2));
end
k=n;
return
end
Mas não está a fazer nada de jeito e não sei como completar o programa, porque para ser sincera, não percebo nada disto, tive duas aulas teoricas e no primeiro projecto o prof tá a pedir isto sem ter explicado quase nada.