function analyticalExampleWithJumps()

k1 = 1; k2 = 10; uhat = 0; tauhat = -1; l = 1;
f1 = @(x)(1); intf1 = @(x)(x); int2f1 = @(x)(0.5*x.^2);
f2 = @(x)(1); intf2 = @(x)(x-1/2); int2f2 = @(x)(0.5*x.^2-0.5*x+1/8);
x = linspace(0,l,100);
u = analyticalSolutionWithJumps(k1,f1,k2,f2,uhat,l,tauhat,intf1,int2f1,intf2,int2f2,x);
plot(x,u)