Matlab - Ayuda iteracion matlab

T=500;

Fa0=2000

Fb0=1000

ya=0.16

thetab=Fb0/Fa0;

r=10;

k1=(1.79E4*(exp(-20800/(1.987*T)

k2=(2.57E16*(exp(-47400/(1.987*T)))

Pa=(ya*(1-x))/(1-2*ya*x)

Pb=(ya*(thetab-3*x))/(1-2*ya*x)

Pc=(2*ya*x)/(1-2*ya*x);

A=(k1)*((Pa*((pb)^1.5)*(200^1.5))/Pc));

B=r+(k2)*(Pc/((pb)^1.5)*(200^1.5)));

clear all

T=500;

Fa0=2000;

Fb0=1000;

ya=0.16;

thetab=Fb0/Fa0;

syms x

r=10;

k1=(1.79E4*(exp(-20800/(1.987*T))));

k2=(2.57E16*(exp(-47400/(1.987*T))));

Pa=(ya*(1-x))/(1-2*ya*x);

Pb=(ya*(thetab-3*x))/(1-2*ya*x);

Pc=(2*ya*x)/(1-2*ya*x);

A=k1*((Pa*(Pb)^1.5*(200^1.5))/Pc);

B=r+k2*(Pc/((Pb)^1.5)*(200^1.5));

D=vpa(A-B,3);  %% aca puede aplicar un método iterativo para resolver por iteraciones.

x=vpa(solve(D,'x'),5)

%d=inline(D);

%Dx=d(x)

>> iteraciones

x =

0.000046279