Matlab - iteración otra vez

Hola, perdón porque no me había explicado bien. Mi problema es el siguiente.
a se calcula con h
h se calcula con I
I se calcula con a

Por eso es un problema iterativo, ¿cómo lo resuelvo?

Muchas gracias

¿tienes algun valor inicia?

Si, tengo un valor inicial de I, le puedes llamar IO. Ese IO lo calculo con unas fórmulas muy parecidas a las que uso pero sin complicarlas tanto, es decir, es una aproximación del valor I, por eso vale como valor inicial

Prueba algo así
I=100*IO; %Valor desproporcionadamente grande para no salir en la
%1ªIteración
I=[ I I0]; %DAdo que no existen problemas de memoria concatenaremos
%todos los valores de I
while(abs(I(length(I)) - (I(length(I))-1) > 0.001)
%Restamos en valor absoluto los dos ultimos elementos de I y comparmos
%con 0.001
a=f1(I);
h=f2(a);
I = [ I f3(h)]; %volvemos a concatenar para que el while compare este ultimo
%con el de la anterior iteracion
end
A ver que tal te va.

Hola, no me ha funcionado, te escribo mis lineas de código para que puedas entender mejor el problema:

epsilon=1.34
H= 2
a2= 3
I0=atand((1+a2)/(l*(1+H))); %Valor inicial que tomo para calcular I
I=100*I0;
I=[I I0];
while abs(I(length(I)) - (I(length(I))-1)) > 0.001
a=2*I*l^2;
h=1+((1-a)*tand(I-epsilon)/l);
f3=atand((1+a)/(l*(1+h)));
I = [I f3];
end

Me dá error:

??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> nuevoeolico at 31
h=1+((1-a)*tand(I-epsilon)/l);

Muchísimas gracias por tu ayuda.

Una par de cosas previas, yo no tengo una funcion atand, pero lo que he puesto atan, si eso tan solo has de cambiarlo; otra cosita, he provao y con este valor inicial de I necesita menos operaciones que con el inicial, puede que sea una tontería pero aún así, de este modo funciona, pruebalo; A partir del comentario es una pichotada mia, puedes borrarla, es solo para ver como funciona; y lo importante para que aprendas de los errores, al ser I un vector en este caso, siempre que quieras referirte a el lo haras a toda la matriz a no ser que especifiques algún valor concreto, como veras te he cambiado un I por un
I(length(I)) ya que lo que realmente necesitas en ese caso es el último elemento de I. Yo me he inventado ele = 2 por ser el unico valro que no tenía y funciona, pero no se si esta ele esta bien o no, repasalo, e codigo es este:

epsilon=1.34;
H= 2;
a2= 3;
I0=atan((1+a2)/(l*(1+H))); %Valor inicial que tomo para calcular I
I=[0 I0];

while abs(I(length(I)) - (I(length(I)-1))) > 0.001
a=2*I*l^2;
h=1+((1-a)*tan(I(length(I))-epsilon)/l);
f3=atan((1+a)/(l*(1+h)));
I = [I f3];
end

%a partir de aqui puedes borrar despues de probar
figure; plot(I); grid on;
title(['Proceso de aproximacion de I en ' int2str(length(I)) ' iteraciones']);