簡略的天星擇日方法真偽

freeye

天機 发表于 2016-12-17 14:49
接著延續鄙人前帖.... (抱歉 ! 周末回家，才有時間好好回帖)
鄙人前帖是以【木星】臨當地子午中天線的【正 ...

感謝天機板主講述其中重點^^

freeye

青山挨星 发表于 2016-12-19 13:36
天文数据直接用西占软件比对是比较可靠的，目前最准的算法就是Jean.Meeus编著的《Astronomical.Algorithms ...

感謝好書
不過現在比起以前很懶得學
查表最快^^

mytncr

freeye 发表于 2016-12-19 15:25
感謝好書
不過現在比起以前很懶得學
查表最快^^

查表？請問用那個表才比較凖確呢？可否推薦一下。

天機

mytncr 发表于 2016-12-19 13:08
我的結果跟你差不多是一樣。相信所用公式都是差不多。
但如果用Moira去引證，2016-12-01 是8:48。所以天 ...

看到 mytncr 先生及 青山挨星 先生，對天星排盤細節孜孜不倦，令人敬佩。
容小弟再囉唆幾句:

1. 鄙人尚未花時間瞭解兩位先生所用的公式，但僅從兩位計算的【觀測地點】、【觀測時間】參數，及計算邏輯來推斷:
    1-1 mytncr 先生似乎是用加拿大 Toronto 為觀測地點，而青山挨星先生應該是以【廣州】為觀測地點。
    1-2 加拿大 Toronto (UTC -5) 與 廣州 (UTC+8) 在時區上相差 13 個小時 ( 8 - -5 = 8 + 5 = 13), 也就是說:
            在廣州是 2016-12-01 08:00 am, 在 加拿大 Toronto 應該是 2016-11-30 19:00pm。
            在加拿大 Toronto是 2016-12-01 08:00 am, 在 廣州 應該是 2016-12-01 21:00pm。
    1-3 如果兩位先生計算出來的【天頂】結果是一樣(或接近)的數值，那代表某人 A 在中國廣州2016-12-01 08:00 am出生的七政四餘星盤
           與某人 B 在加拿大 Toronto是 2016-12-01 08:00 am 出生的七政四餘星盤的【天頂】(第10宮、官祿) 會相同，
          這有可能嗎?   答案不言自明。
    按常理推斷，一定是所使用的參數 (如: 儒略日) 有問题!!

2. 鄙人在第23樓的數據是木星【地平度數】接近【地平正南】的時間。但這不是樓主提問的答案。

3. 鄙人在第44樓的數據是木星【黃經度數】接近【天頂】(即第10宮黃道入宮度) 的時間。這才是樓主提問的答案

4. mytncr 先生用 Moria 軟件重新核算，得出結果是 2016-12-01 08:48 , 是木星【黃經度數】接近【天頂】(即第10宮黃道入宮度) 的時間 , 此結果與鄙人在第 44 樓所計算的是相同的( 2016-12-01 08:47:47, 如下圖所示, 詳細推算數據在第44樓) 。

                               
登录/注册后可看大图

水平有限，鄙見僅供參考 !!

mytncr

天機 发表于 2016-12-20 14:27
看到 mytncr 先生及 青山挨星 先生，對天星排盤細節孜孜不倦，令人敬佩。
容小弟再囉唆幾句:

我用的是來自 astroexcel.wordpress.com 提供的公式。相信它是用New York City 作凖。
我試把它改成廣州，但是有差異。Moira 是用Swiss Ephemeris計算。
Swiss Ephemeris 據說凖確度是到0.00000001度。我已用它在我寫的七政星盤內。
但有些公式，比如星曜的升，落和中天等...計算尚未完全把握得到。

mytncr

我把我的程式加上了地理位置及時區校正。
出來的結果己和天機壇主的差不多是一樣的了。

mytncr

為了提高大家對天星擇日的興趣。我把程式公開給大家下載。
歡迎各位有興趣的壇友多多支持。

freeye

我也覺得用Swiss Ephemeris即可
它有API

青山挨星

东经113度15分与东经120差6.45（差约26分钟）
北纬23度07分

如果大家统一用东经120度的北京时间，相信结论应该都是8:20:51

mytncr

青山挨星 发表于 2016-12-21 14:19
东经113度15分与东经120差6.45（差约26分钟）
北纬23度07分

你說得很對！從我發表的程式可以看到。

青山挨星

加了个选项，两种时间做比对
（1）地平时

（2）广州当地真太阳地方时

（3）北京时间

谦虚好学

学习了~~~~

青山挨星

mytncr 发表于 2016-12-21 22:22
你說得很對！從我發表的程式可以看到。

广州东经113度比东经120度的北京时间晚27分钟左右（即北京时间8:20的时候、广州才7：53），你的程式显示北京时间为8:20，广州时间反而是8:48，这个应该是不对的

木星临天顶的时刻，不管什么城市都是统一的，但是地方时间却是向西逐渐推迟
通俗来说，2016.12.1 北京时间8:20 等于 东七区东经105度四川绵阳的7:20分，广州的经度在这个2个之间，怎么可能是8:48呢？

mytncr

青山挨星 发表于 2016-12-23 16:01
广州东经113度比东经120度的北京时间晚27分钟左右（即北京时间8:20的时候、广州才7：53），你的程式显示 ...

位置越往東移看到木星到天頂的時間就早。越往西移當然要晚一些。
如果以日本的東京計算，(東經139度45分，北緯35度40分)，看到木星到天頂(即正南方)的時間就是6時59分49秒。
亦即是比北京又早一點。所以程式是正確的。

mytncr

freeye 发表于 2016-12-21 09:02
我也覺得用Swiss Ephemeris即可
它有API

請問Swiss Ephemeris的API有沒有詳細的應用資料？如果是中文的就更好。

天機

mytncr 发表于 2016-12-23 23:27
位置越往東移看到木星到天頂的時間就早。越往西移當然要晚一些。
如果以日本的東京計算，(東經139度45分 ...

請問 mytncr 及 青山挨星 兩位大師，對您計算出來的數據有多大的把握?

freeye

mytncr 发表于 2016-12-23 23:42
請問Swiss Ephemeris的API有沒有詳細的應用資料？如果是中文的就更好。

我自己好多年前的資料換幾部電腦之後...
剛才翻找找不著了
以前的文件也非常不成熟
提供出來也貽笑大方


查了網路有些資料
應該就是我自己以前的做法
提供給m大

http://blog.sina.com.cn/s/blog_1360dea890102uy53.html

mytncr

天機 发表于 2016-12-24 09:59
請問 mytncr 及 青山挨星 兩位大師，對您計算出來的數據有多大的把握?

小弟的程式所用的AstroExcel雖然是沒有Swiss Ephemeris 那樣準確。
但差度應是在一度，或一分鐘以下。我已用Stellarium和Moira和結果比較過，結果應是正確的。

天機

按用事時間 (UTC+8h), 從 2016-12-01 08:00am, 每隔一日，計算其真太陽時。
真太陽時計算的觀測地是【廣州】 (東經 113° 15' 00"", 北緯 23° 07' 00")、北京當地(東經 116° 23' 28"", 北緯 39° 54' 20")

用事時間 (UTC+8)                     用事時間真太陽時                    用事時間真太陽時
北京東八區標準時                                廣州當地                                  北京當地
                                        東經 113° 15' 00"", 北緯 23° 07' 00"             東經 116° 23' 28"", 北緯 39° 54' 20"
2016-12-01 08:00:00        2016-12-01 07:43:59.73        2016-12-01 07:56:33.60
2016-12-02 08:00:00        2016-12-02 07:43:36.87        2016-12-02 07:56:10.74
2016-12-03 08:00:00        2016-12-03 07:43:13.41        2016-12-03 07:55:47.27
2016-12-04 08:00:00        2016-12-04 07:42:49.35        2016-12-04 07:55:23.22
2016-12-05 08:00:00        2016-12-05 07:42:24.75        2016-12-05 07:54:58.61
2016-12-06 08:00:00        2016-12-06 07:41:59.61        2016-12-06 07:54:33.47
2016-12-07 08:00:00        2016-12-07 07:41:33.97        2016-12-07 07:54:07.83
2016-12-08 08:00:00        2016-12-08 07:41:07.85        2016-12-08 07:53:41.72
2016-12-09 08:00:00        2016-12-09 07:40:41.29        2016-12-09 07:53:15.16
2016-12-10 08:00:00        2016-12-10 07:40:14.32        2016-12-10 07:52:48.18
2016-12-11 08:00:00        2016-12-11 07:39:46.95        2016-12-11 07:52:20.82
2016-12-12 08:00:00        2016-12-12 07:39:19.22        2016-12-12 07:51:53.09
2016-12-13 08:00:00        2016-12-13 07:38:51.15        2016-12-13 07:51:25.02
2016-12-14 08:00:00        2016-12-14 07:38:22.78        2016-12-14 07:50:56.64
2016-12-15 08:00:00        2016-12-15 07:37:54.11        2016-12-15 07:50:27.98
2016-12-16 08:00:00        2016-12-16 07:37:25.19        2016-12-16 07:49:59.06
2016-12-17 08:00:00        2016-12-17 07:36:56.04        2016-12-17 07:49:29.91
2016-12-18 08:00:00        2016-12-18 07:36:26.69        2016-12-18 07:49:00.55
2016-12-19 08:00:00        2016-12-19 07:35:57.16        2016-12-19 07:48:31.03
2016-12-20 08:00:00        2016-12-20 07:35:27.50        2016-12-20 07:48:01.37
2016-12-21 08:00:00        2016-12-21 07:34:57.74        2016-12-21 07:47:31.60
2016-12-22 08:00:00        2016-12-22 07:34:27.90        2016-12-22 07:47:01.77
2016-12-23 08:00:00        2016-12-23 07:33:58.03        2016-12-23 07:46:31.89
2016-12-24 08:00:00        2016-12-24 07:33:28.15        2016-12-24 07:46:02.02
2016-12-25 08:00:00        2016-12-25 07:32:58.32        2016-12-25 07:45:32.18
2016-12-26 08:00:00        2016-12-26 07:32:28.55        2016-12-26 07:45:02.42
2016-12-27 08:00:00        2016-12-27 07:31:58.90        2016-12-27 07:44:32.76
2016-12-28 08:00:00        2016-12-28 07:31:29.38        2016-12-28 07:44:03.25
2016-12-29 08:00:00        2016-12-29 07:31:00.05        2016-12-29 07:43:33.92
2016-12-30 08:00:00        2016-12-30 07:30:30.94        2016-12-30 07:43:04.81
2016-12-31 08:00:00        2016-12-31 07:30:02.09        2016-12-31 07:42:35.96
2017-01-01 08:00:00        2017-01-01 07:29:33.53        2017-01-01 07:42:07.39
2017-01-02 08:00:00        2017-01-02 07:29:05.29        2017-01-02 07:41:39.16
2017-01-03 08:00:00        2017-01-03 07:28:37.42        2017-01-03 07:41:11.28
2017-01-04 08:00:00        2017-01-04 07:28:09.94        2017-01-04 07:40:43.80
2017-01-05 08:00:00        2017-01-05 07:27:42.88        2017-01-05 07:40:16.75
2017-01-06 08:00:00        2017-01-06 07:27:16.28        2017-01-06 07:39:50.15
2017-01-07 08:00:00        2017-01-07 07:26:50.17        2017-01-07 07:39:24.04
2017-01-08 08:00:00        2017-01-08 07:26:24.57        2017-01-08 07:38:58.44
2017-01-09 08:00:00        2017-01-09 07:25:59.51        2017-01-09 07:38:33.38
2017-01-10 08:00:00        2017-01-10 07:25:35.00        2017-01-10 07:38:08.87

天機

再以前例(在第44樓、第66樓 木星最接近天頂的時間, 從2016-12-01開始，每天計算)，
並列出廣州當地真太陽時及北京的真太陽時, 詳細數據如下，方便易友比對糾錯!

      木星在用事時間之後，                                                                                          北京地方真太陽時
最接近天頂的廣州當地時間(UTC+8h)             廣州【真太陽時】                                 經度=116° 23' 28.00"

2016-12-01 08:47:47.00        2016-12-01 08:31:45.98        2016-12-01 08:44:19.85
2016-12-02 08:44:28.00        2016-12-02 08:28:04.16        2016-12-02 08:40:38.02
2016-12-03 08:41:09.00        2016-12-03 08:24:21.73        2016-12-03 08:36:55.59
2016-12-04 08:37:50.00        2016-12-04 08:20:38.71        2016-12-04 08:33:12.58
2016-12-05 08:34:30.00        2016-12-05 08:16:54.15        2016-12-05 08:29:28.02
2016-12-06 08:31:10.00        2016-12-06 08:13:09.06        2016-12-06 08:25:42.92
2016-12-07 08:27:49.00        2016-12-07 08:09:22.47        2016-12-07 08:21:56.33
2016-12-08 08:24:28.00        2016-12-08 08:05:35.40        2016-12-08 08:18:09.27
2016-12-09 08:21:07.00        2016-12-09 08:01:47.90        2016-12-09 08:14:21.77
2016-12-10 08:17:45.00        2016-12-10 07:57:58.98        2016-12-10 08:10:32.85
2016-12-11 08:14:23.00        2016-12-11 07:54:09.68        2016-12-11 08:06:43.54
2016-12-12 08:11:01.00        2016-12-12 07:50:20.01        2016-12-12 08:02:53.87
2016-12-13 08:07:38.00        2016-12-13 07:46:29.00        2016-12-13 07:59:02.87
2016-12-14 08:04:15.00        2016-12-14 07:42:37.69        2016-12-14 07:55:11.56
2016-12-15 08:00:51.00        2016-12-15 07:38:45.10        2016-12-15 07:51:18.96
2016-12-16 07:57:27.00        2016-12-16 07:34:52.24        2016-12-16 07:47:26.11
2016-12-17 07:54:02.00        2016-12-17 07:30:58.16        2016-12-17 07:43:32.03
2016-12-18 07:50:37.00        2016-12-18 07:27:03.88        2016-12-18 07:39:37.75
2016-12-19 07:47:12.00        2016-12-19 07:23:09.43        2016-12-19 07:35:43.29
2016-12-20 07:43:46.00        2016-12-20 07:19:13.84        2016-12-20 07:31:47.70
2016-12-21 07:40:20.00        2016-12-21 07:15:18.14        2016-12-21 07:27:52.01
2016-12-22 07:36:53.00        2016-12-22 07:11:21.38        2016-12-22 07:23:55.25
2016-12-23 07:33:26.00        2016-12-23 07:07:24.58        2016-12-23 07:19:58.45
2016-12-24 07:29:58.00        2016-12-24 07:03:26.78        2016-12-24 07:16:00.64
2016-12-25 07:26:29.00        2016-12-25 06:59:28.01        2016-12-25 07:12:01.88
2016-12-26 07:23:01.00        2016-12-26 06:55:30.32        2016-12-26 07:08:04.18
2016-12-27 07:19:31.00        2016-12-27 06:51:30.73        2016-12-27 07:04:04.59
2016-12-28 07:16:01.00        2016-12-28 06:47:31.28        2016-12-28 07:00:05.15
2016-12-29 07:12:31.00        2016-12-29 06:43:32.02        2016-12-29 06:56:05.89
2016-12-30 07:09:00.00        2016-12-30 06:39:31.97        2016-12-30 06:52:05.84
2016-12-31 07:05:29.00        2016-12-31 06:35:32.18        2016-12-31 06:48:06.04
2017-01-01 07:01:57.00        2017-01-01 06:31:31.67        2017-01-01 06:44:05.54
2017-01-02 06:58:24.00        2017-01-02 06:27:30.49        2017-01-02 06:40:04.36
2017-01-03 06:54:51.00        2017-01-03 06:23:29.67        2017-01-03 06:36:03.53
2017-01-04 06:51:17.00        2017-01-04 06:19:28.24        2017-01-04 06:32:02.10
2017-01-05 06:47:43.00        2017-01-05 06:15:27.23        2017-01-05 06:28:01.10
2017-01-06 06:44:08.00        2017-01-06 06:11:25.67        2017-01-06 06:23:59.54
2017-01-07 06:40:33.00        2017-01-07 06:07:24.60        2017-01-07 06:19:58.47
2017-01-08 06:36:57.00        2017-01-08 06:03:23.04        2017-01-08 06:15:56.90
2017-01-09 06:33:20.00        2017-01-09 05:59:21.00        2017-01-09 06:11:54.87
2017-01-10 06:29:43.00        2017-01-10 05:55:19.52        2017-01-10 06:07:53.39
2017-01-11 06:26:05.00        2017-01-11 05:51:17.62        2017-01-11 06:03:51.48

青山挨星

天機 发表于 2016-12-24 22:30
再以前例(在第44樓、第66樓 木星最接近天頂的時間, 從2016-12-01開始，每天計算)，
並列出廣州當地真太陽 ...

从天机兄弟的数据来看，有以下1点是共同的：
真太阳时，往西推迟，即广州时间比北京时间晚，或者比东经120度的杭州时间晚；
不同点是：计算出来的数据差些差别

我采用的是许剑伟翻译的算法，即 Jean.Meeus编著的《Astronomical.Algorithms》，计算结果和西洋占星软件astrolog32 非常接近；林国清的Moira和astrolog32 采用的都是瑞士天文星历（其算法也是 Jean.Meeus编著的《Astronomical.Algorithms》），但数据处理可能有问题，有不少和astrolog32对不上。

陈沼州

mytncr 发表于 2016-12-21 06:22
為了提高大家對天星擇日的興趣。我把程式公開給大家下載。
歡迎各位有興趣的壇友多多支持。

可惜网银太少，等以后再来买下

陈沼州

已购买，谢谢分享

mytncr

天機 发表于 2016-12-24 22:30
再以前例(在第44樓、第66樓 木星最接近天頂的時間, 從2016-12-01開始，每天計算)，
並列出廣州當地真太陽 ...

小弟想請教天機壇主，你採用的是什麼公式和軟件計算的呢？看來我用的astroexcel似乎差誤比較是大了些。

青山挨星

mytncr 发表于 2016-12-27 03:51
小弟想請教天機壇主，你採用的是什麼公式和軟件計算的呢？看來我用的astroexcel似乎差誤比較是大了些。

你应该是用法错误，其原始数据应该是对的

mytncr

青山挨星 发表于 2016-12-27 13:41
你应该是用法错误，其原始数据应该是对的

可否請青山兄指正一下是那里出了問題呢？

mytncr

最近在網上收集了不少天文計算公式作為測試。發覺出來的結果都大有差異。
從一到二度或是二到六分鐘不等。主要是採用的公式不同做成。一般公式都是作為學習之用，準確度不高。
以个人所見，專業的當以"Stellarium" "NOAA" 和"Swiss Ephemeris "這幾個為最好。
但其實一點點的差異是否真的對七政四余批命和擇日有那麼大的影響嗎？
我只是初學之人，只是小弟愚見，未知各位大大的看法又如何呢？

天機

mytncr 发表于 2016-12-27 03:51
小弟想請教天機壇主，你採用的是什麼公式和軟件計算的呢？看來我用的astroexcel似乎差誤比較是大了些。

抱歉，鄙人回覆怠慢了。
個人在天文計算也經歷幾個階段。

1. NOAA 的 excel 模板，在某些計算仍然是近似值，例如: EOT 的計算即是。可惜，NOAA 的計算主要在太陽的黃赤、地平、仰角。不能解決中國七政四餘推步、天星計算及回溯古代的天象驗證。

2. Jean Meeus 先生的天文計算書籍 Astronomical Algorithms 是必讀的經典，該書從1982年初版、1985、1988、1990, 已經到第4版。但鄙人發現其中有錯誤，直到 1991 出版，才更正書中的一些錯誤。若您個人想看此書，請選購 1991 年以後版本。鄙人是購買其英文原版書及電子書。國內專家許劍偉先生有把該書做重點翻譯《天文算法》。
      * Jean Meeus 先生的天文計算( Astronomical Algorithms ) 書籍提供行星模型及計算公式，是非常好的進階書籍。鄙人乃逐一按其書中算法，撰寫數萬行的計算機程式，將其整本書的算法全部用軟件完成。這是一個苦功夫，但也非常值得。試想，在中國術數界，有人會花時間去閱讀Jean Meeus 先生英文原版的《天文算法》? 除了許劍偉先生外，在中國術數界是否有人將Jean Meeus 先生英文原版的《天文算法》用軟件實現?  
      * Jean Meeus 先生的天文計算( Astronomical Algorithms ) 書中的行星模型及公式，是採用 VSOP 行星理論。VSOP是"Variations Seculaires des Orbites Planetaires"的縮寫，由巴黎的P.Bretagnon 於1982年,發表了他的行星理論VSOP82。VSOP82由大星行(水星到海王星)的長長的週期項序列組成。給定一個行星及一個時間,對它的序列取和計算,即可獲得密切軌道參數。
Jean Meeus 先生的天文計算( Astronomical Algorithms) 書中採用 1987 年的模型, 即 VSOP87。

  VSOP82 行星理論是採用複雜的微積分及上千個多項式週期項，當不需要完全精度時,可以把計算公式在適當處截斷。幸運的是,1987年 Bretagnon 和 Francou 創建了VSOP87,它提供了直接計算行星日心座標的週期序列項。也就是說可以直接計算出:
  L, 日心黃經
  B, 日心黃緯
  R, 行星到太陽的距離

Jean Meeus 先生的天文計算( Astronomical Algorithms) 書中附錄中有列出七政四餘的 VSOP 87 的簡化週期項數值。







但是如果您按 Jean Meeus 先生的天文計算( Astronomical Algorithms) 書中附錄中列出七政四餘的 VSOP 87 的簡化週期項數值, 去計算七政四餘三王的天文位置，其精度仍然是不夠的。因為原本 VSOP87 的行星模式，是把微積分公式，用數千個週期項展開，而Jean Meeus 先生的天文計算( Astronomical Algorithms) 書中附錄中列出七政四餘的 VSOP 87 的簡化週期項數量, 最多只有 64-100 之間。

以下列出个人所使用 VSOP87 行星模型中【地球】(也就是計算太陽位置) 的完整周期項的函式及公式:
(而這僅是計算太陽行星位置的最基礎的數據。)

天機

/*

   EARTH - VSOP87  

本函式涵蓋的時間跨幅:  西元前2000 BC < Date < 西元 6000 AD
   理論值精確度: ±1 arc second

   L = Longitude in radians
   B = Latitude in radians
   R = Radius vector in AU (Astronomical Units,天文單位)

   t = (JD - 2451545) / 365250

   本原始碼參考資料: :
   Planetary Theories in Rectangular and Spherical Variables, VSOP87 Solutions
   Pierre Bretagnon, Gerard Francou, Journal of Astronomy & Astrophysics
   vol. 202, p309-p315, 1988

*/


   double Earth_L0 (double t) // 623 terms of order 0
{
   double L0 = 0;
   L0 += 1.75347045673;
   L0 += 0.03341656453*cos(4.66925680415 + 6283.0758499914*t);
   L0 += 0.00034894275*cos(4.62610242189 + 12566.1516999828*t);
   L0 += 0.00003417572*cos(2.82886579754 + 3.523118349*t);
   L0 += 0.00003497056*cos(2.74411783405 + 5753.3848848968*t);
   L0 += 0.00003135899*cos(3.62767041756 + 77713.7714681205*t);
   L0 += 0.00002676218*cos(4.41808345438 + 7860.4193924392*t);
   L0 += 0.00002342691*cos(6.13516214446 + 3930.2096962196*t);
   L0 += 0.00001273165*cos(2.03709657878 + 529.6909650946*t);
   L0 += 0.00001324294*cos(0.74246341673 + 11506.7697697936*t);
   L0 += 0.00000901854*cos(2.04505446477 + 26.2983197998*t);
   L0 += 0.00001199167*cos(1.10962946234 + 1577.3435424478*t);
   L0 += 0.00000857223*cos(3.50849152283 + 398.1490034082*t);
   L0 += 0.00000779786*cos(1.17882681962 + 5223.6939198022*t);
   L0 += 0.00000990250*cos(5.23268072088 + 5884.9268465832*t);
   L0 += 0.00000753141*cos(2.53339052847 + 5507.5532386674*t);
   L0 += 0.00000505267*cos(4.58292599973 + 18849.2275499742*t);
   L0 += 0.00000492392*cos(4.20505711826 + 775.522611324*t);
   L0 += 0.00000356672*cos(2.91954114478 + 0.0673103028*t);
   L0 += 0.00000284125*cos(1.89869240932 + 796.2980068164*t);
   L0 += 0.00000242879*cos(0.34481445893 + 5486.777843175*t);
   L0 += 0.00000317087*cos(5.84901948512 + 11790.6290886588*t);
   L0 += 0.00000271112*cos(0.31486255375 + 10977.078804699*t);
   L0 += 0.00000206217*cos(4.80646631478 + 2544.3144198834*t);
   L0 += 0.00000205478*cos(1.86953770281 + 5573.1428014331*t);
   L0 += 0.00000202318*cos(2.45767790232 + 6069.7767545534*t);
   L0 += 0.00000126225*cos(1.08295459501 + 20.7753954924*t);
   L0 += 0.00000155516*cos(0.83306084617 + 213.299095438*t);
   L0 += 0.00000115132*cos(0.64544911683 + 0.9803210682*t);
   L0 += 0.00000102851*cos(0.63599845579 + 4694.0029547076*t);
   L0 += 0.00000101724*cos(4.26679801980 + 7.1135470008*t);
   L0 += 0.00000099206*cos(6.20992926918 + 2146.1654164752*t);
   L0 += 0.00000132212*cos(3.41118292683 + 2942.4634232916*t);
   L0 += 0.00000097607*cos(0.68101342359 + 155.4203994342*t);
   L0 += 0.00000085128*cos(1.29870764804 + 6275.9623029906*t);
   L0 += 0.00000074651*cos(1.75508913300 + 5088.6288397668*t);
   L0 += 0.00000101895*cos(0.97569280312 + 15720.8387848784*t);
   L0 += 0.00000084711*cos(3.67080093031 + 71430.69561812909*t);
   L0 += 0.00000073547*cos(4.67926633877 + 801.8209311238*t);
   L0 += 0.00000073874*cos(3.50319414955 + 3154.6870848956*t);
   L0 += 0.00000078757*cos(3.03697458703 + 12036.4607348882*t);
   L0 += 0.00000079637*cos(1.80791287082 + 17260.1546546904*t);
   L0 += 0.00000085803*cos(5.98322631260 + 161000.6857376741*t);
   L0 += 0.00000056963*cos(2.78430458592 + 6286.5989683404*t);
   L0 += 0.00000061148*cos(1.81839892984 + 7084.8967811152*t);
   L0 += 0.00000069627*cos(0.83297621398 + 9437.762934887*t);
   L0 += 0.00000056116*cos(4.38694865354 + 14143.4952424306*t);
   L0 += 0.00000062449*cos(3.97763912806 + 8827.3902698748*t);
   L0 += 0.00000051145*cos(0.28306832879 + 5856.4776591154*t);
   L0 += 0.00000055577*cos(3.47006059924 + 6279.5527316424*t);
   L0 += 0.00000041036*cos(5.36817592855 + 8429.2412664666*t);
   L0 += 0.00000051605*cos(1.33282739866 + 1748.016413067*t);
   L0 += 0.00000051992*cos(0.18914947184 + 12139.5535091068*t);
   L0 += 0.00000049000*cos(0.48735014197 + 1194.4470102246*t);
   L0 += 0.00000039200*cos(6.16833020996 + 10447.3878396044*t);
   L0 += 0.00000035570*cos(1.77596889200 + 6812.766815086*t);
   L0 += 0.00000036770*cos(6.04133863162 + 10213.285546211*t);
   L0 += 0.00000036596*cos(2.56957481827 + 1059.3819301892*t);
   L0 += 0.00000033296*cos(0.59310278598 + 17789.845619785*t);
   L0 += 0.00000035954*cos(1.70875808777 + 2352.8661537718*t);
   L0 += 0.00000040938*cos(2.39850938714 + 19651.048481098*t);
   L0 += 0.00000030047*cos(2.73975124088 + 1349.8674096588*t);
   L0 += 0.00000030412*cos(0.44294464169 + 83996.84731811189*t);
   L0 += 0.00000023663*cos(0.48473622521 + 8031.0922630584*t);
   L0 += 0.00000023574*cos(2.06528133162 + 3340.6124266998*t);
   L0 += 0.00000021089*cos(4.14825468851 + 951.7184062506*t);
   L0 += 0.00000024738*cos(0.21484762138 + 3.5904286518*t);
   L0 += 0.00000025352*cos(3.16470891653 + 4690.4798363586*t);
   L0 += 0.00000022823*cos(5.22195230819 + 4705.7323075436*t);
   L0 += 0.00000021419*cos(1.42563910473 + 16730.4636895958*t);
   L0 += 0.00000021891*cos(5.55594302779 + 553.5694028424*t);
   L0 += 0.00000017481*cos(4.56052900312 + 135.0650800354*t);
   L0 += 0.00000019927*cos(5.22209149316 + 12168.0026965746*t);
   L0 += 0.00000019860*cos(5.77470242235 + 6309.3741697912*t);
   L0 += 0.00000020300*cos(0.37133792946 + 283.8593188652*t);
   L0 += 0.00000014421*cos(4.19315052005 + 242.728603974*t);
   L0 += 0.00000016225*cos(5.98837767951 + 11769.8536931664*t);
   L0 += 0.00000015077*cos(4.19567163370 + 6256.7775301916*t);
   L0 += 0.00000019124*cos(3.82219958698 + 23581.2581773176*t);
   L0 += 0.00000018888*cos(5.38626892076 + 149854.40013480789*t);
   L0 += 0.00000014346*cos(3.72355084422 + 38.0276726358*t);
   L0 += 0.00000017898*cos(2.21490566029 + 13367.9726311066*t);
   L0 += 0.00000012054*cos(2.62229602614 + 955.5997416086*t);
   L0 += 0.00000011287*cos(0.17739329984 + 4164.311989613*t);
   L0 += 0.00000013973*cos(4.40134615007 + 6681.2248533996*t);
   L0 += 0.00000013621*cos(1.88934516495 + 7632.9432596502*t);
   L0 += 0.00000012503*cos(1.13052412208 + 5.5229243074*t);
   L0 += 0.00000010498*cos(5.35909979317 + 1592.5960136328*t);
   L0 += 0.00000009803*cos(0.99948172646 + 11371.7046897582*t);
   L0 += 0.00000009220*cos(4.57138585348 + 4292.3308329504*t);
   L0 += 0.00000010327*cos(6.19982170609 + 6438.4962494256*t);
   L0 += 0.00000012003*cos(1.00351462266 + 632.7837393132*t);
   L0 += 0.00000010827*cos(0.32734523824 + 103.0927742186*t);
   L0 += 0.00000008356*cos(4.53902748706 + 25132.3033999656*t);
   L0 += 0.00000010005*cos(6.02914963280 + 5746.271337896*t);
   L0 += 0.00000008409*cos(3.29946177848 + 7234.794256242*t);
   L0 += 0.00000008006*cos(5.82145271855 + 28.4491874678*t);
   L0 += 0.00000010523*cos(0.93870455544 + 11926.2544136688*t);
   L0 += 0.00000007686*cos(3.12143640640 + 7238.67559160*t);
   L0 += 0.00000009378*cos(2.62413793196 + 5760.4984318976*t);
   L0 += 0.00000008127*cos(6.11227839253 + 4732.0306273434*t);
   L0 += 0.00000009232*cos(0.48344234496 + 522.5774180938*t);
   L0 += 0.00000009802*cos(5.24413877132 + 27511.4678735372*t);
   L0 += 0.00000007871*cos(0.99590133077 + 5643.1785636774*t);
   L0 += 0.00000008123*cos(6.27053020099 + 426.598190876*t);
   L0 += 0.00000009048*cos(5.33686323585 + 6386.16862421*t);
   L0 += 0.00000008621*cos(4.16537179089 + 7058.5984613154*t);
   L0 += 0.00000006297*cos(4.71723143652 + 6836.6452528338*t);
   L0 += 0.00000007575*cos(3.97381357237 + 11499.6562227928*t);
   L0 += 0.00000007756*cos(2.95728422442 + 23013.5395395872*t);
   L0 += 0.00000007314*cos(0.60652522715 + 11513.8833167944*t);
   L0 += 0.00000005955*cos(2.87641047954 + 6283.14316029419*t);
   L0 += 0.00000006534*cos(5.79046406784 + 18073.7049386502*t);
   L0 += 0.00000007188*cos(3.99831461988 + 74.7815985673*t);
   L0 += 0.00000007346*cos(4.38582423903 + 316.3918696566*t);
   L0 += 0.00000005413*cos(5.39199023275 + 419.4846438752*t);
   L0 += 0.00000005127*cos(2.36059551778 + 10973.55568635*t);
   L0 += 0.00000007056*cos(0.32258442532 + 263.0839233728*t);
   L0 += 0.00000006624*cos(3.66474165840 + 17298.1823273262*t);
   L0 += 0.00000006762*cos(5.91131766896 + 90955.5516944961*t);
   L0 += 0.00000004938*cos(5.73672172371 + 9917.6968745098*t);
   L0 += 0.00000005547*cos(2.45152589382 + 12352.8526045448*t);
   L0 += 0.00000005958*cos(3.32051344660 + 6283.0085396886*t);
   L0 += 0.00000004471*cos(2.06386138131 + 7079.3738568078*t);
   L0 += 0.00000006153*cos(1.45823347458 + 233141.31440436149*t);
   L0 += 0.00000004348*cos(4.42338625285 + 5216.5803728014*t);
   L0 += 0.00000006124*cos(1.07494838623 + 19804.8272915828*t);
   L0 += 0.00000004488*cos(3.65285033073 + 206.1855484372*t);
   L0 += 0.00000004020*cos(0.83995823171 + 20.3553193988*t);
   L0 += 0.00000005188*cos(4.06503864016 + 6208.2942514241*t);
   L0 += 0.00000005307*cos(0.38216728132 + 31441.6775697568*t);
   L0 += 0.00000003785*cos(2.34369213733 + 3.881335358*t);
   L0 += 0.00000004497*cos(3.27230792447 + 11015.1064773348*t);
   L0 += 0.00000004132*cos(0.92129851256 + 3738.761430108*t);
   L0 += 0.00000003521*cos(5.97844803610 + 3894.1818295422*t);
   L0 += 0.00000004215*cos(1.90601721876 + 245.8316462294*t);
   L0 += 0.00000003701*cos(5.03067498875 + 536.8045120954*t);
   L0 += 0.00000003866*cos(1.82632980909 + 11856.2186514245*t);
   L0 += 0.00000003652*cos(1.01840564429 + 16200.7727245012*t);
   L0 += 0.00000003390*cos(0.97784870142 + 8635.9420037632*t);
   L0 += 0.00000003737*cos(2.95378919570 + 3128.3887650958*t);
   L0 += 0.00000003507*cos(3.71291946317 + 6290.1893969922*t);
   L0 += 0.00000003086*cos(3.64646921512 + 10.6366653498*t);
   L0 += 0.00000003397*cos(1.10589356888 + 14712.317116458*t);
   L0 += 0.00000003334*cos(0.83684903082 + 6496.3749454294*t);
   L0 += 0.00000002805*cos(2.58503711584 + 14314.1681130498*t);
   L0 += 0.00000003650*cos(1.08344142571 + 88860.05707098669*t);
   L0 += 0.00000003388*cos(3.20182380957 + 5120.6011455836*t);
   L0 += 0.00000003252*cos(3.47857474229 + 6133.5126528568*t);
   L0 += 0.00000002553*cos(3.94869027260 + 1990.745017041*t);
   L0 += 0.00000003520*cos(2.05559692878 + 244287.60000722769*t);
   L0 += 0.00000002565*cos(1.56072409371 + 23543.23050468179*t);
   L0 += 0.00000002621*cos(3.85639359951 + 266.6070417218*t);
   L0 += 0.00000002954*cos(3.39692614359 + 9225.539273283*t);
   L0 += 0.00000002876*cos(6.02633318445 + 154717.60988768269*t);
   L0 += 0.00000002395*cos(1.16130078696 + 10984.1923516998*t);
   L0 += 0.00000003161*cos(1.32798862116 + 10873.9860304804*t);
   L0 += 0.00000003163*cos(5.08946546862 + 21228.3920235458*t);
   L0 += 0.00000002361*cos(4.27212461943 + 6040.3472460174*t);
   L0 += 0.00000003030*cos(1.80210001168 + 35371.8872659764*t);
   L0 += 0.00000002343*cos(3.57688971514 + 10969.9652576982*t);
   L0 += 0.00000002618*cos(2.57870151918 + 22483.84857449259*t);
   L0 += 0.00000002113*cos(3.71711179417 + 65147.6197681377*t);
   L0 += 0.00000002019*cos(0.81393923319 + 170.6728706192*t);
   L0 += 0.00000002003*cos(0.38091017375 + 6172.869528772*t);
   L0 += 0.00000002506*cos(3.74378169126 + 10575.4066829418*t);
   L0 += 0.00000002381*cos(0.10581361289 + 7.046236698*t);
   L0 += 0.00000001949*cos(4.86892513469 + 36.0278666774*t);
   L0 += 0.00000002074*cos(4.22802468213 + 5650.2921106782*t);
   L0 += 0.00000001924*cos(5.59460549844 + 6282.0955289232*t);
   L0 += 0.00000001949*cos(1.06999605576 + 5230.807466803*t);
   L0 += 0.00000001988*cos(5.19734705445 + 6262.300454499*t);
   L0 += 0.00000001887*cos(3.74365662683 + 23.8784377478*t);
   L0 += 0.00000001787*cos(1.25929659066 + 12559.038152982*t);
   L0 += 0.00000001883*cos(1.90364058477 + 15.252471185*t);
   L0 += 0.00000001816*cos(3.68083794819 + 15110.4661198662*t);
   L0 += 0.00000001701*cos(4.41109562589 + 110.2063212194*t);
   L0 += 0.00000001990*cos(3.93295788548 + 6206.8097787158*t);
   L0 += 0.00000002103*cos(0.75354936681 + 13521.7514415914*t);
   L0 += 0.00000001774*cos(0.48750515837 + 1551.045222648*t);
   L0 += 0.00000001882*cos(0.86685462305 + 22003.9146348698*t);
   L0 += 0.00000001924*cos(1.22901099088 + 709.9330485583*t);
   L0 += 0.00000002073*cos(4.62531597856 + 6037.244203762*t);
   L0 += 0.00000001924*cos(0.60231842492 + 6284.0561710596*t);
   L0 += 0.00000001596*cos(3.98332879712 + 13916.0191096416*t);
   L0 += 0.00000001664*cos(4.41947015623 + 8662.240323563*t);
   L0 += 0.00000001971*cos(1.04560686192 + 18209.33026366019*t);
   L0 += 0.00000001942*cos(4.31335979989 + 6244.9428143536*t);
   L0 += 0.00000001476*cos(0.93274523818 + 2379.1644735716*t);
   L0 += 0.00000001810*cos(0.49112137707 + 1.4844727083*t);
   L0 += 0.00000001346*cos(1.51574753411 + 4136.9104335162*t);
   L0 += 0.00000001528*cos(5.61833568587 + 6127.6554505572*t);
   L0 += 0.00000001791*cos(3.22191142746 + 39302.096962196*t);
   L0 += 0.00000001747*cos(3.05595292589 + 18319.5365848796*t);
   L0 += 0.00000001432*cos(4.51123984264 + 20426.571092422*t);
   L0 += 0.00000001695*cos(0.22049418623 + 25158.6017197654*t);
   L0 += 0.00000001242*cos(4.46665354536 + 17256.6315363414*t);
   L0 += 0.00000001463*cos(4.69248613506 + 14945.3161735544*t);
   L0 += 0.00000001205*cos(1.86911906771 + 4590.910180489*t);
   L0 += 0.00000001190*cos(2.74169967367 + 12569.6748183318*t);
   L0 += 0.00000001222*cos(5.18120087481 + 5333.9002410216*t);
   L0 += 0.00000001390*cos(5.42888623322 + 143571.32428481648*t);
   L0 += 0.00000001473*cos(1.70487100866 + 11712.9553182308*t);
   L0 += 0.00000001362*cos(2.61069503292 + 6062.6632075526*t);
   L0 += 0.00000001148*cos(6.03008430610 + 3634.6210245184*t);
   L0 += 0.00000001198*cos(5.15296117339 + 10177.2576795336*t);
   L0 += 0.00000001266*cos(0.11422490557 + 18422.62935909819*t);
   L0 += 0.00000001411*cos(1.09910890045 + 3496.032826134*t);
   L0 += 0.00000001349*cos(2.99804623019 + 17654.7805397496*t);
   L0 += 0.00000001253*cos(2.79844902576 + 167283.76158766549*t);
   L0 += 0.00000001311*cos(1.60941074300 + 5481.2549188676*t);
   L0 += 0.00000001079*cos(6.20304501787 + 3.2863574178*t);
   L0 += 0.00000001181*cos(1.20653777627 + 131.5419616864*t);
   L0 += 0.00000001254*cos(5.45103277799 + 6076.8903015542*t);
   L0 += 0.00000001036*cos(2.32136959491 + 7342.4577801806*t);
   L0 += 0.00000001117*cos(0.38842340979 + 949.1756089698*t);
   L0 += 0.00000000966*cos(3.18352079941 + 11087.2851259184*t);
   L0 += 0.00000001171*cos(3.39635167732 + 12562.6285816338*t);
   L0 += 0.00000001121*cos(0.72631814699 + 220.4126424388*t);
   L0 += 0.00000001024*cos(2.19381113265 + 11403.676995575*t);
   L0 += 0.00000000888*cos(3.91167196431 + 4686.8894077068*t);
   L0 += 0.00000000910*cos(1.98802695087 + 735.8765135318*t);
   L0 += 0.00000000823*cos(0.48822202854 + 24072.9214697764*t);
   L0 += 0.00000001096*cos(6.17377835617 + 5436.9930152402*t);
   L0 += 0.00000000908*cos(0.44959148878 + 7477.522860216*t);
   L0 += 0.00000000974*cos(1.52996313552 + 9623.6882766912*t);
   L0 += 0.00000000840*cos(1.79540573407 + 5429.8794682394*t);
   L0 += 0.00000000778*cos(6.17703744517 + 38.1330356378*t);
   L0 += 0.00000000776*cos(4.09859968447 + 14.2270940016*t);
   L0 += 0.00000001068*cos(4.64209577648 + 43232.3066584156*t);
   L0 += 0.00000000954*cos(1.49985885818 + 1162.4747044078*t);
   L0 += 0.00000000907*cos(0.86986870809 + 10344.2950653858*t);
   L0 += 0.00000000931*cos(4.06049877517 + 28766.924424484*t);
   L0 += 0.00000000739*cos(5.04368192034 + 639.897286314*t);
   L0 += 0.00000000965*cos(3.44286716197 + 1589.0728952838*t);
   L0 += 0.00000000763*cos(5.86304776787 + 16858.4825329332*t);
   L0 += 0.00000000953*cos(4.20801492835 + 11190.377900137*t);
   L0 += 0.00000000708*cos(1.72432323967 + 13095.8426650774*t);
   L0 += 0.00000000969*cos(1.64437243011 + 29088.811415985*t);
   L0 += 0.00000000717*cos(0.16688678895 + 11.729352836*t);
   L0 += 0.00000000962*cos(3.53101876172 + 12416.5885028482*t);
   L0 += 0.00000000745*cos(5.77741082302 + 12592.4500197826*t);
   L0 += 0.00000000672*cos(1.91091228744 + 3.9321532631*t);
   L0 += 0.00000000671*cos(5.46240758839 + 18052.9295431578*t);
   L0 += 0.00000000675*cos(6.28311649798 + 4535.0594369244*t);
   L0 += 0.00000000684*cos(0.39975011401 + 5849.3641121146*t);
   L0 += 0.00000000799*cos(0.29859056777 + 12132.439962106*t);
   L0 += 0.00000000758*cos(0.96370719224 + 1052.2683831884*t);
   L0 += 0.00000000782*cos(5.33875702541 + 13517.8701062334*t);
   L0 += 0.00000000730*cos(1.70114998543 + 17267.26820169119*t);
   L0 += 0.00000000749*cos(2.59607005624 + 11609.8625440122*t);
   L0 += 0.00000000734*cos(2.78420497780 + 640.8776073822*t);
   L0 += 0.00000000688*cos(5.15097673557 + 16496.3613962024*t);
   L0 += 0.00000000770*cos(1.62459252416 + 4701.1165017084*t);
   L0 += 0.00000000633*cos(2.20588443066 + 25934.1243310894*t);
   L0 += 0.00000000760*cos(4.21317219403 + 377.3736079158*t);
   L0 += 0.00000000584*cos(2.13420121623 + 10557.5941608238*t);
   L0 += 0.00000000572*cos(0.24649745829 + 9779.1086761254*t);
   L0 += 0.00000000573*cos(3.16435264609 + 533.2140834436*t);
   L0 += 0.00000000685*cos(3.19344289472 + 12146.6670561076*t);
   L0 += 0.00000000675*cos(0.96179234176 + 10454.5013866052*t);
   L0 += 0.00000000648*cos(1.46327342554 + 6268.8487559898*t);
   L0 += 0.00000000589*cos(2.50543543638 + 3097.88382272579*t);
   L0 += 0.00000000551*cos(5.28106257475 + 9388.0059094152*t);
   L0 += 0.00000000696*cos(3.65342155550 + 4804.209275927*t);
   L0 += 0.00000000669*cos(2.51030059260 + 2388.8940204492*t);
   L0 += 0.00000000550*cos(0.06883090057 + 20199.094959633*t);
   L0 += 0.00000000629*cos(4.13350997495 + 45892.73043315699*t);
   L0 += 0.00000000678*cos(6.09190163533 + 135.62532501*t);
   L0 += 0.00000000593*cos(1.50129499103 + 226858.23855437008*t);
   L0 += 0.00000000542*cos(3.58582033525 + 6148.010769956*t);
   L0 += 0.00000000599*cos(6.12058050643 + 18875.525869774*t);
   L0 += 0.00000000682*cos(5.02213616830 + 17253.04110768959*t);
   L0 += 0.00000000565*cos(4.29309214275 + 11933.3679606696*t);
   L0 += 0.00000000486*cos(0.77746204893 + 27.4015560968*t);
   L0 += 0.00000000503*cos(0.58974557727 + 15671.0817594066*t);
   L0 += 0.00000000616*cos(4.06539884128 + 227.476132789*t);
   L0 += 0.00000000537*cos(2.15064382406 + 21954.15760939799*t);
   L0 += 0.00000000669*cos(6.06995500278 + 47162.5163546352*t);
   L0 += 0.00000000540*cos(2.83444221432 + 5326.7866940208*t);
   L0 += 0.00000000474*cos(0.40346826846 + 6915.8595893046*t);
   L0 += 0.00000000532*cos(5.26131065063 + 10988.808157535*t);
   L0 += 0.00000000582*cos(3.24533095664 + 153.7788104848*t);
   L0 += 0.00000000641*cos(3.24711790399 + 2107.0345075424*t);
   L0 += 0.00000000619*cos(3.08302108547 + 33019.0211122046*t);
   L0 += 0.00000000466*cos(3.14982369789 + 10440.2742926036*t);
   L0 += 0.00000000466*cos(0.90708835651 + 5966.6839803348*t);
   L0 += 0.00000000528*cos(0.81926454470 + 813.5502839598*t);
   L0 += 0.00000000603*cos(3.81378921927 + 316428.22867391503*t);
   L0 += 0.00000000559*cos(1.81894985730 + 17996.0311682222*t);
   L0 += 0.00000000437*cos(2.28631745987 + 6303.8512454838*t);
   L0 += 0.00000000518*cos(4.86068318058 + 20597.2439630412*t);
   L0 += 0.00000000424*cos(6.23520018697 + 6489.2613984286*t);
   L0 += 0.00000000518*cos(6.17617826756 + 0.2438174835*t);
   L0 += 0.00000000404*cos(5.72804304258 + 5642.1982426092*t);
   L0 += 0.00000000458*cos(1.34117773914 + 6287.0080032545*t);
   L0 += 0.00000000548*cos(5.68454458320 + 155427.54293624099*t);
   L0 += 0.00000000547*cos(1.03391472434 + 3646.3503773544*t);
   L0 += 0.00000000428*cos(4.69809687820 + 846.0828347512*t);
   L0 += 0.00000000413*cos(6.02520699390 + 6279.4854213396*t);
   L0 += 0.00000000534*cos(3.03030638223 + 66567.48586525429*t);
   L0 += 0.00000000383*cos(1.49055588040 + 19800.9459562248*t);
   L0 += 0.00000000411*cos(5.28384176408 + 18451.07854656599*t);
   L0 += 0.00000000352*cos(4.68891600525 + 4907.3020501456*t);
   L0 += 0.00000000480*cos(5.36572651091 + 348.924420448*t);
   L0 += 0.00000000344*cos(5.89157452889 + 6546.1597733642*t);
   L0 += 0.00000000340*cos(0.37557440365 + 13119.72110282519*t);
   L0 += 0.00000000434*cos(4.98417856239 + 6702.5604938666*t);
   L0 += 0.00000000332*cos(2.68909344430 + 29296.6153895786*t);
   L0 += 0.00000000448*cos(2.16478480251 + 5905.7022420756*t);
   L0 += 0.00000000344*cos(2.06546633735 + 49.7570254718*t);
   L0 += 0.00000000315*cos(1.24023810969 + 4061.2192153944*t);
   L0 += 0.00000000324*cos(2.30897526929 + 5017.508371365*t);
   L0 += 0.00000000413*cos(0.17171692945 + 6286.6662786432*t);
   L0 += 0.00000000431*cos(3.86601101380 + 12489.8856287072*t);
   L0 += 0.00000000349*cos(4.55372493131 + 4933.2084403326*t);
   L0 += 0.00000000323*cos(0.41971136084 + 10770.8932562618*t);
   L0 += 0.00000000341*cos(2.68612860807 + 11.0457002639*t);
   L0 += 0.00000000316*cos(3.52966641606 + 17782.7320727842*t);
   L0 += 0.00000000315*cos(5.63357264999 + 568.8218740274*t);
   L0 += 0.00000000340*cos(3.83571212349 + 10660.6869350424*t);
   L0 += 0.00000000296*cos(0.62703270489 + 20995.3929664494*t);
   L0 += 0.00000000405*cos(1.00084965393 + 16460.33352952499*t);
   L0 += 0.00000000414*cos(1.21998752076 + 51092.7260508548*t);
   L0 += 0.00000000336*cos(4.71465945215 + 6179.9830757728*t);
   L0 += 0.00000000361*cos(3.71235613733 + 28237.2334593894*t);
   L0 += 0.00000000327*cos(1.05606504715 + 11919.140866668*t);
   L0 += 0.00000000327*cos(6.14222420980 + 6254.6266625236*t);
   L0 += 0.00000000268*cos(2.47224339737 + 664.75604513*t);
   L0 += 0.00000000269*cos(1.86210872453 + 23141.5583829246*t);
   L0 += 0.00000000345*cos(0.93461290184 + 6058.7310542895*t);
   L0 += 0.00000000353*cos(4.50033650657 + 36949.2308084242*t);
   L0 += 0.00000000344*cos(6.26166140367 + 24356.7807886416*t);
   L0 += 0.00000000300*cos(4.46964001975 + 6418.1409300268*t);
   L0 += 0.00000000260*cos(4.04967464725 + 6525.8044539654*t);
   L0 += 0.00000000298*cos(2.20018811054 + 156137.47598479928*t);
   L0 += 0.00000000253*cos(3.49930797865 + 29864.334027309*t);
   L0 += 0.00000000254*cos(2.44883530154 + 5331.3574437408*t);
   L0 += 0.00000000296*cos(0.84341183907 + 5729.506447149*t);
   L0 += 0.00000000241*cos(2.00721298729 + 16737.5772365966*t);
   L0 += 0.00000000311*cos(1.23668016336 + 6281.5913772831*t);
   L0 += 0.00000000332*cos(3.55576945724 + 7668.6374249425*t);
   L0 += 0.00000000236*cos(2.47437156031 + 6245.0481773556*t);
   L0 += 0.00000000264*cos(4.43924412283 + 12964.300703391*t);
   L0 += 0.00000000257*cos(1.79654471948 + 11080.1715789176*t);
   L0 += 0.00000000260*cos(3.33077596060 + 5888.4499649322*t);
   L0 += 0.00000000285*cos(0.30886361430 + 11823.1616394502*t);
   L0 += 0.00000000290*cos(5.70141882483 + 77.673770428*t);
   L0 += 0.00000000255*cos(4.00939662024 + 5881.4037282342*t);
   L0 += 0.00000000253*cos(4.73318512715 + 16723.350142595*t);
   L0 += 0.00000000228*cos(0.95333661324 + 5540.0857894588*t);
   L0 += 0.00000000281*cos(1.29199646396 + 22805.7355659936*t);
   L0 += 0.00000000319*cos(1.38633229189 + 163096.18036118349*t);
   L0 += 0.00000000224*cos(1.65156322696 + 10027.9031957292*t);
   L0 += 0.00000000226*cos(0.34125379653 + 17796.9591667858*t);
   L0 += 0.00000000236*cos(4.19817431922 + 19.66976089979*t);
   L0 += 0.00000000281*cos(4.14114899916 + 12539.853380183*t);
   L0 += 0.00000000275*cos(5.50306930248 + 32.5325507914*t);
   L0 += 0.00000000223*cos(5.23334210294 + 56.8983749356*t);
   L0 += 0.00000000217*cos(6.08598789777 + 6805.6532680852*t);
   L0 += 0.00000000228*cos(5.17114391778 + 11720.0688652316*t);
   L0 += 0.00000000274*cos(4.50716805713 + 6016.4688082696*t);
   L0 += 0.00000000245*cos(3.96486270306 + 22.7752014508*t);
   L0 += 0.00000000220*cos(4.72078081970 + 6.62855890001*t);
   L0 += 0.00000000207*cos(5.71701403951 + 41.5507909848*t);
   L0 += 0.00000000204*cos(3.91227411250 + 2699.7348193176*t);
   L0 += 0.00000000209*cos(0.86881969024 + 6321.1035226272*t);
   L0 += 0.00000000200*cos(2.11984442601 + 4274.5183108324*t);
   L0 += 0.00000000200*cos(5.39839888151 + 6019.9919266186*t);
   L0 += 0.00000000209*cos(5.67606291663 + 11293.4706743556*t);
   L0 += 0.00000000252*cos(1.64965729350 + 9380.9596727172*t);
   L0 += 0.00000000275*cos(5.04826903506 + 73.297125859*t);
   L0 += 0.00000000208*cos(1.88207277133 + 11300.5842213564*t);
   L0 += 0.00000000206*cos(5.07284889330 + 6277.552925684*t);
   L0 += 0.00000000272*cos(0.74640924904 + 1975.492545856*t);
   L0 += 0.00000000199*cos(3.30813142103 + 22743.4093795164*t);
   L0 += 0.00000000269*cos(4.48560812155 + 64471.99124174489*t);
   L0 += 0.00000000192*cos(2.17463565107 + 5863.5912061162*t);
   L0 += 0.00000000228*cos(5.85373115869 + 128.0188433374*t);
   L0 += 0.00000000261*cos(2.64321183295 + 55022.9357470744*t);
   L0 += 0.00000000196*cos(2.48537622320 + 16062.1845261168*t);
   L0 += 0.00000000187*cos(1.31891760280 + 29826.3063546732*t);
   L0 += 0.00000000220*cos(5.75012110079 + 29.429508536*t);
   L0 += 0.00000000187*cos(4.03230554718 + 467.9649903544*t);
   L0 += 0.00000000200*cos(5.60555262896 + 1066.49547719*t);
   L0 += 0.00000000231*cos(1.09802712785 + 12341.8069042809*t);
   L0 += 0.00000000198*cos(0.29474229005 + 149.5631971346*t);
   L0 += 0.00000000249*cos(5.10477023170 + 7875.6718636242*t);
   L0 += 0.00000000179*cos(0.87066197995 + 12721.572099417*t);
   L0 += 0.00000000203*cos(1.56914310573 + 28286.9904848612*t);
   L0 += 0.00000000198*cos(3.54061588502 + 30.914125635*t);
   L0 += 0.00000000171*cos(3.45366018621 + 5327.4761083828*t);
   L0 += 0.00000000183*cos(0.72325421604 + 6272.0301497275*t);
   L0 += 0.00000000216*cos(2.97175184412 + 19402.7969528166*t);
   L0 += 0.00000000168*cos(2.51559879907 + 23937.856389741*t);
   L0 += 0.00000000195*cos(0.09045393425 + 156.4007205024*t);
   L0 += 0.00000000216*cos(0.42162375972 + 23539.7073863328*t);
   L0 += 0.00000000189*cos(0.37542530191 + 9814.6041002912*t);
   L0 += 0.00000000207*cos(2.01752547259 + 238004.52415723629*t);
   L0 += 0.00000000218*cos(2.36835893645 + 16627.3709153772*t);
   L0 += 0.00000000166*cos(4.23182960518 + 16840.67001081519*t);
   L0 += 0.00000000200*cos(2.02153258098 + 16097.6799502826*t);
   L0 += 0.00000000169*cos(0.91318727000 + 95.9792272178*t);
   L0 += 0.00000000211*cos(5.73370637657 + 151.8972810852*t);
   L0 += 0.00000000204*cos(0.42643085174 + 515.463871093*t);
   L0 += 0.00000000212*cos(3.00223140894 + 12043.574281889*t);
   L0 += 0.00000000192*cos(5.46153589821 + 6379.0550772092*t);
   L0 += 0.00000000160*cos(6.23798383332 + 202.2533951741*t);
   L0 += 0.00000000215*cos(0.20889073407 + 5621.8429232104*t);
   L0 += 0.00000000150*cos(3.12999753018 + 799.8211251654*t);
   L0 += 0.00000000187*cos(2.12345787867 + 491.6632924588*t);
   L0 += 0.00000000192*cos(1.33928820063 + 394.6258850592*t);
   L0 += 0.00000000149*cos(2.65697593276 + 21.335640467*t);
   L0 += 0.00000000146*cos(5.58021191726 + 412.3710968744*t);
   L0 += 0.00000000156*cos(3.75657676380 + 12323.4230960088*t);
   L0 += 0.00000000143*cos(3.28248547724 + 29.8214381488*t);
   L0 += 0.00000000144*cos(1.07862546598 + 1265.5674786264*t);
   L0 += 0.00000000148*cos(0.23386631090 + 10021.8372800994*t);
   L0 += 0.00000000164*cos(0.94288727597 + 14919.0178537546*t);
   L0 += 0.00000000193*cos(5.92751083827 + 40879.4405046438*t);
   L0 += 0.00000000140*cos(4.97612440269 + 158.9435177832*t);
   L0 += 0.00000000148*cos(2.61651818006 + 17157.0618804718*t);
   L0 += 0.00000000140*cos(3.66947933935 + 26084.0218062162*t);
   L0 += 0.00000000147*cos(5.09968173403 + 661.232926781*t);
   L0 += 0.00000000147*cos(1.36976712162 + 4171.4255366138*t);
   L0 += 0.00000000134*cos(4.79432636012 + 111.1866422876*t);
   L0 += 0.00000000140*cos(1.27748013377 + 107.6635239386*t);
   L0 += 0.00000000171*cos(2.77586207403 + 26735.9452622132*t);
   L0 += 0.00000000183*cos(5.43418358741 + 369.6998159404*t);
   L0 += 0.00000000134*cos(3.09132862833 + 17.812522118*t);
   L0 += 0.00000000132*cos(3.05633896779 + 22490.9621214934*t);
   L0 += 0.00000000181*cos(4.22950689891 + 966.9708774356*t);
   L0 += 0.00000000166*cos(3.67660435776 + 15508.6151232744*t);
   L0 += 0.00000000152*cos(5.28885813387 + 12669.2444742014*t);
   L0 += 0.00000000150*cos(5.86819430895 + 97238.62754448749*t);
   L0 += 0.00000000145*cos(5.07330784304 + 87.30820453981*t);
   L0 += 0.00000000133*cos(5.65471067133 + 31.9723058168*t);
   L0 += 0.00000000124*cos(2.83326216907 + 12566.2190102856*t);
   L0 += 0.00000000140*cos(5.84212721453 + 22476.73502749179*t);
   L0 += 0.00000000134*cos(3.12858101887 + 32217.2001810808*t);
   L0 += 0.00000000137*cos(0.86487461904 + 9924.8104215106*t);
   L0 += 0.00000000172*cos(1.98369595114 + 174242.4659640497*t);
   L0 += 0.00000000170*cos(4.41115280254 + 327574.51427678125*t);
   L0 += 0.00000000128*cos(4.49087631612 + 31415.379249957*t);
   L0 += 0.00000000151*cos(0.46542092001 + 39609.6545831656*t);
   L0 += 0.00000000153*cos(3.78801830344 + 17363.24742890899*t);
   L0 += 0.00000000165*cos(5.31654110459 + 16943.7627850338*t);
   L0 += 0.00000000165*cos(4.06747587817 + 58953.145443294*t);
   L0 += 0.00000000118*cos(0.63846333239 + 6.0659156298*t);
   L0 += 0.00000000159*cos(0.86086959274 + 221995.02880149524*t);
   L0 += 0.00000000119*cos(5.96432932413 + 1385.8952763362*t);
   L0 += 0.00000000114*cos(5.16516114595 + 25685.872802808*t);
   L0 += 0.00000000112*cos(4.92889233335 + 56.8032621698*t);
   L0 += 0.00000000119*cos(2.40626699328 + 18635.9284545362*t);
   L0 += 0.00000000115*cos(0.23374479051 + 418.9243989006*t);
   L0 += 0.00000000113*cos(2.79387576740 + 6272.4391846416*t);
   L0 += 0.00000000122*cos(0.93694072400 + 24492.40611365159*t);
   L0 += 0.00000000133*cos(4.87155573413 + 22345.2603761082*t);
   L0 += 0.00000000113*cos(3.80362889046 + 6293.7125153412*t);
   L0 += 0.00000000130*cos(3.72996018683 + 12573.2652469836*t);
   L0 += 0.00000000107*cos(3.40227152756 + 21393.5419698576*t);
   L0 += 0.00000000122*cos(1.00385670948 + 95143.1329209781*t);
   L0 += 0.00000000140*cos(1.09413075700 + 44809.6502008634*t);
   L0 += 0.00000000112*cos(6.05462382871 + 433.7117378768*t);
   L0 += 0.00000000123*cos(4.55640196386 + 239424.39025435288*t);
   L0 += 0.00000000104*cos(1.54931540602 + 127.9515330346*t);
   L0 += 0.00000000111*cos(3.04186517428 + 8982.810669309*t);
   L0 += 0.00000000102*cos(4.12448497391 + 15664.03552270859*t);
   L0 += 0.00000000107*cos(4.67919356465 + 77690.75950573849*t);
   L0 += 0.00000000118*cos(4.52320170120 + 19004.6479494084*t);
   L0 += 0.00000000107*cos(5.71774478555 + 77736.78343050249*t);
   L0 += 0.00000000103*cos(4.79332126649 + 33794.5437235286*t);
   L0 += 0.00000000143*cos(1.81201813018 + 4214.0690150848*t);
   L0 += 0.00000000102*cos(3.75817786570 + 58864.5439181463*t);
   L0 += 0.00000000125*cos(1.14419195615 + 625.6701923124*t);
   L0 += 0.00000000124*cos(3.27736513892 + 12566.08438968*t);
   L0 += 0.00000000110*cos(1.08705709966 + 2787.0430238574*t);
   L0 += 0.00000000102*cos(4.75119578149 + 12242.6462833254*t);
   L0 += 0.00000000101*cos(4.91289409429 + 401.6721217572*t);
   L0 += 0.00000000138*cos(2.89578979744 + 9411.4646150872*t);
   L0 += 0.00000000129*cos(1.23516042371 + 12029.3471878874*t);
   L0 += 0.00000000138*cos(2.45654707076 + 7576.560073574*t);
   L0 += 0.00000000098*cos(5.44771765020 + 29026.48522950779*t);
   L0 += 0.00000000134*cos(1.43105174912 + 86464.61331683119*t);
   L0 += 0.00000000108*cos(0.98989774940 + 5636.0650166766*t);
   L0 += 0.00000000117*cos(5.17362847134 + 34520.3093093808*t);
   L0 += 0.00000000097*cos(3.34717130592 + 16310.9790457206*t);
   L0 += 0.00000000107*cos(2.94547931851 + 24602.61243487099*t);
   L0 += 0.00000000098*cos(4.37041908717 + 34513.2630726828*t);
   L0 += 0.00000000125*cos(2.72182830814 + 24065.80792277559*t);
   L0 += 0.00000000111*cos(0.58899131543 + 6303.4311693902*t);
   L0 += 0.00000000102*cos(0.66938025772 + 10239.5838660108*t);
   L0 += 0.00000000119*cos(1.21666517886 + 1478.8665740644*t);
   L0 += 0.00000000111*cos(1.04321934681 + 16522.6597160022*t);
   L0 += 0.00000000094*cos(4.31076339857 + 26880.3198130326*t);
   L0 += 0.00000000098*cos(4.14248433763 + 6599.467719648*t);
   L0 += 0.00000000095*cos(2.89807657534 + 34911.412076091*t);
   L0 += 0.00000000097*cos(0.89642320201 + 71980.63357473118*t);
   L0 += 0.00000000095*cos(0.65717727948 + 6288.5987742988*t);
   L0 += 0.00000000116*cos(4.19967201116 + 206.7007372966*t);
   L0 += 0.00000000095*cos(1.78315464297 + 18139.2945014159*t);
   L0 += 0.00000000099*cos(1.37437847718 + 1039.0266107904*t);
   L0 += 0.00000000126*cos(3.21642544972 + 305281.94307104882*t);
   L0 += 0.00000000094*cos(0.68995031860 + 7834.1210726394*t);
   L0 += 0.00000000095*cos(5.58111421744 + 3104.9300594238*t);
   L0 += 0.00000000108*cos(0.52696637156 + 276.7457718644*t);
   L0 += 0.00000000124*cos(3.43899862683 + 172146.97134054029*t);
   L0 += 0.00000000087*cos(1.18764938806 + 18842.11400297339*t);
   L0 += 0.00000000087*cos(0.09094166389 + 15141.390794312*t);
   L0 += 0.00000000108*cos(1.03363414379 + 82576.98122099529*t);
   L0 += 0.00000000119*cos(2.86729109648 + 90394.82301305079*t);
   L0 += 0.00000000104*cos(3.39218586218 + 290.972865866*t);
   L0 += 0.00000000094*cos(5.68284937444 + 32367.0976562076*t);
   L0 += 0.00000000117*cos(0.78475956902 + 83286.91426955358*t);
   L0 += 0.00000000105*cos(3.96551057233 + 6357.8574485587*t);
   L0 += 0.00000000094*cos(4.03443174853 + 13341.6743113068*t);
   L0 += 0.00000000096*cos(0.92742567538 + 1062.9050485382*t);
   L0 += 0.00000000089*cos(4.45371820659 + 792.7748884674*t);
   L0 += 0.00000000087*cos(0.40013481685 + 90279.92316810328*t);
   L0 += 0.00000000113*cos(2.48165313368 + 48739.859897083*t);
   L0 += 0.00000000087*cos(3.43122851097 + 27707.5424942948*t);
   L0 += 0.00000000101*cos(5.32081603011 + 2301.58581590939*t);
   L0 += 0.00000000084*cos(0.67020912458 + 28628.3362260996*t);
   L0 += 0.00000000082*cos(0.87060089842 + 10241.2022911672*t);
   L0 += 0.00000000086*cos(4.62064792290 + 36147.4098773004*t);
   L0 += 0.00000000080*cos(4.77649625396 + 6819.8803620868*t);
   L0 += 0.00000000095*cos(2.87032913492 + 23020.65308658799*t);
   L0 += 0.00000000107*cos(5.77864921649 + 34115.1140692746*t);
   L0 += 0.00000000110*cos(3.32898859416 + 72140.62866668739*t);
   L0 += 0.00000000087*cos(4.40657711727 + 142.1786270362*t);
   L0 += 0.00000000103*cos(4.22506726810 + 30666.1549584328*t);
   L0 += 0.00000000082*cos(3.89404392552 + 5547.1993364596*t);
   L0 += 0.00000000109*cos(1.94546065204 + 24279.10701821359*t);
   L0 += 0.00000000087*cos(4.32472045435 + 742.9900605326*t);
   L0 += 0.00000000107*cos(4.91580912547 + 277.0349937414*t);
   L0 += 0.00000000088*cos(2.10180817713 + 26482.1708096244*t);
   L0 += 0.00000000086*cos(4.01895021483 + 12491.3701014155*t);
   L0 += 0.00000000106*cos(5.49092372854 + 62883.3551395136*t);
   L0 += 0.00000000080*cos(6.19785704950 + 6709.6740408674*t);
   L0 += 0.00000000105*cos(2.44166529175 + 6298.3283211764*t);
   L0 += 0.00000000083*cos(4.90662164029 + 51.28033786241*t);
   L0 += 0.00000000074*cos(2.34622575625 + 7018.9523635232*t);
   L0 += 0.00000000078*cos(6.06947270265 + 148434.53403769129*t);
   L0 += 0.00000000079*cos(3.03048221644 + 838.9692877504*t);
   L0 += 0.00000000073*cos(3.05008665738 + 567.7186377304*t);
   L0 += 0.00000000102*cos(3.59223815483 + 22380.755800274*t);
   L0 += 0.00000000084*cos(0.46604373274 + 45.1412196366*t);
   L0 += 0.00000000092*cos(4.12917744733 + 18216.443810661*t);
   L0 += 0.00000000091*cos(0.49382398887 + 6453.7487206106*t);
   L0 += 0.00000000073*cos(1.72446569088 + 21424.4666443034*t);
   L0 += 0.00000000069*cos(1.49372677730 + 21548.9623692918*t);
   L0 += 0.00000000069*cos(4.73181018058 + 8858.3149443206*t);
   L0 += 0.00000000077*cos(0.47683782532 + 11520.9968637952*t);
   L0 += 0.00000000078*cos(5.43942630990 + 15265.8865193004*t);
   L0 += 0.00000000077*cos(2.53773750372 + 76.2660712756*t);
   L0 += 0.00000000069*cos(2.53932635192 + 9910.583327509*t);
   L0 += 0.00000000075*cos(0.15594729210 + 23006.42599258639*t);
   L0 += 0.00000000078*cos(2.41004950269 + 6393.2821712108*t);
   L0 += 0.00000000067*cos(5.05152846816 + 57375.8019008462*t);
   L0 += 0.00000000077*cos(3.91075604860 + 24383.0791084414*t);
   L0 += 0.00000000067*cos(0.78239147387 + 12779.4507954208*t);
   L0 += 0.00000000071*cos(1.24817900687 + 10419.9862835076*t);
   L0 += 0.00000000067*cos(5.77851227760 + 6311.5250374592*t);
   L0 += 0.00000000061*cos(2.66705754411 + 18852.7506683232*t);
   L0 += 0.00000000079*cos(1.89137330427 + 6528.9074962208*t);
   L0 += 0.00000000084*cos(5.11523704221 + 50317.2034395308*t);
   L0 += 0.00000000059*cos(3.38290426621 + 35707.7100829074*t);
   L0 += 0.00000000065*cos(4.66767908854 + 26709.6469424134*t);
   L0 += 0.00000000079*cos(5.59773848156 + 71960.38658322369*t);
   L0 += 0.00000000061*cos(3.30733768968 + 18845.7044316252*t);
   L0 += 0.00000000057*cos(3.90831298022 + 5999.2165311262*t);
   L0 += 0.00000000058*cos(3.92572820333 + 30348.883772767*t);
   L0 += 0.00000000061*cos(0.05695045822 + 7856.89627409019*t);
   L0 += 0.00000000061*cos(5.63297960691 + 7863.9425107882*t);
   L0 += 0.00000000068*cos(2.53986117507 + 20452.8694122218*t);
   L0 += 0.00000000053*cos(5.44021645443 + 32370.9789915656*t);
   L0 += 0.00000000059*cos(4.95695131145 + 11925.2740926006*t);
   L0 += 0.00000000059*cos(5.34668243273 + 137288.24843482509*t);
   L0 += 0.00000000056*cos(4.05779957425 + 17892.93839400359*t);
   L0 += 0.00000000053*cos(3.94182462468 + 18624.8827542723*t);
   L0 += 0.00000000049*cos(2.11782803206 + 22594.05489571199*t);
   L0 += 0.00000000062*cos(2.47154086715 + 12345.739057544*t);
   L0 += 0.00000000049*cos(5.76596937380 + 18606.4989460002*t);
   L0 += 0.00000000052*cos(6.20902099741 + 21947.11137270*t);
   L0 += 0.00000000051*cos(3.13287981215 + 33326.5787331742*t);
   L0 += 0.00000000057*cos(6.17003357597 + 61306.0115970658*t);
   L0 += 0.00000000064*cos(5.66596451842 + 34596.3646546524*t);
   L0 += 0.00000000066*cos(3.92262335500 + 69853.35207568129*t);
   L0 += 0.00000000053*cos(5.51119362049 + 77710.24834977149*t);
   L0 += 0.00000000053*cos(4.88573986965 + 77717.29458646949*t);
   L0 += 0.00000000048*cos(2.71399112516 + 20760.4270331914*t);
   L0 += 0.00000000046*cos(6.08376164442 + 29822.7832363242*t);
   L0 += 0.00000000052*cos(0.31397481380 + 37724.7534197482*t);
   L0 += 0.00000000045*cos(5.53587248663 + 6262.7205305926*t);
   L0 += 0.00000000059*cos(2.45437896854 + 69166.430989505*t);
   L0 += 0.00000000060*cos(5.24261569842 + 56600.2792895222*t);
   L0 += 0.00000000051*cos(6.15794342172 + 11616.976091013*t);
   L0 += 0.00000000060*cos(1.74758109828 + 44034.1275895394*t);
   L0 += 0.00000000049*cos(5.08973141046 + 33990.6183442862*t);
   L0 += 0.00000000044*cos(3.28270864884 + 29424.634232916*t);
   L0 += 0.00000000052*cos(5.34866947943 + 28313.288804661*t);
   L0 += 0.00000000059*cos(2.12384971916 + 36173.7081971002*t);
   L0 += 0.00000000050*cos(0.46480931695 + 25287.7237993998*t);
   L0 += 0.00000000044*cos(2.05503517210 + 63658.8777508376*t);
   L0 += 0.00000000044*cos(1.71009200258 + 13362.4497067992*t);
   L0 += 0.00000000041*cos(3.91721318316 + 6279.1945146334*t);
   L0 += 0.00000000042*cos(1.52106001448 + 23550.34405168259*t);
   L0 += 0.00000000042*cos(3.10938258068 + 84672.47584450469*t);
   L0 += 0.00000000041*cos(3.78863023321 + 35050.00027447539*t);
   L0 += 0.00000000042*cos(5.57543459263 + 38500.2760310722*t);
   L0 += 0.00000000040*cos(5.55145719363 + 12565.1713789146*t);
   L0 += 0.00000000054*cos(0.88935492100 + 50290.905119731*t);
   L0 += 0.00000000055*cos(2.50268487636 + 77828.671313068*t);
   L0 += 0.00000000054*cos(4.16257918787 + 22910.44676536859*t);
   L0 += 0.00000000041*cos(2.32141215064 + 6286.9571853494*t);
   L0 += 0.00000000045*cos(3.18590576311 + 45585.1728121874*t);
   L0 += 0.00000000045*cos(2.24348941683 + 23646.32327890039*t);
   L0 += 0.00000000043*cos(3.87892110950 + 6549.6828917132*t);
   L0 += 0.00000000048*cos(2.11364139445 + 31570.7996493912*t);
   L0 += 0.00000000043*cos(5.55099558987 + 12359.9661515456*t);
   L0 += 0.00000000041*cos(4.88191569433 + 23536.11695768099*t);
   L0 += 0.00000000049*cos(2.44790922235 + 13613.804277336*t);
   L0 += 0.00000000047*cos(3.58181693910 + 4797.0957289262*t);
   return L0;
}



   double Earth_L1 (double t) // 379 terms of order 1
{
   double L1 = 0;
   L1 += 6283.07584999140;
   L1 += 0.00206058863*cos(2.67823455808 + 6283.0758499914*t);
   L1 += 0.00004303419*cos(2.63512233481 + 12566.1516999828*t);
   L1 += 0.00000425264*cos(1.59046982018 + 3.523118349*t);
   L1 += 0.00000109017*cos(2.96631010675 + 1577.3435424478*t);
   L1 += 0.00000093479*cos(2.59211109542 + 18849.2275499742*t);
   L1 += 0.00000119305*cos(5.79555765566 + 26.2983197998*t);
   L1 += 0.00000072121*cos(1.13840581212 + 529.6909650946*t);
   L1 += 0.00000067784*cos(1.87453300345 + 398.1490034082*t);
   L1 += 0.00000067350*cos(4.40932832004 + 5507.5532386674*t);
   L1 += 0.00000059045*cos(2.88815790631 + 5223.6939198022*t);
   L1 += 0.00000055976*cos(2.17471740035 + 155.4203994342*t);
   L1 += 0.00000045411*cos(0.39799502896 + 796.2980068164*t);
   L1 += 0.00000036298*cos(0.46875437227 + 775.522611324*t);
   L1 += 0.00000028962*cos(2.64732254645 + 7.1135470008*t);
   L1 += 0.00000019097*cos(1.84628376049 + 5486.777843175*t);
   L1 += 0.00000020844*cos(5.34138275149 + 0.9803210682*t);
   L1 += 0.00000018508*cos(4.96855179468 + 213.299095438*t);
   L1 += 0.00000016233*cos(0.03216587315 + 2544.3144198834*t);
   L1 += 0.00000017293*cos(2.99116760630 + 6275.9623029906*t);
   L1 += 0.00000015832*cos(1.43049301283 + 2146.1654164752*t);
   L1 += 0.00000014608*cos(1.20469793690 + 10977.078804699*t);
   L1 += 0.00000011877*cos(3.25805082007 + 5088.6288397668*t);
   L1 += 0.00000011514*cos(2.07502080082 + 4694.0029547076*t);
   L1 += 0.00000009721*cos(4.23925865260 + 1349.8674096588*t);
   L1 += 0.00000009969*cos(1.30263423409 + 6286.5989683404*t);
   L1 += 0.00000009452*cos(2.69956827011 + 242.728603974*t);
   L1 += 0.00000012461*cos(2.83432282119 + 1748.016413067*t);
   L1 += 0.00000011808*cos(5.27379760438 + 1194.4470102246*t);
   L1 += 0.00000008577*cos(5.64476085980 + 951.7184062506*t);
   L1 += 0.00000010641*cos(0.76614722966 + 553.5694028424*t);
   L1 += 0.00000007576*cos(5.30056172859 + 2352.8661537718*t);
   L1 += 0.00000005764*cos(1.77228445837 + 1059.3819301892*t);
   L1 += 0.00000006385*cos(2.65034514038 + 9437.762934887*t);
   L1 += 0.00000005223*cos(5.66135782131 + 71430.69561812909*t);
   L1 += 0.00000005315*cos(0.91110018969 + 3154.6870848956*t);
   L1 += 0.00000006101*cos(4.66633726278 + 4690.4798363586*t);
   L1 += 0.00000004335*cos(0.23934560382 + 6812.766815086*t);
   L1 += 0.00000005041*cos(1.42489704722 + 6438.4962494256*t);
   L1 += 0.00000004259*cos(0.77355543889 + 10447.3878396044*t);
   L1 += 0.00000005200*cos(1.85528830215 + 801.8209311238*t);
   L1 += 0.00000003744*cos(2.00119905572 + 8031.0922630584*t);
   L1 += 0.00000003553*cos(2.42789590229 + 14143.4952424306*t);
   L1 += 0.00000003372*cos(3.86210786421 + 1592.5960136328*t);
   L1 += 0.00000003380*cos(0.88545388924 + 12036.4607348882*t);
   L1 += 0.00000003196*cos(3.19713328141 + 4705.7323075436*t);
   L1 += 0.00000003221*cos(0.61601048990 + 8429.2412664666*t);
   L1 += 0.00000004132*cos(5.23992584671 + 7084.8967811152*t);
   L1 += 0.00000002970*cos(6.07029819073 + 4292.3308329504*t);
   L1 += 0.00000002900*cos(2.32464208411 + 20.3553193988*t);
   L1 += 0.00000003504*cos(4.79976712702 + 6279.5527316424*t);
   L1 += 0.00000002950*cos(1.43108874817 + 5746.271337896*t);
   L1 += 0.00000002697*cos(4.80365209201 + 7234.794256242*t);
   L1 += 0.00000002531*cos(6.22289990904 + 6836.6452528338*t);
   L1 += 0.00000002745*cos(0.93466065396 + 5760.4984318976*t);
   L1 += 0.00000003250*cos(3.39951915286 + 7632.9432596502*t);
   L1 += 0.00000002278*cos(5.00339914806 + 17789.845619785*t);
   L1 += 0.00000002076*cos(3.95551309007 + 10213.285546211*t);
   L1 += 0.00000002061*cos(2.22407719190 + 5856.4776591154*t);
   L1 += 0.00000002252*cos(5.67166717686 + 11499.6562227928*t);
   L1 += 0.00000002148*cos(5.20182663314 + 11513.8833167944*t);
   L1 += 0.00000001886*cos(0.53198539077 + 3340.6124266998*t);
   L1 += 0.00000001875*cos(4.73511969924 + 83996.84731811189*t);
   L1 += 0.00000002060*cos(2.54985307819 + 25132.3033999656*t);
   L1 += 0.00000001794*cos(1.47435300254 + 4164.311989613*t);
   L1 += 0.00000001778*cos(3.02473091781 + 5.5229243074*t);
   L1 += 0.00000002036*cos(0.90908165072 + 6256.7775301916*t);
   L1 += 0.00000002064*cos(2.27051326957 + 522.5774180938*t);
   L1 += 0.00000001773*cos(3.03090500693 + 5753.3848848968*t);
   L1 += 0.00000001569*cos(6.12406216872 + 5216.5803728014*t);
   L1 += 0.00000001590*cos(4.63713748247 + 3.2863574178*t);
   L1 += 0.00000001533*cos(4.20305593883 + 13367.9726311066*t);
   L1 += 0.00000001427*cos(1.19087535126 + 3894.1818295422*t);
   L1 += 0.00000001376*cos(4.24955891338 + 426.598190876*t);
   L1 += 0.00000001375*cos(3.09301252193 + 135.0650800354*t);
   L1 += 0.00000001308*cos(3.08492138890 + 5643.1785636774*t);
   L1 += 0.00000001340*cos(5.76513167968 + 6040.3472460174*t);
   L1 += 0.00000001250*cos(3.07748196332 + 11926.2544136688*t);
   L1 += 0.00000001551*cos(3.07664090662 + 6681.2248533996*t);
   L1 += 0.00000001148*cos(3.24144202282 + 12168.0026965746*t);
   L1 += 0.00000001268*cos(2.09201189992 + 6290.1893969922*t);
   L1 += 0.00000001248*cos(3.44506939791 + 536.8045120954*t);
   L1 += 0.00000001118*cos(2.31830078762 + 16730.4636895958*t);
   L1 += 0.00000001105*cos(5.31966001019 + 23.8784377478*t);
   L1 += 0.00000001012*cos(3.74953487087 + 7860.4193924392*t);
   L1 += 0.00000001025*cos(2.44683771610 + 1990.745017041*t);
   L1 += 0.00000000962*cos(0.81771017882 + 3.881335358*t);
   L1 += 0.00000000911*cos(0.41724352112 + 7079.3738568078*t);
   L1 += 0.00000001091*cos(3.98233608618 + 11506.7697697936*t);
   L1 += 0.00000000957*cos(4.07671436153 + 6127.6554505572*t);
   L1 += 0.00000000835*cos(5.28348689371 + 11790.6290886588*t);
   L1 += 0.00000000802*cos(3.88779080089 + 10973.55568635*t);
   L1 += 0.00000000773*cos(2.41044394817 + 1589.0728952838*t);
   L1 += 0.00000000758*cos(1.30034365873 + 103.0927742186*t);
   L1 += 0.00000000749*cos(4.96281442361 + 6496.3749454294*t);
   L1 += 0.00000000765*cos(3.36312388424 + 36.0278666774*t);
   L1 += 0.00000000915*cos(5.41549763095 + 206.1855484372*t);
   L1 += 0.00000000776*cos(2.57589060224 + 11371.7046897582*t);
   L1 += 0.00000000772*cos(3.98363364977 + 955.5997416086*t);
   L1 += 0.00000000749*cos(5.17900231417 + 10969.9652576982*t);
   L1 += 0.00000000806*cos(0.34226117299 + 9917.6968745098*t);
   L1 += 0.00000000728*cos(5.20962563787 + 38.0276726358*t);
   L1 += 0.00000000685*cos(2.77592961854 + 20.7753954924*t);
   L1 += 0.00000000636*cos(4.28242193632 + 28.4491874678*t);
   L1 += 0.00000000608*cos(5.63278510221 + 10984.1923516998*t);
   L1 += 0.00000000704*cos(5.60739437733 + 3738.761430108*t);
   L1 += 0.00000000685*cos(0.38876148682 + 15.252471185*t);
   L1 += 0.00000000601*cos(0.73489602442 + 419.4846438752*t);
   L1 += 0.00000000716*cos(2.65286869987 + 6309.3741697912*t);
   L1 += 0.00000000584*cos(5.54508741381 + 17298.1823273262*t);
   L1 += 0.00000000628*cos(1.11733054796 + 7058.5984613154*t);
   L1 += 0.00000000688*cos(2.59684132401 + 3496.032826134*t);
   L1 += 0.00000000485*cos(0.44470714066 + 12352.8526045448*t);
   L1 += 0.00000000562*cos(2.82510352358 + 3930.2096962196*t);
   L1 += 0.00000000597*cos(5.27675789973 + 10575.4066829418*t);
   L1 += 0.00000000583*cos(3.18931539860 + 4732.0306273434*t);
   L1 += 0.00000000526*cos(5.01737745304 + 5884.9268465832*t);
   L1 += 0.00000000540*cos(1.29182747488 + 640.8776073822*t);
   L1 += 0.00000000481*cos(5.49721461067 + 5230.807466803*t);
   L1 += 0.00000000406*cos(5.21253018484 + 220.4126424388*t);
   L1 += 0.00000000395*cos(1.87489912123 + 16200.7727245012*t);
   L1 += 0.00000000367*cos(0.88533542945 + 6283.14316029419*t);
   L1 += 0.00000000369*cos(3.84778078192 + 18073.7049386502*t);
   L1 += 0.00000000379*cos(0.37991716505 + 10177.2576795336*t);
   L1 += 0.00000000356*cos(3.84152910109 + 11712.9553182308*t);
   L1 += 0.00000000374*cos(5.01577520608 + 7.046236698*t);
   L1 += 0.00000000381*cos(4.30250406634 + 6062.6632075526*t);
   L1 += 0.00000000471*cos(0.86388942467 + 6069.7767545534*t);
   L1 += 0.00000000367*cos(1.32943839930 + 6283.0085396886*t);
   L1 += 0.00000000460*cos(5.19667219582 + 6284.0561710596*t);
   L1 += 0.00000000333*cos(5.54250425107 + 4686.8894077068*t);
   L1 += 0.00000000341*cos(4.36524495363 + 7238.67559160*t);
   L1 += 0.00000000336*cos(4.00205876835 + 3097.88382272579*t);
   L1 += 0.00000000359*cos(6.22679790284 + 245.8316462294*t);
   L1 += 0.00000000307*cos(2.35299010924 + 170.6728706192*t);
   L1 += 0.00000000343*cos(3.77164927142 + 6076.8903015542*t);
   L1 += 0.00000000296*cos(5.44138799494 + 17260.1546546904*t);
   L1 += 0.00000000328*cos(0.13817705132 + 11015.1064773348*t);
   L1 += 0.00000000269*cos(1.13308244952 + 12569.6748183318*t);
   L1 += 0.00000000263*cos(0.00550737510 + 4136.9104335162*t);
   L1 += 0.00000000282*cos(5.04399588559 + 7477.522860216*t);
   L1 += 0.00000000280*cos(3.13703211405 + 12559.038152982*t);
   L1 += 0.00000000259*cos(0.93882269388 + 5642.1982426092*t);
   L1 += 0.00000000292*cos(1.98426314297 + 12132.439962106*t);
   L1 += 0.00000000247*cos(3.84244798673 + 5429.8794682394*t);
   L1 += 0.00000000319*cos(5.04170148790 + 90617.7374312997*t);
   L1 += 0.00000000245*cos(5.70469737024 + 65147.6197681377*t);
   L1 += 0.00000000318*cos(1.35581968834 + 78051.5857313169*t);
   L1 += 0.00000000241*cos(0.99469787369 + 3634.6210245184*t);
   L1 += 0.00000000246*cos(3.06168069393 + 110.2063212194*t);
   L1 += 0.00000000239*cos(6.11854529589 + 11856.2186514245*t);
   L1 += 0.00000000267*cos(0.65297608414 + 21228.3920235458*t);
   L1 += 0.00000000262*cos(1.51070507866 + 12146.6670561076*t);
   L1 += 0.00000000230*cos(1.75923794017 + 9779.1086761254*t);
   L1 += 0.00000000223*cos(2.00967043606 + 6172.869528772*t);
   L1 += 0.00000000246*cos(1.10411690861 + 6282.0955289232*t);
   L1 += 0.00000000214*cos(4.03840492266 + 14314.1681130498*t);
   L1 += 0.00000000212*cos(2.13695623228 + 5849.3641121146*t);
   L1 += 0.00000000207*cos(3.07724246401 + 11.729352836*t);
   L1 += 0.00000000207*cos(6.10303325026 + 23543.23050468179*t);
   L1 += 0.00000000266*cos(1.00720021877 + 2388.8940204492*t);
   L1 += 0.00000000217*cos(6.27840212312 + 17267.26820169119*t);
   L1 += 0.00000000231*cos(5.37372783468 + 13916.0191096416*t);
   L1 += 0.00000000204*cos(2.34615348695 + 266.6070417218*t);
   L1 += 0.00000000195*cos(5.55015549753 + 6133.5126528568*t);
   L1 += 0.00000000203*cos(4.65616806503 + 24072.9214697764*t);
   L1 += 0.00000000188*cos(2.52682282169 + 6525.8044539654*t);
   L1 += 0.00000000177*cos(1.73426919199 + 154717.60988768269*t);
   L1 += 0.00000000187*cos(4.76501318048 + 4535.0594369244*t);
   L1 += 0.00000000186*cos(4.63080493407 + 10440.2742926036*t);
   L1 += 0.00000000183*cos(3.20060840641 + 8635.9420037632*t);
   L1 += 0.00000000172*cos(1.45551703877 + 9225.539273283*t);
   L1 += 0.00000000162*cos(3.30665137166 + 639.897286314*t);
   L1 += 0.00000000168*cos(2.17671416605 + 27.4015560968*t);
   L1 += 0.00000000160*cos(1.68168712750 + 15110.4661198662*t);
   L1 += 0.00000000194*cos(2.79243768345 + 7342.4577801806*t);
   L1 += 0.00000000183*cos(0.56273524797 + 13517.8701062334*t);
   L1 += 0.00000000172*cos(5.97039514134 + 4701.1165017084*t);
   L1 += 0.00000000179*cos(3.58450811616 + 87.30820453981*t);
   L1 += 0.00000000152*cos(2.84070476839 + 5650.2921106782*t);
   L1 += 0.00000000156*cos(1.07156076421 + 18319.5365848796*t);
   L1 += 0.00000000182*cos(0.44053620124 + 17253.04110768959*t);
   L1 += 0.00000000142*cos(1.46290137520 + 11087.2851259184*t);
   L1 += 0.00000000131*cos(5.40912137746 + 2699.7348193176*t);
   L1 += 0.00000000144*cos(2.07312089638 + 25158.6017197654*t);
   L1 += 0.00000000147*cos(6.15107800602 + 9623.6882766912*t);
   L1 += 0.00000000141*cos(5.55739979498 + 10454.5013866052*t);
   L1 += 0.00000000135*cos(0.06098129430 + 16723.350142595*t);
   L1 += 0.00000000123*cos(5.81194797368 + 17256.6315363414*t);
   L1 += 0.00000000124*cos(2.36269386269 + 4933.2084403326*t);
   L1 += 0.00000000126*cos(3.47483886466 + 22483.84857449259*t);
   L1 += 0.00000000159*cos(5.63944722033 + 5729.506447149*t);
   L1 += 0.00000000137*cos(1.93811728826 + 20426.571092422*t);
   L1 += 0.00000000123*cos(3.92815962825 + 17996.0311682222*t);
   L1 += 0.00000000148*cos(3.02542567608 + 1551.045222648*t);
   L1 += 0.00000000121*cos(0.05537321071 + 13095.8426650774*t);
   L1 += 0.00000000120*cos(5.91904349732 + 6206.8097787158*t);
   L1 += 0.00000000134*cos(3.11122937825 + 21954.15760939799*t);
   L1 += 0.00000000119*cos(5.52143897201 + 709.9330485583*t);
   L1 += 0.00000000122*cos(3.00840036775 + 19800.9459562248*t);
   L1 += 0.00000000127*cos(1.37534182407 + 14945.3161735544*t);
   L1 += 0.00000000141*cos(2.56886299638 + 1052.2683831884*t);
   L1 += 0.00000000123*cos(2.83671175442 + 11919.140866668*t);
   L1 += 0.00000000118*cos(0.81918292547 + 5331.3574437408*t);
   L1 += 0.00000000151*cos(2.68728567951 + 11769.8536931664*t);
   L1 += 0.00000000119*cos(5.08654046247 + 5481.2549188676*t);
   L1 += 0.00000000113*cos(4.42675663942 + 18422.62935909819*t);
   L1 += 0.00000000153*cos(2.46021790779 + 11933.3679606696*t);
   L1 += 0.00000000108*cos(1.04936452151 + 11403.676995575*t);
   L1 += 0.00000000128*cos(0.99810456461 + 8827.3902698748*t);
   L1 += 0.00000000144*cos(2.54869747042 + 227.476132789*t);
   L1 += 0.00000000150*cos(4.50631437136 + 2379.1644735716*t);
   L1 += 0.00000000109*cos(0.29269062317 + 16737.5772365966*t);
   L1 += 0.00000000122*cos(4.23040027813 + 29.429508536*t);
   L1 += 0.00000000111*cos(5.16970710025 + 17782.7320727842*t);
   L1 += 0.00000000105*cos(1.61738153441 + 13119.72110282519*t);
   L1 += 0.00000000100*cos(3.52204690579 + 18052.9295431578*t);
   L1 += 0.00000000108*cos(1.08493117155 + 16858.4825329332*t);
   L1 += 0.00000000135*cos(3.20160616970 + 6262.300454499*t);
   L1 += 0.00000000106*cos(1.96085069786 + 74.7815985673*t);
   L1 += 0.00000000129*cos(4.85949366504 + 16496.3613962024*t);
   L1 += 0.00000000110*cos(2.30605777952 + 16460.33352952499*t);
   L1 += 0.00000000097*cos(3.50918940210 + 5333.9002410216*t);
   L1 += 0.00000000099*cos(3.56417337974 + 735.8765135318*t);
   L1 += 0.00000000096*cos(3.40918487598 + 15720.8387848784*t);
   L1 += 0.00000000094*cos(5.01601027363 + 3128.3887650958*t);
   L1 += 0.00000000097*cos(1.65579893894 + 533.2140834436*t);
   L1 += 0.00000000092*cos(0.89219199493 + 29296.6153895786*t);
   L1 += 0.00000000123*cos(3.16062062663 + 9380.9596727172*t);
   L1 += 0.00000000102*cos(1.20493500565 + 23020.65308658799*t);
   L1 += 0.00000000088*cos(2.21265504437 + 12721.572099417*t);
   L1 += 0.00000000089*cos(1.54322669570 + 20199.094959633*t);
   L1 += 0.00000000121*cos(6.19860353182 + 9388.0059094152*t);
   L1 += 0.00000000089*cos(4.08082274765 + 22805.7355659936*t);
   L1 += 0.00000000098*cos(1.09176668094 + 12043.574281889*t);
   L1 += 0.00000000086*cos(1.13649001466 + 143571.32428481648*t);
   L1 += 0.00000000088*cos(5.96980472191 + 107.6635239386*t);
   L1 += 0.00000000082*cos(5.01561173481 + 22003.9146348698*t);
   L1 += 0.00000000094*cos(1.69615700473 + 23006.42599258639*t);
   L1 += 0.00000000081*cos(3.00664741995 + 2118.7638603784*t);
   L1 += 0.00000000098*cos(1.39215287161 + 8662.240323563*t);
   L1 += 0.00000000080*cos(5.16340988714 + 17796.9591667858*t);
   L1 += 0.00000000082*cos(5.86893959287 + 2787.0430238574*t);
   L1 += 0.00000000077*cos(0.80723694712 + 167283.76158766549*t);
   L1 += 0.00000000091*cos(5.74902425304 + 21424.4666443034*t);
   L1 += 0.00000000076*cos(5.67183650604 + 14.2270940016*t);
   L1 += 0.00000000081*cos(6.16619455699 + 1039.0266107904*t);
   L1 += 0.00000000076*cos(3.21449884756 + 111.1866422876*t);
   L1 += 0.00000000078*cos(1.37531518377 + 21947.11137270*t);
   L1 += 0.00000000074*cos(3.58814195051 + 11609.8625440122*t);
   L1 += 0.00000000077*cos(4.84846488388 + 22743.4093795164*t);
   L1 += 0.00000000090*cos(1.48869013606 + 15671.0817594066*t);
   L1 += 0.00000000082*cos(3.48618398216 + 29088.811415985*t);
   L1 += 0.00000000071*cos(2.21591561190 + 12029.3471878874*t);
   L1 += 0.00000000069*cos(1.93625656075 + 135.62532501*t);
   L1 += 0.00000000070*cos(2.66552760898 + 18875.525869774*t);
   L1 += 0.00000000069*cos(5.41478093731 + 26735.9452622132*t);
   L1 += 0.00000000079*cos(5.15158156951 + 12323.4230960088*t);
   L1 += 0.00000000078*cos(4.17014063638 + 1066.49547719*t);
   L1 += 0.00000000065*cos(5.64584720343 + 12139.5535091068*t);
   L1 += 0.00000000071*cos(3.89804774037 + 22779.4372461938*t);
   L1 += 0.00000000063*cos(4.53968787714 + 8982.810669309*t);
   L1 += 0.00000000076*cos(3.29088891716 + 2942.4634232916*t);
   L1 += 0.00000000069*cos(0.94232113005 + 14919.0178537546*t);
   L1 += 0.00000000063*cos(4.09186535490 + 16062.1845261168*t);
   L1 += 0.00000000065*cos(3.34580407184 + 51.28033786241*t);
   L1 += 0.00000000065*cos(5.75757544877 + 52670.0695933026*t);
   L1 += 0.00000000057*cos(5.25050277534 + 20995.3929664494*t);
   L1 += 0.00000000061*cos(1.92290673861 + 6805.6532680852*t);
   L1 += 0.00000000061*cos(0.08878901558 + 13362.4497067992*t);
   L1 += 0.00000000076*cos(1.86947679415 + 25287.7237993998*t);
   L1 += 0.00000000056*cos(4.25396542622 + 6709.6740408674*t);
   L1 += 0.00000000058*cos(4.79429715781 + 6286.3622074092*t);
   L1 += 0.00000000073*cos(0.53299090807 + 2301.58581590939*t);
   L1 += 0.00000000070*cos(4.31243357502 + 19402.7969528166*t);
   L1 += 0.00000000067*cos(2.53852336668 + 377.3736079158*t);
   L1 += 0.00000000056*cos(3.20816844695 + 24889.5747959916*t);
   L1 += 0.00000000054*cos(5.17336469511 + 26084.0218062162*t);
   L1 += 0.00000000053*cos(3.17675406016 + 18451.07854656599*t);
   L1 += 0.00000000053*cos(3.61529270216 + 77.673770428*t);
   L1 += 0.00000000053*cos(0.45467549335 + 30666.1549584328*t);
   L1 += 0.00000000053*cos(2.97761644192 + 21548.9623692918*t);
   L1 += 0.00000000061*cos(0.14805728543 + 23013.5395395872*t);
   L1 += 0.00000000051*cos(3.32803972907 + 56.8983749356*t);
   L1 += 0.00000000052*cos(3.41304011355 + 23141.5583829246*t);
   L1 += 0.00000000058*cos(3.13638677202 + 309.2783226558*t);
   L1 += 0.00000000054*cos(1.60896548545 + 13341.6743113068*t);
   L1 += 0.00000000053*cos(5.81426394852 + 16193.65917750039*t);
   L1 += 0.00000000067*cos(6.27917920454 + 22345.2603761082*t);
   L1 += 0.00000000050*cos(0.42577644151 + 25685.872802808*t);
   L1 += 0.00000000048*cos(0.70204553352 + 1162.4747044078*t);
   L1 += 0.00000000066*cos(3.64350022359 + 15265.8865193004*t);
   L1 += 0.00000000050*cos(5.74382917440 + 19.66976089979*t);
   L1 += 0.00000000054*cos(1.97277370837 + 23581.2581773176*t);
   L1 += 0.00000000051*cos(1.23713196525 + 12539.853380183*t);
   L1 += 0.00000000046*cos(5.41431704639 + 33019.0211122046*t);
   L1 += 0.00000000046*cos(4.80640843261 + 19651.048481098*t);
   L1 += 0.00000000050*cos(1.23847511223 + 22476.73502749179*t);
   L1 += 0.00000000051*cos(4.91913434178 + 12592.4500197826*t);
   L1 += 0.00000000046*cos(2.41369976086 + 98068.53671630539*t);
   L1 += 0.00000000045*cos(3.45227074337 + 30774.5016425748*t);
   L1 += 0.00000000045*cos(4.39659083856 + 433.7117378768*t);
   L1 += 0.00000000045*cos(3.71921056816 + 18209.33026366019*t);
   L1 += 0.00000000044*cos(2.47683925106 + 24356.7807886416*t);
   L1 += 0.00000000049*cos(2.17835058609 + 13521.7514415914*t);
   L1 += 0.00000000046*cos(0.26142733448 + 11.0457002639*t);
   L1 += 0.00000000045*cos(2.46230645202 + 51868.2486621788*t);
   L1 += 0.00000000043*cos(4.29458463014 + 28230.18722269139*t);
   L1 += 0.00000000048*cos(0.89551707131 + 56600.2792895222*t);
   L1 += 0.00000000042*cos(3.63410684699 + 4590.910180489*t);
   L1 += 0.00000000049*cos(3.17757670967 + 6303.8512454838*t);
   L1 += 0.00000000043*cos(4.93350349236 + 10021.8372800994*t);
   L1 += 0.00000000052*cos(3.65410195699 + 7872.1487452752*t);
   L1 += 0.00000000041*cos(4.82166756935 + 10988.808157535*t);
   L1 += 0.00000000040*cos(1.81891629936 + 34596.3646546524*t);
   L1 += 0.00000000043*cos(1.94164978061 + 1903.4368125012*t);
   L1 += 0.00000000041*cos(0.74461854136 + 23937.856389741*t);
   L1 += 0.00000000045*cos(5.45575017530 + 60530.4889857418*t);
   L1 += 0.00000000050*cos(5.67355640472 + 18216.443810661*t);
   L1 += 0.00000000040*cos(0.04502010161 + 38526.574350872*t);
   L1 += 0.00000000053*cos(3.64807615995 + 11925.2740926006*t);
   L1 += 0.00000000042*cos(5.19292937193 + 19004.6479494084*t);
   L1 += 0.00000000041*cos(0.94309683296 + 9924.8104215106*t);
   L1 += 0.00000000039*cos(4.61184303844 + 95.9792272178*t);
   L1 += 0.00000000049*cos(2.05532526216 + 12573.2652469836*t);
   L1 += 0.00000000045*cos(3.73717824543 + 7875.6718636242*t);
   L1 += 0.00000000043*cos(1.14078465002 + 49.7570254718*t);
   L1 += 0.00000000039*cos(1.70539366023 + 32217.2001810808*t);
   L1 += 0.00000000037*cos(1.29390383811 + 310.8407988684*t);
   L1 += 0.00000000038*cos(0.95970925950 + 664.75604513*t);
   L1 += 0.00000000039*cos(0.85957361635 + 16522.6597160022*t);
   L1 += 0.00000000040*cos(1.00170796001 + 36949.2308084242*t);
   L1 += 0.00000000040*cos(3.78164718776 + 55798.4583583984*t);
   L1 += 0.00000000037*cos(5.42237070904 + 6286.6662786432*t);
   L1 += 0.00000000036*cos(1.68167662194 + 10344.2950653858*t);
   L1 += 0.00000000040*cos(5.13217319067 + 15664.03552270859*t);
   L1 += 0.00000000049*cos(3.62741283878 + 77713.7714681205*t);
   L1 += 0.00000000045*cos(6.15877872538 + 28286.9904848612*t);
   L1 += 0.00000000036*cos(3.32158458257 + 16207.886271502*t);
   L1 += 0.00000000035*cos(5.83917764292 + 6321.1035226272*t);
   L1 += 0.00000000036*cos(1.80784164083 + 6279.7894925736*t);
   L1 += 0.00000000035*cos(4.60279245362 + 28237.2334593894*t);
   L1 += 0.00000000033*cos(0.47301775923 + 18635.9284545362*t);
   L1 += 0.00000000035*cos(4.36571027474 + 48739.859897083*t);
   L1 += 0.00000000031*cos(3.06828028412 + 6819.8803620868*t);
   L1 += 0.00000000031*cos(2.40203198480 + 28628.3362260996*t);
   L1 += 0.00000000034*cos(1.90096411242 + 12964.300703391*t);
   L1 += 0.00000000029*cos(6.09291010354 + 18606.4989460002*t);
   L1 += 0.00000000028*cos(3.42046112698 + 6288.5987742988*t);
   L1 += 0.00000000028*cos(3.43736140600 + 34115.1140692746*t);
   L1 += 0.00000000029*cos(1.48920816078 + 6489.2613984286*t);
   L1 += 0.00000000028*cos(3.07474749886 + 29822.7832363242*t);
   L1 += 0.00000000038*cos(2.44608264663 + 31415.379249957*t);
   L1 += 0.00000000028*cos(2.98392582088 + 6277.552925684*t);
   L1 += 0.00000000027*cos(5.03556015623 + 12779.4507954208*t);
   L1 += 0.00000000027*cos(5.40812977287 + 26087.9031415742*t);
   L1 += 0.00000000038*cos(5.56439937893 + 27832.0382192832*t);
   L1 += 0.00000000026*cos(3.87685883153 + 6262.7205305926*t);
   L1 += 0.00000000027*cos(6.15565392840 + 28759.81087748319*t);
   L1 += 0.00000000025*cos(5.17122153205 + 6915.8595893046*t);
   L1 += 0.00000000027*cos(4.03132006944 + 9910.583327509*t);
   L1 += 0.00000000033*cos(3.97763407373 + 12410.7313005486*t);
   L1 += 0.00000000024*cos(5.31307120044 + 29026.48522950779*t);
   L1 += 0.00000000025*cos(0.81055213297 + 36173.7081971002*t);
   L1 += 0.00000000024*cos(1.28706101310 + 24491.4257925834*t);
   L1 += 0.00000000021*cos(2.02548478742 + 28766.924424484*t);
   L1 += 0.00000000021*cos(6.07545114034 + 18139.2945014159*t);
   L1 += 0.00000000022*cos(2.31199937131 + 6303.4311693902*t);
   L1 += 0.00000000023*cos(1.35235057478 + 49515.382508407*t);
   L1 += 0.00000000023*cos(2.92765926961 + 65236.2212932854*t);
   L1 += 0.00000000019*cos(0.03636659763 + 29864.334027309*t);
   L1 += 0.00000000020*cos(0.91374066194 + 45585.1728121874*t);
   L1 += 0.00000000021*cos(4.45003013294 + 22490.9621214934*t);
   L1 += 0.00000000022*cos(1.97119365888 + 34513.2630726828*t);
   L1 += 0.00000000020*cos(4.11682669951 + 17157.0618804718*t);
   L1 += 0.00000000024*cos(1.40243942415 + 14712.317116458*t);
   L1 += 0.00000000022*cos(2.19759737115 + 31570.7996493912*t);
   L1 += 0.00000000021*cos(1.48739821208 + 61306.0115970658*t);
   L1 += 0.00000000025*cos(5.71465573409 + 25934.1243310894*t);
   return L1*t;
}


   double Earth_L2 (double t) // 144 terms of order 2
{
   double L2 = 0;
   L2 += 0.00008721859*cos(1.07253635559 + 6283.0758499914*t);
   L2 -= 0.00000990990;
   L2 += 0.00000294833*cos(0.43717350256 + 12566.1516999828*t);
   L2 += 0.00000027338*cos(0.05295636147 + 3.523118349*t);
   L2 += 0.00000016333*cos(5.18820215724 + 26.2983197998*t);
   L2 += 0.00000015745*cos(3.68504712183 + 155.4203994342*t);
   L2 += 0.00000009425*cos(0.29667114694 + 18849.2275499742*t);
   L2 += 0.00000008938*cos(2.05706319592 + 77713.7714681205*t);
   L2 += 0.00000006940*cos(0.82691541038 + 775.522611324*t);
   L2 += 0.00000005061*cos(4.66243231680 + 1577.3435424478*t);
   L2 += 0.00000004060*cos(1.03067032318 + 7.1135470008*t);
   L2 += 0.00000003464*cos(5.14021224609 + 796.2980068164*t);
   L2 += 0.00000003172*cos(6.05479318507 + 5507.5532386674*t);
   L2 += 0.00000003020*cos(1.19240008524 + 242.728603974*t);
   L2 += 0.00000002885*cos(6.11705865396 + 529.6909650946*t);
   L2 += 0.00000003809*cos(3.44043369494 + 5573.1428014331*t);
   L2 += 0.00000002719*cos(0.30363248164 + 398.1490034082*t);
   L2 += 0.00000002365*cos(4.37666117992 + 5223.6939198022*t);
   L2 += 0.00000002538*cos(2.27966434314 + 553.5694028424*t);
   L2 += 0.00000002078*cos(3.75435095487 + 0.9803210682*t);
   L2 += 0.00000001675*cos(0.90149951436 + 951.7184062506*t);
   L2 += 0.00000001534*cos(5.75895831192 + 1349.8674096588*t);
   L2 += 0.00000001224*cos(2.97285792195 + 2146.1654164752*t);
   L2 += 0.00000001449*cos(4.36401639552 + 1748.016413067*t);
   L2 += 0.00000001341*cos(3.72019379666 + 1194.4470102246*t);
   L2 += 0.00000001253*cos(2.94888726310 + 6438.4962494256*t);
   L2 += 0.00000000999*cos(5.98665341008 + 6286.5989683404*t);
   L2 += 0.00000000917*cos(4.79722195184 + 5088.6288397668*t);
   L2 += 0.00000000829*cos(3.31021398862 + 213.299095438*t);
   L2 += 0.00000001102*cos(1.27094359244 + 161000.6857376741*t);
   L2 += 0.00000000764*cos(3.41231607038 + 5486.777843175*t);
   L2 += 0.00000001046*cos(0.60374190029 + 3154.6870848956*t);
   L2 += 0.00000000887*cos(5.23364022595 + 7084.8967811152*t);
   L2 += 0.00000000644*cos(1.59974355582 + 2544.3144198834*t);
   L2 += 0.00000000681*cos(3.42742947469 + 4694.0029547076*t);
   L2 += 0.00000000606*cos(2.47688996663 + 10977.078804699*t);
   L2 += 0.00000000706*cos(6.19369692903 + 4690.4798363586*t);
   L2 += 0.00000000643*cos(1.98119869589 + 801.8209311238*t);
   L2 += 0.00000000502*cos(1.44415463861 + 6836.6452528338*t);
   L2 += 0.00000000490*cos(2.33889753806 + 1592.5960136328*t);
   L2 += 0.00000000458*cos(1.30867922972 + 4292.3308329504*t);
   L2 += 0.00000000431*cos(0.03542536476 + 7234.794256242*t);
   L2 += 0.00000000349*cos(0.98779272263 + 6040.3472460174*t);
   L2 += 0.00000000385*cos(1.57065592218 + 71430.69561812909*t);
   L2 += 0.00000000371*cos(3.16149051601 + 6309.3741697912*t);
   L2 += 0.00000000348*cos(0.64980950594 + 1059.3819301892*t);
   L2 += 0.00000000458*cos(3.81505682017 + 149854.40013480789*t);
   L2 += 0.00000000302*cos(1.91723873447 + 10447.3878396044*t);
   L2 += 0.00000000306*cos(3.55405857884 + 8031.0922630584*t);
   L2 += 0.00000000395*cos(4.93742673052 + 7632.9432596502*t);
   L2 += 0.00000000314*cos(3.18058352846 + 2352.8661537718*t);
   L2 += 0.00000000281*cos(4.41751404023 + 9437.762934887*t);
   L2 += 0.00000000276*cos(2.71075791682 + 3894.1818295422*t);
   L2 += 0.00000000298*cos(2.52045757001 + 6127.6554505572*t);
   L2 += 0.00000000272*cos(0.24370191144 + 25132.3033999656*t);
   L2 += 0.00000000251*cos(0.55449375765 + 6279.5527316424*t);
   L2 += 0.00000000224*cos(1.40790249012 + 4705.7323075436*t);
   L2 += 0.00000000258*cos(5.29510765044 + 6812.766815086*t);
   L2 += 0.00000000178*cos(0.92636669742 + 1990.745017041*t);
   L2 += 0.00000000217*cos(0.68543630022 + 6256.7775301916*t);
   L2 += 0.00000000154*cos(0.77808632062 + 14143.4952424306*t);
   L2 += 0.00000000150*cos(2.40636982736 + 426.598190876*t);
   L2 += 0.00000000196*cos(6.06877865012 + 640.8776073822*t);
   L2 += 0.00000000137*cos(2.21947617717 + 8429.2412664666*t);
   L2 += 0.00000000127*cos(5.47380312768 + 12036.4607348882*t);
   L2 += 0.00000000121*cos(3.32740512021 + 17789.845619785*t);
   L2 += 0.00000000123*cos(2.16004509785 + 10213.285546211*t);
   L2 += 0.00000000116*cos(0.49705139709 + 7058.5984613154*t);
   L2 += 0.00000000138*cos(2.36181661472 + 11506.7697697936*t);
   L2 += 0.00000000101*cos(0.86299995919 + 6290.1893969922*t);
   L2 += 0.00000000118*cos(5.82317768355 + 7860.4193924392*t);
   L2 += 0.00000000092*cos(5.11639978593 + 7079.3738568078*t);
   L2 += 0.00000000125*cos(2.65424538513 + 88860.05707098669*t);
   L2 += 0.00000000088*cos(4.31188236160 + 83996.84731811189*t);
   L2 += 0.00000000084*cos(3.57682769713 + 16730.4636895958*t);
   L2 += 0.00000000097*cos(5.58011309774 + 13367.9726311066*t);
   L2 += 0.00000000102*cos(2.05853060226 + 87.30820453981*t);
   L2 += 0.00000000080*cos(4.73827128421 + 11926.2544136688*t);
   L2 += 0.00000000080*cos(5.41344057121 + 10973.55568635*t);
   L2 += 0.00000000106*cos(4.10978997399 + 3496.032826134*t);
   L2 += 0.00000000075*cos(4.89166898876 + 5643.1785636774*t);
   L2 += 0.00000000100*cos(3.62645659087 + 244287.60000722769*t);
   L2 += 0.00000000096*cos(1.39443577787 + 6681.2248533996*t);
   L2 += 0.00000000069*cos(1.88399189965 + 10177.2576795336*t);
   L2 += 0.00000000087*cos(0.40842153208 + 11015.1064773348*t);
   L2 += 0.00000000066*cos(0.99444831932 + 6525.8044539654*t);
   L2 += 0.00000000066*cos(1.42471816453 + 9917.6968745098*t);
   L2 += 0.00000000067*cos(5.51240997070 + 3097.88382272579*t);
   L2 += 0.00000000076*cos(2.71564693510 + 4164.311989613*t);
   L2 += 0.00000000085*cos(0.49655523670 + 10575.4066829418*t);
   L2 += 0.00000000077*cos(3.51693861509 + 11856.2186514245*t);
   L2 += 0.00000000062*cos(3.62587628690 + 16496.3613962024*t);
   L2 += 0.00000000054*cos(5.25957420065 + 3340.6124266998*t);
   L2 += 0.00000000053*cos(1.10902178415 + 8635.9420037632*t);
   L2 += 0.00000000049*cos(5.65761054625 + 20426.571092422*t);
   L2 += 0.00000000064*cos(5.79211164779 + 2388.8940204492*t);
   L2 += 0.00000000046*cos(5.45092696155 + 6275.9623029906*t);
   L2 += 0.00000000057*cos(4.97077155798 + 14945.3161735544*t);
   L2 += 0.00000000043*cos(3.30685683359 + 9779.1086761254*t);
   L2 += 0.00000000042*cos(1.61412785248 + 12168.0026965746*t);
   L2 += 0.00000000049*cos(3.92715473768 + 5729.506447149*t);
   L2 += 0.00000000042*cos(0.63481258930 + 2699.7348193176*t);
   L2 += 0.00000000056*cos(4.34038639086 + 90955.5516944961*t);
   L2 += 0.00000000040*cos(5.66871428338 + 11712.9553182308*t);
   L2 += 0.00000000039*cos(3.10911294009 + 16200.7727245012*t);
   L2 += 0.00000000040*cos(5.71338386146 + 709.9330485583*t);
   L2 += 0.00000000053*cos(6.17067257683 + 233141.31440436149*t);
   L2 += 0.00000000037*cos(0.32095173508 + 24356.7807886416*t);
   L2 += 0.00000000035*cos(0.95557073457 + 17298.1823273262*t);
   L2 += 0.00000000035*cos(0.64913397996 + 25158.6017197654*t);
   L2 += 0.00000000041*cos(1.53850422484 + 65147.6197681377*t);
   L2 += 0.00000000035*cos(0.77655626359 + 13916.0191096416*t);
   L2 += 0.00000000031*cos(5.35897350775 + 5331.3574437408*t);
   L2 += 0.00000000030*cos(4.48114682755 + 23543.23050468179*t);
   L2 += 0.00000000032*cos(3.45976963453 + 7477.522860216*t);
   L2 += 0.00000000029*cos(3.46648040769 + 13119.72110282519*t);
   L2 += 0.00000000038*cos(2.90863974625 + 12721.572099417*t);
   L2 += 0.00000000029*cos(3.13390968321 + 4136.9104335162*t);
   L2 += 0.00000000035*cos(3.79717126309 + 143571.32428481648*t);
   L2 += 0.00000000027*cos(0.95726093828 + 12559.038152982*t);
   L2 += 0.00000000026*cos(4.95326877030 + 5753.3848848968*t);
   L2 += 0.00000000032*cos(3.49943896928 + 6284.0561710596*t);
   L2 += 0.00000000026*cos(4.59276256636 + 5884.9268465832*t);
   L2 += 0.00000000026*cos(1.53958920253 + 154717.60988768269*t);
   L2 += 0.00000000023*cos(4.88012908735 + 13095.8426650774*t);
   L2 += 0.00000000023*cos(0.35935706511 + 31415.379249957*t);
   L2 += 0.00000000024*cos(5.11515116629 + 18319.5365848796*t);
   L2 += 0.00000000021*cos(5.73872879912 + 12569.6748183318*t);
   L2 += 0.00000000021*cos(1.32901200081 + 10988.808157535*t);
   L2 += 0.00000000019*cos(2.24263229491 + 18073.7049386502*t);
   L2 += 0.00000000019*cos(3.14253175605 + 6496.3749454294*t);
   L2 += 0.00000000021*cos(2.80122025076 + 6282.0955289232*t);
   L2 += 0.00000000023*cos(0.14288760398 + 6283.0085396886*t);
   L2 += 0.00000000019*cos(4.28209473754 + 3930.2096962196*t);
   L2 += 0.00000000016*cos(0.25933207663 + 11790.6290886588*t);
   L2 += 0.00000000015*cos(6.09623501240 + 13517.8701062334*t);
   L2 += 0.00000000020*cos(5.06358906224 + 6283.14316029419*t);
   L2 += 0.00000000015*cos(1.07680663150 + 4933.2084403326*t);
   L2 += 0.00000000019*cos(5.74000581249 + 3128.3887650958*t);
   L2 += 0.00000000016*cos(6.18924226747 + 7342.4577801806*t);
   L2 += 0.00000000013*cos(1.69105044945 + 4535.0594369244*t);
   L2 += 0.00000000015*cos(3.36968394452 + 17260.1546546904*t);
   L2 += 0.00000000010*cos(3.78885035015 + 22003.9146348698*t);
   L2 += 0.00000000011*cos(2.12851973876 + 7875.6718636242*t);
   return L2*t*t;
}



   double Earth_L3 (double t) // 23 terms of order 3
{
   double L3 = 0;
   L3 += 0.00000289058*cos(5.84173149732 + 6283.0758499914*t);
   L3 += 0.00000020712*cos(6.04983939020 + 12566.1516999828*t);
   L3 += 0.00000002962*cos(5.19560579570 + 155.4203994342*t);
   L3 -= 0.00000002527;
   L3 += 0.00000001288*cos(4.72197611970 + 3.523118349*t);
   L3 += 0.00000000635*cos(5.96904899168 + 242.728603974*t);
   L3 += 0.00000000570*cos(5.54182903238 + 18849.2275499742*t);
   L3 += 0.00000000402*cos(3.78606612895 + 553.5694028424*t);
   L3 += 0.00000000072*cos(4.37131884946 + 6286.5989683404*t);
   L3 += 0.00000000067*cos(0.91133898967 + 6127.6554505572*t);
   L3 += 0.00000000037*cos(5.28611190997 + 6438.4962494256*t);
   L3 += 0.00000000021*cos(2.94917211527 + 6309.3741697912*t);
   L3 += 0.00000000015*cos(3.63037493932 + 71430.69561812909*t);
   L3 += 0.00000000011*cos(4.83261533939 + 25132.3033999656*t);
   L3 += 0.00000000011*cos(5.84259014283 + 6525.8044539654*t);
   L3 += 0.00000000012*cos(3.82296977522 + 7058.5984613154*t);
   L3 += 0.00000000013*cos(2.39991715131 + 5729.506447149*t);
   L3 += 0.00000000008*cos(0.55390332094 + 6040.3472460174*t);
   L3 += 0.00000000008*cos(1.46298993048 + 11856.2186514245*t);
   L3 += 0.00000000008*cos(5.07535888338 + 6256.7775301916*t);
   L3 += 0.00000000006*cos(2.88803526743 + 5507.5532386674*t);
   L3 += 0.00000000005*cos(3.87019253131 + 12036.4607348882*t);
   L3 += 0.00000000005*cos(2.70838853362 + 83996.84731811189*t);
   return L3*t*t*t;
}



   double Earth_L4 (double t) // 11 terms of order 4
{
   double L4 = 0;
   L4 += 0.00000007714*cos(4.14117321449 + 6283.0758499914*t);
   L4 += 0.00000001016*cos(3.27573644241 + 12566.1516999828*t);
   L4 += 0.00000000420*cos(0.41892851415 + 155.4203994342*t);
   L4 += 0.00000000047*cos(3.50591071186 + 18849.2275499742*t);
   L4 += 0.00000000041*cos(3.14032562331 + 3.523118349*t);
   L4 += 0.00000000035*cos(5.01110770000 + 5573.1428014331*t);
   L4 += 0.00000000010*cos(5.64816633449 + 6127.6554505572*t);
   L4 += 0.00000000013*cos(0.48609240774 + 77713.7714681205*t);
   L4 += 0.00000000007*cos(2.84139222289 + 161000.6857376741*t);
   L4 += 0.00000000004*cos(3.65509047070 + 25132.3033999656*t);
   L4 += 0.00000000002*cos(0.54880603487 + 6438.4962494256*t);
   return L4*t*t*t*t;
}



   double Earth_L5 (double t) // 4 terms of order 5
{
   double L5 = 0;
   L5 += 0.00000000172*cos(2.74854172392 + 6283.0758499914*t);
   L5 += 0.00000000050*cos(2.01352986713 + 155.4203994342*t);
   L5 += 0.00000000028*cos(2.93369985477 + 12566.1516999828*t);
   L5 += 0.00000000005*cos(1.93829214518 + 18849.2275499742*t);
   return L5*t*t*t*t*t;
}



   double Earth_B0 (double t) // 184 terms of order 0
{
   double B0 = 0;
   B0 += 0.00000279620*cos(3.19870156017 + 84334.66158130829*t);
   B0 += 0.00000101643*cos(5.42248619256 + 5507.5532386674*t);
   B0 += 0.00000080445*cos(3.88013204458 + 5223.6939198022*t);
   B0 += 0.00000043806*cos(3.70444689759 + 2352.8661537718*t);
   B0 += 0.00000031933*cos(4.00026369781 + 1577.3435424478*t);
   B0 += 0.00000022724*cos(3.98473831560 + 1047.7473117547*t);
   B0 += 0.00000016392*cos(3.56456119782 + 5856.4776591154*t);
   B0 += 0.00000018141*cos(4.98367470262 + 6283.0758499914*t);
   B0 += 0.00000014443*cos(3.70275614915 + 9437.762934887*t);
   B0 += 0.00000014304*cos(3.41117857526 + 10213.285546211*t);
   B0 += 0.00000011246*cos(4.82820690527 + 14143.4952424306*t);
   B0 += 0.00000010900*cos(2.08574562329 + 6812.766815086*t);
   B0 += 0.00000009714*cos(3.47303947751 + 4694.0029547076*t);
   B0 += 0.00000010367*cos(4.05663927945 + 71092.88135493269*t);
   B0 += 0.00000008775*cos(4.44016515666 + 5753.3848848968*t);
   B0 += 0.00000008366*cos(4.99251512183 + 7084.8967811152*t);
   B0 += 0.00000006921*cos(4.32559054073 + 6275.9623029906*t);
   B0 += 0.00000009145*cos(1.14182646613 + 6620.8901131878*t);
   B0 += 0.00000007194*cos(3.60193205744 + 529.6909650946*t);
   B0 += 0.00000007698*cos(5.55425745881 + 167621.57585086189*t);
   B0 += 0.00000005285*cos(2.48446991536 + 4705.7323075436*t);
   B0 += 0.00000005208*cos(6.24992674532 + 18073.7049386502*t);
   B0 += 0.00000004529*cos(2.33827747356 + 6309.3741697912*t);
   B0 += 0.00000005579*cos(4.41023653719 + 7860.4193924392*t);
   B0 += 0.00000004743*cos(0.70995680136 + 5884.9268465832*t);
   B0 += 0.00000004301*cos(1.10255777773 + 6681.2248533996*t);
   B0 += 0.00000003849*cos(1.82229412532 + 5486.777843175*t);
   B0 += 0.00000004093*cos(5.11700141197 + 13367.9726311066*t);
   B0 += 0.00000003681*cos(0.43793170336 + 3154.6870848956*t);
   B0 += 0.00000003420*cos(5.42034800952 + 6069.7767545534*t);
   B0 += 0.00000003617*cos(6.04641937568 + 3930.2096962196*t);
   B0 += 0.00000003670*cos(4.58210192227 + 12194.0329146209*t);
   B0 += 0.00000002918*cos(1.95463881136 + 10977.078804699*t);
   B0 += 0.00000002797*cos(5.61259274877 + 11790.6290886588*t);
   B0 += 0.00000002502*cos(0.60499729368 + 6496.3749454294*t);
   B0 += 0.00000002319*cos(5.01648216088 + 1059.3819301892*t);
   B0 += 0.00000002684*cos(1.39470396487 + 22003.9146348698*t);
   B0 += 0.00000002428*cos(3.24183056545 + 78051.5857313169*t);
   B0 += 0.00000002120*cos(4.30691000285 + 5643.1785636774*t);
   B0 += 0.00000002257*cos(3.15557225087 + 90617.7374312997*t);
   B0 += 0.00000001813*cos(3.75574218286 + 3340.6124266998*t);
   B0 += 0.00000002226*cos(2.79699346673 + 12036.4607348882*t);
   B0 += 0.00000001888*cos(0.86991545944 + 8635.9420037632*t);
   B0 += 0.00000001517*cos(1.95852055710 + 398.1490034082*t);
   B0 += 0.00000001581*cos(3.19976230948 + 5088.6288397668*t);
   B0 += 0.00000001421*cos(6.25530883828 + 2544.3144198834*t);
   B0 += 0.00000001595*cos(0.25619915132 + 17298.1823273262*t);
   B0 += 0.00000001391*cos(4.69964175574 + 7058.5984613154*t);
   B0 += 0.00000001478*cos(2.81808207569 + 25934.1243310894*t);
   B0 += 0.00000001481*cos(3.65823554610 + 11506.7697697936*t);
   B0 += 0.00000001693*cos(4.95689385293 + 156475.2902479957*t);
   B0 += 0.00000001183*cos(1.29343060777 + 775.522611324*t);
   B0 += 0.00000001114*cos(2.37889311847 + 3738.761430108*t);
   B0 += 0.00000000994*cos(4.30088900425 + 9225.539273283*t);
   B0 += 0.00000000924*cos(3.06451026809 + 4164.311989613*t);
   B0 += 0.00000000867*cos(0.55606931068 + 8429.2412664666*t);
   B0 += 0.00000000988*cos(5.97286104208 + 7079.3738568078*t);
   B0 += 0.00000000824*cos(1.50984806177 + 10447.3878396044*t);
   B0 += 0.00000000915*cos(0.12635654605 + 11015.1064773348*t);
   B0 += 0.00000000742*cos(1.99159139281 + 26087.9031415742*t);
   B0 -= 0.00000001039;
   B0 += 0.00000000850*cos(4.24120016095 + 29864.334027309*t);
   B0 += 0.00000000755*cos(2.89631873314 + 4732.0306273434*t);
   B0 += 0.00000000714*cos(1.37548118605 + 2146.1654164752*t);
   B0 += 0.00000000708*cos(1.91406542362 + 8031.0922630584*t);
   B0 += 0.00000000746*cos(0.57893808622 + 796.2980068164*t);
   B0 += 0.00000000802*cos(5.12339137235 + 2942.4634232916*t);
   B0 += 0.00000000751*cos(1.67479850166 + 21228.3920235458*t);
   B0 += 0.00000000602*cos(4.09976538900 + 64809.80550494129*t);
   B0 += 0.00000000594*cos(3.49580704849 + 16496.3613962024*t);
   B0 += 0.00000000592*cos(4.59481504319 + 4690.4798363586*t);
   B0 += 0.00000000530*cos(5.73979295194 + 8827.3902698748*t);
   B0 += 0.00000000503*cos(5.66433137112 + 33794.5437235286*t);
   B0 += 0.00000000483*cos(1.57106522322 + 801.8209311238*t);
   B0 += 0.00000000438*cos(0.06707733720 + 3128.3887650958*t);
   B0 += 0.00000000423*cos(2.86944596145 + 12566.1516999828*t);
   B0 += 0.00000000504*cos(3.26207669160 + 7632.9432596502*t);
   B0 += 0.00000000552*cos(1.02926440457 + 239762.20451754928*t);
   B0 += 0.00000000427*cos(3.67434378208 + 213.299095438*t);
   B0 += 0.00000000404*cos(1.46193297360 + 15720.8387848784*t);
   B0 += 0.00000000503*cos(4.85802444134 + 6290.1893969922*t);
   B0 += 0.00000000417*cos(0.81920713533 + 5216.5803728014*t);
   B0 += 0.00000000365*cos(0.01002966145 + 12168.0026965746*t);
   B0 += 0.00000000363*cos(1.28376436579 + 6206.8097787158*t);
   B0 += 0.00000000353*cos(4.70059133110 + 7234.794256242*t);
   B0 += 0.00000000415*cos(0.96862624176 + 4136.9104335162*t);
   B0 += 0.00000000387*cos(3.09145061418 + 25158.6017197654*t);
   B0 += 0.00000000373*cos(2.65119262808 + 7342.4577801806*t);
   B0 += 0.00000000361*cos(2.97762937735 + 9623.6882766912*t);
   B0 += 0.00000000418*cos(3.75759994446 + 5230.807466803*t);
   B0 += 0.00000000396*cos(1.22507712354 + 6438.4962494256*t);
   B0 += 0.00000000322*cos(1.21162178805 + 8662.240323563*t);
   B0 += 0.00000000284*cos(5.64170320179 + 1589.0728952838*t);
   B0 += 0.00000000379*cos(1.72248432756 + 14945.3161735544*t);
   B0 += 0.00000000320*cos(3.94161159658 + 7330.8231617461*t);
   B0 += 0.00000000313*cos(5.47602376451 + 1194.4470102246*t);
   B0 += 0.00000000292*cos(1.38971327568 + 11769.8536931664*t);
   B0 += 0.00000000305*cos(0.80429352049 + 37724.7534197482*t);
   B0 += 0.00000000257*cos(5.81382810029 + 426.598190876*t);
   B0 += 0.00000000265*cos(6.10358507671 + 6836.6452528338*t);
   B0 += 0.00000000250*cos(4.56452895547 + 7477.522860216*t);
   B0 += 0.00000000266*cos(2.62926282354 + 7238.67559160*t);
   B0 += 0.00000000263*cos(6.22089501237 + 6133.5126528568*t);
   B0 += 0.00000000306*cos(2.79682380532 + 1748.016413067*t);
   B0 += 0.00000000236*cos(2.46093023707 + 11371.7046897582*t);
   B0 += 0.00000000316*cos(1.62662805006 + 250908.49012041549*t);
   B0 += 0.00000000216*cos(3.68721275185 + 5849.3641121146*t);
   B0 += 0.00000000230*cos(0.36165162947 + 5863.5912061162*t);
   B0 += 0.00000000233*cos(5.03509933618 + 20426.571092422*t);
   B0 += 0.00000000200*cos(5.86073159059 + 4535.0594369244*t);
   B0 += 0.00000000277*cos(4.65400292395 + 82239.16695779889*t);
   B0 += 0.00000000209*cos(3.72323200803 + 10973.55568635*t);
   B0 += 0.00000000199*cos(5.05186622555 + 5429.8794682394*t);
   B0 += 0.00000000256*cos(2.40923279886 + 19651.048481098*t);
   B0 += 0.00000000210*cos(4.50691909144 + 29088.811415985*t);
   B0 += 0.00000000181*cos(6.00294783127 + 4292.3308329504*t);
   B0 += 0.00000000249*cos(0.12900984422 + 154379.79562448629*t);
   B0 += 0.00000000209*cos(3.87759458541 + 17789.845619785*t);
   B0 += 0.00000000225*cos(3.18339652605 + 18875.525869774*t);
   B0 += 0.00000000191*cos(4.53897489216 + 18477.1087646123*t);
   B0 += 0.00000000172*cos(2.09694182730 + 13095.8426650774*t);
   B0 += 0.00000000182*cos(3.16107943487 + 16730.4636895958*t);
   B0 += 0.00000000188*cos(2.22746128596 + 41654.9631159678*t);
   B0 += 0.00000000164*cos(5.18686274999 + 5481.2549188676*t);
   B0 += 0.00000000160*cos(2.49298855020 + 12592.4500197826*t);
   B0 += 0.00000000155*cos(1.59595438224 + 10021.8372800994*t);
   B0 += 0.00000000135*cos(0.21349051305 + 10988.808157535*t);
   B0 += 0.00000000178*cos(3.80375178044 + 23581.2581773176*t);
   B0 += 0.00000000123*cos(1.66800739151 + 15110.4661198662*t);
   B0 += 0.00000000122*cos(2.72678272224 + 18849.2275499742*t);
   B0 += 0.00000000126*cos(1.17675512872 + 14919.0178537546*t);
   B0 += 0.00000000142*cos(3.95053440992 + 337.8142631964*t);
   B0 += 0.00000000116*cos(6.06340906212 + 6709.6740408674*t);
   B0 += 0.00000000137*cos(3.52143246244 + 12139.5535091068*t);
   B0 += 0.00000000136*cos(2.92179113491 + 32217.2001810808*t);
   B0 += 0.00000000110*cos(3.51203379263 + 18052.9295431578*t);
   B0 += 0.00000000147*cos(4.63371971408 + 22805.7355659936*t);
   B0 += 0.00000000108*cos(5.45280815225 + 7.1135470008*t);
   B0 += 0.00000000148*cos(0.65447253687 + 95480.9471841745*t);
   B0 += 0.00000000119*cos(5.92110458985 + 33019.0211122046*t);
   B0 += 0.00000000110*cos(5.34824206403 + 639.897286314*t);
   B0 += 0.00000000106*cos(3.71081682614 + 14314.1681130498*t);
   B0 += 0.00000000139*cos(6.17607198262 + 24356.7807886416*t);
   B0 += 0.00000000118*cos(5.59738712949 + 161338.5000008705*t);
   B0 += 0.00000000117*cos(3.65065271640 + 45585.1728121874*t);
   B0 += 0.00000000127*cos(4.74596574209 + 49515.382508407*t);
   B0 += 0.00000000120*cos(1.04211499785 + 6915.8595893046*t);
   B0 += 0.00000000120*cos(5.60638811846 + 5650.2921106782*t);
   B0 += 0.00000000115*cos(3.10668213303 + 14712.317116458*t);
   B0 += 0.00000000099*cos(0.69018939930 + 12779.4507954208*t);
   B0 += 0.00000000097*cos(1.07908724794 + 9917.6968745098*t);
   B0 += 0.00000000093*cos(2.62295197146 + 17260.1546546904*t);
   B0 += 0.00000000099*cos(4.45774681732 + 4933.2084403326*t);
   B0 += 0.00000000123*cos(1.37488921994 + 28286.9904848612*t);
   B0 += 0.00000000121*cos(5.19767249869 + 27511.4678735372*t);
   B0 += 0.00000000105*cos(0.87192268229 + 77375.95720492408*t);
   B0 += 0.00000000087*cos(3.93637812950 + 17654.7805397496*t);
   B0 += 0.00000000122*cos(2.23956068680 + 83997.09113559539*t);
   B0 += 0.00000000087*cos(4.18201600921 + 22779.4372461938*t);
   B0 += 0.00000000104*cos(4.59580877285 + 1349.8674096588*t);
   B0 += 0.00000000102*cos(2.83545248595 + 12352.8526045448*t);
   B0 += 0.00000000102*cos(3.97386522171 + 10818.1352869158*t);
   B0 += 0.00000000101*cos(4.32892825818 + 36147.4098773004*t);
   B0 += 0.00000000094*cos(5.00001709338 + 150192.21439800429*t);
   B0 += 0.00000000077*cos(3.97199369300 + 1592.5960136328*t);
   B0 += 0.00000000100*cos(6.07733097102 + 26735.9452622132*t);
   B0 += 0.00000000086*cos(5.26029638250 + 28313.288804661*t);
   B0 += 0.00000000093*cos(4.31900620254 + 44809.6502008634*t);
   B0 += 0.00000000076*cos(6.22743405935 + 13521.7514415914*t);
   B0 += 0.00000000072*cos(1.55820597747 + 6256.7775301916*t);
   B0 += 0.00000000082*cos(4.95202664555 + 10575.4066829418*t);
   B0 += 0.00000000082*cos(1.69647647075 + 1990.745017041*t);
   B0 += 0.00000000075*cos(2.29836095644 + 3634.6210245184*t);
   B0 += 0.00000000075*cos(2.66367876557 + 16200.7727245012*t);
   B0 += 0.00000000087*cos(0.26630214764 + 31441.6775697568*t);
   B0 += 0.00000000077*cos(2.25530952876 + 5235.3285382367*t);
   B0 += 0.00000000076*cos(1.09869730334 + 12903.9659631792*t);
   B0 += 0.00000000058*cos(4.28246137794 + 12559.038152982*t);
   B0 += 0.00000000064*cos(5.51112829602 + 173904.65170085328*t);
   B0 += 0.00000000056*cos(2.60133794851 + 73188.3759784421*t);
   B0 += 0.00000000055*cos(5.81483150022 + 143233.51002162008*t);
   B0 += 0.00000000054*cos(3.38482031504 + 323049.11878710288*t);
   B0 += 0.00000000039*cos(3.28500401937 + 71768.50988132549*t);
   B0 += 0.00000000039*cos(3.11239910096 + 96900.81328129109*t);
   return B0;
}



   double Earth_B1 (double t) // 134 terms of order 1
{
   double B1 = 0;
   B1 += 0.00227777722*cos(3.41376620530 + 6283.0758499914*t);
   B1 += 0.00003805678*cos(3.37063423795 + 12566.1516999828*t);
   B1 += 0.00003619589;
   B1 += 0.00000071542*cos(3.32777549735 + 18849.2275499742*t);
   B1 += 0.00000007655*cos(1.79489607186 + 5223.6939198022*t);
   B1 += 0.00000008107*cos(3.89190403643 + 5507.5532386674*t);
   B1 += 0.00000006456*cos(5.19789424750 + 2352.8661537718*t);
   B1 += 0.00000003894*cos(2.15568517178 + 6279.5527316424*t);
   B1 += 0.00000003892*cos(1.53021064904 + 6286.5989683404*t);
   B1 += 0.00000003897*cos(4.87293945629 + 10213.285546211*t);
   B1 += 0.00000003812*cos(1.43523182316 + 12036.4607348882*t);
   B1 += 0.00000003577*cos(2.32913869227 + 83996.84731811189*t);
   B1 += 0.00000003570*cos(4.92637739003 + 71430.69561812909*t);
   B1 += 0.00000003494*cos(2.20864641831 + 529.6909650946*t);
   B1 += 0.00000002421*cos(6.22876183393 + 7860.4193924392*t);
   B1 += 0.00000002056*cos(3.06747139741 + 14143.4952424306*t);
   B1 += 0.00000001399*cos(0.50107877909 + 6309.3741697912*t);
   B1 += 0.00000001417*cos(3.28454570977 + 25132.3033999656*t);
   B1 += 0.00000001544*cos(1.82062047625 + 5856.4776591154*t);
   B1 += 0.00000001457*cos(1.75339303307 + 5884.9268465832*t);
   B1 += 0.00000001497*cos(2.19673914456 + 9437.762934887*t);
   B1 += 0.00000001549*cos(5.73650061398 + 17789.845619785*t);
   B1 += 0.00000001277*cos(3.96721717330 + 4705.7323075436*t);
   B1 += 0.00000001038*cos(2.94818701170 + 6256.7775301916*t);
   B1 += 0.00000001018*cos(2.24114547164 + 6681.2248533996*t);
   B1 += 0.00000001021*cos(1.46795026130 + 11790.6290886588*t);
   B1 += 0.00000000916*cos(3.72965830745 + 1059.3819301892*t);
   B1 += 0.00000001156*cos(6.04591336669 + 398.1490034082*t);
   B1 += 0.00000001240*cos(0.77195902957 + 6812.766815086*t);
   B1 += 0.00000001115*cos(3.92255876225 + 12168.0026965746*t);
   B1 += 0.00000000908*cos(1.78447918237 + 3930.2096962196*t);
   B1 += 0.00000000833*cos(0.09941579828 + 11506.7697697936*t);
   B1 += 0.00000000790*cos(2.45555993228 + 775.522611324*t);
   B1 += 0.00000000666*cos(4.40761401665 + 5753.3848848968*t);
   B1 += 0.00000000520*cos(1.60179602491 + 4694.0029547076*t);
   B1 += 0.00000000382*cos(0.26754532042 + 1577.3435424478*t);
   B1 += 0.00000000405*cos(1.97558286872 + 6283.0085396886*t);
   B1 += 0.00000000405*cos(1.53147989887 + 6283.14316029419*t);
   B1 += 0.00000000388*cos(2.59563818411 + 7058.5984613154*t);
   B1 += 0.00000000341*cos(3.61275156842 + 13367.9726311066*t);
   B1 += 0.00000000280*cos(4.63052251735 + 796.2980068164*t);
   B1 += 0.00000000276*cos(0.43350778219 + 7079.3738568078*t);
   B1 += 0.00000000311*cos(5.28219636656 + 17260.1546546904*t);
   B1 += 0.00000000275*cos(5.27597553634 + 11769.8536931664*t);
   B1 += 0.00000000286*cos(2.98639716345 + 6275.9623029906*t);
   B1 += 0.00000000251*cos(2.81315684448 + 6290.1893969922*t);
   B1 += 0.00000000239*cos(5.77837903893 + 10977.078804699*t);
   B1 += 0.00000000228*cos(0.14375973844 + 3738.761430108*t);
   B1 += 0.00000000227*cos(2.51020991853 + 7084.8967811152*t);
   B1 += 0.00000000211*cos(4.72029465470 + 6496.3749454294*t);
   B1 += 0.00000000234*cos(3.11535274010 + 709.9330485583*t);
   B1 += 0.00000000234*cos(0.57075817620 + 11856.2186514245*t);
   B1 += 0.00000000237*cos(2.38791907394 + 213.299095438*t);
   B1 += 0.00000000227*cos(1.12059781634 + 12352.8526045448*t);
   B1 += 0.00000000169*cos(3.20148089605 + 5486.777843175*t);
   B1 += 0.00000000159*cos(3.14604135756 + 8827.3902698748*t);
   B1 += 0.00000000160*cos(1.14784478002 + 5643.1785636774*t);
   B1 += 0.00000000144*cos(5.23285657431 + 78051.5857313169*t);
   B1 += 0.00000000144*cos(1.16454654602 + 90617.7374312997*t);
   B1 += 0.00000000144*cos(3.81203756929 + 6262.300454499*t);
   B1 += 0.00000000143*cos(5.97809021355 + 6303.8512454838*t);
   B1 += 0.00000000138*cos(4.29829933273 + 1589.0728952838*t);
   B1 += 0.00000000161*cos(3.03298851492 + 20426.571092422*t);
   B1 += 0.00000000165*cos(5.09134896587 + 17298.1823273262*t);
   B1 += 0.00000000128*cos(4.41590143012 + 6127.6554505572*t);
   B1 += 0.00000000158*cos(2.00984506334 + 5230.807466803*t);
   B1 += 0.00000000130*cos(1.44170683802 + 12569.6748183318*t);
   B1 += 0.00000000125*cos(1.69950379498 + 3128.3887650958*t);
   B1 += 0.00000000152*cos(4.96946111415 + 65147.6197681377*t);
   B1 += 0.00000000131*cos(4.24961399629 + 6282.0955289232*t);
   B1 += 0.00000000131*cos(5.54051222995 + 6284.0561710596*t);
   B1 += 0.00000000161*cos(3.32421999691 + 6283.3196674749*t);
   B1 += 0.00000000123*cos(2.93221463795 + 15720.8387848784*t);
   B1 += 0.00000000152*cos(1.56488157456 + 18319.5365848796*t);
   B1 += 0.00000000114*cos(1.89110005546 + 18073.7049386502*t);
   B1 += 0.00000000113*cos(4.95214866677 + 4136.9104335162*t);
   B1 += 0.00000000111*cos(4.79699611405 + 8429.2412664666*t);
   B1 += 0.00000000118*cos(3.06184958762 + 22003.9146348698*t);
   B1 += 0.00000000128*cos(2.09693027395 + 12562.6285816338*t);
   B1 += 0.00000000138*cos(0.84648544836 + 6069.7767545534*t);
   B1 += 0.00000000100*cos(0.41938790104 + 5481.2549188676*t);
   B1 += 0.00000000101*cos(4.43791289246 + 19651.048481098*t);
   B1 += 0.00000000116*cos(1.87002428935 + 77713.7714681205*t);
   B1 += 0.00000000095*cos(0.68638183645 + 3340.6124266998*t);
   B1 += 0.00000000104*cos(4.90822646457 + 8635.9420037632*t);
   B1 += 0.00000000114*cos(0.36008214928 + 23543.23050468179*t);
   B1 += 0.00000000091*cos(1.17611219310 + 9225.539273283*t);
   B1 += 0.00000000091*cos(2.74682631104 + 5216.5803728014*t);
   B1 += 0.00000000086*cos(3.09315282195 + 1194.4470102246*t);
   B1 += 0.00000000093*cos(6.27643513190 + 12559.038152982*t);
   B1 += 0.00000000087*cos(3.64926989547 + 7.1135470008*t);
   B1 += 0.00000000083*cos(0.37008971174 + 11371.7046897582*t);
   B1 += 0.00000000098*cos(4.68473306376 + 167283.76158766549*t);
   B1 += 0.00000000098*cos(0.99865886857 + 154717.60988768269*t);
   B1 += 0.00000000086*cos(4.41152307486 + 25934.1243310894*t);
   B1 += 0.00000000098*cos(4.98209568932 + 16496.3613962024*t);
   B1 += 0.00000000085*cos(2.35438163823 + 3154.6870848956*t);
   B1 += 0.00000000083*cos(0.24859477658 + 21228.3920235458*t);
   B1 += 0.00000000076*cos(5.46661600296 + 29864.334027309*t);
   B1 += 0.00000000069*cos(3.22045334237 + 801.8209311238*t);
   B1 += 0.00000000068*cos(2.19928242745 + 37724.7534197482*t);
   B1 += 0.00000000070*cos(5.33711014698 + 2544.3144198834*t);
   B1 += 0.00000000055*cos(5.37872238211 + 11015.1064773348*t);
   B1 += 0.00000000051*cos(5.03600618727 + 5863.5912061162*t);
   B1 += 0.00000000059*cos(3.54984809612 + 4535.0594369244*t);
   B1 += 0.00000000053*cos(1.50812064137 + 7342.4577801806*t);
   B1 += 0.00000000044*cos(1.33712241647 + 4164.311989613*t);
   B1 += 0.00000000055*cos(5.38460045253 + 7477.522860216*t);
   B1 += 0.00000000049*cos(5.61844340512 + 24072.9214697764*t);
   B1 += 0.00000000057*cos(6.17973522898 + 8031.0922630584*t);
   B1 += 0.00000000052*cos(4.24379064407 + 5088.6288397668*t);
   B1 += 0.00000000041*cos(3.41960196246 + 26087.9031415742*t);
   B1 += 0.00000000052*cos(1.68150981131 + 25158.6017197654*t);
   B1 += 0.00000000053*cos(0.11041408311 + 29088.811415985*t);
   B1 += 0.00000000043*cos(0.53354396923 + 2146.1654164752*t);
   B1 += 0.00000000050*cos(0.60270799844 + 33794.5437235286*t);
   B1 += 0.00000000045*cos(1.69577010458 + 6702.5604938666*t);
   B1 += 0.00000000037*cos(6.09033460795 + 64809.80550494129*t);
   B1 += 0.00000000044*cos(2.67993061247 + 15110.4661198662*t);
   B1 += 0.00000000047*cos(3.62555676035 + 41654.9631159678*t);
   B1 += 0.00000000036*cos(0.47219666269 + 13095.8426650774*t);
   B1 += 0.00000000035*cos(5.99520523215 + 23581.2581773176*t);
   B1 += 0.00000000043*cos(4.08726331262 + 156137.47598479928*t);
   B1 += 0.00000000043*cos(5.17376496602 + 18422.62935909819*t);
   B1 += 0.00000000034*cos(2.14460100274 + 4933.2084403326*t);
   B1 += 0.00000000030*cos(4.12992793541 + 426.598190876*t);
   B1 += 0.00000000036*cos(2.28596930002 + 90279.92316810328*t);
   B1 += 0.00000000032*cos(0.81117747619 + 26.2983197998*t);
   B1 += 0.00000000028*cos(4.07036560467 + 14712.317116458*t);
   B1 += 0.00000000025*cos(1.92905243842 + 1748.016413067*t);
   B1 += 0.00000000024*cos(5.06152880842 + 16730.4636895958*t);
   B1 += 0.00000000023*cos(3.24047012195 + 31415.379249957*t);
   B1 += 0.00000000027*cos(3.90672018305 + 18451.07854656599*t);
   B1 += 0.00000000021*cos(3.51419811826 + 10447.3878396044*t);
   return B1*t;
}



   double Earth_B2 (double t) // 62 terms of order 2
{
   double B2 = 0;
   B2 += 0.00009721424*cos(5.15192809920 + 6283.0758499914*t);
   B2 -= 0.00000233002;
   B2 += 0.00000134188*cos(0.64406212977 + 12566.1516999828*t);
   B2 += 0.00000006504*cos(1.07333397797 + 18849.2275499742*t);
   B2 += 0.00000001662*cos(1.62746869551 + 84334.66158130829*t);
   B2 += 0.00000000635*cos(3.51985338656 + 6279.5527316424*t);
   B2 += 0.00000000492*cos(2.41382223971 + 1047.7473117547*t);
   B2 += 0.00000000307*cos(6.10181422085 + 5223.6939198022*t);
   B2 += 0.00000000322*cos(0.37660897389 + 6286.5989683404*t);
   B2 += 0.00000000326*cos(2.35727931602 + 5507.5532386674*t);
   B2 += 0.00000000274*cos(1.65307581765 + 7860.4193924392*t);
   B2 += 0.00000000228*cos(1.14082932988 + 25132.3033999656*t);
   B2 += 0.00000000202*cos(0.49836682530 + 2352.8661537718*t);
   B2 += 0.00000000201*cos(0.15552765600 + 10213.285546211*t);
   B2 += 0.00000000167*cos(3.98005254015 + 529.6909650946*t);
   B2 += 0.00000000170*cos(5.28668290523 + 6256.7775301916*t);
   B2 += 0.00000000166*cos(3.04613930284 + 12036.4607348882*t);
   B2 += 0.00000000153*cos(4.06779216239 + 83996.84731811189*t);
   B2 += 0.00000000150*cos(3.18772213951 + 71430.69561812909*t);
   B2 += 0.00000000120*cos(3.13558669517 + 5884.9268465832*t);
   B2 += 0.00000000123*cos(4.17102530625 + 6309.3741697912*t);
   B2 += 0.00000000100*cos(1.46356761368 + 11506.7697697936*t);
   B2 += 0.00000000100*cos(3.00322421365 + 11790.6290886588*t);
   B2 += 0.00000000077*cos(1.65643898948 + 4705.7323075436*t);
   B2 += 0.00000000068*cos(6.03791904123 + 6812.766815086*t);
   B2 += 0.00000000075*cos(4.85191600582 + 14143.4952424306*t);
   B2 += 0.00000000079*cos(4.12628805658 + 5753.3848848968*t);
   B2 += 0.00000000064*cos(2.71944160460 + 6127.6554505572*t);
   B2 += 0.00000000063*cos(0.78655326011 + 6438.4962494256*t);
   B2 += 0.00000000075*cos(0.84213523741 + 167621.57585086189*t);
   B2 += 0.00000000051*cos(4.55059044701 + 5486.777843175*t);
   B2 += 0.00000000050*cos(5.29314320585 + 7079.3738568078*t);
   B2 += 0.00000000045*cos(0.33147576416 + 775.522611324*t);
   B2 += 0.00000000043*cos(3.61793711640 + 1577.3435424478*t);
   B2 += 0.00000000057*cos(6.15295833679 + 12194.0329146209*t);
   B2 += 0.00000000045*cos(1.18274698508 + 17789.845619785*t);
   B2 += 0.00000000047*cos(1.11643162773 + 398.1490034082*t);
   B2 += 0.00000000047*cos(5.52984238430 + 12168.0026965746*t);
   B2 += 0.00000000042*cos(2.64375300620 + 10988.808157535*t);
   B2 += 0.00000000037*cos(0.31226891972 + 6681.2248533996*t);
   B2 += 0.00000000046*cos(3.38617099014 + 156475.2902479957*t);
   B2 += 0.00000000033*cos(1.26266496002 + 6290.1893969922*t);
   B2 += 0.00000000033*cos(2.25067065498 + 6275.9623029906*t);
   B2 += 0.00000000027*cos(1.72859626293 + 1059.3819301892*t);
   B2 += 0.00000000028*cos(2.88681054153 + 796.2980068164*t);
   B2 += 0.00000000030*cos(5.52021264181 + 7058.5984613154*t);
   B2 += 0.00000000021*cos(6.22275008403 + 6069.7767545534*t);
   B2 += 0.00000000024*cos(5.83690442827 + 6282.0955289232*t);
   B2 += 0.00000000024*cos(3.95322179797 + 6284.0561710596*t);
   B2 += 0.00000000021*cos(4.87022458966 + 3738.761430108*t);
   B2 += 0.00000000021*cos(3.58048145762 + 6496.3749454294*t);
   B2 += 0.00000000022*cos(1.97770318395 + 7.1135470008*t);
   B2 += 0.00000000019*cos(3.45979888770 + 4136.9104335162*t);
   B2 += 0.00000000015*cos(2.81379880771 + 1589.0728952838*t);
   B2 += 0.00000000015*cos(1.61313151838 + 1194.4470102246*t);
   B2 += 0.00000000015*cos(0.63083434831 + 11769.8536931664*t);
   B2 += 0.00000000019*cos(0.28082047767 + 17260.1546546904*t);
   B2 += 0.00000000017*cos(1.35165136146 + 78051.5857313169*t);
   B2 += 0.00000000015*cos(2.30941599560 + 11856.2186514245*t);
   B2 += 0.00000000012*cos(1.12997370513 + 90617.7374312997*t);
   B2 += 0.00000000012*cos(0.00611601597 + 18073.7049386502*t);
   B2 += 0.00000000012*cos(1.70058636365 + 12559.038152982*t);
   return B2*t*t;
}



   double Earth_B3 (double t) // 14 terms of order 3
{
   double B3 = 0;
   B3 += 0.00000275993*cos(0.59480097092 + 6283.0758499914*t);
   B3 -= 0.00000017034;
   B3 += 0.00000003617*cos(0.11750575325 + 12566.1516999828*t);
   B3 += 0.00000000339*cos(5.66087461682 + 18849.2275499742*t);
   B3 += 0.00000000056*cos(5.02765554835 + 6279.5527316424*t);
   B3 += 0.00000000019*cos(5.99007646261 + 6256.7775301916*t);
   B3 += 0.00000000018*cos(3.80004734567 + 6309.3741697912*t);
   B3 += 0.00000000018*cos(1.21049250774 + 6127.6554505572*t);
   B3 += 0.00000000018*cos(2.29734567137 + 6438.4962494256*t);
   B3 += 0.00000000015*cos(4.72881467263 + 6286.5989683404*t);
   B3 += 0.00000000012*cos(0.41481671808 + 83996.84731811189*t);
   B3 += 0.00000000013*cos(5.54637369296 + 25132.3033999656*t);
   B3 += 0.00000000010*cos(2.91937214232 + 71430.69561812909*t);
   B3 += 0.00000000006*cos(2.14173241210 + 11856.2186514245*t);
   return B3*t*t*t;
}



   double Earth_B4 (double t) // 6 terms of order 4
{
   double B4 = 0;
   B4 += 0.00000005745*cos(2.26734029843 + 6283.0758499914*t);
   B4 += 0.00000000870;
   B4 += 0.00000000119*cos(4.26807972611 + 12566.1516999828*t);
   B4 += 0.00000000017*cos(4.07422620440 + 18849.2275499742*t);
   B4 += 0.00000000005*cos(0.84308705203 + 1047.7473117547*t);
   B4 += 0.00000000005*cos(0.05711572303 + 84334.66158130829*t);
   return B4*t*t*t*t;
}



   double Earth_B5 (double t) // 2 terms of order 5
{
   double B5 = 0;
   B5 += 0.00000000114*cos(4.31455980099 + 6283.0758499914*t);
   B5 += 0.00000000024;
   return B5*t*t*t*t*t;
}



   double Earth_R0 (double t) // 523 terms of order 0
{
   double R0 = 0;
   R0 += 1.00013988784;
   R0 += 0.01670699632*cos(3.09846350258 + 6283.0758499914*t);
   R0 += 0.00013956024*cos(3.05524609456 + 12566.1516999828*t);
   R0 += 0.00003083720*cos(5.19846674381 + 77713.7714681205*t);
   R0 += 0.00001628463*cos(1.17387558054 + 5753.3848848968*t);
   R0 += 0.00001575572*cos(2.84685214877 + 7860.4193924392*t);
   R0 += 0.00000924799*cos(5.45292236722 + 11506.7697697936*t);
   R0 += 0.00000542439*cos(4.56409151453 + 3930.2096962196*t);
   R0 += 0.00000472110*cos(3.66100022149 + 5884.9268465832*t);
   R0 += 0.00000328780*cos(5.89983686142 + 5223.6939198022*t);
   R0 += 0.00000345969*cos(0.96368627272 + 5507.5532386674*t);
   R0 += 0.00000306784*cos(0.29867139512 + 5573.1428014331*t);
   R0 += 0.00000174844*cos(3.01193636733 + 18849.2275499742*t);
   R0 += 0.00000243181*cos(4.27349530790 + 11790.6290886588*t);
   R0 += 0.00000211836*cos(5.84714461348 + 1577.3435424478*t);
   R0 += 0.00000185740*cos(5.02199710705 + 10977.078804699*t);
   R0 += 0.00000109835*cos(5.05510635860 + 5486.777843175*t);
   R0 += 0.00000098316*cos(0.88681311278 + 6069.7767545534*t);
   R0 += 0.00000086500*cos(5.68956418946 + 15720.8387848784*t);
   R0 += 0.00000085831*cos(1.27079125277 + 161000.6857376741*t);
   R0 += 0.00000062917*cos(0.92177053978 + 529.6909650946*t);
   R0 += 0.00000057056*cos(2.01374292245 + 83996.84731811189*t);
   R0 += 0.00000064908*cos(0.27251341435 + 17260.1546546904*t);
   R0 += 0.00000049384*cos(3.24501240359 + 2544.3144198834*t);
   R0 += 0.00000055736*cos(5.24159799170 + 71430.69561812909*t);
   R0 += 0.00000042520*cos(6.01110257982 + 6275.9623029906*t);
   R0 += 0.00000046966*cos(2.57799853213 + 775.522611324*t);
   R0 += 0.00000038963*cos(5.36063832897 + 4694.0029547076*t);
   R0 += 0.00000044666*cos(5.53715663816 + 9437.762934887*t);
   R0 += 0.00000035661*cos(1.67447135798 + 12036.4607348882*t);
   R0 += 0.00000031922*cos(0.18368299942 + 5088.6288397668*t);
   R0 += 0.00000031846*cos(1.77775642078 + 398.1490034082*t);
   R0 += 0.00000033193*cos(0.24370221704 + 7084.8967811152*t);
   R0 += 0.00000038245*cos(2.39255343973 + 8827.3902698748*t);
   R0 += 0.00000028468*cos(1.21344887533 + 6286.5989683404*t);
   R0 += 0.00000037486*cos(0.82961281844 + 19651.048481098*t);
   R0 += 0.00000036957*cos(4.90107587287 + 12139.5535091068*t);
   R0 += 0.00000034537*cos(1.84270693281 + 2942.4634232916*t);
   R0 += 0.00000026275*cos(4.58896863104 + 10447.3878396044*t);
   R0 += 0.00000024596*cos(3.78660838036 + 8429.2412664666*t);
   R0 += 0.00000023587*cos(0.26866098169 + 796.2980068164*t);
   R0 += 0.00000027795*cos(1.89934427832 + 6279.5527316424*t);
   R0 += 0.00000023927*cos(4.99598548145 + 5856.4776591154*t);
   R0 += 0.00000020345*cos(4.65282190725 + 2146.1654164752*t);
   R0 += 0.00000023287*cos(2.80783632869 + 14143.4952424306*t);
   R0 += 0.00000022099*cos(1.95002636847 + 3154.6870848956*t);
   R0 += 0.00000019509*cos(5.38233922479 + 2352.8661537718*t);
   R0 += 0.00000017958*cos(0.19871369960 + 6812.766815086*t);
   R0 += 0.00000017178*cos(4.43322156854 + 10213.285546211*t);
   R0 += 0.00000016190*cos(5.23159323213 + 17789.845619785*t);
   R0 += 0.00000017315*cos(6.15224075188 + 16730.4636895958*t);
   R0 += 0.00000013814*cos(5.18962074032 + 8031.0922630584*t);
   R0 += 0.00000018834*cos(0.67280058021 + 149854.40013480789*t);
   R0 += 0.00000018330*cos(2.25348717053 + 23581.2581773176*t);
   R0 += 0.00000013639*cos(3.68511810757 + 4705.7323075436*t);
   R0 += 0.00000013142*cos(0.65267698994 + 13367.9726311066*t);
   R0 += 0.00000010414*cos(4.33285688501 + 11769.8536931664*t);
   R0 += 0.00000009978*cos(4.20126336356 + 6309.3741697912*t);
   R0 += 0.00000010170*cos(1.59366684542 + 4690.4798363586*t);
   R0 += 0.00000007564*cos(2.62560597391 + 6256.7775301916*t);
   R0 += 0.00000009654*cos(3.67583728703 + 27511.4678735372*t);
   R0 += 0.00000006743*cos(0.56269927047 + 3340.6124266998*t);
   R0 += 0.00000008743*cos(6.06359123461 + 1748.016413067*t);
   R0 += 0.00000007786*cos(3.67371235367 + 12168.0026965746*t);
   R0 += 0.00000006633*cos(5.66149277789 + 11371.7046897582*t);
   R0 += 0.00000007712*cos(0.31242577788 + 7632.9432596502*t);
   R0 += 0.00000006586*cos(3.13580054586 + 801.8209311238*t);
   R0 += 0.00000007460*cos(5.64758066660 + 11926.2544136688*t);
   R0 += 0.00000006933*cos(2.92384586372 + 6681.2248533996*t);
   R0 += 0.00000006805*cos(1.42327153767 + 23013.5395395872*t);
   R0 += 0.00000006118*cos(5.13395999022 + 1194.4470102246*t);
   R0 += 0.00000006477*cos(2.64986648493 + 19804.8272915828*t);
   R0 += 0.00000005233*cos(4.62432817299 + 6438.4962494256*t);
   R0 += 0.00000006147*cos(3.02863936662 + 233141.31440436149*t);
   R0 += 0.00000004608*cos(1.72194702724 + 7234.794256242*t);
   R0 += 0.00000004221*cos(1.55697533726 + 7238.67559160*t);
   R0 += 0.00000005310*cos(2.40821524293 + 11499.6562227928*t);
   R0 += 0.00000005128*cos(5.32398965690 + 11513.8833167944*t);
   R0 += 0.00000004770*cos(0.25554311730 + 11856.2186514245*t);
   R0 += 0.00000005519*cos(2.09089153789 + 17298.1823273262*t);
   R0 += 0.00000005625*cos(4.34052903053 + 90955.5516944961*t);
   R0 += 0.00000004578*cos(4.46569641570 + 5746.271337896*t);
   R0 += 0.00000003788*cos(4.90728294810 + 4164.311989613*t);
   R0 += 0.00000005337*cos(5.09957905103 + 31441.6775697568*t);
   R0 += 0.00000003967*cos(1.20054555175 + 1349.8674096588*t);
   R0 += 0.00000004005*cos(3.02853885902 + 1059.3819301892*t);
   R0 += 0.00000003480*cos(0.76066308841 + 10973.55568635*t);
   R0 += 0.00000004232*cos(1.05485713117 + 5760.4984318976*t);
   R0 += 0.00000004582*cos(3.76570026763 + 6386.16862421*t);
   R0 += 0.00000003335*cos(3.13829943354 + 6836.6452528338*t);
   R0 += 0.00000003420*cos(3.00043974511 + 4292.3308329504*t);
   R0 += 0.00000003595*cos(5.70703236079 + 5643.1785636774*t);
   R0 += 0.00000003236*cos(4.16387400645 + 9917.6968745098*t);
   R0 += 0.00000004154*cos(2.59940749519 + 7058.5984613154*t);
   R0 += 0.00000003362*cos(4.54577164994 + 4732.0306273434*t);
   R0 += 0.00000002978*cos(1.30561268820 + 6283.14316029419*t);
   R0 += 0.00000002765*cos(0.51311975671 + 26.2983197998*t);
   R0 += 0.00000002807*cos(5.66230537649 + 8635.9420037632*t);
   R0 += 0.00000002927*cos(5.73787834080 + 16200.7727245012*t);
   R0 += 0.00000003167*cos(1.69181759900 + 11015.1064773348*t);
   R0 += 0.00000002598*cos(2.96244118358 + 25132.3033999656*t);
   R0 += 0.00000003519*cos(3.62639325753 + 244287.60000722769*t);
   R0 += 0.00000002676*cos(4.20727719487 + 18073.7049386502*t);
   R0 += 0.00000002978*cos(1.74971565805 + 6283.0085396886*t);
   R0 += 0.00000002287*cos(1.06976449088 + 14314.1681130498*t);
   R0 += 0.00000002863*cos(5.92838917309 + 14712.317116458*t);
   R0 += 0.00000003071*cos(0.23793217000 + 35371.8872659764*t);
   R0 += 0.00000002656*cos(0.89959301615 + 12352.8526045448*t);
   R0 += 0.00000002415*cos(2.79975176800 + 709.9330485583*t);
   R0 += 0.00000002811*cos(3.51513864541 + 21228.3920235458*t);
   R0 += 0.00000001977*cos(2.61358297551 + 951.7184062506*t);
   R0 += 0.00000002548*cos(2.47684686575 + 6208.2942514241*t);
   R0 += 0.00000001999*cos(0.56090396506 + 7079.3738568078*t);
   R0 += 0.00000002305*cos(1.05376463592 + 22483.84857449259*t);
   R0 += 0.00000001855*cos(2.86093570752 + 5216.5803728014*t);
   R0 += 0.00000002157*cos(1.31395211105 + 154717.60988768269*t);
   R0 += 0.00000001970*cos(4.36931551625 + 167283.76158766549*t);
   R0 += 0.00000001754*cos(2.14452400686 + 6290.1893969922*t);
   R0 += 0.00000001628*cos(5.85704450617 + 10984.1923516998*t);
   R0 += 0.00000002154*cos(6.03828353794 + 10873.9860304804*t);
   R0 += 0.00000001714*cos(3.70158195222 + 1592.5960136328*t);
   R0 += 0.00000001541*cos(6.21599512982 + 23543.23050468179*t);
   R0 += 0.00000001602*cos(1.99860679677 + 10969.9652576982*t);
   R0 += 0.00000001712*cos(1.34295218697 + 3128.3887650958*t);
   R0 += 0.00000001647*cos(5.54948299069 + 6496.3749454294*t);
   R0 += 0.00000001495*cos(5.43980459648 + 155.4203994342*t);
   R0 += 0.00000001827*cos(5.91227480351 + 3738.761430108*t);
   R0 += 0.00000001726*cos(2.16765465036 + 10575.4066829418*t);
   R0 += 0.00000001532*cos(5.35683107063 + 13521.7514415914*t);
   R0 += 0.00000001824*cos(1.66056145084 + 39302.096962196*t);
   R0 += 0.00000001605*cos(1.90930973224 + 6133.5126528568*t);
   R0 += 0.00000001282*cos(2.46013372544 + 13916.0191096416*t);
   R0 += 0.00000001211*cos(4.41360631550 + 3894.1818295422*t);
   R0 += 0.00000001394*cos(1.77801929250 + 9225.539273283*t);
   R0 += 0.00000001571*cos(4.95512957606 + 25158.6017197654*t);
   R0 += 0.00000001205*cos(1.19212756308 + 3.523118349*t);
   R0 += 0.00000001132*cos(2.69830084955 + 6040.3472460174*t);
   R0 += 0.00000001504*cos(5.77577388271 + 18209.33026366019*t);
   R0 += 0.00000001393*cos(1.62625077326 + 5120.6011455836*t);
   R0 += 0.00000001081*cos(2.93726744446 + 17256.6315363414*t);
   R0 += 0.00000001232*cos(0.71651766504 + 143571.32428481648*t);
   R0 += 0.00000001087*cos(0.99769687961 + 955.5997416086*t);
   R0 += 0.00000001068*cos(5.28472576591 + 65147.6197681377*t);
   R0 += 0.00000001169*cos(3.11663802316 + 14945.3161735544*t);
   R0 += 0.00000000975*cos(5.10887260780 + 6172.869528772*t);
   R0 += 0.00000001202*cos(4.02992510403 + 553.5694028424*t);
   R0 += 0.00000000979*cos(2.00000879106 + 15110.4661198662*t);
   R0 += 0.00000000962*cos(4.02380771400 + 6282.0955289232*t);
   R0 += 0.00000000999*cos(3.62643002790 + 6262.300454499*t);
   R0 += 0.00000001030*cos(5.84987815239 + 213.299095438*t);
   R0 += 0.00000001014*cos(2.84227679965 + 8662.240323563*t);
   R0 += 0.00000001185*cos(1.51330629149 + 17654.7805397496*t);
   R0 += 0.00000000967*cos(2.67081017562 + 5650.2921106782*t);
   R0 += 0.00000001222*cos(2.65423784904 + 88860.05707098669*t);
   R0 += 0.00000000986*cos(2.36212814824 + 6206.8097787158*t);
   R0 += 0.00000001034*cos(0.13634950642 + 11712.9553182308*t);
   R0 += 0.00000001103*cos(3.08477302937 + 43232.3066584156*t);
   R0 += 0.00000000781*cos(2.53374971725 + 16496.3613962024*t);
   R0 += 0.00000001019*cos(3.04569392376 + 6037.244203762*t);
   R0 += 0.00000000795*cos(5.80662989126 + 5230.807466803*t);
   R0 += 0.00000000813*cos(3.57702871938 + 10177.2576795336*t);
   R0 += 0.00000000962*cos(5.31470594766 + 6284.0561710596*t);
   R0 += 0.00000000717*cos(5.95797471837 + 12559.038152982*t);
   R0 += 0.00000000967*cos(2.74413738053 + 6244.9428143536*t);
   R0 += 0.00000000921*cos(0.10160160830 + 29088.811415985*t);
   R0 += 0.00000000719*cos(5.91788189939 + 4136.9104335162*t);
   R0 += 0.00000000688*cos(3.89489045092 + 1589.0728952838*t);
   R0 += 0.00000000772*cos(4.05505380285 + 6127.6554505572*t);
   R0 += 0.00000000706*cos(5.49323197725 + 22003.9146348698*t);
   R0 += 0.00000000665*cos(1.60002747134 + 11087.2851259184*t);
   R0 += 0.00000000690*cos(4.50539825729 + 426.598190876*t);
   R0 += 0.00000000854*cos(3.26104645060 + 20426.571092422*t);
   R0 += 0.00000000656*cos(4.32410182940 + 16858.4825329332*t);
   R0 += 0.00000000840*cos(2.59572585212 + 28766.924424484*t);
   R0 += 0.00000000686*cos(0.61944033771 + 11403.676995575*t);
   R0 += 0.00000000700*cos(3.40901412473 + 7.1135470008*t);
   R0 += 0.00000000728*cos(0.04050185963 + 5481.2549188676*t);
   R0 += 0.00000000653*cos(1.03869451230 + 6062.6632075526*t);
   R0 += 0.00000000559*cos(4.79221805695 + 20199.094959633*t);
   R0 += 0.00000000633*cos(5.70229959167 + 45892.73043315699*t);
   R0 += 0.00000000591*cos(6.10986487621 + 9623.6882766912*t);
   R0 += 0.00000000520*cos(3.62310356479 + 5333.9002410216*t);
   R0 += 0.00000000602*cos(5.58381898589 + 10344.2950653858*t);
   R0 += 0.00000000496*cos(2.21027756314 + 1990.745017041*t);
   R0 += 0.00000000691*cos(1.96733114988 + 12416.5885028482*t);
   R0 += 0.00000000640*cos(1.59062417043 + 18319.5365848796*t);
   R0 += 0.00000000625*cos(3.82358168221 + 13517.8701062334*t);
   R0 += 0.00000000475*cos(1.17025904180 + 12569.6748183318*t);
   R0 += 0.00000000660*cos(5.08498512995 + 283.8593188652*t);
   R0 += 0.00000000664*cos(4.50029469969 + 47162.5163546352*t);
   R0 += 0.00000000569*cos(0.16318535463 + 17267.26820169119*t);
   R0 += 0.00000000568*cos(3.86100969474 + 6076.8903015542*t);
   R0 += 0.00000000462*cos(0.26368763517 + 4590.910180489*t);
   R0 += 0.00000000535*cos(4.83225423196 + 18422.62935909819*t);
   R0 += 0.00000000466*cos(0.75873879417 + 7342.4577801806*t);
   R0 += 0.00000000541*cos(3.07212190556 + 226858.23855437008*t);
   R0 += 0.00000000610*cos(1.53597089605 + 33019.0211122046*t);
   R0 += 0.00000000617*cos(2.62356328726 + 11190.377900137*t);
   R0 += 0.00000000548*cos(4.55798855803 + 18875.525869774*t);
   R0 += 0.00000000633*cos(4.60110281228 + 66567.48586525429*t);
   R0 += 0.00000000587*cos(5.78087907808 + 632.7837393132*t);
   R0 += 0.00000000603*cos(5.38458554802 + 316428.22867391503*t);
   R0 += 0.00000000525*cos(5.01522072363 + 12132.439962106*t);
   R0 += 0.00000000469*cos(0.59975173763 + 21954.15760939799*t);
   R0 += 0.00000000548*cos(3.50627043672 + 17253.04110768959*t);
   R0 += 0.00000000502*cos(0.98804327589 + 11609.8625440122*t);
   R0 += 0.00000000568*cos(1.98497313089 + 7668.6374249425*t);
   R0 += 0.00000000482*cos(1.62460405687 + 12146.6670561076*t);
   R0 += 0.00000000391*cos(3.68718382972 + 18052.9295431578*t);
   R0 += 0.00000000457*cos(3.77214896610 + 156137.47598479928*t);
   R0 += 0.00000000401*cos(5.29221540240 + 15671.0817594066*t);
   R0 += 0.00000000469*cos(1.80963351735 + 12562.6285816338*t);
   R0 += 0.00000000514*cos(3.37031288919 + 20597.2439630412*t);
   R0 += 0.00000000452*cos(5.66811219778 + 10454.5013866052*t);
   R0 += 0.00000000375*cos(4.98528185039 + 9779.1086761254*t);
   R0 += 0.00000000523*cos(0.97215560834 + 155427.54293624099*t);
   R0 += 0.00000000403*cos(5.13948189770 + 1551.045222648*t);
   R0 += 0.00000000372*cos(3.69883738807 + 9388.0059094152*t);
   R0 += 0.00000000367*cos(4.43875659833 + 4535.0594369244*t);
   R0 += 0.00000000406*cos(4.20863156497 + 12592.4500197826*t);
   R0 += 0.00000000362*cos(2.55099560446 + 242.728603974*t);
   R0 += 0.00000000471*cos(4.61907324819 + 5436.9930152402*t);
   R0 += 0.00000000388*cos(4.96020928400 + 24356.7807886416*t);
   R0 += 0.00000000441*cos(5.83872966262 + 3496.032826134*t);
   R0 += 0.00000000349*cos(6.16307810648 + 19800.9459562248*t);
   R0 += 0.00000000356*cos(0.23819081240 + 5429.8794682394*t);
   R0 += 0.00000000346*cos(5.60809622572 + 2379.1644735716*t);
   R0 += 0.00000000380*cos(2.72105213132 + 11933.3679606696*t);
   R0 += 0.00000000432*cos(0.24215988572 + 17996.0311682222*t);
   R0 += 0.00000000378*cos(5.22516848076 + 7477.522860216*t);
   R0 += 0.00000000337*cos(5.10885555836 + 5849.3641121146*t);
   R0 += 0.00000000315*cos(0.57827745123 + 10557.5941608238*t);
   R0 += 0.00000000318*cos(4.49949007320 + 3634.6210245184*t);
   R0 += 0.00000000323*cos(1.55850824803 + 10440.2742926036*t);
   R0 += 0.00000000314*cos(5.77154773334 + 20.7753954924*t);
   R0 += 0.00000000303*cos(2.34615580398 + 4686.8894077068*t);
   R0 += 0.00000000414*cos(5.93237602310 + 51092.7260508548*t);
   R0 += 0.00000000362*cos(2.17561997119 + 28237.2334593894*t);
   R0 += 0.00000000288*cos(0.18377405421 + 13095.8426650774*t);
   R0 += 0.00000000277*cos(5.12952205030 + 13119.72110282519*t);
   R0 += 0.00000000325*cos(6.18608287927 + 6268.8487559898*t);
   R0 += 0.00000000273*cos(0.30522428863 + 23141.5583829246*t);
   R0 += 0.00000000267*cos(5.76152585786 + 5966.6839803348*t);
   R0 += 0.00000000345*cos(2.94246040875 + 36949.2308084242*t);
   R0 += 0.00000000253*cos(5.20994580359 + 24072.9214697764*t);
   R0 += 0.00000000342*cos(5.76212804329 + 16460.33352952499*t);
   R0 += 0.00000000307*cos(6.01039067183 + 22805.7355659936*t);
   R0 += 0.00000000261*cos(2.00304796059 + 6148.010769956*t);
   R0 += 0.00000000238*cos(5.08241964961 + 6915.8595893046*t);
   R0 += 0.00000000249*cos(2.94762789744 + 135.0650800354*t);
   R0 += 0.00000000306*cos(3.89765478921 + 10988.808157535*t);
   R0 += 0.00000000308*cos(0.05451027736 + 4701.1165017084*t);
   R0 += 0.00000000319*cos(2.95712862064 + 163096.18036118349*t);
   R0 += 0.00000000272*cos(2.07967681309 + 4804.209275927*t);
   R0 += 0.00000000209*cos(4.43768461442 + 6546.1597733642*t);
   R0 += 0.00000000217*cos(0.73691592312 + 6303.8512454838*t);
   R0 += 0.00000000203*cos(0.32033085531 + 25934.1243310894*t);
   R0 += 0.00000000205*cos(5.22936478995 + 20995.3929664494*t);
   R0 += 0.00000000213*cos(0.20671418919 + 28286.9904848612*t);
   R0 += 0.00000000197*cos(0.48286131290 + 16737.5772365966*t);
   R0 += 0.00000000230*cos(6.06567392849 + 6287.0080032545*t);
   R0 += 0.00000000219*cos(1.29194216300 + 5326.7866940208*t);
   R0 += 0.00000000201*cos(1.74700937253 + 22743.4093795164*t);
   R0 += 0.00000000207*cos(4.45440927276 + 6279.4854213396*t);
   R0 += 0.00000000269*cos(6.05640445030 + 64471.99124174489*t);
   R0 += 0.00000000190*cos(0.99261116842 + 29296.6153895786*t);
   R0 += 0.00000000194*cos(3.82656562755 + 419.4846438752*t);
   R0 += 0.00000000262*cos(5.26961924126 + 522.5774180938*t);
   R0 += 0.00000000210*cos(4.68618183158 + 6254.6266625236*t);
   R0 += 0.00000000197*cos(2.80624554186 + 4933.2084403326*t);
   R0 += 0.00000000252*cos(4.36220154620 + 40879.4405046438*t);
   R0 += 0.00000000261*cos(1.07241516738 + 55022.9357470744*t);
   R0 += 0.00000000233*cos(5.41751014958 + 39609.6545831656*t);
   R0 += 0.00000000185*cos(4.14324541379 + 5642.1982426092*t);
   R0 += 0.00000000247*cos(3.44855612987 + 6702.5604938666*t);
   R0 += 0.00000000205*cos(4.04424043226 + 536.8045120954*t);
   R0 += 0.00000000191*cos(3.15807087926 + 16723.350142595*t);
   R0 += 0.00000000222*cos(5.16259496507 + 23539.7073863328*t);
   R0 += 0.00000000180*cos(4.56214752149 + 6489.2613984286*t);
   R0 += 0.00000000227*cos(0.60156339452 + 5905.7022420756*t);
   R0 += 0.00000000170*cos(0.93185903228 + 16062.1845261168*t);
   R0 += 0.00000000159*cos(0.92751013112 + 23937.856389741*t);
   R0 += 0.00000000157*cos(4.69607868164 + 6805.6532680852*t);
   R0 += 0.00000000218*cos(0.85533373430 + 16627.3709153772*t);
   R0 += 0.00000000169*cos(0.94641052064 + 3097.88382272579*t);
   R0 += 0.00000000207*cos(4.88410451334 + 6286.6662786432*t);
   R0 += 0.00000000160*cos(4.95943826819 + 10021.8372800994*t);
   R0 += 0.00000000175*cos(6.12762824563 + 239424.39025435288*t);
   R0 += 0.00000000173*cos(3.13887234973 + 6179.9830757728*t);
   R0 += 0.00000000157*cos(3.62822057807 + 18451.07854656599*t);
   R0 += 0.00000000206*cos(5.74617821138 + 3646.3503773544*t);
   R0 += 0.00000000157*cos(4.67695912207 + 6709.6740408674*t);
   R0 += 0.00000000146*cos(3.09506069745 + 4907.3020501456*t);
   R0 += 0.00000000165*cos(2.27139128760 + 10660.6869350424*t);
   R0 += 0.00000000144*cos(3.96947747592 + 6019.9919266186*t);
   R0 += 0.00000000171*cos(5.91302216729 + 6058.7310542895*t);
   R0 += 0.00000000144*cos(2.13155655120 + 26084.0218062162*t);
   R0 += 0.00000000151*cos(0.67417383565 + 2388.8940204492*t);
   R0 += 0.00000000196*cos(1.67718461229 + 2107.0345075424*t);
   R0 += 0.00000000146*cos(5.10373877968 + 10770.8932562618*t);
   R0 += 0.00000000187*cos(1.23915444627 + 19402.7969528166*t);
   R0 += 0.00000000137*cos(1.26247412216 + 12566.2190102856*t);
   R0 += 0.00000000191*cos(5.03547476279 + 263.0839233728*t);
   R0 += 0.00000000137*cos(3.52825454595 + 639.897286314*t);
   R0 += 0.00000000135*cos(0.73840670927 + 5017.508371365*t);
   R0 += 0.00000000164*cos(2.39195095081 + 6357.8574485587*t);
   R0 += 0.00000000168*cos(0.05515907462 + 9380.9596727172*t);
   R0 += 0.00000000161*cos(1.15721259392 + 26735.9452622132*t);
   R0 += 0.00000000144*cos(1.76097645199 + 5888.4499649322*t);
   R0 += 0.00000000131*cos(2.51859277344 + 6599.467719648*t);
   R0 += 0.00000000142*cos(2.43802911123 + 5881.4037282342*t);
   R0 += 0.00000000159*cos(5.90325893762 + 6281.5913772831*t);
   R0 += 0.00000000151*cos(3.72338532519 + 12669.2444742014*t);
   R0 += 0.00000000132*cos(2.38417741883 + 6525.8044539654*t);
   R0 += 0.00000000127*cos(0.00254936441 + 10027.9031957292*t);
   R0 += 0.00000000148*cos(2.85102145528 + 6418.1409300268*t);
   R0 += 0.00000000143*cos(5.74460279560 + 26087.9031415742*t);
   R0 += 0.00000000172*cos(0.41289962240 + 174242.4659640497*t);
   R0 += 0.00000000136*cos(4.15497742275 + 6311.5250374592*t);
   R0 += 0.00000000170*cos(5.98194913129 + 327574.51427678125*t);
   R0 += 0.00000000136*cos(2.48430537541 + 13341.6743113068*t);
   R0 += 0.00000000149*cos(0.33002271275 + 245.8316462294*t);
   R0 += 0.00000000165*cos(2.49667924600 + 58953.145443294*t);
   R0 += 0.00000000123*cos(1.67328384813 + 32217.2001810808*t);
   R0 += 0.00000000123*cos(3.45660563754 + 6277.552925684*t);
   R0 += 0.00000000117*cos(0.86065134175 + 6245.0481773556*t);
   R0 += 0.00000000149*cos(5.61358281003 + 5729.506447149*t);
   R0 += 0.00000000128*cos(0.71204006448 + 103.0927742186*t);
   R0 += 0.00000000159*cos(2.43166592149 + 221995.02880149524*t);
   R0 += 0.00000000137*cos(1.70657709200 + 12566.08438968*t);
   R0 += 0.00000000129*cos(2.80667872683 + 6016.4688082696*t);
   R0 += 0.00000000113*cos(3.58302904101 + 25685.872802808*t);
   R0 += 0.00000000109*cos(3.26403795962 + 6819.8803620868*t);
   R0 += 0.00000000122*cos(0.34120688202 + 1162.4747044078*t);
   R0 += 0.00000000106*cos(1.59721172719 + 17782.7320727842*t);
   R0 += 0.00000000144*cos(2.28891651774 + 12489.8856287072*t);
   R0 += 0.00000000137*cos(5.82029768354 + 44809.6502008634*t);
   R0 += 0.00000000134*cos(1.26539983018 + 5331.3574437408*t);
   R0 += 0.00000000103*cos(5.96518130595 + 6321.1035226272*t);
   R0 += 0.00000000109*cos(0.33808549034 + 11300.5842213564*t);
   R0 += 0.00000000129*cos(5.89187277190 + 12029.3471878874*t);
   R0 += 0.00000000122*cos(5.77325634636 + 11919.140866668*t);
   R0 += 0.00000000107*cos(6.24998989350 + 77690.75950573849*t);
   R0 += 0.00000000107*cos(1.00535580713 + 77736.78343050249*t);
   R0 += 0.00000000115*cos(5.86963518266 + 12721.572099417*t);
   R0 += 0.00000000102*cos(5.66283467269 + 5540.0857894588*t);
   R0 += 0.00000000143*cos(0.24122178432 + 4214.0690150848*t);
   R0 += 0.00000000143*cos(0.88529649733 + 7576.560073574*t);
   R0 += 0.00000000107*cos(2.92124030351 + 31415.379249957*t);
   R0 += 0.00000000100*cos(5.99485644501 + 4061.2192153944*t);
   R0 += 0.00000000103*cos(2.41941934525 + 5547.1993364596*t);
   R0 += 0.00000000104*cos(4.44106051277 + 2118.7638603784*t);
   R0 += 0.00000000110*cos(0.37559635174 + 5863.5912061162*t);
   R0 += 0.00000000124*cos(2.55619029611 + 12539.853380183*t);
   R0 += 0.00000000110*cos(3.66952094465 + 238004.52415723629*t);
   R0 += 0.00000000112*cos(4.32512422724 + 97238.62754448749*t);
   R0 += 0.00000000120*cos(1.26895630075 + 12043.574281889*t);
   R0 += 0.00000000097*cos(5.42612959752 + 7834.1210726394*t);
   R0 += 0.00000000094*cos(2.56461130309 + 19004.6479494084*t);
   R0 += 0.00000000105*cos(5.68272475301 + 16522.6597160022*t);
   R0 += 0.00000000117*cos(3.65425622684 + 34520.3093093808*t);
   R0 += 0.00000000108*cos(1.24206843948 + 84672.47584450469*t);
   R0 += 0.00000000098*cos(0.13589994287 + 11080.1715789176*t);
   R0 += 0.00000000097*cos(2.46722096722 + 71980.63357473118*t);
   R0 += 0.00000000095*cos(5.36958330451 + 6288.5987742988*t);
   R0 += 0.00000000096*cos(0.20796618776 + 18139.2945014159*t);
   R0 += 0.00000000111*cos(5.01961920313 + 11823.1616394502*t);
   R0 += 0.00000000090*cos(2.72355843779 + 26880.3198130326*t);
   R0 += 0.00000000099*cos(0.90164266199 + 18635.9284545362*t);
   R0 += 0.00000000126*cos(4.78722177847 + 305281.94307104882*t);
   R0 += 0.00000000124*cos(5.00979495566 + 172146.97134054029*t);
   R0 += 0.00000000090*cos(4.50544881196 + 40077.61957352*t);
   R0 += 0.00000000104*cos(5.63679680710 + 2787.0430238574*t);
   R0 += 0.00000000091*cos(5.43564326147 + 6272.0301497275*t);
   R0 += 0.00000000100*cos(2.00639461597 + 12323.4230960088*t);
   R0 += 0.00000000117*cos(2.35555589778 + 83286.91426955358*t);
   R0 += 0.00000000105*cos(2.59824000109 + 30666.1549584328*t);
   R0 += 0.00000000090*cos(2.35779490026 + 12491.3701014155*t);
   R0 += 0.00000000089*cos(3.57152453732 + 11720.0688652316*t);
   R0 += 0.00000000095*cos(5.67015349858 + 14919.0178537546*t);
   R0 += 0.00000000087*cos(1.86043406047 + 27707.5424942948*t);
   R0 += 0.00000000106*cos(3.04150600352 + 22345.2603761082*t);
   R0 += 0.00000000082*cos(5.58298993353 + 10241.2022911672*t);
   R0 += 0.00000000083*cos(3.10607039533 + 36147.4098773004*t);
   R0 += 0.00000000094*cos(5.47749711149 + 9924.8104215106*t);
   R0 += 0.00000000082*cos(4.71988314145 + 15141.390794312*t);
   R0 += 0.00000000096*cos(3.89073946348 + 6379.0550772092*t);
   R0 += 0.00000000110*cos(4.92131611151 + 5621.8429232104*t);
   R0 += 0.00000000110*cos(4.89978492291 + 72140.62866668739*t);
   R0 += 0.00000000097*cos(5.20764563059 + 6303.4311693902*t);
   R0 += 0.00000000085*cos(1.61269222311 + 33326.5787331742*t);
   R0 += 0.00000000093*cos(1.32651591333 + 23020.65308658799*t);
   R0 += 0.00000000090*cos(0.57733016380 + 26482.1708096244*t);
   R0 += 0.00000000078*cos(3.99588630754 + 11293.4706743556*t);
   R0 += 0.00000000106*cos(3.92012705073 + 62883.3551395136*t);
   R0 += 0.00000000098*cos(2.94397773524 + 316.3918696566*t);
   R0 += 0.00000000076*cos(3.96310417608 + 29026.48522950779*t);
   R0 += 0.00000000098*cos(0.95914722366 + 48739.859897083*t);
   R0 += 0.00000000078*cos(1.97068528043 + 90279.92316810328*t);
   R0 += 0.00000000076*cos(0.23027966596 + 21424.4666443034*t);
   R0 += 0.00000000079*cos(1.46227790922 + 8982.810669309*t);
   R0 += 0.00000000078*cos(2.28840998832 + 266.6070417218*t);
   R0 += 0.00000000071*cos(1.51940765590 + 33794.5437235286*t);
   R0 += 0.00000000076*cos(0.22880641443 + 57375.8019008462*t);
   R0 += 0.00000000097*cos(0.39449562097 + 24279.10701821359*t);
   R0 += 0.00000000075*cos(2.77638584795 + 12964.300703391*t);
   R0 += 0.00000000077*cos(5.18846946344 + 11520.9968637952*t);
   R0 += 0.00000000068*cos(0.50006599129 + 4274.5183108324*t);
   R0 += 0.00000000075*cos(2.07323762803 + 15664.03552270859*t);
   R0 += 0.00000000077*cos(0.46665178780 + 16207.886271502*t);
   R0 += 0.00000000081*cos(4.10452219483 + 161710.61878623239*t);
   R0 += 0.00000000071*cos(3.91415328513 + 7875.6718636242*t);
   R0 += 0.00000000081*cos(0.91938383406 + 74.7815985673*t);
   R0 += 0.00000000083*cos(4.69916218791 + 23006.42599258639*t);
   R0 += 0.00000000069*cos(0.98999300277 + 6393.2821712108*t);
   R0 += 0.00000000065*cos(5.41938745446 + 28628.3362260996*t);
   R0 += 0.00000000073*cos(2.45564765251 + 15508.6151232744*t);
   R0 += 0.00000000065*cos(3.02336771694 + 5959.570433334*t);
   R0 += 0.00000000064*cos(0.18375587635 + 1066.49547719*t);
   R0 += 0.00000000080*cos(5.81239171612 + 12341.8069042809*t);
   R0 += 0.00000000066*cos(2.15105504851 + 38.0276726358*t);
   R0 += 0.00000000067*cos(5.14047250153 + 9814.6041002912*t);
   R0 += 0.00000000062*cos(2.43313614978 + 10138.1095169486*t);
   R0 += 0.00000000068*cos(2.24442548639 + 24383.0791084414*t);
   R0 += 0.00000000078*cos(1.39649333997 + 9411.4646150872*t);
   R0 += 0.00000000059*cos(4.95362151577 + 35707.7100829074*t);
   R0 += 0.00000000073*cos(1.35229143121 + 5327.4761083828*t);
   R0 += 0.00000000057*cos(3.16018882154 + 5490.300961524*t);
   R0 += 0.00000000072*cos(5.91833527334 + 10881.0995774812*t);
   R0 += 0.00000000067*cos(0.66414713064 + 29864.334027309*t);
   R0 += 0.00000000065*cos(0.30352816135 + 7018.9523635232*t);
   R0 += 0.00000000059*cos(5.36231868425 + 10239.5838660108*t);
   R0 += 0.00000000056*cos(3.22196331515 + 2636.725472637*t);
   R0 += 0.00000000068*cos(5.32086226658 + 3116.6594122598*t);
   R0 += 0.00000000059*cos(1.63156134967 + 61306.0115970658*t);
   R0 += 0.00000000054*cos(4.29491690425 + 21947.11137270*t);
   R0 += 0.00000000070*cos(0.29271565928 + 6528.9074962208*t);
   R0 += 0.00000000057*cos(5.89190132575 + 34513.2630726828*t);
   R0 += 0.00000000054*cos(2.51856815404 + 6279.1945146334*t);
   R0 += 0.00000000074*cos(1.38235845304 + 9967.4538999816*t);
   R0 += 0.00000000054*cos(0.92276712152 + 6286.9571853494*t);
   R0 += 0.00000000070*cos(5.00933012248 + 6453.7487206106*t);
   R0 += 0.00000000053*cos(3.86543309344 + 32370.9789915656*t);
   R0 += 0.00000000055*cos(4.51794544854 + 34911.412076091*t);
   R0 += 0.00000000063*cos(5.41479412056 + 11502.8376165305*t);
   R0 += 0.00000000063*cos(2.34416220742 + 11510.7019230567*t);
   R0 += 0.00000000056*cos(0.91310629913 + 9910.583327509*t);
   R0 += 0.00000000067*cos(4.03308763854 + 34596.3646546524*t);
   R0 += 0.00000000060*cos(5.57024703495 + 5756.9080032458*t);
   R0 += 0.00000000072*cos(2.80863088166 + 10866.8724834796*t);
   R0 += 0.00000000066*cos(6.12047940728 + 12074.488407524*t);
   R0 += 0.00000000051*cos(2.59519527563 + 11396.5634485742*t);
   R0 += 0.00000000062*cos(5.14746754396 + 25287.7237993998*t);
   R0 += 0.00000000054*cos(2.50994032776 + 5999.2165311262*t);
   R0 += 0.00000000051*cos(4.51195861837 + 29822.7832363242*t);
   R0 += 0.00000000059*cos(0.44167237876 + 250570.67585721909*t);
   R0 += 0.00000000051*cos(3.68849066760 + 6262.7205305926*t);
   R0 += 0.00000000049*cos(0.54704693048 + 22594.05489571199*t);
   R0 += 0.00000000065*cos(2.38423614501 + 52670.0695933026*t);
   R0 += 0.00000000069*cos(5.34363738671 + 66813.5648357332*t);
   R0 += 0.00000000056*cos(2.67216180349 + 17892.93839400359*t);
   R0 += 0.00000000049*cos(4.18361320516 + 18606.4989460002*t);
   R0 += 0.00000000055*cos(0.83886167974 + 20452.8694122218*t);
   R0 += 0.00000000050*cos(1.46327331958 + 37455.7264959744*t);
   R0 += 0.00000000058*cos(3.34847975377 + 33990.6183442862*t);
   R0 += 0.00000000065*cos(1.45522693982 + 76251.32777062019*t);
   R0 += 0.00000000056*cos(2.35650664200 + 37724.7534197482*t);
   R0 += 0.00000000048*cos(1.80689447612 + 206.1855484372*t);
   R0 += 0.00000000056*cos(3.84224878744 + 5483.254724826*t);
   R0 += 0.00000000053*cos(0.17334326094 + 77717.29458646949*t);
   R0 += 0.00000000053*cos(0.79879700631 + 77710.24834977149*t);
   R0 += 0.00000000047*cos(0.43240779709 + 735.8765135318*t);
   R0 += 0.00000000053*cos(4.58786566028 + 11616.976091013*t);
   R0 += 0.00000000048*cos(6.20230111054 + 4171.4255366138*t);
   R0 += 0.00000000052*cos(2.91719053030 + 6993.0088985497*t);
   R0 += 0.00000000057*cos(3.42008310383 + 50317.2034395308*t);
   R0 += 0.00000000048*cos(0.12356889012 + 13362.4497067992*t);
   R0 += 0.00000000060*cos(5.52056066934 + 949.1756089698*t);
   R0 += 0.00000000045*cos(3.37963782356 + 10763.779709261*t);
   R0 += 0.00000000047*cos(5.50958184902 + 12779.4507954208*t);
   R0 += 0.00000000052*cos(5.42770349015 + 310145.15282392364*t);
   R0 += 0.00000000061*cos(2.93237974631 + 5791.4125575326*t);
   R0 += 0.00000000044*cos(2.87440620802 + 8584.6616659008*t);
   R0 += 0.00000000046*cos(4.03141796560 + 10667.8004820432*t);
   R0 += 0.00000000044*cos(1.21579107625 + 6272.4391846416*t);
   R0 += 0.00000000047*cos(2.57670800912 + 11492.542675792*t);
   R0 += 0.00000000044*cos(3.62570223167 + 63658.8777508376*t);
   R0 += 0.00000000051*cos(0.84531181151 + 12345.739057544*t);
   R0 += 0.00000000046*cos(1.17584556517 + 149.5631971346*t);
   R0 += 0.00000000043*cos(0.01524970172 + 37853.8754993826*t);
   R0 += 0.00000000043*cos(0.79038834934 + 640.8776073822*t);
   R0 += 0.00000000044*cos(2.22554419931 + 6293.7125153412*t);
   R0 += 0.00000000049*cos(1.01528394907 + 149144.46708624958*t);
   R0 += 0.00000000041*cos(3.27146326065 + 8858.3149443206*t);
   R0 += 0.00000000045*cos(3.03765521215 + 65236.2212932854*t);
   R0 += 0.00000000058*cos(5.45843180927 + 1975.492545856*t);
   R0 += 0.00000000041*cos(1.32190847146 + 2547.8375382324*t);
   R0 += 0.00000000047*cos(3.67626039848 + 28313.288804661*t);
   R0 += 0.00000000047*cos(6.21438985953 + 10991.3058987006*t);
   R0 += 0.00000000040*cos(2.37237751212 + 8273.8208670324*t);
   R0 += 0.00000000056*cos(1.09773690181 + 77376.20102240759*t);
   R0 += 0.00000000040*cos(2.35698541041 + 2699.7348193176*t);
   R0 += 0.00000000043*cos(5.28030897946 + 17796.9591667858*t);
   R0 += 0.00000000054*cos(2.59175932091 + 22910.44676536859*t);
   R0 += 0.00000000055*cos(0.07988985505 + 83467.15635301729*t);
   R0 += 0.00000000041*cos(4.47510694062 + 5618.3198048614*t);
   R0 += 0.00000000040*cos(1.35670430524 + 27177.8515292002*t);
   R0 += 0.00000000041*cos(2.48011323946 + 6549.6828917132*t);
   R0 += 0.00000000050*cos(2.56387920528 + 82576.98122099529*t);
   R0 += 0.00000000042*cos(4.78798367468 + 7856.89627409019*t);
   R0 += 0.00000000047*cos(2.75482175292 + 18202.21671665939*t);
   R0 += 0.00000000039*cos(1.97008298629 + 24491.4257925834*t);
   R0 += 0.00000000042*cos(4.04346599946 + 7863.9425107882*t);
   R0 += 0.00000000039*cos(3.01033936420 + 853.196381752*t);
   R0 += 0.00000000038*cos(0.49178679251 + 38650.173506199*t);
   R0 += 0.00000000044*cos(1.35931241699 + 21393.5419698576*t);
   R0 += 0.00000000036*cos(4.86047906533 + 4157.1984426122*t);
   R0 += 0.00000000043*cos(5.64354880978 + 1062.9050485382*t);
   R0 += 0.00000000039*cos(3.92736779879 + 3903.9113764198*t);
   R0 += 0.00000000040*cos(5.39694918320 + 9498.2122306346*t);
   R0 += 0.00000000043*cos(2.40863861919 + 29424.634232916*t);
   R0 += 0.00000000046*cos(2.08022244271 + 12573.2652469836*t);
   R0 += 0.00000000050*cos(6.15760345261 + 78051.34191383338*t);
   return R0;
}



   double Earth_R1 (double t) // 290 terms of order 1
{
   double R1 = 0;
   R1 += 0.00103018607*cos(1.10748968172 + 6283.0758499914*t);
   R1 += 0.00001721238*cos(1.06442300386 + 12566.1516999828*t);
   R1 -= 0.00000702217;
   R1 += 0.00000032345*cos(1.02168583254 + 18849.2275499742*t);
   R1 += 0.00000030801*cos(2.84358443952 + 5507.5532386674*t);
   R1 += 0.00000024978*cos(1.31906570344 + 5223.6939198022*t);
   R1 += 0.00000018487*cos(1.42428709076 + 1577.3435424478*t);
   R1 += 0.00000010077*cos(5.91385248388 + 10977.078804699*t);
   R1 += 0.00000008635*cos(0.27158192945 + 5486.777843175*t);
   R1 += 0.00000008654*cos(1.42046854427 + 6275.9623029906*t);
   R1 += 0.00000005069*cos(1.68613408916 + 5088.6288397668*t);
   R1 += 0.00000004985*cos(6.01402338185 + 6286.5989683404*t);
   R1 += 0.00000004667*cos(5.98749245692 + 529.6909650946*t);
   R1 += 0.00000004395*cos(0.51800423445 + 4694.0029547076*t);
   R1 += 0.00000003870*cos(4.74932206877 + 2544.3144198834*t);
   R1 += 0.00000003755*cos(5.07053801166 + 796.2980068164*t);
   R1 += 0.00000004100*cos(1.08424801084 + 9437.762934887*t);
   R1 += 0.00000003518*cos(0.02290216978 + 83996.84731811189*t);
   R1 += 0.00000003436*cos(0.94937503872 + 71430.69561812909*t);
   R1 += 0.00000003221*cos(6.15628775321 + 2146.1654164752*t);
   R1 += 0.00000003418*cos(5.41151581880 + 775.522611324*t);
   R1 += 0.00000002863*cos(5.48433323746 + 10447.3878396044*t);
   R1 += 0.00000002525*cos(0.24296913555 + 398.1490034082*t);
   R1 += 0.00000002205*cos(4.94892172085 + 6812.766815086*t);
   R1 += 0.00000002186*cos(0.41991932164 + 8031.0922630584*t);
   R1 += 0.00000002828*cos(3.41986300734 + 2352.8661537718*t);
   R1 += 0.00000002554*cos(6.13241770582 + 6438.4962494256*t);
   R1 += 0.00000001932*cos(5.31374587091 + 8429.2412664666*t);
   R1 += 0.00000002427*cos(3.09118902115 + 4690.4798363586*t);
   R1 += 0.00000001730*cos(1.53685999718 + 4705.7323075436*t);
   R1 += 0.00000002250*cos(3.68836395620 + 7084.8967811152*t);
   R1 += 0.00000002094*cos(1.28169060400 + 1748.016413067*t);
   R1 += 0.00000001483*cos(3.22226346483 + 7234.794256242*t);
   R1 += 0.00000001434*cos(0.81293662216 + 14143.4952424306*t);
   R1 += 0.00000001754*cos(3.22883705112 + 6279.5527316424*t);
   R1 += 0.00000001583*cos(4.09815978783 + 11499.6562227928*t);
   R1 += 0.00000001575*cos(5.53890314149 + 3154.6870848956*t);
   R1 += 0.00000001847*cos(1.82041234937 + 7632.9432596502*t);
   R1 += 0.00000001499*cos(3.63177937611 + 11513.8833167944*t);
   R1 += 0.00000001337*cos(4.64442556061 + 6836.6452528338*t);
   R1 += 0.00000001275*cos(2.69329661394 + 1349.8674096588*t);
   R1 += 0.00000001348*cos(6.15284035323 + 5746.271337896*t);
   R1 += 0.00000001126*cos(3.35676107739 + 17789.845619785*t);
   R1 += 0.00000001470*cos(3.65282991735 + 1194.4470102246*t);
   R1 += 0.00000001101*cos(4.49747427670 + 4292.3308329504*t);
   R1 += 0.00000001168*cos(2.58033028504 + 13367.9726311066*t);
   R1 += 0.00000001236*cos(5.64980098028 + 5760.4984318976*t);
   R1 += 0.00000000985*cos(0.65326301914 + 5856.4776591154*t);
   R1 += 0.00000000928*cos(2.32555018290 + 10213.285546211*t);
   R1 += 0.00000001073*cos(5.82672338169 + 12036.4607348882*t);
   R1 += 0.00000000918*cos(0.76907130762 + 16730.4636895958*t);
   R1 += 0.00000000876*cos(1.50335727807 + 11926.2544136688*t);
   R1 += 0.00000001023*cos(5.62071200879 + 6256.7775301916*t);
   R1 += 0.00000000853*cos(0.65678134630 + 155.4203994342*t);
   R1 += 0.00000000802*cos(4.10519132094 + 951.7184062506*t);
   R1 += 0.00000000859*cos(1.42880883564 + 5753.3848848968*t);
   R1 += 0.00000000992*cos(1.14238001610 + 1059.3819301892*t);
   R1 += 0.00000000814*cos(1.63584008733 + 6681.2248533996*t);
   R1 += 0.00000000664*cos(4.55039663226 + 5216.5803728014*t);
   R1 += 0.00000000627*cos(1.50782904323 + 5643.1785636774*t);
   R1 += 0.00000000644*cos(4.19480024859 + 6040.3472460174*t);
   R1 += 0.00000000590*cos(6.18371704849 + 4164.311989613*t);
   R1 += 0.00000000635*cos(0.52423584770 + 6290.1893969922*t);
   R1 += 0.00000000650*cos(0.97935492869 + 25132.3033999656*t);
   R1 += 0.00000000568*cos(2.30121525349 + 10973.55568635*t);
   R1 += 0.00000000549*cos(5.26737827342 + 3340.6124266998*t);
   R1 += 0.00000000547*cos(2.20143332641 + 1592.5960136328*t);
   R1 += 0.00000000526*cos(0.92464258271 + 11371.7046897582*t);
   R1 += 0.00000000493*cos(5.91036281399 + 3894.1818295422*t);
   R1 += 0.00000000483*cos(1.66005711540 + 12168.0026965746*t);
   R1 += 0.00000000514*cos(3.59683072524 + 10969.9652576982*t);
   R1 += 0.00000000516*cos(3.97164781773 + 17298.1823273262*t);
   R1 += 0.00000000529*cos(5.03538677680 + 9917.6968745098*t);
   R1 += 0.00000000487*cos(2.50544745305 + 6127.6554505572*t);
   R1 += 0.00000000419*cos(4.05235655996 + 10984.1923516998*t);
   R1 += 0.00000000538*cos(5.54081539813 + 553.5694028424*t);
   R1 += 0.00000000402*cos(2.16859478359 + 7860.4193924392*t);
   R1 += 0.00000000552*cos(2.32219865498 + 11506.7697697936*t);
   R1 += 0.00000000367*cos(3.39145698451 + 6496.3749454294*t);
   R1 += 0.00000000360*cos(5.34467204596 + 7079.3738568078*t);
   R1 += 0.00000000334*cos(3.61346365667 + 11790.6290886588*t);
   R1 += 0.00000000454*cos(0.28755421898 + 801.8209311238*t);
   R1 += 0.00000000419*cos(3.69613970002 + 10575.4066829418*t);
   R1 += 0.00000000319*cos(0.30793759304 + 16200.7727245012*t);
   R1 += 0.00000000376*cos(5.81560210508 + 7058.5984613154*t);
   R1 += 0.00000000364*cos(1.08425056923 + 6309.3741697912*t);
   R1 += 0.00000000294*cos(4.54798604178 + 11856.2186514245*t);
   R1 += 0.00000000290*cos(1.26451946335 + 8635.9420037632*t);
   R1 += 0.00000000394*cos(4.15683669084 + 26.2983197998*t);
   R1 += 0.00000000260*cos(5.09424572996 + 10177.2576795336*t);
   R1 += 0.00000000241*cos(2.25766000302 + 11712.9553182308*t);
   R1 += 0.00000000239*cos(1.06936978753 + 242.728603974*t);
   R1 += 0.00000000276*cos(3.44260568764 + 5884.9268465832*t);
   R1 += 0.00000000255*cos(5.38496803122 + 21228.3920235458*t);
   R1 += 0.00000000307*cos(4.24313885601 + 3738.761430108*t);
   R1 += 0.00000000213*cos(3.44661200485 + 213.299095438*t);
   R1 += 0.00000000198*cos(0.69427265195 + 1990.745017041*t);
   R1 += 0.00000000195*cos(5.16563409007 + 12352.8526045448*t);
   R1 += 0.00000000213*cos(3.89937836808 + 13916.0191096416*t);
   R1 += 0.00000000214*cos(4.00445200772 + 5230.807466803*t);
   R1 += 0.00000000184*cos(5.59805976614 + 6283.14316029419*t);
   R1 += 0.00000000184*cos(2.85275392124 + 7238.67559160*t);
   R1 += 0.00000000179*cos(2.54259058252 + 14314.1681130498*t);
   R1 += 0.00000000236*cos(5.58826125715 + 6069.7767545534*t);
   R1 += 0.00000000189*cos(2.72689937708 + 6062.6632075526*t);
   R1 += 0.00000000184*cos(6.04216273598 + 6283.0085396886*t);
   R1 += 0.00000000225*cos(1.66128561344 + 4732.0306273434*t);
   R1 += 0.00000000230*cos(3.62591335086 + 6284.0561710596*t);
   R1 += 0.00000000172*cos(0.97566476085 + 3930.2096962196*t);
   R1 += 0.00000000162*cos(2.19467339429 + 18073.7049386502*t);
   R1 += 0.00000000215*cos(1.04672844028 + 3496.032826134*t);
   R1 += 0.00000000182*cos(5.17782354566 + 17253.04110768959*t);
   R1 += 0.00000000167*cos(2.17754938066 + 6076.8903015542*t);
   R1 += 0.00000000167*cos(4.75672473773 + 17267.26820169119*t);
   R1 += 0.00000000149*cos(0.80944185798 + 709.9330485583*t);
   R1 += 0.00000000149*cos(0.17584214812 + 9779.1086761254*t);
   R1 += 0.00000000192*cos(5.00680790235 + 11015.1064773348*t);
   R1 += 0.00000000141*cos(4.38420380014 + 4136.9104335162*t);
   R1 += 0.00000000158*cos(4.60969054283 + 9623.6882766912*t);
   R1 += 0.00000000133*cos(3.30507062245 + 154717.60988768269*t);
   R1 += 0.00000000166*cos(6.13191098325 + 3.523118349*t);
   R1 += 0.00000000181*cos(1.60715321141 + 7.1135470008*t);
   R1 += 0.00000000150*cos(5.28136702046 + 13517.8701062334*t);
   R1 += 0.00000000142*cos(0.49788089569 + 25158.6017197654*t);
   R1 += 0.00000000124*cos(6.03440459813 + 9225.539273283*t);
   R1 += 0.00000000124*cos(0.99251562639 + 65147.6197681377*t);
   R1 += 0.00000000128*cos(1.92032744711 + 22483.84857449259*t);
   R1 += 0.00000000124*cos(3.99739675184 + 4686.8894077068*t);
   R1 += 0.00000000121*cos(2.37814805239 + 167283.76158766549*t);
   R1 += 0.00000000123*cos(5.62315112940 + 5642.1982426092*t);
   R1 += 0.00000000117*cos(5.81755956156 + 12569.6748183318*t);
   R1 += 0.00000000157*cos(3.40236948518 + 16496.3613962024*t);
   R1 += 0.00000000130*cos(2.10499918142 + 1589.0728952838*t);
   R1 += 0.00000000116*cos(0.55839966736 + 5849.3641121146*t);
   R1 += 0.00000000123*cos(5.81645568991 + 6282.0955289232*t);
   R1 += 0.00000000110*cos(0.42176497674 + 6172.869528772*t);
   R1 += 0.00000000150*cos(4.26279600865 + 3128.3887650958*t);
   R1 += 0.00000000106*cos(2.27436561182 + 5429.8794682394*t);
   R1 += 0.00000000114*cos(1.52894564202 + 12559.038152982*t);
   R1 += 0.00000000121*cos(0.39459045915 + 12132.439962106*t);
   R1 += 0.00000000104*cos(2.41845930933 + 426.598190876*t);
   R1 += 0.00000000109*cos(5.82786999856 + 16858.4825329332*t);
   R1 += 0.00000000102*cos(4.46626484910 + 23543.23050468179*t);
   R1 += 0.00000000100*cos(2.93812275274 + 4535.0594369244*t);
   R1 += 0.00000000097*cos(3.97935904984 + 6133.5126528568*t);
   R1 += 0.00000000098*cos(0.87616810121 + 6525.8044539654*t);
   R1 += 0.00000000110*cos(6.22339014386 + 12146.6670561076*t);
   R1 += 0.00000000098*cos(3.17344332543 + 10440.2742926036*t);
   R1 += 0.00000000096*cos(2.44128701699 + 3097.88382272579*t);
   R1 += 0.00000000099*cos(5.75642493267 + 7342.4577801806*t);
   R1 += 0.00000000090*cos(0.18984343165 + 13119.72110282519*t);
   R1 += 0.00000000099*cos(5.58884724219 + 2388.8940204492*t);
   R1 += 0.00000000091*cos(6.04278320182 + 20426.571092422*t);
   R1 += 0.00000000080*cos(1.29028142103 + 5650.2921106782*t);
   R1 += 0.00000000086*cos(3.94529200528 + 10454.5013866052*t);
   R1 += 0.00000000085*cos(1.92836879835 + 29088.811415985*t);
   R1 += 0.00000000076*cos(2.70726317966 + 143571.32428481648*t);
   R1 += 0.00000000091*cos(5.63859073351 + 8827.3902698748*t);
   R1 += 0.00000000076*cos(1.80783856698 + 28286.9904848612*t);
   R1 += 0.00000000075*cos(3.40858032804 + 5481.2549188676*t);
   R1 += 0.00000000070*cos(4.53719487231 + 17256.6315363414*t);
   R1 += 0.00000000089*cos(1.10064490942 + 11769.8536931664*t);
   R1 += 0.00000000066*cos(2.78384937771 + 536.8045120954*t);
   R1 += 0.00000000068*cos(3.88199295043 + 17260.1546546904*t);
   R1 += 0.00000000088*cos(3.88075269535 + 7477.522860216*t);
   R1 += 0.00000000061*cos(6.17558202197 + 11087.2851259184*t);
   R1 += 0.00000000060*cos(4.34824715818 + 6206.8097787158*t);
   R1 += 0.00000000082*cos(4.59843208943 + 9388.0059094152*t);
   R1 += 0.00000000079*cos(1.63139280394 + 4933.2084403326*t);
   R1 += 0.00000000081*cos(1.55550779371 + 9380.9596727172*t);
   R1 += 0.00000000078*cos(4.20905757519 + 5729.506447149*t);
   R1 += 0.00000000058*cos(5.76889633224 + 3634.6210245184*t);
   R1 += 0.00000000060*cos(0.93813100594 + 12721.572099417*t);
   R1 += 0.00000000071*cos(6.11408885148 + 8662.240323563*t);
   R1 += 0.00000000057*cos(5.48112524468 + 18319.5365848796*t);
   R1 += 0.00000000070*cos(0.01749174864 + 14945.3161735544*t);
   R1 += 0.00000000074*cos(1.09976045820 + 16460.33352952499*t);
   R1 += 0.00000000056*cos(1.63036186739 + 15720.8387848784*t);
   R1 += 0.00000000055*cos(4.86788348404 + 13095.8426650774*t);
   R1 += 0.00000000060*cos(5.93729841267 + 12539.853380183*t);
   R1 += 0.00000000054*cos(0.22608242982 + 15110.4661198662*t);
   R1 += 0.00000000054*cos(2.30250047594 + 16062.1845261168*t);
   R1 += 0.00000000064*cos(2.13513754101 + 7875.6718636242*t);
   R1 += 0.00000000059*cos(5.87963500139 + 5331.3574437408*t);
   R1 += 0.00000000058*cos(2.30546168615 + 955.5997416086*t);
   R1 += 0.00000000049*cos(1.93839278478 + 5333.9002410216*t);
   R1 += 0.00000000054*cos(5.80331607119 + 12043.574281889*t);
   R1 += 0.00000000054*cos(4.44671053809 + 4701.1165017084*t);
   R1 += 0.00000000049*cos(0.30241161485 + 6805.6532680852*t);
   R1 += 0.00000000046*cos(2.76898193028 + 6709.6740408674*t);
   R1 += 0.00000000046*cos(3.98449608961 + 98068.53671630539*t);
   R1 += 0.00000000049*cos(3.72022009896 + 12323.4230960088*t);
   R1 += 0.00000000045*cos(3.30065998328 + 22003.9146348698*t);
   R1 += 0.00000000048*cos(0.71071357303 + 6303.4311693902*t);
   R1 += 0.00000000061*cos(1.66030429494 + 6262.300454499*t);
   R1 += 0.00000000047*cos(1.26317154881 + 11919.140866668*t);
   R1 += 0.00000000051*cos(1.08020906825 + 10988.808157535*t);
   R1 += 0.00000000045*cos(0.89150445122 + 51868.2486621788*t);
   R1 += 0.00000000043*cos(0.57756724285 + 24356.7807886416*t);
   R1 += 0.00000000043*cos(1.61526242998 + 6277.552925684*t);
   R1 += 0.00000000045*cos(2.96132920534 + 8982.810669309*t);
   R1 += 0.00000000043*cos(5.74295325645 + 11403.676995575*t);
   R1 += 0.00000000055*cos(3.14274403422 + 33019.0211122046*t);
   R1 += 0.00000000057*cos(0.06379726305 + 15671.0817594066*t);
   R1 += 0.00000000041*cos(2.53761820726 + 6262.7205305926*t);
   R1 += 0.00000000040*cos(1.53130436944 + 18451.07854656599*t);
   R1 += 0.00000000052*cos(1.71451922581 + 1551.045222648*t);
   R1 += 0.00000000055*cos(0.89439119424 + 11933.3679606696*t);
   R1 += 0.00000000045*cos(3.88495384656 + 60530.4889857418*t);
   R1 += 0.00000000040*cos(4.75740908001 + 38526.574350872*t);
   R1 += 0.00000000040*cos(3.77498297348 + 26087.9031415742*t);
   R1 += 0.00000000039*cos(2.97113832621 + 2118.7638603784*t);
   R1 += 0.00000000040*cos(3.36050962605 + 10021.8372800994*t);
   R1 += 0.00000000047*cos(1.67051113434 + 6303.8512454838*t);
   R1 += 0.00000000052*cos(5.21827368711 + 77713.7714681205*t);
   R1 += 0.00000000047*cos(4.26356628717 + 21424.4666443034*t);
   R1 += 0.00000000037*cos(1.66712389942 + 6819.8803620868*t);
   R1 += 0.00000000037*cos(0.65746800933 + 12029.3471878874*t);
   R1 += 0.00000000035*cos(3.36255650927 + 24072.9214697764*t);
   R1 += 0.00000000036*cos(0.11087914947 + 10344.2950653858*t);
   R1 += 0.00000000040*cos(4.14725582115 + 2787.0430238574*t);
   R1 += 0.00000000035*cos(5.93650887012 + 31570.7996493912*t);
   R1 += 0.00000000036*cos(2.15108874765 + 30774.5016425748*t);
   R1 += 0.00000000036*cos(1.75078825382 + 16207.886271502*t);
   R1 += 0.00000000034*cos(2.75708224536 + 12139.5535091068*t);
   R1 += 0.00000000034*cos(6.16891378800 + 24491.4257925834*t);
   R1 += 0.00000000034*cos(2.31528650443 + 55798.4583583984*t);
   R1 += 0.00000000032*cos(4.21446357042 + 15664.03552270859*t);
   R1 += 0.00000000034*cos(3.19783054699 + 32217.2001810808*t);
   R1 += 0.00000000039*cos(1.24979117796 + 6418.1409300268*t);
   R1 += 0.00000000038*cos(5.89832942685 + 640.8776073822*t);
   R1 += 0.00000000033*cos(4.80200120107 + 16723.350142595*t);
   R1 += 0.00000000032*cos(1.72442327688 + 27433.88921587499*t);
   R1 += 0.00000000035*cos(4.44608896525 + 18202.21671665939*t);
   R1 += 0.00000000031*cos(4.52790731280 + 6702.5604938666*t);
   R1 += 0.00000000034*cos(3.96287980676 + 18216.443810661*t);
   R1 += 0.00000000030*cos(5.06259854444 + 226858.23855437008*t);
   R1 += 0.00000000034*cos(1.43910280005 + 49515.382508407*t);
   R1 += 0.00000000030*cos(0.29303163371 + 13521.7514415914*t);
   R1 += 0.00000000029*cos(2.02633840220 + 11609.8625440122*t);
   R1 += 0.00000000030*cos(2.54923230240 + 9924.8104215106*t);
   R1 += 0.00000000032*cos(4.91793198558 + 11300.5842213564*t);
   R1 += 0.00000000030*cos(0.23284423547 + 23581.2581773176*t);
   R1 += 0.00000000029*cos(1.62807736495 + 639.897286314*t);
   R1 += 0.00000000028*cos(3.84568936822 + 2699.7348193176*t);
   R1 += 0.00000000029*cos(1.83149729794 + 29822.7832363242*t);
   R1 += 0.00000000033*cos(4.60320094415 + 19004.6479494084*t);
   R1 += 0.00000000027*cos(1.86151121799 + 6288.5987742988*t);
   R1 += 0.00000000030*cos(4.46494072240 + 36147.4098773004*t);
   R1 += 0.00000000028*cos(5.19684492912 + 5863.5912061162*t);
   R1 += 0.00000000035*cos(4.52695674113 + 36949.2308084242*t);
   R1 += 0.00000000027*cos(3.52528177609 + 10770.8932562618*t);
   R1 += 0.00000000026*cos(1.48499438453 + 11080.1715789176*t);
   R1 += 0.00000000035*cos(2.82154380962 + 19402.7969528166*t);
   R1 += 0.00000000025*cos(2.46339998836 + 6279.4854213396*t);
   R1 += 0.00000000026*cos(4.97688894643 + 16737.5772365966*t);
   R1 += 0.00000000027*cos(0.40827112500 + 12964.300703391*t);
   R1 += 0.00000000029*cos(4.15148654061 + 45892.73043315699*t);
   R1 += 0.00000000026*cos(4.56404104286 + 17796.9591667858*t);
   R1 += 0.00000000025*cos(2.89309528854 + 6286.6662786432*t);
   R1 += 0.00000000026*cos(4.82914580957 + 1066.49547719*t);
   R1 += 0.00000000031*cos(3.93096113738 + 29864.334027309*t);
   R1 += 0.00000000024*cos(6.14987193584 + 18606.4989460002*t);
   R1 += 0.00000000024*cos(3.74225964547 + 29026.48522950779*t);
   R1 += 0.00000000025*cos(5.70460621565 + 27707.5424942948*t);
   R1 += 0.00000000025*cos(5.33928840652 + 15141.390794312*t);
   R1 += 0.00000000023*cos(2.37624087345 + 17996.0311682222*t);
   R1 += 0.00000000026*cos(1.34231351782 + 18875.525869774*t);
   R1 += 0.00000000022*cos(5.50791626120 + 6245.0481773556*t);
   R1 += 0.00000000024*cos(1.33998410121 + 19800.9459562248*t);
   R1 += 0.00000000023*cos(0.22512280890 + 6279.7894925736*t);
   R1 += 0.00000000022*cos(1.17576471775 + 11925.2740926006*t);
   R1 += 0.00000000022*cos(3.58603606640 + 6915.8595893046*t);
   R1 += 0.00000000023*cos(3.21621246666 + 6286.3622074092*t);
   R1 += 0.00000000029*cos(2.09564449439 + 15265.8865193004*t);
   R1 += 0.00000000022*cos(4.74660932338 + 28230.18722269139*t);
   R1 += 0.00000000021*cos(2.30688751432 + 5999.2165311262*t);
   R1 += 0.00000000028*cos(3.92087592807 + 18208.349942592*t);
   R1 += 0.00000000021*cos(3.22643339385 + 25934.1243310894*t);
   R1 += 0.00000000021*cos(3.04956726238 + 6566.9351688566*t);
   R1 += 0.00000000027*cos(5.35645770522 + 33794.5437235286*t);
   R1 += 0.00000000025*cos(5.91542362188 + 6489.2613984286*t);
   R1 += 0.00000000020*cos(1.52296293311 + 135.0650800354*t);
   R1 += 0.00000000019*cos(1.78134428631 + 156137.47598479928*t);
   R1 += 0.00000000019*cos(0.34388684087 + 5327.4761083828*t);
   R1 += 0.00000000026*cos(3.41701003233 + 25287.7237993998*t);
   R1 += 0.00000000019*cos(2.86664271911 + 18422.62935909819*t);
   R1 += 0.00000000019*cos(4.71432851499 + 77690.75950573849*t);
   R1 += 0.00000000019*cos(2.54227398241 + 77736.78343050249*t);
   R1 += 0.00000000020*cos(5.91915117116 + 48739.859897083*t);
   return R1*t;
}



   double Earth_R2 (double t) // 134 terms of order 2
{
   double R2 = 0;
   R2 += 0.00004359385*cos(5.78455133808 + 6283.0758499914*t);
   R2 += 0.00000123633*cos(5.57935427994 + 12566.1516999828*t);
   R2 -= 0.00000012342;
   R2 += 0.00000008792*cos(3.62777893099 + 77713.7714681205*t);
   R2 += 0.00000005689*cos(1.86958905084 + 5573.1428014331*t);
   R2 += 0.00000003302*cos(5.47034879713 + 18849.2275499742*t);
   R2 += 0.00000001471*cos(4.47964125007 + 5507.5532386674*t);
   R2 += 0.00000001013*cos(2.81323115556 + 5223.6939198022*t);
   R2 += 0.00000000854*cos(3.10776566900 + 1577.3435424478*t);
   R2 += 0.00000001102*cos(2.84173992403 + 161000.6857376741*t);
   R2 += 0.00000000648*cos(5.47348203398 + 775.522611324*t);
   R2 += 0.00000000608*cos(1.37894173533 + 6438.4962494256*t);
   R2 += 0.00000000499*cos(4.41649242250 + 6286.5989683404*t);
   R2 += 0.00000000416*cos(0.90332697974 + 10977.078804699*t);
   R2 += 0.00000000404*cos(3.20567269530 + 5088.6288397668*t);
   R2 += 0.00000000351*cos(1.81081728907 + 5486.777843175*t);
   R2 += 0.00000000466*cos(3.65086758149 + 7084.8967811152*t);
   R2 += 0.00000000458*cos(5.38585314743 + 149854.40013480789*t);
   R2 += 0.00000000304*cos(3.51015066341 + 796.2980068164*t);
   R2 += 0.00000000266*cos(6.17413982699 + 6836.6452528338*t);
   R2 += 0.00000000281*cos(1.83874672540 + 4694.0029547076*t);
   R2 += 0.00000000262*cos(1.41420110644 + 2146.1654164752*t);
   R2 += 0.00000000264*cos(3.14103683911 + 71430.69561812909*t);
   R2 += 0.00000000319*cos(5.35037932146 + 3154.6870848956*t);
   R2 += 0.00000000238*cos(2.17695432424 + 155.4203994342*t);
   R2 += 0.00000000229*cos(4.75969588070 + 7234.794256242*t);
   R2 += 0.00000000291*cos(4.61776401638 + 4690.4798363586*t);
   R2 += 0.00000000211*cos(0.21864885298 + 4705.7323075436*t);
   R2 += 0.00000000204*cos(4.22895113488 + 1349.8674096588*t);
   R2 += 0.00000000195*cos(4.58550676556 + 529.6909650946*t);
   R2 += 0.00000000255*cos(2.81442711144 + 1748.016413067*t);
   R2 += 0.00000000182*cos(5.70454011389 + 6040.3472460174*t);
   R2 += 0.00000000180*cos(6.02147727878 + 4292.3308329504*t);
   R2 += 0.00000000186*cos(1.58690991244 + 6309.3741697912*t);
   R2 += 0.00000000167*cos(2.88802733052 + 9437.762934887*t);
   R2 += 0.00000000166*cos(1.99990574734 + 8031.0922630584*t);
   R2 += 0.00000000160*cos(0.04412738495 + 2544.3144198834*t);
   R2 += 0.00000000197*cos(2.01089431842 + 1194.4470102246*t);
   R2 += 0.00000000165*cos(5.78372596774 + 83996.84731811189*t);
   R2 += 0.00000000214*cos(3.38300910371 + 7632.9432596502*t);
   R2 += 0.00000000140*cos(0.36669664351 + 10447.3878396044*t);
   R2 += 0.00000000151*cos(0.95519595275 + 6127.6554505572*t);
   R2 += 0.00000000136*cos(1.48417295645 + 2352.8661537718*t);
   R2 += 0.00000000128*cos(5.48057748834 + 951.7184062506*t);
   R2 += 0.00000000126*cos(5.26866506592 + 6279.5527316424*t);
   R2 += 0.00000000127*cos(3.77552907014 + 6812.766815086*t);
   R2 += 0.00000000103*cos(4.95897533789 + 398.1490034082*t);
   R2 += 0.00000000104*cos(0.70183576826 + 1592.5960136328*t);
   R2 += 0.00000000101*cos(1.14481598642 + 3894.1818295422*t);
   R2 += 0.00000000131*cos(0.76624310306 + 553.5694028424*t);
   R2 += 0.00000000109*cos(5.41063597567 + 6256.7775301916*t);
   R2 += 0.00000000078*cos(5.84775340741 + 242.728603974*t);
   R2 += 0.00000000097*cos(1.94685257714 + 11856.2186514245*t);
   R2 += 0.00000000100*cos(5.19725292131 + 244287.60000722769*t);
   R2 += 0.00000000076*cos(0.70480774041 + 8429.2412664666*t);
   R2 += 0.00000000080*cos(6.18430772683 + 1059.3819301892*t);
   R2 += 0.00000000068*cos(5.29561709093 + 14143.4952424306*t);
   R2 += 0.00000000085*cos(5.39487308005 + 25132.3033999656*t);
   R2 += 0.00000000055*cos(5.16874637579 + 7058.5984613154*t);
   R2 += 0.00000000063*cos(0.48494730699 + 801.8209311238*t);
   R2 += 0.00000000058*cos(4.07254840265 + 13367.9726311066*t);
   R2 += 0.00000000051*cos(3.89696552232 + 12036.4607348882*t);
   R2 += 0.00000000051*cos(5.56335232286 + 1990.745017041*t);
   R2 += 0.00000000060*cos(2.25046596710 + 8635.9420037632*t);
   R2 += 0.00000000049*cos(5.58163417371 + 6290.1893969922*t);
   R2 += 0.00000000051*cos(3.87240194908 + 26.2983197998*t);
   R2 += 0.00000000051*cos(4.19300909995 + 7860.4193924392*t);
   R2 += 0.00000000041*cos(3.97169191582 + 10973.55568635*t);
   R2 += 0.00000000041*cos(3.57080919230 + 7079.3738568078*t);
   R2 += 0.00000000056*cos(2.76959005761 + 90955.5516944961*t);
   R2 += 0.00000000042*cos(1.91461189163 + 7477.522860216*t);
   R2 += 0.00000000042*cos(0.42775891995 + 10213.285546211*t);
   R2 += 0.00000000042*cos(1.06925480488 + 709.9330485583*t);
   R2 += 0.00000000038*cos(6.17935925345 + 9917.6968745098*t);
   R2 += 0.00000000050*cos(0.81691517401 + 11506.7697697936*t);
   R2 += 0.00000000053*cos(1.45828359397 + 233141.31440436149*t);
   R2 += 0.00000000038*cos(3.32444534628 + 5643.1785636774*t);
   R2 += 0.00000000047*cos(6.21543665927 + 6681.2248533996*t);
   R2 += 0.00000000037*cos(0.36359309980 + 10177.2576795336*t);
   R2 += 0.00000000045*cos(5.29587706357 + 10575.4066829418*t);
   R2 += 0.00000000034*cos(5.63446915337 + 6525.8044539654*t);
   R2 += 0.00000000034*cos(5.36385158519 + 4933.2084403326*t);
   R2 += 0.00000000035*cos(5.36152295839 + 25158.6017197654*t);
   R2 += 0.00000000042*cos(5.08837645072 + 11015.1064773348*t);
   R2 += 0.00000000042*cos(4.22496037505 + 88860.05707098669*t);
   R2 += 0.00000000039*cos(1.99171699618 + 6284.0561710596*t);
   R2 += 0.00000000029*cos(3.19088628170 + 11926.2544136688*t);
   R2 += 0.00000000029*cos(0.14996158324 + 12168.0026965746*t);
   R2 += 0.00000000030*cos(1.58346276808 + 9779.1086761254*t);
   R2 += 0.00000000026*cos(4.16210340581 + 12569.6748183318*t);
   R2 += 0.00000000036*cos(2.74684637873 + 3738.761430108*t);
   R2 += 0.00000000026*cos(0.72824915320 + 1589.0728952838*t);
   R2 += 0.00000000031*cos(5.34906371821 + 143571.32428481648*t);
   R2 += 0.00000000025*cos(0.10240267494 + 22483.84857449259*t);
   R2 += 0.00000000030*cos(3.47110495524 + 14945.3161735544*t);
   R2 += 0.00000000026*cos(3.89359701125 + 5753.3848848968*t);
   R2 += 0.00000000024*cos(1.18744224678 + 4535.0594369244*t);
   R2 += 0.00000000033*cos(2.99317143244 + 3930.2096962196*t);
   R2 += 0.00000000024*cos(1.57253767584 + 6496.3749454294*t);
   R2 += 0.00000000024*cos(3.47434797542 + 4136.9104335162*t);
   R2 += 0.00000000022*cos(3.91230073719 + 6275.9623029906*t);
   R2 += 0.00000000025*cos(4.02978941287 + 3128.3887650958*t);
   R2 += 0.00000000023*cos(1.07724492065 + 12721.572099417*t);
   R2 += 0.00000000021*cos(1.89591807148 + 16730.4636895958*t);
   R2 += 0.00000000025*cos(2.42198937013 + 5729.506447149*t);
   R2 += 0.00000000020*cos(1.78163489101 + 17789.845619785*t);
   R2 += 0.00000000021*cos(0.49258939822 + 29088.811415985*t);
   R2 += 0.00000000026*cos(4.14947806747 + 2388.8940204492*t);
   R2 += 0.00000000027*cos(2.54785812264 + 3496.032826134*t);
   R2 += 0.00000000020*cos(4.29944129273 + 16858.4825329332*t);
   R2 += 0.00000000021*cos(5.97796936723 + 7.1135470008*t);
   R2 += 0.00000000019*cos(0.80292033311 + 16062.1845261168*t);
   R2 += 0.00000000024*cos(4.89894141052 + 17260.1546546904*t);
   R2 += 0.00000000025*cos(1.37003752175 + 6282.0955289232*t);
   R2 += 0.00000000022*cos(4.92663152168 + 18875.525869774*t);
   R2 += 0.00000000023*cos(5.68902059771 + 16460.33352952499*t);
   R2 += 0.00000000023*cos(3.03021283729 + 66567.48586525429*t);
   R2 += 0.00000000016*cos(3.89713736666 + 5331.3574437408*t);
   R2 += 0.00000000016*cos(5.68562539832 + 12559.038152982*t);
   R2 += 0.00000000016*cos(3.95085099736 + 3097.88382272579*t);
   R2 += 0.00000000016*cos(3.99041783945 + 6283.14316029419*t);
   R2 += 0.00000000020*cos(6.10643919100 + 167283.76158766549*t);
   R2 += 0.00000000015*cos(4.09775914607 + 11712.9553182308*t);
   R2 += 0.00000000016*cos(5.71769940700 + 17298.1823273262*t);
   R2 += 0.00000000016*cos(3.28894009404 + 5884.9268465832*t);
   R2 += 0.00000000015*cos(4.42564243680 + 13517.8701062334*t);
   R2 += 0.00000000016*cos(4.43452080930 + 6283.0085396886*t);
   R2 += 0.00000000014*cos(1.44384279999 + 4164.311989613*t);
   R2 += 0.00000000014*cos(4.47380919159 + 11790.6290886588*t);
   R2 += 0.00000000014*cos(4.77646531825 + 7342.4577801806*t);
   R2 += 0.00000000011*cos(2.56768522896 + 5481.2549188676*t);
   R2 += 0.00000000011*cos(1.51443332200 + 16200.7727245012*t);
   R2 += 0.00000000011*cos(0.88708889185 + 21228.3920235458*t);
   R2 += 0.00000000014*cos(4.50116508534 + 640.8776073822*t);
   return R2*t*t;
}



   double Earth_R3 (double t) // 20 terms of order 3
{
   double R3 = 0;
   R3 += 0.00000144595*cos(4.27319433901 + 6283.0758499914*t);
   R3 += 0.00000006729*cos(3.91706261708 + 12566.1516999828*t);
   R3 += 0.00000000774;
   R3 += 0.00000000247*cos(3.73021571217 + 18849.2275499742*t);
   R3 += 0.00000000036*cos(2.80081409050 + 6286.5989683404*t);
   R3 += 0.00000000033*cos(5.62990083112 + 6127.6554505572*t);
   R3 += 0.00000000018*cos(3.72826142555 + 6438.4962494256*t);
   R3 += 0.00000000016*cos(4.26011484232 + 6525.8044539654*t);
   R3 += 0.00000000014*cos(3.47817116396 + 6256.7775301916*t);
   R3 += 0.00000000012*cos(3.55747379482 + 25132.3033999656*t);
   R3 += 0.00000000010*cos(4.43995693209 + 4705.7323075436*t);
   R3 += 0.00000000010*cos(4.28045255470 + 83996.84731811189*t);
   R3 += 0.00000000009*cos(5.36457057335 + 6040.3472460174*t);
   R3 += 0.00000000008*cos(1.78458957263 + 5507.5532386674*t);
   R3 += 0.00000000009*cos(0.47275199930 + 6279.5527316424*t);
   R3 += 0.00000000009*cos(1.34741231639 + 6309.3741697912*t);
   R3 += 0.00000000009*cos(0.77092900708 + 5729.506447149*t);
   R3 += 0.00000000007*cos(3.50146897332 + 7058.5984613154*t);
   R3 += 0.00000000005*cos(2.89071061700 + 775.522611324*t);
   R3 += 0.00000000006*cos(2.36514111314 + 6836.6452528338*t);
   return R3*t*t*t;
}



   double Earth_R4 (double t) // 9 terms of order 4
{
   double R4 = 0;
   R4 += 0.00000003858*cos(2.56389016346 + 6283.0758499914*t);
   R4 += 0.00000000306*cos(2.26911740541 + 12566.1516999828*t);
   R4 += 0.00000000053*cos(3.44031471924 + 5573.1428014331*t);
   R4 += 0.00000000015*cos(2.03136359366 + 18849.2275499742*t);
   R4 += 0.00000000013*cos(2.05688873673 + 77713.7714681205*t);
   R4 += 0.00000000007*cos(4.41218854480 + 161000.6857376741*t);
   R4 += 0.00000000004*cos(5.33854414781 + 6438.4962494256*t);
   R4 += 0.00000000006*cos(3.81514213664 + 149854.40013480789*t);
   R4 += 0.00000000004*cos(4.26602478239 + 6127.6554505572*t);
   return R4*t*t*t*t;
}



   double Earth_R5 (double t) // 2 terms of order 5
{
   double R5 = 0;
   R5 += 0.00000000086*cos(1.21805304895 + 6283.0758499914*t);
   R5 += 0.00000000012*cos(0.65572878044 + 12566.1516999828*t);
   return R5*t*t*t*t*t;
}

mytncr

天機 发表于 2017-1-15 09:24
/*

   EARTH - VSOP87  

多謝天機壇主給了在下和大家提供了這寶貴的一課。
小弟著實得益良多。很多以前的疑問都一一得以解釋。
這本天文計算書籍 Astronomical Algorithms ，我己買了1998年那個版本。
連同另一本 Practical_Astronomy_with_your_Calculator，我相信己足夠我去學習了。
我希望也能多抽出時間像你一様，逐一按其書中算法，將其整本書的算法全部用計算機程式完成。
同時又可以把我已完成的一些七政四余論命程式，孌得再凖確些。
希望各位有心學習天文計算的朋友。請勿錯過這篇寶貴的文章！