t = linspace(0,60,10000);
x0 = 5;
w_z = 4.88E-5;
w0 = 0.34;
w_d = sqrt(w0^2+w_z^2);

q = x0*exp(-i*w_z*t);
q = q.*(cos(w_d*t)+i*w_z/w_d*sin(w_d*t));

plot(q);

xlabel('x');
ylabel('y');
title('8.09 PS#8, P2: Foucalt pendulum in the first minute');
grid on;