function s = an1( tt,T,omega1,n)
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here
a=n*omega1;
s=tt/a*sin(a*tt)+1/a/a*cos(a*tt);

end

