function w=q(t,l)
w=zeros(3,1);
w(3)=sin(pi*t)*pi^3+ 2* l^3 *cos(pi*t) + l^2*sin(pi*t)*pi+2*l*cos(pi*t)* pi^2;

