utilización de algoritmos

Carlos Enrique Del Río Montes 200040184
A) clc;clear all;
n=input('Digite la dimension de la matriz cuadrada: ');
L=[];
for i=1:n
for j=1:i
L(i,j)=-1;

end
enddisp(L)
disp(L')
A=L*L';
disp (A)
utilizando este algoritmo se obtiene A apartir de una matriz L triangular inferior nxn, esta dimensión es introducida por el usuario. Entonces con A=L*L´QUE ES ELALGORITMO DE CHOLESKY se llega al resultado. Por ejemplo si n=4:
Digite la dimension de la matriz cuadrada: 4
-1 0 0 0
-1 -1 0 0
-1 -1 -1 0
-1 -1-1 -1

-1 -1 -1 -1
0 -1 -1 -1
0 0 -1 -1
0 0 0 -1

1 1 1 1
1 2 2 2
1 2 3 31 2 3 4
Donde se muestra L,L´y A.
1 1 1
1 2 2
1 2 3

-3 -5 -6
-2 -3 -3
-1 -1 -1
clc;clear all;n=input('Digite la dimension de la matriz cuadrada: ');
L=[];
P=zeros(n);
for i=1:n
for j=1:i
L(i,j)=-1;
end
end
disp(L)
disp(L')
for i=1:n
for j=i:-1:1
P(j,j)=1;end
end
P(1,1)=0;
P(1,n)=1;
P(n,1)=1;
P(n,n)=0;
disp(P)
A=L*L';
disp (A)
%PA=LU
%L'PA=L'LU
U=L'*P*A;
disp(U)
Para P se intercambian la primera con la ultima fila:P(1,1)=0;P(1,n)=1;P(n,1)=1;P(n,n)=0;

Para hallar U se usa PA=LU DONDE U=L’*P*A

2. a)
function [C,D]=newpoly(X, Y)
n=length(X);
D=zeros(n,n);
D(:,1)=Y';
for j=2:n
for k=j:nD(k,j)=(D(k,j-1)-D(k-1,j-1))/(X(k)-X(k-j+1));
end
end
C=D(n,n);
for k=(n-1):-1:1
C=conv(C,poly(X(k)));
m=length(C);
C(m)=C(m)+D(k,k);

end
disp(C)utilizando el algoritmo del polinomio de newton se halla el polinomio que queramos con n+1 puntos:
Por ejemplo si metemos los siguientes puntos:
clc;clear all;
newpoly([2 5 3 6 1 0 9 10 23 16], [5...