// 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)')
a=gca()
a.data_bounds=[0,0.45;30,0.75]

// la serie converge al valor 0.6180

pause

clf  // cierra la gráfica anterior

x(1) = 1 // nuevas condiciones iniciales

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

plot(x)
xlabel('k')
ylabel('x(k)')
a=gca()
a.data_bounds=[0,-0.1;30,1.1]
