Matlab - ¿Cómo hacer iteraciones ?

Hola que tal, buenas tardes a todos.

Espero me puedan ayudar con un problema.

Tengo que encontrar el valor de Z mediante iteraciones pero no se como poner la condición para que me realice las iteraciones.

Les dejo lo que llevo

» P= 70

P =

70

» T= (273.15+200)

T =

473.1500

» Pc= 40.43

Pc =

40.4300

» Tc=420

Tc =

420

» Pr= P/Pc

Pr =

1.7314

» Tr= T/Tc

Tr =

1.1265

» w= .191

w =

0.1910

» Qo= .627

Qo =

0.6270

» Q1= 1.096

Q1 =

1.0960

» Q=(Qo*Q1^w)

Q =

0.6381

» f=P*Q

f =

44.6652

» R= .08314


R =

0.0831
A= .42748*((R^2*Tc^2.5)/Pc)

A =

264.2137

» B= .08664*((R*Tc)/Pc)

B =

0.0748

» Zsup= 1

Zsup =

1

» H=(B*P)/(Zsup*R*T)

H =

0.1332

» Zcal= (Zsup/(Zsup-H))-(A/(B*R*T^1.5))*(H/(Zsup+H))

Zcal =

0.6687


Ahora aquí quiero hacer que Zcal sea igual a Zsup mediante iteraciones.

P= 70;
T= (273.15+200);
Pc= 40.43;
Tc=420;
Pr= P/Pc;
Tr= T/Tc;
w= .191;
Qo= .627;
Q1= 1.096;
Q=(Qo*Q1^w);
f=P*Q;
R= .08314;
A= .42748*((R^2*Tc^2.5)/Pc);
B= .08664*((R*Tc)/Pc);

for Zsup= 1:300
H=(B*P)/(Zsup*R*T);
Zcal= (Zsup/(Zsup-H))-(A/(B*R*T^1.5))*(H/(Zsup+H));
disp(Zcal)
end

Hola Brian;

Es necesario que determines cual serán las variables dependientes e independientes?

Para que apartir de ellas puedas hacer un proceso iterativo.

Saludos
Dave C
[email protected]
http://fismatlab.blogspot.com