Matlab - Romper un bucle

% condiciones iniciales, suponemos x2=0

y1=input('Razón inicial de soluto en el gas: ');

y2=input('Razón de soluto a la salida en el gas: ');

G=input('Caudal o flujo del gas l.s. : ');

L=input('Caudal o flujo del líquido l.s. : ');

A=input('AEPT: '); % altura de relleno por plato teórico

x1=linsolve(L/G,y1-y2)

re=input('Introduzca polinomio de equilibrio en vector: '); % recta de equilibrio en razones, polinimio como [a b c]

yN=y2;

i=1;

yn=y2;

for i=1:100,

    yN=yn;

    RE=re-[0 0 yN];

    raiz=roots(RE);

    xN=raiz(2:2);

    yn=(L/G)*xN+y2;

    if yn==y1,

        break;

    end;

end;

z=i*A;

fprintf('El número de platos teóricos es de %d, por lo que la columna mide %f', i, z)

for n=1:20

   disp(n)

   if n ==10

      break

   end

end

disp(['Termino, n= ',num2str(n)])

1

     2

     3

     4

     5

     6

     7

     8

     9

    10

Termino, n= 10

>> soluto

Razón inicial de soluto en el gas: 2

Razón de soluto a la salida en el gas: 3

Caudal o flujo del gas l.s. : 4

Caudal o flujo del líquido l.s. : 5

AEPT: 6

Introduzca polinomio de equilibrio en vector: [4 5 6]

i=  1 yn=2.21875000000000 y1=2.00000000

i=  2 yn=2.37382322451287 y1=2.00000000

i=  3 yn=2.41912861225509 y1=2.00000000

i=  4 yn=2.41627278927592 y1=2.00000000

i=  5 yn=2.41403460776414 y1=2.00000000

i=  6 yn=2.41396273221502 y1=2.00000000

i=  7 yn=2.41405494656519 y1=2.00000000

i=  8 yn=2.41406622622462 y1=2.00000000

i=  9 yn=2.41406316047161 y1=2.00000000

i= 10 yn=2.41406239413373 y1=2.00000000

i= 11 yn=2.41406246161073 y1=2.00000000

i= 12 yn=2.41406250128027 y1=2.00000000

i= 13 yn=2.41406250180237 y1=2.00000000

i= 14 yn=2.41406250010231 y1=2.00000000

i= 15 yn=2.41406249992968 y1=2.00000000

i= 16 yn=2.41406249998931 y1=2.00000000

i= 17 yn=2.41406250000215 y1=2.00000000

i= 18 yn=2.41406250000066 y1=2.00000000

i= 19 yn=2.41406249999996 y1=2.00000000

i= 20 yn=2.41406249999997 y1=2.00000000

i= 21 yn=2.41406250000000 y1=2.00000000

i= 22 yn=2.41406250000000 y1=2.00000000

i= 23 yn=2.41406250000000 y1=2.00000000

i= 24 yn=2.41406250000000 y1=2.00000000

i= 25 yn=2.41406250000000 y1=2.00000000

i= 26 yn=2.41406250000000 y1=2.00000000

i= 27 yn=2.41406250000000 y1=2.00000000

i= 28 yn=2.41406250000000 y1=2.00000000

i= 29 yn=2.41406250000000 y1=2.00000000

i= 30 yn=2.41406250000000 y1=2.00000000

i= 31 yn=2.41406250000000 y1=2.00000000

i= 32 yn=2.41406250000000 y1=2.00000000

i= 33 yn=2.41406250000000 y1=2.00000000

i= 34 yn=2.41406250000000 y1=2.00000000

i= 35 yn=2.41406250000000 y1=2.00000000

i= 36 yn=2.41406250000000 y1=2.00000000

i= 37 yn=2.41406250000000 y1=2.00000000

i= 38 yn=2.41406250000000 y1=2.00000000

i= 39 yn=2.41406250000000 y1=2.00000000

i= 40 yn=2.41406250000000 y1=2.00000000

i= 41 yn=2.41406250000000 y1=2.00000000

i= 42 yn=2.41406250000000 y1=2.00000000

i= 43 yn=2.41406250000000 y1=2.00000000

i= 44 yn=2.41406250000000 y1=2.00000000

i= 45 yn=2.41406250000000 y1=2.00000000

i= 46 yn=2.41406250000000 y1=2.00000000

i= 47 yn=2.41406250000000 y1=2.00000000

i= 48 yn=2.41406250000000 y1=2.00000000

i= 49 yn=2.41406250000000 y1=2.00000000

i= 50 yn=2.41406250000000 y1=2.00000000

i= 51 yn=2.41406250000000 y1=2.00000000

i= 52 yn=2.41406250000000 y1=2.00000000

i= 53 yn=2.41406250000000 y1=2.00000000

i= 54 yn=2.41406250000000 y1=2.00000000

i= 55 yn=2.41406250000000 y1=2.00000000

i= 56 yn=2.41406250000000 y1=2.00000000

i= 57 yn=2.41406250000000 y1=2.00000000

i= 58 yn=2.41406250000000 y1=2.00000000

i= 59 yn=2.41406250000000 y1=2.00000000

i= 60 yn=2.41406250000000 y1=2.00000000

i= 61 yn=2.41406250000000 y1=2.00000000

i= 62 yn=2.41406250000000 y1=2.00000000

i= 63 yn=2.41406250000000 y1=2.00000000

i= 64 yn=2.41406250000000 y1=2.00000000

i= 65 yn=2.41406250000000 y1=2.00000000

i= 66 yn=2.41406250000000 y1=2.00000000

i= 67 yn=2.41406250000000 y1=2.00000000

i= 68 yn=2.41406250000000 y1=2.00000000

i= 69 yn=2.41406250000000 y1=2.00000000

i= 70 yn=2.41406250000000 y1=2.00000000

i= 71 yn=2.41406250000000 y1=2.00000000

i= 72 yn=2.41406250000000 y1=2.00000000

i= 73 yn=2.41406250000000 y1=2.00000000

i= 74 yn=2.41406250000000 y1=2.00000000

i= 75 yn=2.41406250000000 y1=2.00000000

i= 76 yn=2.41406250000000 y1=2.00000000

i= 77 yn=2.41406250000000 y1=2.00000000

i= 78 yn=2.41406250000000 y1=2.00000000

i= 79 yn=2.41406250000000 y1=2.00000000

i= 80 yn=2.41406250000000 y1=2.00000000

i= 81 yn=2.41406250000000 y1=2.00000000

i= 82 yn=2.41406250000000 y1=2.00000000

i= 83 yn=2.41406250000000 y1=2.00000000

i= 84 yn=2.41406250000000 y1=2.00000000

i= 85 yn=2.41406250000000 y1=2.00000000

i= 86 yn=2.41406250000000 y1=2.00000000

i= 87 yn=2.41406250000000 y1=2.00000000

i= 88 yn=2.41406250000000 y1=2.00000000

i= 89 yn=2.41406250000000 y1=2.00000000

i= 90 yn=2.41406250000000 y1=2.00000000

i= 91 yn=2.41406250000000 y1=2.00000000

i= 92 yn=2.41406250000000 y1=2.00000000

i= 93 yn=2.41406250000000 y1=2.00000000

i= 94 yn=2.41406250000000 y1=2.00000000

i= 95 yn=2.41406250000000 y1=2.00000000

i= 96 yn=2.41406250000000 y1=2.00000000

i= 97 yn=2.41406250000000 y1=2.00000000

i= 98 yn=2.41406250000000 y1=2.00000000

i= 99 yn=2.41406250000000 y1=2.00000000

i=100 yn=2.41406250000000 y1=2.00000000

Razón inicial de soluto en el gas: 0.111

Razón de soluto a la salida en el gas: 0.0011

Caudal o flujo del gas l.s. : 0.034

Caudal o flujo del líquido l.s. : 0.0361

AEPT: 0.31

Introduzca polinomio de equilibrio en vector: [2.7 0.685 0]

for n=1:20

disp(n)

if n >=10

break

end

end

disp(['Termino, n= ',num2str(n)])

yN=y2;

yn=1;

for i=1:200,

    RE=re-[0 0 yN];

    raiz=roots(RE);

    xN=raiz(2:2);

    yn=(L/G)*xN+y2;

    if yn>=y1,

        break;

    end;

    yN=yn;

end;

z=i*A;