Diffusion Computational physics (master) Browse

compphys/exo3/forward_sub.m
33bba04159d6master

forward_sub.m
View Options

			function y = forward_sub(L, b)
			[Ni, Nj] = size(L);
			if length(b) ~= Ni
			return; % TODO, throw error
			end
			y = zeros(Nj, 1);
			y(1) = b(1) / L(1,1);
			for i=2:Nj
			y(i) = b(i);
			for j=1:(i-1)
			y(i) = y(i) - L(i,j) * y(j);
			end
			y(i) = y(i) / L(i,i);
			end
			end