fork download
copy 
  1. Program AddNums(output);
  2. // GEO->MAG Еclipse 25.10.2022,
  3. // [LAT, LON]=[RAD]
  4. procedure MagCoord(LAT, LON: real; var mLAT, mLON: real);
  5. var
  6.   Teta0, LAT0, LON0, Teta, mTeta, cos_mTeta, cos_mLON: real;
  7. begin
  8.   LAT0:=86.275*PI/180;  // Nord Pole
  9.   LON0:=151.272*PI/180; // Ost
  10.   Teta0:=PI/2-LAT0;
  11.   Teta:=PI/2-LAT;
  12.  
  13.   cos_mTeta:=cos(Teta0)*cos(Teta)+sin(Teta0)*sin(Teta)*cos(LON-LON0);
  14.   mTeta:=1/(cos_mTeta);
  15.   cos_mLON:=(-sin(Teta0)*cos(Teta)+cos(Teta0)*sin(Teta)*cos(LON-LON0))/sin(mTeta);
  16.   mLAT:=PI/2-mTeta;
  17.   mLON:=1/(cos_mLON);
  18. end;
  19. begin
  20. end.
Success #stdin #stdout 0s 5460KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
Standard output is empty
https://ideone.com/ETwkpH
language:
Pascal (gpc 20070904)
created:
1 year ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!