backward_sub.m
33bba04159d6master

backward_sub.m
View Options

			function x = backward_sub(U, y)
			[Ni, Nj] = size(U);
			if length(y) ~= Ni
			return; % TODO, throw error
			end
			x = zeros(Nj, 1);
			x(Nj) = y(Nj) / L(Nj,Nj);
			for i=(Nj-1):1 % this should go backwards eheh xD
			x(i) = y(i);
			for j=(i+1):Nj
			x(i) = x(i) - U(i,j) * x(j);
			end
			x(i) = x(i) / U(i,i);
			end
			end