Dudas propagador orbital

Me llamo Jorge Ariza. Soy estudiante de ingeniería de telecomunicaciones de la universidad de Málaga (España). Os escribo porque estoy implementando un programa para poder realizar el cálculo de la potencia emitida por un satélite en órbita. Básicamente es un propagador orbital que, a partir de los parámetros orbitales del satélite, tomados en una fecha concreta, debe proporcionar el punto subsatélite del satélite en una fecha posterior. Me he basado en el ejemplo que aparece en el libro Satellite communications de Timothy Pratt, Charles W. Bostian (no sé si sabréis cuál es). Las fórmulas y matrices que he utilizado son las que aparecen en él. El problema es que en la fecha de medición de los parámetros el programa me da el valor correcto del punto subsatélite. Sin embargo, si quiero calcular la posición media hora más tarde me da un valor erróneo. Se trata de una órbita geoestacionaria, por lo que el punto subsatélite debería ser constante (sin tener en cuenta las perturbaiones orbitales). Simplemente quería saber si crees que pueda haber algún error de concepto o que aparentemente se me haya olvidado algo, no es mi intención que hagáis el trabajo por mí. Sé que este tipo de cuestiones no son fáciles de entender por email, pero os agradeceré cualquier ayuda que me podáis prestar. Muchas gracias por adelantado. Un saludo, Jorge.
.