ang_head=function(Ta,La,Tb,Lb){
Ta=rad(Ta)
Tb=rad(Tb)
La=rad(La)
Lb=rad(Lb)
dL=abs(Lb-La)
x=cos(Tb)*sin(dL)
y=cos(Ta)*sin(Tb)-sin(Ta)*cos(Tb)*cos(dL)
b=atan2(x,y)
b=deg(b)
return(b)}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.