t=0:1:24*30;
M2=sin(2*pi*t/12.421);
K1=sin(2*pi*t/23.93);
S2=sin(2*pi*t/12);
O1=sin(2*pi*t/25.82);
plot(t,24.233*M2);
pause
hold
plot(t,24.233*M2+14.156*K1,'r');
pause
plot(t,24.233*M2+14.156*K1+11.284*S2,'g');
pause
plot(t,24.233*M2+14.156*K1+11.284*S2+10.051*O1,'m');
pause
hold off
plot(t,24.233*M2+14.156*K1+11.284*S2+10.051*O1);

