(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 34209, 859] NotebookOptionsPosition[ 32500, 797] NotebookOutlinePosition[ 33131, 822] CellTagsIndexPosition[ 33046, 817] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{"ClearAll", "[", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"Off", "[", StyleBox[ RowBox[{"General", "::", "obspkg"}], "MessageName"], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"SetDirectory", "[", RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"<<", "PhysicalConstants`"}], ";"}]}], "Input", CellChangeTimes->{{3.590755855381109*^9, 3.5907559447380657`*^9}, { 3.590756007949377*^9, 3.5907560326598206`*^9}, {3.591355415798644*^9, 3.5913554168438463`*^9}, {3.5913555040479994`*^9, 3.591355509180408*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"kappa", "=", RowBox[{"GravitationalConstant", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.591355517588823*^9, 3.5913555317380476`*^9}, { 3.5913560759445114`*^9, 3.5913560789085164`*^9}, {3.59135610956257*^9, 3.591356202335933*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"t0", "=", RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[CapitalDelta]t", " ", "=", " ", "0.01"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r0Sl", "=", RowBox[{"AstronomicalData", "[", RowBox[{"\"\\"", ",", " ", RowBox[{"{", RowBox[{"\"\\"", ",", " ", RowBox[{"DateList", "[", "t0", "]"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r0Z", ":=", RowBox[{"AstronomicalData", "[", RowBox[{"\"\\"", ",", " ", RowBox[{"{", RowBox[{"\"\\"", ",", " ", RowBox[{"DateList", "[", "t0", "]"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r1Z", ":=", " ", RowBox[{"AstronomicalData", "[", RowBox[{"\"\\"", ",", " ", RowBox[{"{", RowBox[{"\"\\"", ",", " ", RowBox[{"DateList", "[", RowBox[{"t0", "+", "\[CapitalDelta]t"}], "]"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"v0Z", "=", FractionBox[ RowBox[{"r1Z", "-", "r0Z"}], "\[CapitalDelta]t"]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"v0S", "=", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox["m", "sl"], "=", RowBox[{"AstronomicalData", "[", RowBox[{"\"\\"", ",", " ", "\"\\""}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ SubscriptBox["m", "z"], "=", RowBox[{"AstronomicalData", "[", RowBox[{"\"\\"", ",", " ", "\"\\""}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.590757380316388*^9, 3.5907573927964096`*^9}, { 3.5907575070198107`*^9, 3.5907575081898127`*^9}, {3.590757562243907*^9, 3.5907575674387164`*^9}, {3.5907575980147705`*^9, 3.5907576213924108`*^9}, {3.5907577311698036`*^9, 3.5907577316066046`*^9}, {3.5907579047201085`*^9, 3.590757909010116*^9}, { 3.5907579557165985`*^9, 3.590757972502228*^9}, {3.590758027320724*^9, 3.590758164444965*^9}, {3.5907582861563787`*^9, 3.590758455985077*^9}, 3.590758499353153*^9, {3.5907586141537547`*^9, 3.5907586328113875`*^9}, { 3.590758681951474*^9, 3.5907587371755705`*^9}, {3.590758781011648*^9, 3.590758912098678*^9}, {3.590759201416786*^9, 3.5907594987689085`*^9}, { 3.5907595303121643`*^9, 3.5907596141467113`*^9}, {3.5907597141740866`*^9, 3.5907597385257297`*^9}, 3.5913551789683604`*^9, {3.5913554463746977`*^9, 3.5913554869815693`*^9}, {3.5913556110636945`*^9, 3.591355685304225*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmax", "=", RowBox[{"t0", "+", RowBox[{"365", "*", "24", "*", "3600"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eq1", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["r", "sl"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["r", "z"], "[", "t", "]"}]}], ")"}], "*", FractionBox[ RowBox[{ SubscriptBox["m", "z"], SubscriptBox["m", "sl"]}], SuperscriptBox[ RowBox[{"Norm", "[", RowBox[{ RowBox[{ SubscriptBox["r", "sl"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["r", "z"], "[", "t", "]"}]}], "]"}], "3"]], "*", "kappa"}], "\[Equal]", RowBox[{ SubscriptBox["m", "z"], RowBox[{ RowBox[{ SubscriptBox["r", "z"], "''"}], "[", "t", "]"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eq2", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["r", "z"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["r", "sl"], "[", "t", "]"}]}], ")"}], "*", FractionBox[ RowBox[{ SubscriptBox["m", "z"], SubscriptBox["m", "sl"]}], SuperscriptBox[ RowBox[{"Norm", "[", RowBox[{ RowBox[{ SubscriptBox["r", "z"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["r", "sl"], "[", "t", "]"}]}], "]"}], "3"]], "*", "kappa"}], "\[Equal]", RowBox[{ SubscriptBox["m", "sl"], RowBox[{ RowBox[{ SubscriptBox["r", "sl"], "''"}], "[", "t", "]"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"pocs", "=", RowBox[{"List", "[", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "sl"], "[", "t0", "]"}], "\[Equal]", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], ",", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "sl"], "'"}], "[", "t0", "]"}], "\[Equal]", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], ",", RowBox[{ RowBox[{ SubscriptBox["r", "z"], "[", "t0", "]"}], "==", "r0Z"}], ",", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "z"], "'"}], "[", "t0", "]"}], "==", "v0Z"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"res", "=", RowBox[{ RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"Union", "[", RowBox[{ RowBox[{"{", RowBox[{"eq1", ",", "eq2"}], "}"}], ",", "pocs"}], "]"}], ",", " ", RowBox[{"{", RowBox[{ SubscriptBox["r", "z"], ",", SubscriptBox["r", "sl"]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "t0", ",", " ", "tmax"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.591355187008642*^9, 3.591355227303513*^9}, { 3.591356395339472*^9, 3.5913564084941573`*^9}, {3.591356762739579*^9, 3.591356802129648*^9}, {3.591356961827129*^9, 3.5913571193874054`*^9}, { 3.5913573748066545`*^9, 3.5913574171763287`*^9}, {3.59135746353961*^9, 3.591357503662881*^9}, {3.5913575373121395`*^9, 3.5913575566405735`*^9}, { 3.591358418417287*^9, 3.591358429914508*^9}, {3.5913600542452216`*^9, 3.5913600968020964`*^9}, {3.591360214551103*^9, 3.5913602541595726`*^9}, { 3.5913606468278627`*^9, 3.5913606575918813`*^9}, {3.591360734968017*^9, 3.5913607393360248`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{ SubscriptBox["r", "z"], "\[Rule]", TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"3.5913634196987634`*^9", ",", "3.6228994196987634`*^9"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False]}], ",", RowBox[{ SubscriptBox["r", "sl"], "\[Rule]", TagBox[ RowBox[{"InterpolatingFunction", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{"3.5913634196987634`*^9", ",", "3.6228994196987634`*^9"}], "}"}], "}"}], ",", "\<\"<>\"\>"}], "]"}], False, Editable->False]}]}], "}"}]], "Output", CellChangeTimes->{{3.591357403370304*^9, 3.5913574191263323`*^9}, 3.591357464647212*^9, 3.5913575065332856`*^9, {3.5913575433493505`*^9, 3.5913575571709747`*^9}, 3.5913584310065093`*^9, 3.591360098377699*^9, 3.5913602217271156`*^9, 3.5913602555323753`*^9, 3.5913606588710833`*^9, 3.5913607401940265`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "z"], "[", "t", "]"}], "/.", "res"}], ",", RowBox[{"{", RowBox[{"t", ",", " ", "t0", ",", " ", "tmax"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.5913600282087755`*^9, 3.5913600469600086`*^9}, 3.5913601110137215`*^9, {3.5913606152690067`*^9, 3.5913606200894156`*^9}, { 3.5913606816471233`*^9, 3.591360711801976*^9}, 3.59136079334332*^9}], Cell[BoxData[ Graphics3DBox[{{}, {}, {Hue[0.67, 0.6, 0.6], Line3DBox[CompressedData[" 1:eJwV13c8ld8fAHAre5a9995ku/c5x1bCvbKVRCRRspJRKqF8ZVSIQiqkMlIS 5Zyk7JVSUkIoiURF1O/5/XVfr/freT3n3M96zpELiKAFMTEwMDQwMjAwk78e jx6t1R2nQDbPTIO/R3Whw3Bbw8534qBZ8Wq/eRwFuo8dv82Trgs1Std2GMZL gIGSMYWWCArUTxQQFMvXhU+HfOS5RSXBm28ijpW7KNAjdfCSwn1dOOLnNXqD LgWqPW90qVhQIDU3iY9tSRd+EbZ1OdUpAwpednk/WbGEK/K1ESz79WCNqdrZ XT7yoL6isWffXkto56RQ/piuD0cT6+mc1ipgFpvl8KpZQIMr/sl7JAzh/dSK nqgTuiDT0Ehyu5wprOe5anaqxhg2q7rEpA2ZgfbUbTLRE1thQhcjx+brZuQ+ lcrYmSFg8vRoq/PVh4xPBfchWUvIe3AXumRpC+Ltq7rzDmtCv+feAeXOVMhg ILxjKMIRXAv3L02iKMKxGU2hO6EEbBy8xE+ZcwLF9Q8zFMUV4bp17nrIQQIe qT6AaKtOoO6Dx4G2FQV4ovvqivQhAsorfo3cxbYDhKwNWS5VK0DtAhOxhGgC zso8avdW2AEC8j092SQUIP53mfdPMgGDF8ofTnvvAIt59d9kf8rBdoWB35x5 BKy+EShs/2IH6BFvaW12kYF/J4tcDzYScPvrXZlbipzB+kTHsxpWGSi0EPYh rYmAXfsu6zFUOoPnxnOxl5ulocBLsO9qMwH9hxRokw3OwPKMtsVOVWmoO7jD t6WVgHtlnH6f63MGa1zP6nz/ScLR7Cuvm18QUM/yXIMRowsQjvsx2HtLHC5Y ZSZ2vCZgafVJrGPgAnI2W7q++ysEWcSfHxRaIqDPoY6NmFwXkKE+d1mnQQjK NnqJpv0g4MOB80JGBS6AMh51OeGAEEwN6h5ZWSZgXlOw3+IVF5AUeJCXZUQQ JrCfK8C/CCiSlRLjXuUCWk9Wi3XVbIGR3bt45TYIaGsXeXW21QUMLo4f3b5H ABIf/9v/nRVAsxYDL4V5F6BwzLU3MJEDwt6n7mPCAE4IeMrNAFeQN2t1rmuW HeZKd2/UigDIzbvGz2HnCnSXzfk13Nhhn/+euydEAbT9HrahtN0V3JMrEB5W Z4OXKVp6m8UBvNiO2l3cXYHYN8mXb16zwOArcRW/JAEUldr7fX+oKwh6HzbY pcsAI/faGb6UA9Awaa9fZbYrCK0Ycblh/A8EahWFeMkDWH+9NyT+oivY5fty /hbPXzCdLVL+hvTwPbe9bC+7gonOM27FD/+AoIoxnS4FAPeoKE73lrsC/eNL DvICvwDfgYSLkUrkPt9MppQ9cAWcapt5GJPmgMmlqpODKgByOOrm0965gjVl a/cQ18/AstikgkUVQPukAZVT464AsyCzDsUZ8G6hLFaf9FGLBwa1U65AkWmQ Lb5rAtg6DAufIF1Tg69546sreNmfxtF1ZRyMXTo1XUr6m5a9RlLfXQF/ZkGt ztQ7MIHf+z8mvewGUjVZcQXNNVWKO/ReAoXZ8J1fSH956MBh/w1X8CJCwaPq bh8I9r/Qt0r6P5drzw4y0kB+QWlgzXo7UKu1W2ZRA7A/pTkzZhMNLNbF2e/J bQZn1o4/5SK9k7XmRTwHDZR8i5Nr8L4B0k9aWvKS3iIYFx3PQwOtyVpTjyby kXRT/B5u0leF+XOjBWiAn3154f6VB0jZwFiXlXT921GyYUI00G2QGeub0oZi qiIr/5D72Xn3muQuMRpQFXtu77m1B00wafd8Jb1DtSh1mxQNBGZf8lVeHkS2 tgcuvyHd1WeXn6EcDSgDta7LTq9R/HUtQUz6EJq7KqpEAzz15j9MGcaRBcVt Uyrp/8lEBQ5p0kB8fafaWN0s4hGZm35D5uum4FCfgAW57jPiceyHL6g/iIuz nHSrnUmfh6k0YDPOnc7FPY8evP7xej/pjcGawResaMDoTqXdf3WLSDvbLWpW GcAbPG0ZrNtoIIfZJTmUawVtHharrSXrxFLzrGetFw0cev1pKU/uJ8Jlxppe pFfYfot18SPjaXeg76HxL6TkV0VbVwQwU+Xtjy/+NCCeYln1fe8qOvv+S+lW 0ifO39gQCKGBAGV1k83NG+jA0Svbo8i6VTzTR1WNIde9Y1AWtpUZczF2fpWT BjDmpF9jVDYNSB1m4nAc58LSjHsCqZsBjBfMOnumjQZuBPu/iYjhxjWHtIwd BMh9zrJPDTyngY6Olq2Z3Dw4ZbfPbyd+AM9/vVEh3EUDrPI647dNeLFdB4+d PS+AqcPRkVkDZNwknXiBPT/umDwpzcAJ4GffnnnL96SrnyxVuboZsw0WMyQy ADgmlHcv8ycNlIu/aiqzF8H8+6rM73wmoJrNvpAleTr4zhzsxpwrgo95tZwT niXgMy259DQlOriAxL/5jongvrC2kaPTBCzUz1wTVaWDY+73O9YOi2LL4qq9 epMEVO093q2qRQf0wFeN7wvFsMgefnrgOwKWfz4iymJMBye1i4wN5ySweLsG 39seAkaq0OsrHOiA68+84mMoiz+J3NrEcZuAqdBq//VQOvB5pHp5NVAW/wuq K/lVRcDYP363A8PowFxmdkj3jCxODFf++rGCgMKsR+akw+ngPsv1/TldsnjP /ZyLleUE3CLxjn7qMB1MJGvryLnJYYOaCIdNRQT86qHSuTmODrppURFBgfK4 8lZl1bd0AmbKc69YnqID8WPhR7adUMScL8PqMwPIOb85p68znw5gQarJtzJF fDuJPe3HbgJSrtdcni2gg8h/RhVn2xRx6EsjuZ1+BOR9POvOeJkONi63Bd5n U8KL26XKuTxJr71eoFBMB1WKFf6dmUq4t4mHwdaJgI91Sh00SungY9nxN1P5 yvi4rNi9RSMC0vIL5e0q6GDhTPOrm1Wq+Km/UJACGwFf//kTd7eODlymb5R/ 7lTFztmfWLewELBvlnefaj0dPEhJ2KY0p4pTf5rc/ctAwD/DhRFFpOu0SLFk aKjh+Ed6WnkrVBjFP/Q6+h4dKEWwW7FXq+GnCadvB7yjwkxepj/zDXTw4Tn/ mki3GmYRck9lekWFwSICndT7dJBXX+go91UNW4nJHi3so0IbNqWGc6Szbi3J l9VUx6w5RH0dpkL7F6YfRB7QQYT/INfALXVsMnW64lQFFUbvPZzyg/Rj6uZ1 1V3qOM5+OvlTKRWOb1vVlW4k8zKpXX9iTh0r8DYkUS5T4cELe/9Ykd40ueu6 iIYGDs5NWh3JpMIiy+nJ46QrP/mROuKogTXbvsWJn6FCA+tNTPmkX3/cdi0n VAPrGTUqux2nwq4/WyyqSN/zkWnZJkMDK1QUsqQepUL6jGjOA9KnhT5GLFVq YI0RD8GaSCqs6ZNgaSX9zNY8vvwODTxwqt5j8AAV1j1VynlKejYrvcvoswaO DjvS/zWQCjlFgQUmPZnifrmXXRN3bN99/J8fGYfgo0zNpIcmjsTuVtXE7ZPE fk4PKuzPHp2qIf1koLDXZztNXP1l8hyPCxUOtURPlpDemLHDOCxYE0eICs+z O1Ah9PD4e5Z0lsu1Xb43NDFsDQmdNaNClg+G6XTS6af3RnY+08QGhdpnuw2o 0JfRfVWH9PkqqV+6nzSxpFfvQoUmFfb9Yk9lI134v95jswpaOJ11/ayjNBUm PE5YuUnm5aJtViFDqRb2FJK31mOkwqe+Tjp3yPw+MSqP7S3Qxo82FXn/94wC l0YcVPvJugrUtQtqb9TGWsN7Cn+1UKBd48nwXaRX1EXcv/9aG7d5BVt43afA DBbroRmybj8MrgYeF9LB36l8dzfdpMAOP8bVuVo6eFHAKvPgvA5O+LIjlTmN AsV0PZse3qWDuuw9D4ahLr7k+fRxnj0F1rBMUQ9U0UHvzcHYXXRdrH/6vNMw lQKrSvT/FVbSQddlmZiJvbr4UwWNi8eYAqesA5mfkX301dBO7N0pXfzK0F86 QIkCubfWjrLdJOOcsB51oV0X+74ZnS9jpMBVC5Fh+2vk/6194C/poIclfqK8 E42WsOr1imljIR3s5p5x83bUxxm/dhW3y1hCvdsSV5xPk+9RPVd20cIQV3aO IcFhc9hy5Tx7Ijnfkh0rDNsdDXHvABPzUqc5LNxi+TvYjpwnr5fLFz0NsQBo snjRag5rzt5j22FDB1MnhDcbRRvizIms8t3V5tClICyNC9BBR0R2c8ptQ/wC sGabnTSHO17Iym83IeNpMGMnK7UVj9qHWlzTMYemGY+tpsi5zRJ5T9/t91bc aLCsTz9pBnWCNeaiNmjg6Okms6liY7z36fjHHmFTCNyDLcEVGoieWtwn+s4M h+hW7Bu/tRVWiN3bJnneFWSGBfsQ5y2xXQsl+eq8Hnyh2G8aa+QCmmLy5p8X U7Gvpemil5MmjCs+xSoa6AR2hl5S4rlBYE4Nq8eLVxXg8Aud/xzyHYBIWJ1T WQWBrww4zT9PUoBp84hqmO4Aei7auBrcIrDloJVVoZ8CrKTfuiQa7wCI/JpZ +xoCs547eFqXvB/4tmg69vo4gPRNow6wicD+v+2cFC/Kw5G21yrfpB3A7yHR vQd7CCyiyg51z8rBsTJsylxuD/p7Uy+WLhJ42WmPidqUNJxe+XTlSqUdCNx3 uCFUHeBq4duf+/YLw2OsGt8Tr9mA+ZaI3i2aAGt7UU5PyAjDswI2pVnZNmC8 dmrugRbA/nWJDd+GheC1Oke+y8k2IDjujf53XYA7KzPVl4EQHHIuNbvqYwNO vB2cUTEG2PpQV1uDmCB06f4vyG+LDajJ2lR1HZLH0LGkKuUXApAmtdHRcMIa VPtqtYzvBLhp9fCd25ATzpR+Yd3wsgJ1nZ/efHYHOHSz26+ztzng8fDRvy52 VoARy/yY9wB483T1UIAoB3Ty7tp71dAKcBYniM94AayrJbq8Ns8G56RnMzT4 rEBzxReTYj+AmZ4IZb/J3wTtuSS+MLZBEPBtsUA0EOA35x4rxMwzwCtRT/pv aULwZHt84LaDALeaJa6HO30DT78ml0msEcDsli5rHemKnpsM9jN+BXv/tJvI zxDA2bipf0s4wPrnOK/5NHwG94/aSim+JIB3c+D4C9JD72W3aEh9AurRmdNi dwgw42k49zuCfM+6e7zQm9cg9WinQO8eAhifvfTS5BDA+dGIqCWGwPZR39q7 OwjQ0KU5c4T0kitRe7JcekDNRFLnOXMCaJ80ShslnXD7WHxjrAGUCusmmwoR oO6gQPOhwwDzj8Q+vPeyAzUa7RZ6dJ4KmG1sv6eT/u2g48VnKQPI5xeQNg6l gqf6AS5XSNdt6f2hIfEKnVg+eLjKigr6P554c5v0sJ5tOufq3yLzjn4xISkq +EB/eq6RdNb810yftn1AWp+85GN/UoDdR5vQx6QbdW7PE1/6iIQL1i7091HA exP1yCekR+7fd9uSMoVesNxIlqukgObdOeWPSK/g+qjmkzGNeDe5ftmfQgFc taWM9aQ/UF3YF/lqFukSS+OVPhQwMpaQfp30CxndKEV+DiXop4R/NKSAD18d LHL//7z3k0tnw+eRkeSfYn5eCjiakyLv//91Q5wEIfN3JKeiyePeagkU+oU9 qKQf2z7lbKq2hJh9XQLDCizBiSj4QIz0VzY9JzWcf6Crfh7p8ZGWYIBF8M8T Ms62R0qy/hWuILVQsf1JipagW1Pkw1nSWeukn423/kRWpu9NozYswK9dYQs0 0sXuit55NP0LLcYWsAS8sgC8k68uDJJ5p4h+2+mnv4Zi+cXvKaVZgDO74v5d IutkJvOMtF3HXxQTV16zY84c3FuqiDM/APD6RNjb+2vMmKNwrU+3xgxE7Fa1 aAoF2OFtxJaATBacW9wVcSvNDPgbP9QxJP1LT/EVNtlN+Ffw2ybpPWagd59W hdh+gB/u4uAENqz49nvT5AUBM9DcnjhdtQ/gXSMPu0Ky2LGUrtCq2RFTsN1w YWzaH+BsbYUnIfI8uEchLfqkoQmQ2pBCyqQLTU3+96qKB6sRHXkWfCZg2Dw1 LmA3wJZf91ItDXixD0/m/MJnY8BurxPUQ/YdY4Vn67IVH/6ywXXI5qoxYL3n TDnsDbDfQITXqTp+nMgq11jNYQy8OOyNvtMArvxV+3tocDNOk+t0jhrbCmqq k8rJdsIy5x2uL/KJYJ0z8WUqJwxArOfAJJUf4CErT92iOBl8/8RyHYzRAs/E U1VceQEe/uiQwo9ksNC72KhZeS1g/Yz9mh83wPJDp3+2sMhiq5wmztQ+TfBd U886hB3ge9y73ks4yWKz4V9bKlU1wXHz/dCaAeCD+9+kl4zJ4rj337Sz36oD uwMpPSILBN5v+Ci7bUMOZ/XwGf21UAVGPT6e+h0ENjzPNaFnpojjesI4v6/L g5PBrXp7IwhcKTBlnlCqijtbve163vGDLjad7zraVDx2utSyyVYX99a+iHKq EUbK7ZOcrgKWuGBxSXxtzgBLH/lZvitbGzk0h0nVO5lhlbxSV8ZMYyzkv/zA OdoUNUe5WSCaMX76M0C/Is8ML6gVOHyyB+ifQWpkT4MB1vG5iAbDLLFQ6eq9 8AUbxHrqUFH1LR18dWyPxp0mKr7bozIlqeCI3hzQ5Hj4WwUb2CqwnrhH4BCJ f0nTt5xQWcyXDb1j0vjB3aZcNnK+Nrtxj9BOuqA17aM7RkyksbWhQe8Emc+A Ept0wXQXxDOudSH2pxSWKZzfeOgL8FuZYv/B/1yQfetmq9LDUlhCOcvcl6wj VsXQBFjogmLzPIfKgiVxinu64qFggI/0slgv17ggNa/wXWE0cWywEu45HAUw bWH3UtGYC/rKECOfqiCMWf797ss5B3BKxfrSqoErCv3a4bN/XAjfHy3bmMsk +4Zl64KlqSuCW94l2hUL4TmnZCGQBbCKitvnRIorYv063LQsLITtLlYxfcgG mPfgscHv9q5opiW7V5BTEM99+tX5+SLASxICHpV+rkhR4fqWowsCWD+1/7VZ CcAxVSFy7qmuSF6J71NyOBcm6hS2dtwl978wW3XlpSuydeHbPvOZEyuG5v0r rAF42b1P3/6NK1qRNrprH8SJvbJGh0Nqye+YZVTL/Jgr6vlYMr7qw4HL/9pc /10HcKbk5CvNGVfkl/1EW8yBDaME1a7ZBnLeO8xuxK+6ouC4l4eD5JnxDpee zVyPAC7699vcXYqGJqz/nIovZsL6RV1pxaSf3Tkf80SOhtbKdkacFWXCLGpB 7JrNAFtxVZUrKtPQqdpTRCkvI3505beCZQv5/dcufjmhTUOUC4LqxMoGciLu d+g8AXhUJPmANkFDDOHqQhzlK8i15XeTOwYYtNeUMAbQ0I35RobOgmXUdkP6 SSXpO5U7Lqjto6EacU1aatYP1PDoWNdv0mcEs8OdQmmoueOJ4FL8d7St8+p6 xlOAM2bOPE+LpCHf/t1bTa99QeomUk+S2wDWuVS90ZNCQ5KyX8KT/ptFXF6n mOtJF4l59fNlKg2Vq9E1Hx+dRuIQOU+QHu4Q2/06g4bWq0Z4dVwmyP93+Z/h M/K7apTDMZBDQyaRTAvtq0Po64JkcxXpL80y67LLaEj3+ma6SkkfmvCJoTwm Xe/ei5jYG+T+H5xpfLj9OQpbPPqih3S5t/H8XlU0dDydezhzqQURjMreb0jf UvQ8yegODRF7nghaXapC0g+P/hgnXbx7qJWvjowbx6b7B99mgAGx6Lwp0lPl ng1PNdCQyyshIwume0BxXdTy/y5piJsbHtKQqGudVK4bBgt8wfMfSH+QtBSV 0kJD+XoEG0tFJ2B/t+f6a9I/7knacEQ0NFLYk/ggcQBsf8wd1EW6kEzWTr5n NPTC051v/81hUDZ3QPMR6b6jOxL6XtBQ/xef1ua1McCo7Tv4H+kqu9P0rfpp 6LfipxPamz+CiQ+rNZGkj056t/4coiFHd8rQmvIkCBQOukQjXWImU+zGaxoq bB5STXWZAWxp7cdYSWeo61f//Z6Mj1rf1sKczyAqXPXcKzIv8i0PJj9P0FDG vb60my/ngIjO1sXS/+cxkfv6q2kyDuYh/Ts1F8E3ocx6ddJl2D9ZlH6joS0z sQmKNcug9KyqGgtZJ9mBH8vXNmgoZ9utysr7G8A359DkK7I+s30uxLRI0FGa kyr0fv0XrOxQNfElXZhhzHK3DB1ZvNZ4zPT7H9i+nWPi3WOAn5xp2vgjT0cj UdmixZKMUGV339Z+sv6vsR7cq6xOR62p897FfszwpGbwxySyj7j3hVdTTOgo pMx7/up7Nih2KSJsnezHEcPf+xLodGS71ph5PZwdZjRKKhuTbmh72r7HnY5E nBVR+V92eF7usvHBe2S/UAQExbzpKFHlukOWNCd8H/4+oIvs9/Ud1UGl/nTk +dC6hX8XNwwITBl3JueJ6+rE26iD5D7vKHMOHuSHJr9Vq4/eBPiwupSrciod WfJVsDPNCkLuXPY9tDyA6Yix93Y9Hbn22QnXqAvBK6rtH/hyAe45nn5T8z4d XS+T2+Z1UAiuHnbz7CDnXupS4caNRjqqkz0wX7QkBJmEWE10yTm50cMt9V8L HS39mQqZWReGNxUmH/akk3XoQIgpP6ejR0L+k8f4xeCFvkwu+yTyPHvjUUbC GzriCin3NzGQgm6wydCCPF8XWVPLn6zR0Y5i10GtaCl4/+nPjKMBADc4VDPf Xqej7Mz8U5IPpGAh+7NPdeR3YZOKQN+Fv3TEvzMy/JOZNGSzqawSI88RyGrp pheTG0q2l6hRgDKwZBwlPCDvDfxp+TurONzQ8kt2n7KbsvDGzwtXDG3Jcxb7 HDVCxA05SkQ0vq6Vh72Hic/GimRejuceP6XjhsLVDvFsGZaHse+0vj+TAzj3 04vDf3TdUHT4krPDb3k4bjDq4ywDcEJD11KYvhvKV4iIvUFRgLLVzL5u4uR9 pSizjLrVDaVvtrCV71KAMdQaMRXyPHLUo+j2NTM3JDn6/KLIpCLUYkmSWVwj sBUAEvus3VBiNI/rxiYlKK2bUmb4m7zPUTa9PW3jhvTZjfaMqSrBm/lsuVEr BD7y4Ix8ia0bigmO9jsXrgTfPot4Mk2eOxw3RTu32buh6m3vW8pXleCeH1+8 8qYI3LJ/la9ruxtS7vGfus+rAgNebzYXIu+FMVzJswN0N+TS2ZlrqacCn0pK fOPvJLBHhYt0oZsbejHfzfSYrgI9YnPF2Z+T91dpoze+O92Qwl2TtVv5KnBe eXVitpXAdl1wbMDdDaU98nIwl1eFU4yjdlHkOWGQc0wk3ssNHb3AayVjqAY/ nB4XqS8kMOXOMZWnfm6oLed+pqW7GmzxK8P7LhFYTUAjz2CXG5IYkGPeGacG KyK//RDKI7BoTmXKFdJfpFg5RDSrQbacvPl9mQR+MXpkau9uN8QhLEq4WKtD d8/AnUPJBC7p1XhQ7e+GThw0eG+0Tx1ulD/avjeBwNpnTbS/k96kSvUTTlOH HX0FA9/iCMw+9FtJd48b0sZ+b550qUOCK/nM78ME3vtL9lIx6UvdfG3p8+Tz ivSmmPD/38tp/N2kq+cshm7n04Cp2vZJi6EEVt3azbFCeklg1otHNA2oNWDe 8movgRN5TBONAtwQa6Rpi1akBmQSLP+b7E/gKIkLJ51Ifyo5ZFGQrQFDcquC lPwInHXJ+PFu0t8zfaH+rdGAIv/4f7R7EbjwnYDWQdIFJXme+PZrwPufCm4G uhP4rodwfzTp1Te76+4taMDR++yn/tIIrJ9mcvMo6aLmTVtY+TThKSSXnudM 4Ol/YXVxpKf3h7x10daEJ/WrHihvJ/CPgMrvkaS33wv7l+dEunAIf4M9gb8F TwWGkG6S5XNsMEwT6nvJ51FtCPzxtRC3F+nBe7ttOc5pQslnJVbtgMDu/xnN WJNuM3nxTkCnJrxzmUfthRmBLxVom/KSLn5IyeL0Z014fTnL38qYwNZMG3fn yLglNFz9U8KuBb/T/7Y3GRDY5U2FdxvpEwn9S09ttWBwZr/CFU3y/Pk8Z1sw 6WPr7RJdQVqw2a5PglONwK1n+rP1SO83U/PqOkX68wWbSCUCi3HMcf0i8z55 uPXHfawFxfyHDYylCSxxCF0+Qnrz8+i3phRtmOIvnZnLR+Clt4nqBmRdwdV8 llkDHfg9seaD1iQVG7zl4u72cUOfeYVZEmk6cMnM+e25USoWGCvKdSA94fxU GedhHdizaZk+NUTFl5/VO2NvN/SNd28p3x0dmHLZoiuhjYqndJ8GXSf7Qvao z9yHLbpQzv3aGv06FYu9G+LS93BDp13Vgt+c04Vegbqfi4KoOONBiXuQqxvC Jy66jBzTgy6l4vj6OAXXJt5+ZQjdUOo481mfDD245oSErUYo+MR572fSgJwb 9Uldb/P1oIyGx8E3fRR8LCRXkoUg64Fhm2tPgx5ky72m8v0xBfvnR29Clm6o UV17S8qCHmzv8mHpKaLgQLH+dmZTch7mW2aoBOhDlwsEg4kHBYuYPDW9rU3O DQ1LsM3aANasWgXMtltiXw+h8mxhN/RKtjyO+GMIwZr7gHyOBX4pueNL7Qgd pW+ixUxFmEBN6eLsk0WmuJS93rlNkI7Ch7aJa20zh2Ma1pOr0UZYrT/KeHEL eQ6EN8tzFyyh4sGCfVK39XEh3clVrNYFha6s/3ZkJmCJ0/6NdEctbB+jciBI bAfabOexdfQxAQPjz5hCuhJuCWn+I5/niI4ZTN9KRgR0K3r58JuOEhYr+9Px 66QjcosZSpJpIz3Zw/cCtxKe/NNHPDviiAoL9hm5dhAQnzirPfBMEUunBZs6 0h1RtzXzcMIQAe9qrguxmCjikh7xvnUBRxR96ono8jQBX303Ct0toYAJR2Fv 9v8ckJTnccKPE8DGVOit+14We1Zy97w+bo/kdqf4LkEARde/yrZXi2K7TH0J d5otUpQs2uVsDeD+S2le10JEsbDIsFoFsEVVetPaFTYA3ujgKY9XFMWBprT+ ZV1bZJXjYuJkD+CXv0drRIpE8Mruc8uJfLbozGjBRKATgFfY3kbynBXG3yoU 2V26bZCXhLRppzuAuw3fsfwMFsSbmz8N/LS2QeGvDaxC9wF4S/Eqi440P3YN jfsaqm+NKHzX7IhgAH/eIsbe1fHh1pCzt0TkrFHVqivv5hAAx4r8f9HV+PBp lraUJ3zW6PVng48V+wE0YbY9qSHMi4/K8WUwzFuhsqZredfCAGSYTbt6coEL G/mZRdjftELrm3YvdR8GMOzPDsWVEjbcWfVK74e4FUpK+vlrVySAgov+DOwi bJhb7lAAncMK8Y66P/pKeks8zVQokxVz2h2Zv/MLIm/duPg/RwC0XfhpJ3N0 E07yjFXwfAlRinHJ0KdoALnQx5MrLsx49Qzr6SOZEBk8WnsyHwfg4Jz/w6x/ G6gMifYN/AVIoo2H0+wogPvEVgayFtdRux+tw2oeoI+QxzSFdEW0VSr94x+0 K1Ejs3YUIKVRBRvWeABbVfsdQp6uIs2ZoRPHGwHask9RYIz0pHyGhQ+pK8g+ +lPW18MABQdyJLEfAzBD98D5mthlRDSu2Vn6A3RVOL9Xh3S7Cwp7EkJ+oBsp G5/TdwA05q0gHkk6zpPPYnX8jmLesa8LawAkWP9WMJP0X63tB56ZLaLTi/pR HmIADW2rXbpGunzKzFcGj2+oKTdzOpcNoGy6hccz0mnDJWojYl/QhI5iJ8MU gfT168f6SPfO3XA2+zODzCcyrPQGCfRet9f21f/XvbSv/+K7T2guL/a5XyuB HOYP5YyQ7r5+UP1byyQajWfyTr1DoI9rya3/f56h4PYYcfUjSvps/a+qiEAr E+P9/aSblGQNSTN+QOW7vR51ZRCoXzH+aTvpqkI9RfeevkVbEn1yZuMI9HKr QX4j6aIhUcLWqa8Qx1pwGmMwgS6dWN5xg3QL3hVqr/0g2nynqFR4J4HYQm/O ZJEuGVysVyHcjaYvik6oWBEo3c4oIJr0/iBpcPvyU/S974/TVj0CcQpeXDcg /XhsRb2P1gXkLpDUaM1DPu/s0ZZM5stX9+dLZ9s+4Fhdlx+WQUWDye/225Pu XXSurGR8CORcef+Ato+KaqXqF7lJT1wKuy0oPQIuG6gf0oNU9O2BXFk6WScn mVr9orePA6eEyJ6RVQrq03iV40rWm8XKW2+x57PgXkayYn8IBVEE3K5LkPVp si53NG9sGcj62ReE2Fgi78ETWfkHADw876/xcJoJRtUI1h5bNUP2DXygJxRA s9F7LJeOMkOviF7lPQNmqDc8SmyD7LvYyNzICG4W+Pn5HiVqpRmSTpWj7iD7 9EXJQDCf/ibo+PTH7JSnGQoyWfS/Hwhgr6OzOVsiG0zrPtQy8sAUVWqziF/x BVBsJrqoToAbNswPgANRJsjk1+nUbnLONJ0cjwor3AwdL/NGpM9uRZ8H2Wbr NgPYzfbq7b1qCRj5zIJT9okuuhyg7jDMT84TsVYHpkUJaNlvG2afr4tshliP /+AFUKTNfsbeQBIW7/dHYYd10ZUh90MKXAB6rg4dQA8l4eCyfWilgi7q3nwo yJUZQKszrp1+7VJwP8PksRk+HfRSdGTl0xIBr0n0+2x9LwPLEkZolVc1keVw c9pyLwEvTMr9XWGWh0kxhce3NKkgxeSwuyiegI/F2nZ8k1GGq7sA8G6RQOeL l06lxhLw+FObKGitDC1MPe7WckugtE9qLrZRBLx+VrYjK0QZZn8NHWbwFUfb XES67h8k4KWSqiGhWmVYa/xP9syaKGJbMz7g40+uqzar0wtUoLy29mVTY2FU 2QIl71kT0ENGW216jyp0WxtwD7rLj8YpZwcrOQhIbSvPDSxSh70RKYrbDv9s Hb9Srdi3iYBx3jvfvXiiDjlfxT9Tqf7RutIwl7HASMAEgf55xUl1GBX5xOvv 9GKrqMGqjtUyFR6Z7qa0qWlA013tp00fTrd6JN8uv/CSCntrfgZcua8BZV22 uhfOZ7eaufAcDcqlwtQv8lwPejXhtpxyYurjCrG/XPL6lrOkrx6uy/muCb/o UuyYbv4mfsJyycYUKuwLFo8MEtSCy3b8TaJhfwhfQY9fk4epsN/nKu+ilxaM +2HfInSMAbh4lhQ9d6bCiWLJjPwpLWhuH0fMsrAC5ce1Qap2VJhwugtpsGnD U4qI0SmYDVw9x5aZTKFCbBvU1KCmDbVZWLOqOtiB5w9FBgEtKhzxHsC3w7Vh 4DZjzx3/cYFlw+cu2xSpkN9lulkiWxt+R7EPcxe5wYeBAI4ECSpEw6Hbj9dr w/3h3nr9NF5w5EaY5rXNVLi1/PS+0WFtWFqQTGMme8gf+Da3clAhJ9tuBq3f 2lD/VKvHxbP84H/97Cta "]]}}, AutomaticImageSize->True, Axes->True, ImageSize->{344.2610448066234, 125.03002588122521`}, Method->{}, PlotRange->{{-1.4984563208444427`*^11, 1.501844697298538*^11}, {-1.526690715942017*^11, 1.4740648613902762`*^11}, {-5.34371861282544*^6, 5.169934283554958*^6}}, PlotRangePadding->{ Scaled[0.02], Scaled[0.02], Scaled[0.02]}, ViewPoint->{0.7661617307437147, -3.1833249725374575`, 0.8540716138376191}, ViewVertical->{0.07328844345630683, -0.3690813104598395, 23136.73196513034}]], "Output", CellChangeTimes->{ 3.5913601128077245`*^9, 3.591360225299522*^9, 3.5913602595103817`*^9, 3.5913606212438173`*^9, {3.591360671897106*^9, 3.591360712815978*^9}, 3.591360742830431*^9, 3.5913607942949214`*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"t0", "=", RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmax", "=", RowBox[{"t0", "+", RowBox[{"365", "*", "24", "*", "3600"}]}]}], ";"}]}], "Input", CellChangeTimes->{{3.591361764226625*^9, 3.5913618448787665`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"telesa", "=", RowBox[{"Join", "[", RowBox[{ RowBox[{"{", "\"\\"", "}"}], ",", RowBox[{"AstronomicalData", "[", "\"\\"", "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.591361501646964*^9, 3.5913615405690317`*^9}, { 3.591361582314706*^9, 3.591361584249109*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"Sun\"\>", ",", "\<\"Mercury\"\>", ",", "\<\"Venus\"\>", ",", "\<\"Earth\"\>", ",", "\<\"Mars\"\>", ",", "\<\"Jupiter\"\>", ",", "\<\"Saturn\"\>", ",", "\<\"Uranus\"\>", ",", "\<\"Neptune\"\>"}], "}"}]], "Output", CellChangeTimes->{ 3.591361460135291*^9, {3.5913615183077927`*^9, 3.591361541239833*^9}, 3.591361585309911*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"\[CapitalDelta]t", " ", "=", " ", "0.01"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"hmotnosti", "=", RowBox[{ RowBox[{ RowBox[{"AstronomicalData", "[", RowBox[{"#", ",", " ", "\"\\""}], "]"}], "&"}], "/@", "telesa"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"pocPolohy", "=", RowBox[{ RowBox[{ RowBox[{"AstronomicalData", "[", RowBox[{"#", ",", " ", RowBox[{"{", RowBox[{"\"\\"", ",", " ", RowBox[{"DateList", "[", "t0", "]"}]}], "}"}]}], "]"}], "&"}], "/@", "telesa"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"pocPolohyPlus", "=", RowBox[{ RowBox[{ RowBox[{"AstronomicalData", "[", RowBox[{"#", ",", " ", RowBox[{"{", RowBox[{"\"\\"", ",", " ", RowBox[{"DateList", "[", RowBox[{"t0", "+", "\[CapitalDelta]t"}], "]"}]}], "}"}]}], "]"}], "&"}], "/@", "telesa"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"pocRychlosti", "=", RowBox[{ RowBox[{"(", RowBox[{"pocPolohyPlus", "-", "pocPolohy"}], ")"}], "/", "\[CapitalDelta]t"}]}], ";"}], "\n"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"pocPodmR", "=", RowBox[{"MapThread", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "#1"], "[", "t0", "]"}], "\[Equal]", "#2"}], "&"}], ",", RowBox[{"{", RowBox[{"telesa", ",", "pocPolohy"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"pocPodmV", "=", RowBox[{"MapThread", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{ SubscriptBox["r", "#1"], "'"}], "[", "t0", "]"}], "\[Equal]", "#2"}], "&"}], ",", RowBox[{"{", RowBox[{"telesa", ",", "pocRychlosti"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.591361848825574*^9, 3.591361906592475*^9}, { 3.5913620600655446`*^9, 3.591362307544379*^9}, {3.5913627136910925`*^9, 3.5913627374967346`*^9}, 3.591363054411291*^9, {3.5913632079155607`*^9, 3.5913633285193725`*^9}, {3.5913636693019714`*^9, 3.5913637160864534`*^9}, {3.591363756974125*^9, 3.591363757130125*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Position", "[", RowBox[{"telesa", ",", "\"\\""}], "]"}]], "Input", CellChangeTimes->{{3.59136456227874*^9, 3.591364590483589*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", "1", "}"}], "}"}]], "Output", CellChangeTimes->{{3.591364579267169*^9, 3.5913645908423896`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"hmotnost", "[", "teleso_", "]"}], ":=", RowBox[{ RowBox[{"Cases", "[", RowBox[{ RowBox[{"Transpose", "[", RowBox[{"{", RowBox[{"telesa", ",", "hmotnosti"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"teleso", ",", "_"}], "}"}]}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}]}]], "Input", CellChangeTimes->{{3.591364526149076*^9, 3.591364537412296*^9}, { 3.591364703178187*^9, 3.5913647070313935`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Fna2", "[", RowBox[{"teleso1_", ",", "teleso2_"}], "]"}], ":=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["r", "teleso1"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["r", "teleso2"], "[", "t", "]"}]}], ")"}], "*", FractionBox[ RowBox[{ RowBox[{"hmotnost", "[", "teleso2", "]"}], RowBox[{"hmotnost", "[", "teleso1", "]"}]}], SuperscriptBox[ RowBox[{"Norm", "[", RowBox[{ RowBox[{ SubscriptBox["r", "teleso1"], "[", "t", "]"}], "-", RowBox[{ SubscriptBox["r", "teleso2"], "[", "t", "]"}]}], "]"}], "3"]], "*", "kappa"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Fna2", "[", RowBox[{"teleso_", ",", "teleso_"}], "]"}], ":=", "0"}], ";"}]}], "Input", CellChangeTimes->{{3.591364489723012*^9, 3.591364490783814*^9}, { 3.5913647649698954`*^9, 3.5913648640768695`*^9}, {3.591365019671543*^9, 3.591365033649167*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"??", "Fna2"}]], "Input", CellChangeTimes->{{3.591365038516376*^9, 3.5913650427751837`*^9}}], Cell[CellGroupData[{ Cell["Global`Fna2", "Print", "PrintUsage", CellChangeTimes->{3.591365043820385*^9}, CellTags->"Info3591372243-3521966"], Cell[BoxData[ InterpretationBox[GridBox[{ {GridBox[{ { RowBox[{ RowBox[{"Fna2", "[", RowBox[{"teleso_", ",", "teleso_"}], "]"}], ":=", "0"}]}, {" "}, { RowBox[{ RowBox[{"Fna2", "[", RowBox[{"teleso1_", ",", "teleso2_"}], "]"}], ":=", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"\<\"Subscript\"\>", "[", RowBox[{"r", ",", "teleso1"}], "]"}], "[", "t", "]"}], "-", RowBox[{ RowBox[{"\<\"Subscript\"\>", "[", RowBox[{"r", ",", "teleso2"}], "]"}], "[", "t", "]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"hmotnost", "[", "teleso2", "]"}], " ", RowBox[{"hmotnost", "[", "teleso1", "]"}]}], ")"}], " ", "kappa"}], SuperscriptBox[ RowBox[{"Norm", "[", RowBox[{ RowBox[{ RowBox[{"\<\"Subscript\"\>", "[", RowBox[{"r", ",", "teleso1"}], "]"}], "[", "t", "]"}], "-", RowBox[{ RowBox[{"\<\"Subscript\"\>", "[", RowBox[{"r", ",", "teleso2"}], "]"}], "[", "t", "]"}]}], "]"}], "3"]]}]} }, BaselinePosition->{Baseline, {1, 1}}, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxItemSize->{"Columns" -> {{ Scaled[0.999]}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, "RowsIndexed" -> {}}]} }, BaselinePosition->{Baseline, {1, 1}}, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}], Definition[$CellContext`Fna2], Editable->False]], "Print", CellChangeTimes->{3.5913650438671856`*^9}, CellTags->"Info3591372243-3521966"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Fna2", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]], "Input", CellChangeTimes->{{3.5913648956669254`*^9, 3.591364909036148*^9}}], Cell[BoxData[ FractionBox[ RowBox[{"7.925926899061922`*^44", " ", RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["r", "\<\"Earth\"\>"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["r", "\<\"Sun\"\>"], "[", "t", "]"}]}], ")"}]}], SuperscriptBox[ RowBox[{"Norm", "[", RowBox[{ RowBox[{"-", RowBox[{ SubscriptBox["r", "\<\"Earth\"\>"], "[", "t", "]"}]}], "+", RowBox[{ SubscriptBox["r", "\<\"Sun\"\>"], "[", "t", "]"}]}], "]"}], "3"]]], "Output", CellChangeTimes->{3.59136490976935*^9, 3.591365054506404*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Fna2", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]], "Input", CellChangeTimes->{3.591365061698017*^9}], Cell[BoxData["0"], "Output", CellChangeTimes->{3.5913650626184187`*^9}] }, Open ]] }, WindowSize->{852, 948}, WindowMargins->{{-1, Automatic}, {3, Automatic}}, ShowSelection->True, Magnification:>1.25 Inherited, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{ "Info3591372243-3521966"->{ Cell[29311, 695, 121, 2, 50, "Print", CellTags->"Info3591372243-3521966"], Cell[29435, 699, 1966, 54, 100, "Print", CellTags->"Info3591372243-3521966"]} } *) (*CellTagsIndex CellTagsIndex->{ {"Info3591372243-3521966", 32859, 809} } *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 648, 16, 106, "Input"], Cell[1208, 38, 356, 9, 60, "Input"], Cell[1567, 49, 2764, 68, 310, "Input"], Cell[CellGroupData[{ Cell[4356, 121, 3546, 107, 217, "Input"], Cell[7905, 230, 1050, 29, 68, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8992, 264, 468, 10, 36, "Input"], Cell[9463, 276, 14356, 243, 174, "Output"] }, Open ]], Cell[23834, 522, 310, 8, 60, "Input"], Cell[CellGroupData[{ Cell[24169, 534, 323, 8, 36, "Input"], Cell[24495, 544, 386, 8, 36, "Output"] }, Open ]], Cell[24896, 555, 2357, 69, 267, "Input"], Cell[CellGroupData[{ Cell[27278, 628, 164, 3, 36, "Input"], Cell[27445, 633, 141, 3, 36, "Output"] }, Open ]], Cell[27601, 639, 508, 14, 60, "Input"], Cell[28112, 655, 1033, 30, 107, "Input"], Cell[CellGroupData[{ Cell[29170, 689, 116, 2, 36, "Input"], Cell[CellGroupData[{ Cell[29311, 695, 121, 2, 50, "Print", CellTags->"Info3591372243-3521966"], Cell[29435, 699, 1966, 54, 100, "Print", CellTags->"Info3591372243-3521966"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[31450, 759, 170, 3, 36, "Input"], Cell[31623, 764, 607, 19, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[32267, 788, 142, 3, 36, "Input"], Cell[32412, 793, 72, 1, 36, "Output"] }, Open ]] } ] *) (* End of internal cache information *)