function analyticalExample1()

k = 1; uhat = 0; tauhat = 1; l = 1;
f = @(x)(1); intf = @(x)(x); int2f = @(x)(0.5*x.^2);
% 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)