function k=force(t)

     l=mod(t,2*pi);
     k=lght(pi/2,l);