Originally Posted by

**jucabala789** I need some help with the following code.

%%%%%%%%%%%%%%%%%%

function [XY] = Test

V=[18.068;55.918];

K=[1005;213.2];

KAB=50;

phi1=[0.121];

phi2=1-phi1;

F = @(C) [-phi1+(C(1)/(1-K(1,1)*C(1))^2)*(1+((V(1,1)*KAB*C(2))/(V(2,1)*(1-K(2,1)*C(2)))));

-phi2+(C(2)/(1-K(2,1)*C(2))^2)*(1+((KAB*C(1))/(1-K(1,1)*C(1))))];

InitialGuess = [0.1;0.1];

Options = optimset('Display','iter');

XY = fsolve(F, InitialGuess, Options)

%%%%%%%%%%%%%%%%%%

When I use the phi1=[0.121] it's ok, but when I try to use

phi1=[0.121;0.214;0.341;0.423;0.512;0.611;0.721;0.834;92 5];

This code give an error. I need to solve this with more than one value for phi1 and I don't know how to do it.

Thanks in advance for any help.