http://www.buenastareas.com/ensayos/Fisica/47716329.html

// Programa para hallar numéricamente la posición en función del tiempo
// para un movimiento circular uniforme - Movimiento circular uniforme - Actividad 2
// Posición de encuentro de las agujasde un reloj analógico

// Condiciones iniciales
// posicion Consideremos el radio de la aguja igual a 1
// y ambas empiezan juntas en la posición correspondiente a la hora 12
// resolveremos laposición en función del tiempo para el minutero
// para la hora será el mismo movimiento pero multiplicando el tiempo por 60
// Por otro lado resolveremos una sola vuelta ya que las posiciones luegose repetirán.
// Lo hacemos asi ya que sino los pequeños errores del calculo se van acumulando y se llega a soluciones incorrectas.

xm(1)=0;
ym(1)=1;

// velocidad
// aquí sabemos que el laaguja de los minutos da una vuelta por minutos
// luego la velocidad angular w será de 1/min = 1 / 60 s
// por lo que la velocidad lineal tendrá un módulo de
// v = 2 Pi R . w = 2 Pi . 1m / 60 s = Pi/ 30 (m/s)
// otra forma de calcular es que da una vuelta de perímetro 2 Pi R en 1 minuto
// por lo que la velocidad será 2 Pi R / 1 min
vmx(1)=3.1415926535/30;
vmy(1)=0;

// aceleración
//modulo constante e igual a v*v/R
// y el módulo del vector posición es igual al radio

amx(1)= - vmx(1)*vmx(1) * xm(1);
amy(1)= - vmx(1)*vmx(1) * ym(1);

tiempo(1)=0;

// calculo una vueltacompleta del minutero en 30000 pasos
for i=1:30000
tiempo(i+1)= tiempo(i) + (60/30000);

// hora
//xh(i+1)= xh(i) + vhx(i) * 2;
//yh(i+1)= yh(i) + vhy(i) * 2;

//vhx(i+1)= vhx(i) +ahx(i) * 2;
//vhy(i+1)= vhy(i) + ahy(i) * 2;

//ahx(i+1)= - vhx(1)*vhx(1) * xh(i+1);
//ahy(i+1)= - vhx(1)*vhx(1) * yh(i+1);

// minutero
xm(i+1)= xm(i) + vmx(i) * 0.002;
ym(i+1)=ym(i) + vmy(i) * 0.002;

vmx(i+1)= vmx(i) + amx(i) * 0.002;
vmy(i+1)= vmy(i) + amy(i) * 0.002;

amx(i+1)= - vmx(1)*vmx(1) * xm(i+1);
amy(i+1)= - vmx(1)*vmx(1) * ym(i+1);

end...