На самой последней странице дельта Vx и дельта Vy надо взять 
значительно меньше.   Не 12 м/с, а можно 3 м/с

Графики получились правильные. 
Только нужно в канале широты подобрать для функции Фz(t) 
такой масштаб, чтобы она была как можно больше похожа на прямую. 
А дельта омега икс со штрихом от t не должна иметь значительных 
колебаний. Когда будешь строить, поймешь. ))

Успехов )