Páginas

procedures

CREATE
@IDC
PROCEDURE SP_INSERTACLIENTE VARCHAR(20) ,@IDN VARCHAR(20)AS
BEGIN
INSERT
INTO CLIENTE(CLIENTEID,CLIENTENOMBRE)VALUES(@IDC,@IDN)END
GO
EXECUTE
SP_INSERTACLIENTE '5','RUPEREZ' /*PARA INSERTAR A CUENTA*/CREATE
@CUID
PROCEDURE SP_INSERTACUENTA VARCHAR(20) ,@SALCU FLOAT,@SUCID VARCHAR(20)AS

ejercicio 07 newton

clc
clear all
syms x y z
f1=x*y*z-x^2+y^2-1.34
f2=x*y-z^2-0.09

metodo de jacobi para 2 incognitas

clc
clear all
syms x y
A=[5 2;1 -4]
b=[1;0]
D=[A(1,1) 0;0 A(2,2)]
L=[0 0;A(2,1) 0]
U=[0 A(1,2);0 0]
sol=inv(D)*(b-(L+U)*[x;y])
x=1;
y=2;
error=20;
tolerancia =0.001;
fprintf('%10.2f %10.2f %10.2f\n',x,y,tolerancia)


while error >tolerancia
x1=eval(sol(1));
y1=eval(sol(2));
error=((x1-x)^2+(y1-y)^2)^0.5;
fprintf('%10.5f %10.5f %10.5f\n',x,y,error)
x=x1;
y=y1;


end

algoritmo de newton

clc
clear all
syms x y
f1=x*x+y*y-1
f2=x^2-y^2+0.5
df1x=diff(f1,x);
df1y=diff(f1,y);
df2x=diff(f2,x);
df2y=diff(f2,y);


A=[df1x df1y;df2x df2y];
B=[f1;f2];
sol=[x;y]-inv(A)*B;
x=1;
y=3;
error=20;
tolerancia=0.01
fprintf('%15.3f ,%15.3f\n ',x,y)
disp(' n x y error')
disp('------------------------------------------------------------------')
n=-1;
while error>tolerancia
n=n+1;
x1=eval(sol(1));
y1=eval(sol(2));
error=((x1-x)^2+(y1-y)^2)^0.5;
fprintf('%15.0f %15.7f %15.7f %15.7f \n ',n,x,y,error)
x=x1;
y=y1;
end

METODOS DE JACOBI

%METODO JACOBI
clc
clear
tol=0.0001;
x(1)=0;
y(1)=0;
z(1)=0;
w(1)=0;
for i=2:50
x(i)=(8+2*y(i-1)-4*z(i-1)+3*w(i-1))/17;
y(i)=(10-x(i-1)+z(i-1)-2*w(i-1))/9;
z(i)=(14-2*x(i-1)+y(i-1)+w(i-1))/6;
w(i)=(24-x(i-1)-y(i-1)-z(i-1))/8;
error=sqrt((x(i)-x(i-1))^2+(y(i)-y(i-1))^2+(z(i)-z(i-1))^2+(w(i)-w(i-1))^2);
if(error<=0.0001)
break
end
end
x(i)
y(i)
z(i)
w(i)
error

punto fijo multivariable con syms

clc 
clear all 
syms x10 x20 x30
y1=cos(x20*x30)/3+1/6
y2=sqrt(x10*x10+sin(x30) + 1.06)/9-0.1
y3=(1-exp(-x10*x20))-pi/6
x10=0.1;
x20=0.1;
x30=-0.1;
tolerancia=0.00001;
error= 20;
disp(' x11 x21 x31')
while( tolerancia<error)
x11=eval(y1);  
x21=eval(y2);
x31=eval(y3);
error= sqrt(power(x11-x10,2)+power(x21-x20,2)+power(x31-x30,2));
fprintf(' %2d %15.2f %15.2f \n',x11,x21, x31)
x10=x11;
x20=x21;
x30=x31;
end
fprintf(' %2d %15.2f %15.2f \n',x11,x21, x31)

newton raphson optimizacion nro 10

clc
clear all
syms x y
f=-((x-2)^2)-x-y^2
dfx=diff(f,x);
dfy=diff(f,y);
dfxx=diff(dfx,x);
dfxy=diff(dfx,y);
dfyx=diff(dfy,x);
dfyy=diff(dfy,y);
grad=[dfx; dfy];
hess=[dfxx dfxy;dfyx dfyy];
sol=[x y]'-inv(hess)*grad;
tol=0.001;
error =20;
x=1;
y=1;
fprintf(' x0= %15.5f y0 = %15.5f \n',x,y)


disp(' x y error ')
disp('-----------------------------------------------------')
while error>tol
x1=eval(sol(1));
y1=eval(sol(2));
f1=-((x1-2)^2)-x1-y1^2
error=abs(f1)-abs(eval(f);
fprintf(' %15.5f %15.5f %15.5f\n',x1,y1,error)
x=x1;
y=y1;
end

newton raphson ejercicio nro 08

clc
clear all
syms x0 y0
f1=((x0-y0)*(x0))/((2-x0-y0)*(1-x0))-2.6
f2=((2*y0^2)/((2-x0-y0)*(x0-y0)))-3.1
df1x=diff(f1,x0);
df1y=diff(f1,y0);
df2x=diff(f2,x0);
df2y=diff(f2,y0);
error =20;
A=[df1x df1y;df2x df2y];
B=[f1;f2];
sol=[x0;y0]-inv(A)*B;
tol=0.001;
x0=0.8;
y0=0.4;
fprintf(' x0= %15.5f y0 = %15.5f \n',x0,y0)


disp(' x y error ')
disp('-----------------------------------------------------')
while error>tol
x=eval(sol(1));
y=eval(sol(2));
error=(((x-x0)^2)+(y-y0)^2)^0.5;
fprintf(' %15.5f %15.5f %15.5f\n',x,y,error)
x0=x;
y0=y;
end