function analyticalExample2()

k = 1; uhat = 0; tauhat = 0; l = 1;
f = @(x)(x.*(1-x)); intf = @(y)(-(y.^2.*(2*y - 3))/6); int2f = @(z)(-(z.^3.*(z - 2))/12);
% f = @(x)(sin(pi/l*x)); intf = @(x)(-cos(pi/l*x)/(pi/l)); int2f = @(x)(-sin(pi/l*x)/(pi/l)^2);
x = linspace(0,l,100);
u = analyticalSolution(k,f,uhat,l,tauhat,intf,int2f,x);
plot(x,u)