From 1060ad7f8a1919082228af4a0636340c6e557c7a Mon Sep 17 00:00:00 2001 From: LinlyBoi Date: Tue, 21 Mar 2023 07:57:45 +0200 Subject: [PATCH] AI TIME --- LU.m | 26 +++++++++++++------------- 1 file changed, 13 insertions(+), 13 deletions(-) diff --git a/LU.m b/LU.m index f3cfbdd..7992fd8 100644 --- a/LU.m +++ b/LU.m @@ -1,32 +1,32 @@ -# upper matrix ballocks. +# upper matrix. clear,clc #U2 = -A(2,1)/A(1,1) * A(1, :) + A(2, :) #U3 = -A(3,1)/A(1,1) * A(1, :) + A(3, :) -#U_temp = [A(1,:); U2; U3] -#U3 = -A(3,2)/U_temp(2,2) * U_temp(2, :) + U_temp(3, :) #U = [A(1,:); U2; U3] -#L = [1, 0, 0; A(2,1)/A(1,1),1,0;(A(3,1)/A(1,1)),(A(3,2)/U_temp(2,2)),1] +#U3 = -A(3,2)/U(2,2) * U(2, :) + U(3, :) +#U = [A(1,:); U2; U3] +#L = [1, 0, 0; A(2,1)/A(1,1),1,0;(A(3,1)/A(1,1)),(A(3,2)/U(2,2)),1] -# U_temp var = A (matrix) +# U var = A (matrix) # for i < n # j = i + 1 # for j <= n - # U_temp(j, :) = -U_temp(j,i)/U_temp(j,i) * U_temp(i, :) + U_temp(j, :) + # U(j, :) = -U(j,i)/U(j,i) * U(i, :) + U(j, :) A = [5,2,4;2,1,1;1,-5,-3]; b = [17;0;13;]; - U_temp = A - L_temp = eye(3); + U = A + L = eye(3); for i = [1:length(A)-1] for j = [i+1: length(A)] - L_temp(j,i) = U_temp(j,i)/U_temp(i,i); - U_temp(j, :) = -U_temp(j,i)/U_temp(i,i) * U_temp(i, :) + U_temp(j, :); + L(j,i) = U(j,i)/U(i,i); + U(j, :) = -U(j,i)/U(i,i) * U(i, :) + U(j, :); endfor endfor - U_temp - L_temp + U + L #another for loop that's like we create another matrix for Z -L_temp*U_temp\b +L*U\b