function numericalExampleWithJumps(k1,f1,k2,f2,uhat,l,tauhat,n)

A = zeros(n,n);
b = zeros(n,1);
for i=1:n
    for j=1:n
%         A(i,j) = k1*i*j/(i+j-1);
        A(i,j) = k1*i*j/(i+j-1)*(1/2)^(i+j-1) + k2*i*j/(i+j-1)*(l^(i+j-1)-(1/2)^(i+j-1));
    end
%     b(i) = f1/(i+1) + tauhat;
    b(i) = f1/(i+1)*(1/2)^(i+1) + f2/(i+1)*(l^(i+1)-(1/2)^(i+1)) + tauhat*l^i;
end

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)