From 1b5effaece00cacf460fd79ec84cfa950c6d9316 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Renato=20Rodr=C3=ADguez=20Llanos?= <51269109+renatoseb@users.noreply.github.com> Date: Sun, 19 Apr 2020 07:57:04 -0500 Subject: [PATCH] Nuevo Pull (#6) * agregeminombre * Agrege el nombre de Gabo * Gabo * dos2unix Co-authored-by: otreblan --- Crout.sci | 52 +++++++++++++++++++++++++------------------------- Doolittle.sci | 53 ++++++++++++++++++++++++++------------------------- README.md | 3 ++- d.sci | 1 - gauss.sci | 1 - 5 files changed, 55 insertions(+), 55 deletions(-) diff --git a/Crout.sci b/Crout.sci index aa08fe6..f8ac573 100644 --- a/Crout.sci +++ b/Crout.sci @@ -1,26 +1,26 @@ -function ans = Crout(A) - - // Obtener dimensiones de la matriz - [m, n] = size(A); - - // Comprobar que sea cuadrada - if( m ~= n ) then - error("La matriz debe ser cuadrada") - end - - // L se inicializa como una copia de A - L = A - - // U se inicializa como una matriz identidad de mxn - U = eye(m,n) - - for i = 1:n-1 - for j = i+1:n - U(i,j) = L(i,j)/L(i,i) // Se llena la matriz U con los multiplicadores - L(:,j) = L(:,j) - U(i,j)*L(:,i) // Operaciones elementales por columnas - end - end - - ans = struct("L", L, "U", U) - -endfunction +function ans = Crout(A) + + // Obtener dimensiones de la matriz + [m, n] = size(A); + + // Comprobar que sea cuadrada + if( m ~= n ) then + error("La matriz debe ser cuadrada") + end + + // L se inicializa como una copia de A + L = A + + // U se inicializa como una matriz identidad de mxn + U = eye(m,n) + + for i = 1:n-1 + for j = i+1:n + U(i,j) = L(i,j)/L(i,i) // Se llena la matriz U con los multiplicadores + L(:,j) = L(:,j) - U(i,j)*L(:,i) // Operaciones elementales por columnas + end + end + + ans = struct("L", L, "U", U) + +endfunction diff --git a/Doolittle.sci b/Doolittle.sci index 1487595..1570b70 100644 --- a/Doolittle.sci +++ b/Doolittle.sci @@ -1,26 +1,27 @@ -function ans = Doolittle(A) - - // Obtener dimensiones de la matriz - [m, n] = size(A); - - // Comprobar que sea cuadrada - if( m ~= n ) then - error("La matriz debe ser cuadrada") - end - - // U se inicializa como una copia de A - U = A - - // L se inicializa como una matriz identidad de mxn - L = eye(m,n) - - for i = 1:n-1 - for j = i+1:n - L(j,i) = U(j,i)/U(i,i) // Se llena la matriz L con los multiplicadores - U(j,:) = U(j,:) - L(j,i)*U(i,:) // Operaciones elementales por filas - end - end - - ans = struct("L", L, "U", U) - -endfunction +function ans = Doolittle(A) + + // Obtener dimensiones de la matriz + [m, n] = size(A); + + // Comprobar que sea cuadrada + if( m ~= n ) then + error("La matriz debe ser cuadrada") + end + + // U se inicializa como una copia de A + U = A + + // L se inicializa como una matriz identidad de mxn + L = eye(m,n) + + for i = 1:n-1 + for j = i+1:n + L(j,i) = U(j,i)/U(i,i) // Se llena la matriz L con los multiplicadores + U(j,:) = U(j,:) - L(j,i)*U(i,:) // Operaciones elementales por filas + end + end + + ans = struct("L", L, "U", U) + + +endfunction diff --git a/README.md b/README.md index 7d3519d..46292c0 100644 --- a/README.md +++ b/README.md @@ -5,4 +5,5 @@ - Frings Barrueta - Alberto Oporto - Luis Carbajal -- . +- Renato Rodríguez +- Gabriel Jimenez diff --git a/d.sci b/d.sci index 463cb94..5e63dc8 100644 --- a/d.sci +++ b/d.sci @@ -39,4 +39,3 @@ if(rA == rAb) else printf('Sin solución\n') end - diff --git a/gauss.sci b/gauss.sci index ad1316c..fe2ec52 100644 --- a/gauss.sci +++ b/gauss.sci @@ -29,4 +29,3 @@ function sol = gauss(A, b) sol = x endfunction -