a = 1; k = 1;
w = sqrt(a*k); s10 = 1; Q = 1;
c = 1; lambda = 1;
t = linspace(0, 2*pi/w, 100);

S = @(t) s10*cos(w*t) + Q*sin(w*t);
F = @(t) (-c*w*s10/(lambda*a))*sin(w*t) + (c*w*Q/(lambda*a))*cos(w*t);

comet(S(t),F(t));