clear all; close all; clc;
% Supongamos que tenemos un RCAI con reacción de 2º orden 

% dCa/dt= v/V*Cain-v/V*Ca-k*Ca^2

% sup. v/V = 1 min-1; Cain = 1 mol/L; k=1 L/mol.min en estado estacionario
% llamemos x = Ca

% f(x) = -x^2 - x + 1

% g(x) = (-x + 1)^1/2

x(1) = 0.5 ; % condiciones iniciales

for k = 1:29
  x(k+1) = sqrt( -x(k)+1 );
end

plot(x)
xlabel('k')
ylabel('x(k)')
axis([0,30, 0.45,0.75])

% la serie converge al valor 0.6180

pause


x(1) = 1; % nuevas condiciones iniciales

for k = 1:29
  x(k+1) = sqrt( -x(k)+1 );
end

plot(x) % al llamar nuevamente a plot, se borra el gráfico anterior y se hace la curva nueva
xlabel('k')
ylabel('x(k)')
axis([0,30,-0.1,1.1])