function u = numericalExample2(n)

% -u''(x) = x*(1-x) na (0,1)
% u(0) = 0
% -u'(1) = 0

A = zeros(n,n);
b = zeros(n,1);
for i=1:n-1
    for j=1:n
        if j>=2
            A(i,j) = -j*(j-1)/(j-2+i) * 1;
        end
    end
    b(i) = -1/(i+2) + 1/(i+1);
end
for j=1:n
    A(n,j) = j;
end
b(n) = 0;

alpha = A\b;
x = linspace(0,1,100);
u = 0*x;
for j=1:n
    u = u + alpha(j)*x.^j;
end
plot(x,u);