function x = F3(u,f,b)
m = length(b); x = zeros(m,1);
x(m) = b(m)/u(m);
for i=m-1:-1:1
    x(i) = (b(i) - f(i)*x(i+1))/u(i);
end
