anetarymotionsarestochastic.Sincethepurposeofthispaperistoexhibitandoverviewtheresultsofourlong-termnumericalintegrations，weshowtypicalexamplefiguresasevidenceoftheverylong-termstabilityofSolarsystemplanetarymotion.Forreaderswhohavemorespecificanddeeperinterestsinournumericalresults，wehavepreparedawebpage(access)，whereweshowraworbitalelements，theirlow-passfilteredresults，variationofDelaunayelementsandangularmomentumdeficit，andresultsofoursimpletime–frequencyanalysisonallofourintegrations.

    InSection2webrieflyexplainourdynamicalmodel，numericalmethodandinitialconditionsusedinourintegrations.Section3isdevotedtoadescriptionofthequickresultsofthenumericalintegrations.Verylong-termstabilityofSolarsystemplanetarymotionisapparentbothinplanetarypositionsandorbitalelements.Aroughestimationofnumericalerrorsisalsogiven.Section4goesontoadiscussionofthelongest-termvariationofplanetaryorbitsusingalow-passfilterandincludesadiscussionofangularmomentumdeficit.InSection5，wepresentasetofnumericalintegrationsfortheouterfiveplanetsthatspans±5×1010yr.InSection6wealsodiscussthelong-termstabilityoftheplanetarymotionanditspossiblecause.

    2Descriptionofthenumericalintegrations

    （本部分涉及比较复杂的积分计算，作者君就不贴上来了，贴上来了起点也不一定能成功显示。）

    2.3Numericalmethod

    Weutilizeasecond-orderWisdom–Holmansymplecticmapasourmainintegrationmethod(Wisdom&Holman1991;Kinoshita，Yoshida&Nakai1991)withaspecialstart-upproceduretoreducethetruncationerrorofanglevariables，‘warmstart’(Saha&T

    -->>(第5/18页)（本章未完，请点击下一页继续阅读）