[Matlab] Skriva en enkel kod för solens upp/nedgång

Tänkte som projektarbete nu i 3an skriva en liten kod som beräknar solens upp och nedgång 1 Januari genom longitud/latiutud men känner att jag tagit lite vatten över huvudet.
Det är ingen som har några förslag på hur man kan gå tillväga?

Formler jag fått via handledare:

Deklinatiosvinkel: X=23.45sin[(360/365)*(284+(dag räknat från 1 jan))]
Timvinkel: Y=(+-)arccos[-tan(latitud)*tan(X)]
Jordens hastighet: Z=9.87sin(2b) - 7.53cos(b) - 1.5sin(b) b= (360/364)*(dag räknat från första jan - 81)
Lokala förskjutningen: K= -Z - 4(lst - longitud)

Någon som har en aning?