Practica de metodos

ELIMINACION GAUSIANA CON PIVOTEO PARCIAL
function X=gauss(A,B)
% Datos
% A matriz invertible de orden nxn
% B matriz de orden nx1
% X es la matriz solucion del sistema AX=B
% Inicializamos X y unamatriz C que sirve de almacen temporal
% Ejemplo: A=[4 8 4 0;1 5 4 -3;1 4 7 2;1 3 0 -2],b=[8;-4;10;-4]
% Ejemplo: A=[1 1 1;1 2 4;1 3 9] , b=[4;7;14]
% Llamada al programa: gauss(A,B)
format long[N,N]=size(A);
X=zeros(N,1);
C=zeros(1,N+1);
%Construyendo la matriz aumentada Aum=[A|B]
Aum=[A B];
%pivoteo parcial en la columna q-esima
for q=1:N-1
[Y,j]=max(abs(Aum(q:N,q)));
%intercambiandofilas q-esima y (j+q-1)esima
C=Aum(q,:);
Aum(q,:)=Aum(j+q-1,:);
Aum(j+q-1,:)=C;
if Aum(q,q)==0
'A es singular. No hay solucion o no es unica'
break
end
%Proceso deeliminacion en la columna q-esima
for k=q+1:N
m=Aum(k,q)/Aum(q,q);
Aum(k,q:N+1)=Aum(k,q:N+1)-m*Aum(q,q:N+1);
end
end
%sustitucion regresiva en [U|Y]X=backsup(Aum(1:N,1:N),Aum(1:N,N+1));


function X=backsup(A,B)
%Datos
% A matriz triangular superior de orden nxn
% B matriz de orden nx1
% X es la solucion del sistema lineal AX=B
format long
n=length(B);X=zeros(n,1);X(n)=B(n)/A(n,n);
for k=n-1:-1:1
X(k)=(B(k)-A(k,k+1:n)*X(k+1:n))/A(k,k);
end

Ejemplo:
Halle la cubica que pasa por los puntos (0,0), (1,1), (2,2) y (3,2)

PRACTICA:
1. Elabore un archivo scriptque grafique las curvas seno , coseno, tangente y logaritmo natural con leyenda y titulo :” Curvas trascendentes”, dominio de x=-2*pi,2*pi












2. Elabore un archivo function que permitaobtener la presión para :
donde R= la constante universal de los gases =0.518KJ/Kg.°K. T=temperatura absoluta (°K). p=presión absoluta (KPa) y v=volumen de un Kg de gas (.
Losparámetros a y b se calculan mediante:
,
Donde , en funcion de v, T
3. En un proceso de ingeniería química el vapor de agua se calienta a temperaturas lo suficientemente altas para que una...