Untitled

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function problema
clc,clear,syms CA1 CA2 CA3
disp(' Problema de Aplicación 3.4 ')
disp(' -------------------------- ')
v=sym([CA1 CA2 CA3]);
g=sym([(10+5*CA3-5*CA1^0.5)/15;
 (15*CA1-2.5*CA2^0.5)/15;
 (15*CA2-2.5*CA3^0.5)/15]);
xo=[0.5 0.4 0.3];I=0;N=0;m=1;t=0;k=0;l=1;
format long
disp(' ')
disp(' Iter.(k) Variables Distancia')
while N==0;
 while k==0;
 t=g(m);
 while l<=3;
 t=subs(t,sym(v(l)),xo(l));l=l+1;
 end
 x(m)=single(t);m=m+1;l=1;t=0;
 if m>3
 k=2;m=1;
 end
 end
 d=norm(x-single(xo));
 if d>10^-4
 xo=sym(x);I=I+1;disp([I single(xo) d]);k=0;
 else
 disp([I+1 x d]) ;N=1;
 end
end
disp('Solución: ')
for l=1:3;
fprintf(' El valor de la %dº variable es: %f \n',l,x(l));
end 
end