Program AddNums(output);
// GEO->MAG Еclipse 25.10.2022,
// [LAT, LON]=[RAD]
procedure MagCoord(LAT, LON: real; var mLAT, mLON: real);
var
Teta0, LAT0, LON0, Teta, mTeta, cos_mTeta, cos_mLON: real;
begin
LAT0:=86.275*PI/180; // Nord Pole
LON0:=151.272*PI/180; // Ost
Teta0:=PI/2-LAT0;
Teta:=PI/2-LAT;
cos_mTeta:=cos(Teta0)*cos(Teta)+sin(Teta0)*sin(Teta)*cos(LON-LON0);
mTeta:=1/(cos_mTeta);
cos_mLON:=(-sin(Teta0)*cos(Teta)+cos(Teta0)*sin(Teta)*cos(LON-LON0))/sin(mTeta);
mLAT:=PI/2-mTeta;
mLON:=1/(cos_mLON);
end;
begin
end.
UHJvZ3JhbSBBZGROdW1zKG91dHB1dCk7Ci8vIEdFTy0+TUFHINCVY2xpcHNlIDI1LjEwLjIwMjIsCi8vIFtMQVQsIExPTl09W1JBRF0KcHJvY2VkdXJlIE1hZ0Nvb3JkKExBVCwgTE9OOiByZWFsOyB2YXIgbUxBVCwgbUxPTjogcmVhbCk7CnZhcgogIFRldGEwLCBMQVQwLCBMT04wLCBUZXRhLCBtVGV0YSwgY29zX21UZXRhLCBjb3NfbUxPTjogcmVhbDsKYmVnaW4KICBMQVQwOj04Ni4yNzUqUEkvMTgwOyAgLy8gTm9yZCBQb2xlCiAgTE9OMDo9MTUxLjI3MipQSS8xODA7IC8vIE9zdAogIFRldGEwOj1QSS8yLUxBVDA7CiAgVGV0YTo9UEkvMi1MQVQ7CiAgCiAgY29zX21UZXRhOj1jb3MoVGV0YTApKmNvcyhUZXRhKStzaW4oVGV0YTApKnNpbihUZXRhKSpjb3MoTE9OLUxPTjApOwogIG1UZXRhOj0xLyhjb3NfbVRldGEpOwogIGNvc19tTE9OOj0oLXNpbihUZXRhMCkqY29zKFRldGEpK2NvcyhUZXRhMCkqc2luKFRldGEpKmNvcyhMT04tTE9OMCkpL3NpbihtVGV0YSk7CiAgbUxBVDo9UEkvMi1tVGV0YTsKICBtTE9OOj0xLyhjb3NfbUxPTik7CmVuZDsKYmVnaW4KZW5kLg==