diff --git a/advalgo/exams/Final2017.pdf b/advalgo/exams/Final2017.pdf new file mode 100644 index 0000000..1fc2875 Binary files /dev/null and b/advalgo/exams/Final2017.pdf differ diff --git a/pdes/session03_solutions.pdf b/pdes/session03_solutions.pdf new file mode 100644 index 0000000..7af8c22 Binary files /dev/null and b/pdes/session03_solutions.pdf differ diff --git a/pdes/session04.pdf b/pdes/session04.pdf new file mode 100644 index 0000000..550a026 Binary files /dev/null and b/pdes/session04.pdf differ diff --git a/pdes/session05.pdf b/pdes/session05.pdf new file mode 100644 index 0000000..1297f10 Binary files /dev/null and b/pdes/session05.pdf differ diff --git a/pdes/session06_corr.pdf b/pdes/session06_corr.pdf new file mode 100644 index 0000000..af253a8 Binary files /dev/null and b/pdes/session06_corr.pdf differ diff --git a/pdes/session07_corr.pdf b/pdes/session07_corr.pdf new file mode 100644 index 0000000..aac1e49 Binary files /dev/null and b/pdes/session07_corr.pdf differ diff --git a/pdes/session09.pdf b/pdes/session09.pdf new file mode 100644 index 0000000..a6954b4 Binary files /dev/null and b/pdes/session09.pdf differ diff --git a/pdes/session10.pdf b/pdes/session10.pdf new file mode 100644 index 0000000..bffc115 Binary files /dev/null and b/pdes/session10.pdf differ diff --git a/pdes/session8.pdf b/pdes/session8.pdf new file mode 100644 index 0000000..4cdb912 Binary files /dev/null and b/pdes/session8.pdf differ diff --git a/project/6050a761f5eacc94eb4d3cfc b/project/6050a761f5eacc94eb4d3cfc new file mode 160000 index 0000000..180d862 --- /dev/null +++ b/project/6050a761f5eacc94eb4d3cfc @@ -0,0 +1 @@ +Subproject commit 180d8623470f2b8d062de299e8c37ead41189693 diff --git a/project/MONLI1D b/project/MONLI1D index 2c5a77f..b14e695 160000 --- a/project/MONLI1D +++ b/project/MONLI1D @@ -1 +1 @@ -Subproject commit 2c5a77f780346be4518335204b2b0d35d74618e0 +Subproject commit b14e695f9964f8187320d08c0e3df0fcc2a2adb6 diff --git a/project/MONLI1D_Sem2_Raffaele_Ancarola.pdf b/project/MONLI1D_Sem2_Raffaele_Ancarola.pdf new file mode 100644 index 0000000..f665015 Binary files /dev/null and b/project/MONLI1D_Sem2_Raffaele_Ancarola.pdf differ diff --git a/project/collmom.out b/project/collmom.out new file mode 100644 index 0000000..df7473f --- /dev/null +++ b/project/collmom.out @@ -0,0 +1,199 @@ +p nu1 nu0.5 nu0.1 nu0.01 nu0.001 + 3 0.0003576 0.000419121264713 0.000484 0.0005011 0.0005028 + 4 0.0007381 0.000847284216434 0.0009724 0.001007 0.001011 + 5 0.0007635 0.000859574065544 0.0009927 0.001029 0.001032 + 6 0.0007792 0.000877529557329 0.001004 0.001038 0.001042 + 7 0.0007847 0.000883671105839 0.001002 0.001034 0.001037 + 8 0.0007888 0.000887378526386 0.0009969 0.001027 0.00103 + 9 0.0007984 0.000888742972165 0.0009895 0.001017 0.00102 + 10 0.0008058 0.000887800357305 0.0009794 0.001004 0.001007 + 11 0.0008114 0.000884791836143 0.0009671 0.0009897 0.0009921 + 12 0.000815 0.000879737024661 0.0009527 0.0009727 0.0009748 + 13 0.0008167 0.000872684060596 0.0009361 0.0009536 0.0009554 + 14 0.0008166 0.000863675726578 0.0009175 0.0009324 0.0009339 + 15 0.0008145 0.00085275596939 0.0008968 0.0009091 0.0009104 + 16 0.0008107 0.000840020074975 0.0008746 0.0008842 0.0008852 + 17 0.000805 0.000825622759294 0.0008505 0.0008575 0.0008583 + 18 0.0007976 0.000809451565146 0.0008246 0.000829 0.0008295 + 19 0.0007884 0.000791560159996 0.000797 0.0007988 0.000799 + 20 0.0007776 0.000772004423197 0.0007678 0.000767 0.0007669 + 21 0.0007652 0.000750841922127 0.000737 0.0007336 0.0007333 + 22 0.0007511 0.000728132086806 0.0007048 0.0006988 0.0006982 + 23 0.0007355 0.000703936035279 0.0006711 0.0006626 0.0006617 + 24 0.0007183 0.000678316922858 0.0006362 0.0006251 0.000624 + 25 0.0006996 0.00065133778844 0.0006 0.0005864 0.000585 + 26 0.0006795 0.00062306295149 0.0005626 0.0005467 0.000545 + 27 0.000658 0.000593558186665 0.0005243 0.0005059 0.000504 + 28 0.0006351 0.000562890025321 0.0004849 0.0004642 0.000462 + 29 0.000611 0.000531140598468 0.0004447 0.0004216 0.0004192 + 30 0.0005857 0.000498390756547 0.0004036 0.0003782 0.0003756 + 31 0.0005592 0.000464678945718 0.0003617 0.0003342 0.0003313 + 32 0.0005315 0.00043007437489 0.0003192 0.0002895 0.0002864 + 33 0.0005029 0.000394646660425 0.0002761 0.0002443 0.000241 + 34 0.0004732 0.000358465767931 0.0002325 0.0001986 0.0001951 + 35 0.0004426 0.000321603787597 0.0001885 0.0001526 0.0001489 + 36 0.0004112 0.000284130888758 0.0001441 0.0001872 0.0001929 + 37 0.000379 0.000246117124334 0.0001756 0.0002318 0.0002377 + 38 0.000346 0.000207634016988 0.0002189 0.0002772 0.0002834 + 39 0.0003123 0.000168753438629 0.0002623 0.0003223 0.0003286 + 40 0.0002781 0.000129547319375 0.0003061 0.0003671 0.0003735 + 41 0.0002434 0.000125934879179 0.0003496 0.0004113 0.0004177 + 42 0.0002082 0.000165540535818 0.0003924 0.0004547 0.0004611 + 43 0.0001726 0.000205097705475 0.0004345 0.0004972 0.0005037 + 44 0.0001366 0.000244476745138 0.0004758 0.0005389 0.0005454 + 45 0.0001004 0.000283530942397 0.0005163 0.0005795 0.000586 + 46 0.0001125 0.000322174921166 0.0005558 0.000619 0.0006255 + 47 0.0001507 0.000360323349014 0.0005943 0.0006573 0.0006638 + 48 0.0001889 0.000397891213652 0.0006316 0.0006944 0.0007009 + 49 0.000227 0.000434903573478 0.0006678 0.0007301 0.0007365 + 50 0.0002648 0.000471616862342 0.0007026 0.0007643 0.0007707 + 51 0.0003024 0.000507531105541 0.0007361 0.000797 0.0008034 + 52 0.0003394 0.000542523281183 0.0007681 0.0008282 0.0008344 + 53 0.000376 0.00057650799863 0.0007985 0.0008576 0.0008637 + 54 0.000412 0.000609399110544 0.0008273 0.0008853 0.0008913 + 55 0.0004472 0.000641108956188 0.0008545 0.0009111 0.000917 + 56 0.0004816 0.000671551621053 0.0008798 0.0009351 0.0009408 + 57 0.0005152 0.000700641656294 0.0009034 0.0009571 0.0009627 + 58 0.0005477 0.000728294893634 0.0009251 0.0009771 0.0009825 + 59 0.0005791 0.0007544279797 0.0009447 0.0009949 0.001 + 60 0.0006093 0.000778958958108 0.0009624 0.001011 0.001016 + 61 0.0006382 0.000801807793323 0.0009781 0.001024 0.001029 + 62 0.0006657 0.000822900678031 0.0009915 0.001035 0.00104 + 63 0.0006918 0.000842202163767 0.001003 0.001044 0.001049 + 64 0.0007162 0.000859598279931 0.001012 0.001051 0.001055 + 65 0.000739 0.000875019351952 0.001018 0.001055 0.001059 + 66 0.00076 0.000888399430551 0.001022 0.001056 0.00106 + 67 0.0007792 0.000899677514099 0.001024 0.001055 0.001059 + 68 0.0007965 0.000908796035219 0.001023 0.001052 0.001055 + 69 0.0008118 0.000915701850317 0.001019 0.001045 0.001048 + 70 0.0008251 0.000920347461943 0.001013 0.001037 0.001039 + 71 0.0008362 0.000922691193409 0.001005 0.001028 0.001032 + 72 0.0008461 0.000922697770875 0.0009937 0.00103 0.001033 + 73 0.0008541 0.000920337974094 0.0009939 0.001029 0.001032 + 74 0.0008599 0.00091558980057 0.0009926 0.001025 0.001028 + 75 0.0008635 0.000908438232727 0.0009888 0.001019 0.001022 + 76 0.0008648 0.000898876285646 0.0009826 0.00101 0.001013 + 77 0.0008639 0.000886904832441 0.0009739 0.000999 0.001002 + 78 0.0008607 0.000881054031197 0.000963 0.0009855 0.0009878 + 79 0.0008552 0.000877235666849 0.0009495 0.0009693 0.0009714 + 80 0.0008473 0.000871325377375 0.0009336 0.0009506 0.0009524 + 81 0.0008372 0.00086314312648 0.0009151 0.0009294 0.0009309 + 82 0.0008247 0.000852683035191 0.0008942 0.0009056 0.0009068 + 83 0.0008128 0.00083995051682 0.0008708 0.0008793 0.0008802 + 84 0.0008074 0.000824940914754 0.000845 0.0008505 0.0008511 + 85 0.0008 0.000807667034678 0.0008168 0.0008194 0.0008197 + 86 0.0007903 0.000788150005974 0.0007863 0.0007859 0.0007858 + 87 0.0007785 0.000766418525018 0.0007535 0.0007501 0.0007497 + 88 0.0007646 0.000742509844713 0.0007185 0.0007121 0.0007114 + 89 0.0007485 0.000716469599865 0.0006813 0.0006719 0.000671 + 90 0.0007304 0.000688352622092 0.0006422 0.0006298 0.0006285 + 91 0.0007101 0.000658222008497 0.0006011 0.0005858 0.0005842 + 92 0.0006878 0.000626149587333 0.0005582 0.0005399 0.000538 + 93 0.0006636 0.000592216209043 0.0005136 0.0004924 0.0004902 + 94 0.0006374 0.000556512561161 0.0004674 0.0004434 0.0004409 + 95 0.0006094 0.000519133347552 0.0004198 0.000393 0.0003902 + 96 0.0005795 0.000480183400214 0.0003708 0.0003414 0.0003383 + 97 0.000548 0.000439774506958 0.0003208 0.0002887 0.0002854 + 98 0.0005148 0.000398025673348 0.0002697 0.0002351 0.0002315 + 99 0.0004802 0.000355098512955 0.0002177 0.0001807 0.0001769 +100 0.0004441 0.000311105803121 0.0001651 0.0002216 0.0002278 +101 0.0004066 0.000266159244347 0.0002124 0.0002728 0.0002791 +102 0.000368 0.000220398447709 0.0002623 0.0003264 0.0003331 +103 0.0003283 0.00017396689509 0.0003142 0.0003793 0.0003861 +104 0.0002877 0.000127011095174 0.0003652 0.0004313 0.0004381 +105 0.0002463 0.00016911138664 0.0004151 0.0004818 0.0004887 +106 0.0002041 0.000215410880628 0.0004637 0.0005309 0.0005379 +107 0.0001614 0.000260993867414 0.0005108 0.0005783 0.0005854 +108 0.0001182 0.000305612687953 0.0005563 0.000624 0.000631 +109 0.000121 0.000349166890373 0.0006001 0.0006677 0.0006748 +110 0.0001631 0.0003915122943 0.000642 0.0007095 0.0007165 +111 0.0002046 0.000432536762673 0.000682 0.000749 0.000756 +112 0.0002452 0.000472135201562 0.0007199 0.0007864 0.0007933 +113 0.0002849 0.000510209705681 0.0007556 0.0008214 0.0008283 +114 0.0003236 0.000546670577023 0.0007892 0.0008542 0.0008609 +115 0.0003611 0.000581480795518 0.0008206 0.0008845 0.0008911 +116 0.0003974 0.000614544376731 0.0008495 0.0009123 0.0009188 +117 0.0004324 0.000645776628517 0.0008761 0.0009375 0.0009439 +118 0.000466 0.000675118702929 0.0009003 0.0009602 0.0009664 +119 0.0004982 0.000702519551851 0.0009219 0.0009802 0.0009863 +120 0.0005288 0.00072793540312 0.0009411 0.0009976 0.001003 +121 0.000558 0.000751329876948 0.0009577 0.001012 0.001018 +122 0.0005855 0.000772673694883 0.0009718 0.001024 0.00103 +123 0.0006113 0.000791948870756 0.0009834 0.001034 0.001039 +124 0.0006355 0.000809136836324 0.0009925 0.001041 0.001046 +125 0.0006579 0.000824225717224 0.0009991 0.001045 0.00105 +126 0.0006786 0.000837211555336 0.001003 0.001047 0.001051 +127 0.0006976 0.000848095514812 0.001005 0.001046 0.00105 +128 0.0007149 0.000856884173118 0.001004 0.001042 0.001046 +129 0.0007303 0.000863644294441 0.001001 0.001036 0.00104 +130 0.0003946 0.000453585438663 0.0005113 0.0005346 0.0005372 +131 3.004e-05 2.7630423574e-05 2.952e-05 3.346e-05 3.389e-05 +132 1.073e-05 1.21954917631e-05 1.427e-05 1.495e-05 1.502e-05 +133 1.365e-06 1.55092220666e-06 1.976e-06 2.143e-06 2.162e-06 +134 5.324e-07 7.38075073059e-07 1.077e-06 1.198e-06 1.211e-06 +135 2.076e-07 2.91000787911e-07 4.409e-07 4.943e-07 5.001e-07 +136 3.209e-08 4.57999327352e-08 8.088e-08 9.269e-08 9.398e-08 +137 1.841e-08 2.85304935232e-08 4.535e-08 5.12e-08 5.184e-08 +138 4.618e-09 8.02485544682e-09 1.452e-08 1.695e-08 1.722e-08 +139 1.358e-09 2.53870102718e-09 4.771e-09 5.596e-09 5.688e-09 +140 5.39e-10 9.65674873399e-10 1.752e-09 2.045e-09 2.077e-09 +141 1.177e-10 2.27674157305e-10 4.593e-10 5.522e-10 5.627e-10 +142 5.473e-11 1.15327553063e-10 2.434e-10 2.941e-10 2.999e-10 +143 1.457e-11 3.20400789233e-11 7.076e-11 8.651e-11 8.831e-11 +144 4.326e-12 1.11098109878e-11 2.734e-11 3.415e-11 3.493e-11 +145 2.036e-12 4.77614735955e-12 1.128e-11 1.404e-11 1.436e-11 +146 4.891e-13 1.20253123865e-12 3.026e-12 3.833e-12 3.927e-12 +147 2.201e-13 6.27210739942e-13 1.753e-12 2.261e-12 2.321e-12 +148 6.153e-14 1.84986466249e-13 5.412e-13 7.06e-13 7.254e-13 +149 1.913e-14 7.10181154336e-14 2.321e-13 3.074e-13 3.162e-13 +150 9.709e-15 2.96852115372e-14 8.512e-14 1.086e-13 1.112e-13 +151 2.48e-15 7.7984146579e-15 2.364e-14 3.317e-14 3.433e-14 +152 1.147e-15 4.44754638933e-15 1.501e-14 1.939e-14 1.986e-14 +153 3.27e-16 1.33335437282e-15 4.346e-15 5.377e-15 5.487e-15 +154 1.11e-16 5.74616669385e-16 2.035e-15 2.683e-15 2.778e-15 +155 5.644e-17 2.16915476322e-16 6.214e-16 8.821e-16 9.217e-16 +156 1.499e-17 5.78894476936e-17 2.276e-16 3.961e-16 4.317e-16 +157 7.373e-18 3.71387220552e-17 9.741e-17 2.125e-16 2.293e-16 +158 2.164e-18 1.17506779734e-17 3.454e-17 1.401e-16 1.694e-16 +159 8.727e-19 6.0505073027e-18 2.881e-17 1.133e-16 1.242e-16 +160 4.027e-19 2.37846562076e-18 2.238e-17 8.549e-17 1.035e-16 +161 1.616e-19 1.28519816914e-18 1.826e-17 7.115e-17 7.88e-17 +162 9.207e-20 8.52087382689e-19 1.586e-17 5.521e-17 7.062e-17 +163 7.099e-20 6.43470571046e-19 1.249e-17 4.601e-17 5.303e-17 +164 4.892e-20 6.95572403084e-19 1.225e-17 3.999e-17 5.641e-17 +165 4.608e-20 5.37860596531e-19 1.157e-17 3.896e-17 4.41e-17 +166 3.474e-20 5.82497743926e-19 1.006e-17 3.049e-17 4.986e-17 +167 2.838e-20 4.80045738195e-19 1.111e-17 3.409e-17 4.191e-17 +168 2.475e-20 4.5577382918e-19 9.511e-18 2.41e-17 4.625e-17 +169 1.713e-20 4.11059366247e-19 1.059e-17 3.042e-17 4.164e-17 +170 1.678e-20 3.36413217752e-19 9.345e-18 2.086e-17 4.34e-17 +171 1.017e-20 3.47723723557e-19 1.001e-17 2.741e-17 4.131e-17 +172 1.089e-20 2.4585431147e-19 9.064e-18 1.868e-17 4.19e-17 +173 5.967e-21 2.85573301437e-19 9.511e-18 2.559e-17 4.068e-17 +174 6.847e-21 1.87217700894e-19 8.618e-18 1.856e-17 4.076e-17 +175 3.517e-21 2.33626979354e-19 9.045e-18 2.44e-17 3.93e-17 +176 4.221e-21 1.41841054031e-19 8.013e-18 1.951e-17 4e-17 +177 2.07e-21 1.8733821077e-19 8.373e-18 2.313e-17 3.713e-17 +178 2.552e-21 1.11233885027e-19 7.344e-18 2.131e-17 3.982e-17 +179 1.235e-21 1.46508924701e-19 7.512e-18 2.185e-17 3.435e-17 +180 1.509e-21 8.67694432616e-20 6.703e-18 2.265e-17 3.943e-17 +181 7.895e-22 1.11851905267e-19 6.594e-18 2.082e-17 3.135e-17 +182 8.669e-22 6.5751255648e-20 6.109e-18 2.33e-17 3.865e-17 +183 5.583e-22 8.44253555773e-20 5.847e-18 2.039e-17 2.875e-17 +184 4.789e-22 4.95787114499e-20 5.507e-18 2.34e-17 3.739e-17 +185 4.073e-22 6.25702938284e-20 5.478e-18 2.015e-17 2.665e-17 +186 2.487e-22 3.78170206265e-20 4.893e-18 2.313e-17 3.53e-17 +187 3.007e-22 4.52877572587e-20 5.117e-18 2.005e-17 2.475e-17 +188 1.385e-22 2.90956030639e-20 4.274e-18 2.236e-17 3.264e-17 +189 2.188e-22 3.3770104103e-20 4.763e-18 1.969e-17 2.328e-17 +190 8.693e-23 2.29024711216e-20 3.661e-18 2.127e-17 2.964e-17 +191 1.554e-22 2.54138625471e-20 4.405e-18 1.923e-17 2.188e-17 +192 5.721e-23 1.88968990885e-20 3.148e-18 1.999e-17 2.623e-17 +193 1.085e-22 2.00017781267e-20 4.139e-18 1.879e-17 2.072e-17 +194 3.891e-23 1.55302675331e-20 2.942e-18 1.867e-17 2.284e-17 +195 7.463e-23 1.65146189244e-20 4.041e-18 1.819e-17 2.003e-17 +196 2.644e-23 1.26501280377e-20 2.896e-18 1.726e-17 1.966e-17 +197 5.054e-23 1.51047735808e-20 3.976e-18 1.744e-17 1.95e-17 +198 1.833e-23 1.03050928116e-20 2.823e-18 1.606e-17 1.703e-17 +199 3.377e-23 1.39956245458e-20 3.896e-18 1.688e-17 1.913e-17 +200 1.229e-23 8.15328180663e-21 2.706e-18 1.53e-17 1.542e-17 diff --git a/project/collmom_new.out b/project/collmom_new.out new file mode 100644 index 0000000..4ec7db2 --- /dev/null +++ b/project/collmom_new.out @@ -0,0 +1,149 @@ +p nu5e-1 nu1e-1 nu1e-2 nu1e-3 + 3 8.318e-05 0.0001011 0.0001058 0.0001063 + 4 0.0001667 0.0002027 0.0002121 0.0002131 + 5 0.0001728 0.0002119 0.0002223 0.0002234 + 6 0.000178 0.00022 0.0002312 0.0002324 + 7 0.0001786 0.000221 0.0002324 0.0002336 + 8 0.0001786 0.0002214 0.0002329 0.0002341 + 9 0.000178 0.000221 0.0002325 0.0002337 + 10 0.000177 0.0002199 0.0002315 0.0002327 + 11 0.0001755 0.0002186 0.0002302 0.0002314 + 12 0.0001738 0.0002167 0.0002284 0.0002296 + 13 0.0001716 0.0002144 0.000226 0.0002272 + 14 0.0001691 0.0002115 0.000223 0.0002242 + 15 0.0001661 0.0002082 0.0002196 0.0002208 + 16 0.0001628 0.0002043 0.0002156 0.0002168 + 17 0.0001591 0.0002 0.0002111 0.0002123 + 18 0.000155 0.0001952 0.0002061 0.0002073 + 19 0.0001505 0.0001899 0.0002007 0.0002018 + 20 0.0001457 0.0001842 0.0001947 0.0001958 + 21 0.0001406 0.0001781 0.0001884 0.0001894 + 22 0.0001351 0.0001716 0.0001816 0.0001826 + 23 0.0001293 0.0001647 0.0001744 0.0001754 + 24 0.0001233 0.0001574 0.0001668 0.0001677 + 25 0.0001169 0.0001497 0.0001588 0.0001597 + 26 0.0001103 0.0001418 0.0001504 0.0001513 + 27 0.0001034 0.0001334 0.0001417 0.0001426 + 28 9.635e-05 0.0001248 0.0001327 0.0001336 + 29 8.903e-05 0.0001159 0.0001234 0.0001242 + 30 8.151e-05 0.0001068 0.0001138 0.0001146 + 31 7.38e-05 9.737e-05 0.000104 0.0001047 + 32 6.592e-05 8.774e-05 9.388e-05 9.452e-05 + 33 5.79e-05 7.792e-05 8.357e-05 8.417e-05 + 34 4.974e-05 6.791e-05 7.308e-05 7.362e-05 + 35 4.147e-05 5.776e-05 6.242e-05 6.291e-05 + 36 4.916e-05 6.698e-05 7.204e-05 7.257e-05 + 37 5.674e-05 7.606e-05 8.15e-05 8.207e-05 + 38 6.493e-05 8.609e-05 9.202e-05 9.264e-05 + 39 7.297e-05 9.593e-05 0.0001023 0.000103 + 40 8.089e-05 0.0001056 0.0001125 0.0001132 + 41 8.863e-05 0.0001151 0.0001225 0.0001232 + 42 9.618e-05 0.0001243 0.0001321 0.0001329 + 43 0.0001035 0.0001333 0.0001415 0.0001423 + 44 0.0001106 0.0001419 0.0001505 0.0001514 + 45 0.0001174 0.0001502 0.0001592 0.0001602 + 46 0.000124 0.0001582 0.0001676 0.0001685 + 47 0.0001303 0.0001658 0.0001755 0.0001765 + 48 0.0001363 0.000173 0.0001831 0.0001841 + 49 0.0001419 0.0001798 0.0001902 0.0001913 + 50 0.0001473 0.0001862 0.0001969 0.000198 + 51 0.0001523 0.0001922 0.0002031 0.0002042 + 52 0.0001569 0.0001977 0.0002089 0.00021 + 53 0.0001612 0.0002028 0.0002141 0.0002153 + 54 0.0001651 0.0002074 0.0002189 0.0002201 + 55 0.0001686 0.0002115 0.0002231 0.0002244 + 56 0.0001717 0.0002151 0.0002269 0.0002281 + 57 0.0001744 0.0002182 0.0002301 0.0002313 + 58 0.0001767 0.0002208 0.0002327 0.000234 + 59 0.0001786 0.0002229 0.0002348 0.0002361 + 60 0.0001801 0.0002244 0.0002364 0.0002376 + 61 0.0001812 0.0002255 0.0002374 0.0002387 + 62 0.0001818 0.000226 0.0002379 0.0002391 + 63 0.0001819 0.0002259 0.0002377 0.000239 + 64 0.0001817 0.0002253 0.000237 0.0002383 + 65 0.000181 0.0002242 0.0002358 0.000237 + 66 0.00018 0.0002225 0.0002339 0.0002351 + 67 0.0001786 0.0002205 0.0002318 0.000233 + 68 0.0001768 0.0002181 0.0002291 0.0002303 + 69 0.0001757 0.0002167 0.0002277 0.0002288 + 70 0.0001769 0.0002185 0.0002297 0.0002308 + 71 0.0001777 0.0002198 0.0002311 0.0002323 + 72 0.000178 0.0002206 0.000232 0.0002332 + 73 0.000178 0.0002208 0.0002324 0.0002336 + 74 0.0001775 0.0002206 0.0002322 0.0002334 + 75 0.0001765 0.0002197 0.0002314 0.0002326 + 76 0.0001752 0.0002184 0.0002301 0.0002313 + 77 0.0001734 0.0002165 0.0002282 0.0002294 + 78 0.0001712 0.0002141 0.0002257 0.0002269 + 79 0.0001685 0.0002112 0.0002227 0.0002239 + 80 0.0001655 0.0002077 0.0002192 0.0002204 + 81 0.000162 0.0002037 0.0002151 0.0002163 + 82 0.0001582 0.0001993 0.0002105 0.0002116 + 83 0.0001539 0.0001943 0.0002053 0.0002065 + 84 0.0001493 0.0001889 0.0001997 0.0002008 + 85 0.0001443 0.0001829 0.0001935 0.0001946 + 86 0.0001389 0.0001766 0.0001869 0.0001879 + 87 0.0001332 0.0001697 0.0001798 0.0001808 + 88 0.0001272 0.0001625 0.0001722 0.0001732 + 89 0.0001208 0.0001548 0.0001642 0.0001652 + 90 0.0001141 0.0001468 0.0001558 0.0001567 + 91 0.0001071 0.0001383 0.000147 0.0001479 + 92 9.989e-05 0.0001295 0.0001378 0.0001386 + 93 9.24e-05 0.0001204 0.0001282 0.000129 + 94 8.467e-05 0.000111 0.0001184 0.0001191 + 95 7.672e-05 0.0001013 0.0001082 0.0001089 + 96 6.858e-05 9.135e-05 9.775e-05 9.842e-05 + 97 6.026e-05 8.115e-05 8.706e-05 8.768e-05 + 98 5.179e-05 7.075e-05 7.614e-05 7.671e-05 + 99 4.318e-05 6.017e-05 6.503e-05 6.554e-05 +100 5.107e-05 6.96e-05 7.486e-05 7.541e-05 +101 5.882e-05 7.886e-05 8.45e-05 8.509e-05 +102 6.737e-05 8.936e-05 9.553e-05 9.617e-05 +103 7.574e-05 9.964e-05 0.0001063 0.000107 +104 8.393e-05 0.0001097 0.0001168 0.0001176 +105 9.19e-05 0.0001194 0.0001271 0.0001279 +106 9.965e-05 0.0001289 0.000137 0.0001379 +107 0.0001071 0.0001381 0.0001466 0.0001475 +108 0.0001144 0.0001468 0.0001558 0.0001567 +109 0.0001213 0.0001553 0.0001646 0.0001656 +110 0.0001279 0.0001633 0.000173 0.000174 +111 0.0001342 0.0001709 0.000181 0.000182 +112 0.0001402 0.0001781 0.0001885 0.0001895 +113 0.0001458 0.0001848 0.0001955 0.0001966 +114 0.000151 0.0001911 0.000202 0.0002032 +115 0.0001559 0.0001969 0.000208 0.0002092 +116 0.0001603 0.0002022 0.0002136 0.0002147 +117 0.0001644 0.0002069 0.0002185 0.0002197 +118 0.0001681 0.0002112 0.000223 0.0002242 +119 0.0001713 0.000215 0.0002268 0.0002281 +120 0.0001742 0.0002182 0.0002302 0.0002314 +121 0.0001766 0.0002209 0.0002329 0.0002342 +122 0.0001785 0.000223 0.0002351 0.0002364 +123 0.0001801 0.0002246 0.0002367 0.000238 +124 0.0001812 0.0002257 0.0002377 0.000239 +125 0.0001818 0.0002262 0.0002382 0.0002394 +126 0.000182 0.0002262 0.0002381 0.0002393 +127 0.0001818 0.0002256 0.0002374 0.0002386 +128 0.0001811 0.0002245 0.0002361 0.0002373 +129 0.00018 0.0002228 0.0002343 0.0002355 +130 9.533e-05 0.0001195 0.0001261 0.0001268 +131 9.872e-06 1.53e-05 1.694e-05 1.712e-05 +132 5.229e-06 8.162e-06 9.058e-06 9.154e-06 +133 6.418e-07 1.102e-06 1.255e-06 1.272e-06 +134 2.491e-07 4.786e-07 5.573e-07 5.659e-07 +135 5.796e-08 1.144e-07 1.344e-07 1.366e-07 +136 1.389e-08 3.33e-08 4.077e-08 4.16e-08 +137 6.215e-09 1.455e-08 1.775e-08 1.81e-08 +138 1.156e-09 3.096e-09 3.904e-09 3.997e-09 +139 5.714e-10 1.66e-09 2.129e-09 2.183e-09 +140 1.268e-10 3.937e-10 5.155e-10 5.297e-10 +141 4.685e-11 1.698e-10 2.292e-10 2.362e-10 +142 1.37e-11 5.138e-11 7.001e-11 7.222e-11 +143 4.037e-12 1.825e-11 2.583e-11 2.675e-11 +144 1.854e-12 7.746e-12 1.08e-11 1.117e-11 +145 4.309e-13 2.027e-12 2.914e-12 3.023e-12 +146 2.34e-13 1.229e-12 1.807e-12 1.878e-12 +147 6.14e-14 3.494e-13 5.269e-13 5.492e-13 +148 2.614e-14 1.744e-13 2.707e-13 2.829e-13 +149 8.166e-15 5.648e-14 8.869e-14 9.281e-14 +150 2.962e-15 2.394e-14 3.872e-14 4.064e-14 diff --git a/project/dispersion-mom.dat b/project/dispersion-mom.dat new file mode 100644 index 0000000..d39fdc3 --- /dev/null +++ b/project/dispersion-mom.dat @@ -0,0 +1,8 @@ +p damp +0.004 -0.0338488 +0.005 -0.0338373 +0.0066 -0.0337049 +0.01 -0.0336916 +0.0166 -0.0438038 +0.0333 -0.0527543 +0.05 -0.0713638 diff --git a/project/dispersion-nu.dat b/project/dispersion-nu.dat new file mode 100644 index 0000000..e56ea5a --- /dev/null +++ b/project/dispersion-nu.dat @@ -0,0 +1,6 @@ +nu damp +1e-3 -0.0341514 +1e-2 -0.0359197 +1e-1 -0.0552398 +5e-1 -0.154032 +1 -0.21962 diff --git a/project/dispersion.out b/project/dispersion.out new file mode 100644 index 0000000..0ad6fbc --- /dev/null +++ b/project/dispersion.out @@ -0,0 +1,8 @@ +p pdamp nu nudamp +250 -0.0338488 0 -0.0338488 +200 -0.0338373 1e-3 -0.0341514 +150 -0.0337049 1e-2 -0.0359197 +100 -0.0336916 1e-1 -0.0552398 +60 -0.0438038 5e-1 -0.154032 +30 -0.0527543 1 -0.21962 +20 -0.0713638 diff --git a/project/dopri5-estim.out b/project/dopri5-estim.out new file mode 100644 index 0000000..e59c4b7 --- /dev/null +++ b/project/dopri5-estim.out @@ -0,0 +1,547 @@ +t l2 inf + 7.7000 1.5015 1.5119 + 7.8000 1.4827 1.5048 + 7.9000 1.4612 1.4956 + 8.0000 1.4392 1.4824 + 8.1000 1.4184 1.4628 + 8.2000 1.4005 1.4333 + 8.3000 1.3869 1.3025 + 8.4000 1.3784 1.0340 + 8.5000 1.3756 0.9410 + 8.6000 1.3787 0.8741 + 8.7000 1.3877 0.8890 + 8.8000 1.4022 0.9044 + 8.9000 1.4214 0.9208 + 9.0000 1.4445 0.9386 + 9.1000 1.4705 0.9592 + 9.2000 1.4983 0.9830 + 9.3000 1.5266 1.0111 + 9.4000 1.5544 1.0449 + 9.5000 1.5804 1.1612 + 9.6000 1.6036 1.4631 + 9.7000 1.6227 1.8088 + 9.8000 1.6371 1.9172 + 9.9000 1.6460 1.8000 +10.0000 1.6494 1.7026 +10.1000 1.6474 1.6215 +10.2000 1.6407 1.5528 +10.3000 1.6303 1.5295 +10.4000 1.6175 1.5146 +10.5000 1.6038 1.4926 +10.6000 1.5905 1.4628 +10.7000 1.5789 1.4258 +10.8000 1.5700 1.3851 +10.9000 1.5648 1.2108 +11.0000 1.5638 1.1717 +11.1000 1.5673 1.0692 +11.2000 1.5752 0.9186 +11.3000 1.5872 0.8836 +11.4000 1.6028 0.8937 +11.5000 1.6216 0.9531 +11.6000 1.6428 0.9599 +11.7000 1.6657 0.9707 +11.8000 1.6891 1.2004 +11.9000 1.7121 1.4070 +12.0000 1.7340 1.5766 +12.1000 1.7541 1.7646 +12.2000 1.7720 1.8326 +12.3000 1.7870 1.9403 +12.4000 1.7985 1.7920 +12.5000 1.8057 1.6223 +12.6000 1.8085 1.5803 +12.7000 1.8071 1.5459 +12.8000 1.8024 1.5000 +12.9000 1.7953 1.4420 +13.0000 1.7868 1.3780 +13.1000 1.7776 1.3764 +13.2000 1.7686 1.4014 +13.3000 1.7608 1.4282 +13.4000 1.7554 1.6713 +13.5000 1.7529 1.5765 +13.6000 1.7537 1.3913 +13.7000 1.7571 1.2781 +13.8000 1.7630 1.1455 +13.9000 1.7714 1.2462 +14.0000 1.7833 1.2788 +14.1000 1.7992 1.2552 +14.2000 1.8189 1.2233 +14.3000 1.8405 1.3001 +14.4000 1.8615 1.6692 +14.5000 1.8803 2.0763 +14.6000 1.8966 2.4133 +14.7000 1.9113 2.5923 +14.8000 1.9257 2.6619 +14.9000 1.9402 2.3940 +15.0000 1.9538 1.9942 +15.1000 1.9645 2.0013 +15.2000 1.9706 1.9709 +15.3000 1.9714 1.8735 +15.4000 1.9671 1.7036 +15.5000 1.9587 1.4861 +15.6000 1.9481 1.4788 +15.7000 1.9374 1.6645 +15.8000 1.9288 1.8521 +15.9000 1.9240 1.9544 +16.0000 1.9232 2.1856 +16.1000 1.9246 1.9083 +16.2000 1.9259 1.7066 +16.3000 1.9252 1.4780 +16.4000 1.9231 1.4021 +16.5000 1.9231 1.3489 +16.6000 1.9292 1.2649 +16.7000 1.9442 1.3105 +16.8000 1.9669 1.3438 +16.9000 1.9928 1.4738 +17.0000 2.0158 2.2546 +17.1000 2.0318 2.8899 +17.2000 2.0405 3.2704 +17.3000 2.0450 2.9923 +17.4000 2.0500 2.5638 +17.5000 2.0590 2.1203 +17.6000 2.0726 1.9911 +17.7000 2.0883 1.9527 +17.8000 2.1017 1.9658 +17.9000 2.1088 1.9766 +18.0000 2.1075 2.0562 +18.1000 2.0985 2.0744 +18.2000 2.0852 1.9856 +18.3000 2.0722 2.0207 +18.4000 2.0635 2.3623 +18.5000 2.0607 2.5015 +18.6000 2.0620 2.4885 +18.7000 2.0634 2.5374 +18.8000 2.0611 2.2937 +18.9000 2.0538 2.0001 +19.0000 2.0445 1.5702 +19.1000 2.0390 1.4040 +19.2000 2.0437 1.5602 +19.3000 2.0617 1.6465 +19.4000 2.0907 1.6476 +19.5000 2.1241 1.9038 +19.6000 2.1540 2.9918 +19.7000 2.1741 3.4475 +19.8000 2.1821 3.2617 +19.9000 2.1801 2.9793 +20.0000 2.1734 2.6058 +20.1000 2.1678 2.2785 +20.2000 2.1673 2.1965 +20.3000 2.1724 2.0876 +20.4000 2.1806 1.9413 +20.5000 2.1880 1.8784 +20.6000 2.1924 2.2891 +20.7000 2.1937 2.4232 +20.8000 2.1939 2.2915 +20.9000 2.1957 2.3605 +21.0000 2.2001 2.5441 +21.1000 2.2059 2.6422 +21.2000 2.2100 2.9599 +21.3000 2.2085 3.2511 +21.4000 2.1997 3.2663 +21.5000 2.1855 2.4682 +21.6000 2.1714 1.6980 +21.7000 2.1647 1.7223 +21.8000 2.1712 1.8048 +21.9000 2.1923 1.8043 +22.0000 2.2239 1.7318 +22.1000 2.2593 2.8193 +22.2000 2.2921 4.0594 +22.3000 2.3178 3.7198 +22.4000 2.3342 3.2374 +22.5000 2.3408 2.6234 +22.6000 2.3389 2.2775 +22.7000 2.3317 2.3072 +22.8000 2.3236 2.2279 +22.9000 2.3180 2.2134 +23.0000 2.3167 2.2711 +23.1000 2.3193 2.4265 +23.2000 2.3247 2.3720 +23.3000 2.3323 2.3862 +23.4000 2.3421 2.5617 +23.5000 2.3537 2.8767 +23.6000 2.3655 2.9752 +23.7000 2.3748 3.0322 +23.8000 2.3782 3.0795 +23.9000 2.3733 3.6255 +24.0000 2.3605 4.1062 +24.1000 2.3437 2.8782 +24.2000 2.3301 2.0300 +24.3000 2.3277 2.1005 +24.4000 2.3408 2.1909 +24.5000 2.3668 2.2427 +24.6000 2.3975 3.0401 +24.7000 2.4242 4.3496 +24.8000 2.4421 4.1198 +24.9000 2.4521 3.7196 +25.0000 2.4571 3.1780 +25.1000 2.4597 2.4749 +25.2000 2.4598 2.4007 +25.3000 2.4566 2.3499 +25.4000 2.4499 2.3652 +25.5000 2.4421 2.4762 +25.6000 2.4369 2.5281 +25.7000 2.4379 2.3592 +25.8000 2.4466 2.2717 +25.9000 2.4615 2.4415 +26.0000 2.4788 3.0443 +26.1000 2.4939 3.4151 +26.2000 2.5030 3.3302 +26.3000 2.5036 3.2279 +26.4000 2.4955 3.1444 +26.5000 2.4799 3.5998 +26.6000 2.4600 3.7000 +26.7000 2.4410 2.6084 +26.8000 2.4302 2.1979 +26.9000 2.4349 2.2403 +27.0000 2.4584 2.3048 +27.1000 2.4971 2.3490 +27.2000 2.5403 4.0328 +27.3000 2.5748 4.4218 +27.4000 2.5914 4.1770 +27.5000 2.5880 3.7851 +27.6000 2.5707 3.1564 +27.7000 2.5495 2.5329 +27.8000 2.5338 2.3956 +27.9000 2.5287 2.2778 +28.0000 2.5337 2.1734 +28.1000 2.5455 2.2597 +28.2000 2.5601 2.4356 +28.3000 2.5749 2.6004 +28.4000 2.5892 2.7216 +28.5000 2.6028 2.9480 +28.6000 2.6157 3.3442 +28.7000 2.6269 3.6019 +28.8000 2.6348 3.4684 +28.9000 2.6374 3.3438 +29.0000 2.6330 3.2292 +29.1000 2.6211 3.6573 +29.2000 2.6035 3.6178 +29.3000 2.5854 2.4238 +29.4000 2.5747 2.1192 +29.5000 2.5788 2.2074 +29.6000 2.5993 2.3029 +29.7000 2.6304 2.7533 +29.8000 2.6624 4.6279 +29.9000 2.6876 4.4573 +30.0000 2.7023 4.0206 +30.1000 2.7040 3.3887 +30.2000 2.6919 2.6120 +30.3000 2.6684 2.4617 +30.4000 2.6409 2.3811 +30.5000 2.6193 2.4090 +30.6000 2.6109 2.4456 +30.7000 2.6173 2.4093 +30.8000 2.6344 2.6850 +30.9000 2.6555 3.3718 +31.0000 2.6747 3.7141 +31.1000 2.6891 3.8379 +31.2000 2.6980 3.8379 +31.3000 2.7026 3.8498 +31.4000 2.7035 3.8525 +31.5000 2.7015 3.8188 +31.6000 2.6965 3.7627 +31.7000 2.6880 3.5238 +31.8000 2.6754 3.3233 +31.9000 2.6608 2.3533 +32.0000 2.6519 2.2135 +32.1000 2.6599 2.3328 +32.2000 2.6906 2.4492 +32.3000 2.7358 2.9743 +32.4000 2.7767 4.6367 +32.5000 2.7967 4.5109 +32.6000 2.7926 4.1460 +32.7000 2.7749 3.5690 +32.8000 2.7564 2.9362 +32.9000 2.7431 2.6207 +33.0000 2.7327 2.3537 +33.1000 2.7220 2.2644 +33.2000 2.7123 2.4427 +33.3000 2.7088 2.8041 +33.4000 2.7157 2.8352 +33.5000 2.7325 3.1614 +33.6000 2.7533 3.6566 +33.7000 2.7706 3.8634 +33.8000 2.7798 3.7547 +33.9000 2.7813 3.4525 +34.0000 2.7788 3.1332 +34.1000 2.7768 2.9483 +34.2000 2.7772 2.8637 +34.3000 2.7783 2.9959 +34.4000 2.7765 3.1705 +34.5000 2.7708 2.3148 +34.6000 2.7678 2.0203 +34.7000 2.7783 2.2010 +34.8000 2.8065 2.2869 +34.9000 2.8447 3.5941 +35.0000 2.8823 5.0066 +35.1000 2.9126 4.3017 +35.2000 2.9288 3.3774 +35.3000 2.9233 2.5198 +35.4000 2.8971 2.2547 +35.5000 2.8634 2.9940 +35.6000 2.8373 3.0853 +35.7000 2.8244 2.9329 +35.8000 2.8226 2.7436 +35.9000 2.8297 2.8585 +36.0000 2.8458 2.7612 +36.1000 2.8701 2.7349 +36.2000 2.8974 2.9219 +36.3000 2.9196 3.1853 +36.4000 2.9306 3.3904 +36.5000 2.9297 3.3779 +36.6000 2.9209 3.2696 +36.7000 2.9093 3.1969 +36.8000 2.8979 3.1285 +36.9000 2.8862 3.1629 +37.0000 2.8718 3.3483 +37.1000 2.8557 2.4876 +37.2000 2.8475 2.0736 +37.3000 2.8630 2.2261 +37.4000 2.9069 2.3077 +37.5000 2.9631 3.5398 +37.6000 3.0072 5.0939 +37.7000 3.0266 4.4539 +37.8000 3.0265 3.5901 +37.9000 3.0204 2.5960 +38.0000 3.0155 2.5675 +38.1000 3.0084 2.4722 +38.2000 2.9956 2.6370 +38.3000 2.9808 3.1408 +38.4000 2.9699 3.4208 +38.5000 2.9645 3.1711 +38.6000 2.9634 2.9513 +38.7000 2.9676 2.9384 +38.8000 2.9805 2.7868 +38.9000 3.0050 2.7581 +39.0000 3.0399 2.9162 +39.1000 3.0793 3.1506 +39.2000 3.1151 3.3143 +39.3000 3.1399 3.2292 +39.4000 3.1486 3.2043 +39.5000 3.1394 3.7550 +39.6000 3.1142 3.9449 +39.7000 3.0822 2.7047 +39.8000 3.0634 2.1865 +39.9000 3.0794 2.3224 +40.0000 3.1326 2.3293 +40.1000 3.2029 4.1282 +40.2000 3.2710 5.3143 +40.3000 3.3290 4.1919 +40.4000 3.3738 2.9702 +40.5000 3.4053 2.5566 +40.6000 3.4308 2.4254 +40.7000 3.4585 3.3623 +40.8000 3.4868 2.7704 +40.9000 3.5077 3.8206 +41.0000 3.5206 3.2631 +41.1000 3.5359 3.3764 +41.2000 3.5617 3.4951 +41.3000 3.5953 3.6522 +41.4000 3.6289 3.7698 +41.5000 3.6608 3.8985 +41.6000 3.6994 3.7541 +41.7000 3.7570 3.6687 +41.8000 3.8418 3.5482 +41.9000 3.9540 3.7144 +42.0000 4.0861 3.5874 +42.1000 4.2262 4.1079 +42.2000 4.3648 5.0002 +42.3000 4.5049 3.2342 +42.4000 4.6683 2.5143 +42.5000 4.8823 2.4176 +42.6000 5.1485 2.8050 +42.7000 5.4370 4.4757 +42.8000 5.7193 5.7495 +42.9000 5.9926 5.4658 +43.0000 6.2676 5.6304 +43.1000 6.5447 7.4736 +43.2000 6.8063 7.1153 +43.3000 7.0289 9.0537 +43.4000 7.1991 10.2142 +43.5000 7.3171 11.5665 +43.6000 7.3898 13.1804 +43.7000 7.4302 12.9262 +43.8000 7.4649 11.7906 +43.9000 7.5353 11.0636 +44.0000 7.6844 10.6178 +44.1000 7.9443 10.8463 +44.2000 8.3355 10.6177 +44.3000 8.8723 11.5468 +44.4000 9.5662 13.9297 +44.5000 10.4216 15.3499 +44.6000 11.4289 15.3551 +44.7000 12.5594 17.1899 +44.8000 13.7707 17.9040 +44.9000 15.0301 13.4360 +45.0000 16.3457 11.6864 +45.1000 17.7600 11.7695 +45.2000 19.2945 14.6262 +45.3000 20.9100 24.1558 +45.4000 22.5302 25.7212 +45.5000 24.0879 27.1021 +45.6000 25.5376 29.6476 +45.7000 26.8486 36.7177 +45.8000 27.9900 43.3418 +45.9000 28.9209 49.5892 +46.0000 29.6049 57.6418 +46.1000 30.0481 64.0814 +46.2000 30.3336 62.1373 +46.3000 30.6085 56.7349 +46.4000 31.0416 52.0246 +46.5000 31.7863 48.9694 +46.6000 32.9861 46.9365 +46.7000 34.7705 52.0996 +46.8000 37.2333 58.6447 +46.9000 40.3964 70.1562 +47.0000 44.2107 77.0848 +47.1000 48.5773 82.0185 +47.2000 53.3520 82.2928 +47.3000 58.3272 90.4458 +47.4000 63.2546 97.4713 +47.5000 67.9477 80.3523 +47.6000 72.3674 76.5010 +47.7000 76.5499 82.4504 +47.8000 80.4531 108.6111 +47.9000 83.9517 94.8132 +48.0000 86.9624 99.1434 +48.1000 89.4891 88.7061 +48.2000 91.5638 97.5522 +48.3000 93.2099 93.2620 +48.4000 94.4529 112.8862 +48.5000 95.3374 106.2050 +48.6000 95.9343 109.4269 +48.7000 96.3139 111.9826 +48.8000 96.5594 112.0209 +48.9000 96.7679 107.7305 +49.0000 97.0816 109.4482 +49.1000 97.6256 110.7368 +49.2000 98.4688 110.5140 +49.3000 99.5817 108.6136 +49.4000 100.8864 108.9902 +49.5000 102.2840 105.5908 +49.6000 103.6652 107.5762 +49.7000 104.9275 104.5201 +49.8000 106.0048 107.1801 +49.9000 106.8826 111.8003 +50.0000 107.5845 109.8322 +50.1000 108.1499 112.7777 +50.2000 108.6153 112.7711 +50.3000 109.0001 105.6903 +50.4000 109.3086 107.7021 +50.5000 109.5459 110.2181 +50.6000 109.7255 103.3274 +50.7000 109.8640 112.2041 +50.8000 109.9720 103.8979 +50.9000 110.0509 115.9229 +51.0000 110.1005 106.8103 +51.1000 110.1140 109.7323 +51.2000 110.0952 110.5089 +51.3000 110.0407 109.7485 +51.4000 109.9648 111.9342 +51.5000 109.8721 110.7906 +51.6000 109.7825 112.5287 +51.7000 109.7071 112.4406 +51.8000 109.6611 112.7281 +51.9000 109.6503 108.1335 +52.0000 109.6736 109.4097 +52.1000 109.7264 110.2488 +52.2000 109.8014 112.9146 +52.3000 109.8931 106.3314 +52.4000 109.9962 109.4330 +52.5000 110.1057 112.8982 +52.6000 110.2181 107.3449 +52.7000 110.3297 110.7347 +52.8000 110.4373 113.3563 +52.9000 110.5375 106.2983 +53.0000 110.6273 114.1591 +53.1000 110.7054 110.0901 +53.2000 110.7709 109.8047 +53.3000 110.8234 111.9676 +53.4000 110.8623 108.3027 +53.5000 110.8843 113.6404 +53.6000 110.8876 108.3901 +53.7000 110.8654 110.9662 +53.8000 110.8170 110.2344 +53.9000 110.7373 110.2682 +54.0000 110.6303 112.0771 +54.1000 110.5002 109.9154 +54.2000 110.3575 111.3108 +54.3000 110.2173 111.8900 +54.4000 110.0914 111.9099 +54.5000 109.9935 109.1618 +54.6000 109.9282 108.9591 +54.7000 109.8984 115.0331 +54.8000 109.9009 112.1358 +54.9000 109.9305 106.6837 +55.0000 109.9813 107.2941 +55.1000 110.0463 114.6216 +55.2000 110.1201 106.2452 +55.3000 110.1974 107.6057 +55.4000 110.2747 111.6553 +55.5000 110.3500 114.5522 +55.6000 110.4214 104.6730 +55.7000 110.4886 112.4311 +55.8000 110.5490 108.2776 +55.9000 110.6000 109.8214 +56.0000 110.6372 109.8468 +56.1000 110.6547 109.6932 +56.2000 110.6483 111.2319 +56.3000 110.6130 109.0113 +56.4000 110.5479 112.9227 +56.5000 110.4559 109.0443 +56.6000 110.3412 112.3414 +56.7000 110.2181 109.6553 +56.8000 110.0946 111.6544 +56.9000 109.9885 113.5577 +57.0000 109.9057 110.4023 +57.1000 109.8565 106.8017 +57.2000 109.8416 108.5771 +57.3000 109.8613 114.3187 +57.4000 109.9113 110.9765 +57.5000 109.9857 104.6741 +57.6000 110.0782 113.3852 +57.7000 110.1815 110.1415 +57.8000 110.2907 108.7656 +57.9000 110.4003 112.8906 +58.0000 110.5079 106.8445 +58.1000 110.6113 107.9533 +58.2000 110.7100 111.7745 +58.3000 110.8039 114.3637 +58.4000 110.8929 107.3912 +58.5000 110.9762 112.2256 +58.6000 111.0517 108.8292 +58.7000 111.1172 111.4456 +58.8000 111.1689 110.5295 +58.9000 111.2050 110.1860 +59.0000 111.2218 111.2979 +59.1000 111.2200 111.2750 +59.2000 111.1972 110.3359 +59.3000 111.1575 112.3448 +59.4000 111.1017 110.4275 +59.5000 111.0368 112.7702 +59.6000 110.9684 111.2722 +59.7000 110.9041 112.3445 +59.8000 110.8506 110.6671 +59.9000 110.8129 112.3284 +60.0000 110.7925 111.8962 +60.1000 110.7899 109.7376 +60.2000 110.8009 112.1868 +60.3000 110.8222 111.1264 +60.4000 110.8462 112.8006 +60.5000 110.8664 108.7430 +60.6000 110.8724 113.8865 +60.7000 110.8548 108.4519 +60.8000 110.7980 112.9855 +60.9000 110.6908 110.8250 +61.0000 110.5119 109.9846 +61.1000 110.2580 108.9705 +61.2000 109.9264 110.9136 +61.3000 109.5604 108.3379 +61.4000 109.2098 113.6157 +61.5000 108.9356 110.5189 +61.6000 108.7565 100.3124 +61.7000 108.6684 108.5466 +61.8000 108.6516 116.8288 +61.9000 108.6920 110.8823 +62.0000 108.7780 99.0097 +62.1000 108.9041 111.9406 +62.2000 109.0627 113.1317 diff --git a/project/dopri5-p.out b/project/dopri5-p.out new file mode 100644 index 0000000..7d7fac4 --- /dev/null +++ b/project/dopri5-p.out @@ -0,0 +1,8 @@ +p eig err +80 160 6 +70 150 5 +60 138 4 +50 125 4 +40 110 3 +30 94 2 +20 75 2 diff --git a/project/dopri5-p.py b/project/dopri5-p.py new file mode 100755 index 0000000..8c1723d --- /dev/null +++ b/project/dopri5-p.py @@ -0,0 +1,43 @@ +#!/usr/bin/python + +import matplotlib +import numpy as np +import matplotlib.pyplot as plt + +p = [ 80, + 70, + 60, + 50, + 40, + 30, + 20 ] + +eig = [ 160, + 150, + 138, + 125, + 110, + 94, + 75 ] + +p = np.array(p) +eig = np.array(eig) + +fit = np.polyfit(np.log(p), np.log(eig), 1) + +fig, axis = plt.subplots() +axis.set_title("DOPRI5 estim") +axis.set_xlabel("p") +axis.set_ylabel("lambda") + + +axis.loglog(p, eig, 'r*', label='data') +axis.loglog(p, np.exp(fit[0] * np.log(p) + fit[1]), 'k-', label='fit') + +axis.legend() + +plt.show() + +print("Fit: %g * p^%g" % (np.exp(fit[1]), fit[0])) + + diff --git a/project/eig_lambda.out b/project/eig_lambda.out new file mode 100644 index 0000000..2b2a4d5 --- /dev/null +++ b/project/eig_lambda.out @@ -0,0 +1,74 @@ +t l2 ana +1.0000 1.1305 18.6327 +1.1000 1.1164 18.7627 +1.2000 1.1152 18.8802 +1.3000 1.1305 18.9846 +1.4000 1.1681 19.0751 +1.5000 1.2265 19.1510 +1.6000 1.3025 19.2112 +1.7000 1.3846 19.2548 +1.8000 1.4612 19.2810 +1.9000 1.5230 19.2893 +2.0000 1.5597 19.2794 +2.1000 1.5687 19.2517 +2.2000 1.5514 19.2067 +2.3000 1.5170 19.1458 +2.4000 1.4754 19.0704 +2.5000 1.4350 18.9827 +2.6000 1.4065 18.8851 +2.7000 1.3927 18.7815 +2.8000 1.3928 18.6777 +2.9000 1.4009 18.5831 +3.0000 1.4031 18.5173 +3.1000 1.3821 18.4812 +3.2000 1.3246 18.4711 +3.3000 1.2285 18.4779 +3.4000 1.1095 18.4865 +3.5000 0.9983 18.4888 +3.6000 0.9385 18.5189 +3.7000 0.9792 18.5469 +3.8000 1.0510 18.5621 +3.9000 1.2738 18.5659 +4.0000 1.4545 18.5634 +4.1000 1.5612 18.5641 +4.2000 2.0124 18.5772 +4.3000 2.0539 18.6030 +4.4000 2.1134 18.6322 +4.5000 2.6045 18.6553 +4.6000 2.5942 18.6659 +4.7000 2.4522 18.6708 +4.8000 2.6404 18.6758 +4.9000 2.7359 18.6646 +5.0000 2.6325 18.6312 +5.1000 2.5482 18.5873 +5.2000 2.5650 18.5388 +5.3000 2.6198 18.4871 +5.4000 2.8205 18.4382 +5.5000 3.1982 18.3984 +5.6000 3.6838 18.3647 +5.7000 4.9797 18.3284 +5.8000 4.9051 18.3013 +5.9000 5.9687 18.2956 +6.0000 8.1196 18.3140 +6.1000 5.4013 18.3442 +6.2000 5.9701 18.3738 +6.3000 12.2895 18.3949 +6.4000 7.3903 18.4024 +6.5000 2.9802 18.4139 +6.6000 9.5110 18.4370 +6.7000 18.3143 18.4528 +6.8000 5.8830 18.4650 +6.9000 3.6886 18.4811 +7.0000 8.8508 18.5056 +7.1000 14.5809 18.5319 +7.2000 7.8262 18.5449 +7.3000 6.5433 18.5349 +7.4000 7.8123 18.5108 +7.5000 9.4798 18.4887 +7.6000 10.6287 18.4674 +7.7000 10.6861 18.4292 +7.8000 10.5240 18.3770 +7.9000 11.2621 18.3233 +8.0000 11.9804 18.2705 +8.1000 12.0554 18.2371 +8.2000 12.0242 18.2272 diff --git a/project/energy_distribution.pdf b/project/energy_distribution.pdf new file mode 100644 index 0000000..ee90600 Binary files /dev/null and b/project/energy_distribution.pdf differ diff --git a/project/imex-erk-stability.pdf b/project/imex-erk-stability.pdf new file mode 100644 index 0000000..a02e96f Binary files /dev/null and b/project/imex-erk-stability.pdf differ diff --git a/project/imex-imrk-stability.pdf b/project/imex-imrk-stability.pdf new file mode 100644 index 0000000..5e148c1 Binary files /dev/null and b/project/imex-imrk-stability.pdf differ diff --git a/project/phi.fig b/project/phi.fig new file mode 100644 index 0000000..473afb0 Binary files /dev/null and b/project/phi.fig differ diff --git a/project/rkc4-stability.pdf b/project/rkc4-stability.pdf new file mode 100644 index 0000000..61ea377 Binary files /dev/null and b/project/rkc4-stability.pdf differ diff --git a/stocha/NISDE_information_exam_preparation_2021.pdf b/stocha/NISDE_information_exam_preparation_2021.pdf new file mode 100644 index 0000000..1f3fd13 Binary files /dev/null and b/stocha/NISDE_information_exam_preparation_2021.pdf differ