Matlab - bucle, cadena, reescribir...

hola gente, no se como se identifica lo que quiero hacer, asiq no conozco el termino exacto que deberia usar.
el valor de la variable ro, eta, r, vr y fi ya estan determinados anteriormente y no dependen de los valores posteriores
1 cl=1; % valor inicial de cl para poder empezar la iteracion
2 cd=0.001; % valor inicial de cd para poder empezar
3 c=8*pi*r*cl;
4 Re=vr*c*ro/eta;
5 for i=1:20
6 alfa(i,1)=i-5;
7 cl(i,1)=9.91779e-08*alfa(i,1)-1.56464e-06;
8 cd(i,1)=4.47348e-03*cl(i,1)^2-1.22549e-02*cl(i,1)+1.37882e-02;
9 ct(i,1)=cl(i,1)*sin(fi)-cd(i,1)*cos(fi);
10 end
11 ctmax=max(ct);
12 b=find(ct==ctmax);
13 cl=cl(b)

como pueden ver en la linea 1, empece con un valor de cl=1, solo a efectos de poder calcular 'c' y Re. yo no puse por ahora, como van a variar las relaciones de cl-alfa y cd-alfa según el valor de Re pero ya lo voy a hacer. el asunto es q una vez q calculo el cl (q es el q maximiza ct, el de la linea 13) deberia volver arriba (linea 1) y modifiacar el valor inicial de cl (con el nuevo valor de cl obtenido en la linea 13). como puedo hacer para q el valor de cl de la linea 1 este en funcion del q obtengo en la linea 13?
PD: si se fijan, el cl de la linea 13 no depende de nada de lo de arriba (solo de alfa), pero esto es porq no esta completo el codigo, deberia haber distintas curvas de cl-alfa y cd-alfa segun el valor q tenga Re, q esta en funcion de 'cl' a traves de 'c'
bueno, alguien sabe como solucionar esto?
gracias muchas!