fork download
copy 
  1. def inveuler(q,p,dt,t):
  2. 	invd=1/(1+dt**2) 
  3. 	qt = invd*q + invd*p*dt
  4. 	pt = invd*p - invd*(4*α*q**3 - 2*q)*dt
  5. 	return qt, pt
  6. dt = 0.02
  7. istep = 300
  8. q = 3
  9. p = 10
  10. α = 0.02
  11. for i in range(1, istep): 
  12. 	t = i *dt
  13. 	qt,pt = inveuler(q,p,dt,t) 
  14. 	print(t,qt,pt)
  15. 	q=qt
  16. 	p=pt
Success #stdin #stdout 0.1s 14036KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
0.02 3.198720511795282 10.072770891643342
0.04 3.398816403066922 10.144296057953857
0.06 3.600262219338264 10.213342629591256
0.08 3.8030078687825757 10.278575774726118
0.1 4.006976593639642 10.33855682028242
0.12 4.212062904883337 10.39174244959444
0.14 4.418130501674556 10.436485245045294
0.16 4.6250102024944635 10.471035867507307
0.18 4.832497920676339 10.49354718418998
0.2 5.04035272327083 10.50208067090467
0.22 5.248295018681451 10.494615421526913
0.24 5.456004925141932 10.469060094017514
0.26 5.663120878670814 10.423268105940293
0.28 5.869238545371471 10.355056360021738
0.3 6.073910108528494 10.262227728965652
0.32 6.276644005505606 10.142597455692991
0.34 6.476905192542449 9.994023528115331
0.36 6.674116016698076 9.814440964944268
0.38 6.8676577728878065 9.601899800513019
0.4 7.056873019690191 9.354606383310712
0.42 7.241068719868458 9.070967408014667
0.44 7.419520259924782 8.749635889652403
0.46 7.591476387162964 8.389558069079987
0.48 7.756165082511559 7.990020021816355
0.5 7.9128003628027646 7.550692540530989
0.52 8.060589977622335 7.0716726903621305
0.54 8.198743933856036 6.553520312338211
0.56 8.32648374660416 5.99728769038121
0.58 8.443052279499984 5.404540617489489
0.6 8.547724002248874 4.777369209800109
0.62 8.639815460260772 4.118387031905051
0.64 8.71869572260983 3.4307174154620257
0.66 8.783796552298151 2.717966270655154
0.68 8.834622028899695 1.9841811928443365
0.7000000000000001 8.870757349816653 1.2337972326499085
0.72 8.891876543852112 0.4715702968106753
0.74 8.897748850248226 -0.2975002565756794
0.76 8.888243547697632 -1.0682577108094358
0.78 8.86333306025734 -1.8354801815214636
0.8 8.823094218939334 -2.5939759615471996
0.8200000000000001 8.767707616661726 -3.3386788676665007
0.84 8.697455057285481 -4.064740394356194
0.86 8.612715163333023 -4.7676156080360785
0.88 8.513957268264996 -5.443140007453101
0.9 8.401733774606091 -6.087595011760184
0.92 8.276671205888531 -6.697760286524344
0.9400000000000001 8.139460216071615 -7.270951740500276
0.96 7.99084484332428 -7.805044679607127
0.98 7.831611305210052 -8.29848224618585
1.0 7.662576629634481 -8.750269862098774
1.02 7.484577401431933 -9.1599569012747
1.04 7.298458879854497 -9.527607217683263
1.06 7.105064709616984 -9.853760435326636
1.08 6.905227409946472 -10.139386064576335
1.1 6.699759784741049 -10.385832549958275
1.12 6.489447354799964 -10.594773291445495
1.1400000000000001 6.2750418722221655 -10.768151532693258
1.16 6.057255939192623 -10.908125796514181
1.18 5.836758719774429 -11.01701729202642
1.2 5.61417270485196 -11.097260440033331
1.22 5.390071467464307 -11.15135738162378
1.24 5.164978328500431 -11.181837064775058
1.26 4.939365840868582 -11.19121925629391
1.28 4.713655993345366 -11.181983609475193
1.3 4.488221032742765 -11.156543735735765
1.32 4.2633848041064075 -11.117226082632683
1.34 4.0394245126486945 -11.066253310217862
1.36 3.8165728173174105 -11.005731780025057
1.3800000000000001 3.5950201736474505 -10.937642722407194
1.4000000000000001 3.3749173522584033 -10.863836624097354
1.42 3.156378068549037 -10.78603037420272
1.44 2.9394816683976233 -10.705806718915127
1.46 2.7242758236898448 -10.624615598925882
1.48 2.5107792000313145 -10.543776975230926
1.5 2.298984066899936 -10.464484785614124
1.52 2.0888588276565914 -10.387811713072548
1.54 1.8803504532138546 -10.314714486810015
1.56 1.6733868087541528 -10.24603947468079
1.58 1.4678788677134516 -10.18252836205576
1.6 1.2637228113477974 -10.124823745317357
1.62 1.060802015635196 -10.07347449817734
1.6400000000000001 0.8589889300996093 -10.028940795602152
1.6600000000000001 0.6581468554453881 -9.991598703356038
1.68 0.4581316287267766 -9.961744261200044
1.7 0.25879322621229073 -9.939597004878964
1.72 0.05997729519663281 -9.92530288651656
1.74 -0.13847337318442463 -9.918936565290387
1.76 -0.33671741752322304 -9.920503050662086
1.78 -0.5349135131312123 -9.929938689397845
1.8 -0.7332189993194415 -9.947111495528254
1.82 -0.9317885138244768 -9.971820829664763
1.84 -1.1307726213692244 -10.003796441101873
1.86 -1.3303164236218132 -10.042696893271685
1.8800000000000001 -1.530558138231954 -10.08810740075264
1.9000000000000001 -1.7316276351929296 -10.13953711452826
1.92 -1.9336449195156886 -10.196415901892424
1.94 -2.1367185501334838 -10.258090678636437
1.96 -2.3409439861117676 -10.323821364238093
1.98 -2.5464018526554675 -10.392776545980183
2.0 -2.7531561211266204 -10.464028955486489
2.02 -2.9612521993566077 -10.53655088124533
2.04 -3.1707149310091105 -10.609209663377474
2.06 -3.38154650567439 -10.680763442166045
2.08 -3.5937242848037894 -10.749857359523432
2.1 -3.8071985525732286 -10.815020442259163
2.12 -4.021890205336278 -10.874663427136982
2.14 -4.23768839851961 -10.927077819392933
2.16 -4.454448175637213 -10.970436507424084
2.18 -4.671988110541478 -11.002796285169364
2.2 -4.890088001044447 -11.022102658294877
2.22 -5.108486659546526 -11.02619732823267
2.24 -5.326879854169511 -11.01282875653028
2.2600000000000002 -5.54491846191535 -10.979666207566705
2.2800000000000002 -5.762206903305363 -10.924317646864548
2.3000000000000003 -5.978301935468466 -10.844351831165167
2.32 -6.192711887336834 -10.737324861322424
2.34 -6.404896425992885 -10.600811376378838
2.36 -6.614266946741766 -10.432440444019893
2.38 -6.820187680549943 -10.229936047134032
2.4 -7.021977610448444 -9.991161878182162
2.42 -7.218913282699008 -9.714169934347918
2.44 -7.410232588350626 -9.397252161470119
2.46 -7.5951395757497275 -9.038994131102356
2.48 -7.772810334238079 -8.638329463639694
2.5 -7.942399963525462 -8.194593445698391
2.52 -8.103050612194552 -7.7075740494446
2.54 -8.253900532970256 -7.177558365517138
2.56 -8.394094062655537 -6.605372331206498
2.58 -8.522792392322737 -5.9924115931469135
2.6 -8.639184950205594 -5.34066140835374
2.62 -8.74250117790151 -4.652703674020552
2.64 -8.83202244240496 -3.9317094933771584
2.66 -8.907093794754601 -3.181416131519156
2.68 -8.967135263279673 -2.406087780357388
2.7 -9.011652357943642 -1.6104602133074302
2.72 -9.040245464024181 -0.7996701346857575
2.74 -9.05261781959006 0.02082922669853271
2.7600000000000002 -9.048581802335155 0.8453640146822659
2.7800000000000002 -9.02806329672282 1.6681442333767733
2.8000000000000003 -8.991103970467098 2.4833734193683945
2.82 -8.937861357536715 3.2853597580572043
2.84 -8.868606719687694 4.06862477007966
2.86 -8.783720735991704 4.828005775343642
2.88 -8.68368714562658 5.5587486359735525
2.9 -8.569084539091472 6.25658776506968
2.92 -8.44057655316881 6.917811029880207
2.94 -8.298900772262302 7.539307928246669
2.96 -8.144856671028956 8.118600221569679
2.98 -7.979292949417796 8.65385501009155
3.0 -7.8030946113714155 9.143880984802612
3.02 -7.617170123625913 9.588109240795507
3.04 -7.422438963224713 9.986560556945436
3.06 -7.219819824156142 10.339801417103525
3.08 -7.010219707930899 10.648891262620667
3.1 -6.794524073049266 10.91532353116976
3.12 -6.5735881671590075 11.14096296859479
3.14 -6.348229615940736 11.327981522005453
3.16 -6.119222296581994 11.47879486065302
3.18 -5.887291482775824 11.596001254462953
3.2 -5.653110213601125 11.682324195325357
3.22 -5.41729681097023 11.740559797063684
3.24 -5.180413449649097 11.77352967580924
3.2600000000000002 -4.942965669864966 11.784039707853063
3.2800000000000002 -4.705402714622056 11.774844796634634
3.3000000000000003 -4.468118571260859 11.748619559626206
3.3200000000000003 -4.231453598628883 11.707934670976547
3.34 -3.9956966265587277 11.655238465407011
3.36 -3.761087422281675 11.592843319374353
3.38 -3.5278194281229385 11.522916271915904
3.4 -3.2960426856103764 11.447473324148763
3.42 -3.0658668723784497 11.368376857204277
3.44 -2.837364389478573 11.287335627766682
3.46 -2.6105734475442217 11.20590683317382
3.48 -2.3855011104365706 11.125499779713877
3.5 -2.1621262643365586 11.047380734536311
3.52 -1.9404024886503723 10.972678590450444
3.54 -1.720260812516357 10.902391021474626
3.56 -1.5016123471480052 10.837390853615346
3.58 -1.2843507897597943 10.778432418835374
3.6 -1.0683547994633016 10.726157699801815
3.62 -0.8534902493675183 10.681102108441861
3.64 -0.6396123622537795 10.6436997725129
3.66 -0.42656773970763845 10.614288231464858
3.68 -0.2141962965597174 10.593112466140516
3.7 -0.0023331139913105314 10.580328206735492
3.72 0.20918977423370586 10.576004480403999
3.74 0.42054164718291265 10.580125374471283
3.7600000000000002 0.6318913981130931 10.59259100394181
3.7800000000000002 0.8434058558495895 10.613217683423747
3.8000000000000003 1.0552481102739548 10.641737314273138
3.8200000000000003 1.267575826228926 10.677796008234944
3.84 1.4805395305813924 10.72095197965057
3.86 1.6942808578312714 10.770672749924309
3.88 1.908930740533544 10.826331720886602
3.9 2.1246075319385005 10.887204188426816
3.92 2.341415049687162 10.952462884724678
3.94 2.5594405311691877 11.021173156973148
3.96 2.7787524933113263 11.092287912976458
3.98 2.999398492173986 11.164642489656574
4.0 3.221402780854775 11.236949629420108
4.0200000000000005 3.4447638678960186 11.307794781476595
4.04 3.6694519827324577 11.37563198030254
4.0600000000000005 3.8954064597546068 11.438780590996133
4.08 4.1225330583511886 11.49542325042401
4.1 4.350701242862523 11.543605372593346
4.12 4.579741453732896 11.581236624903053
4.14 4.80944240926725 11.606094816631801
4.16 5.039548486205404 11.615832669436532
4.18 5.269757236699455 11.607987958399104
4.2 5.499717109023828 11.579997517325554
4.22 5.7290254491906625 11.52921558910792
4.24 5.957226870224731 11.45293696614915
4.26 6.1838120847138285 11.348425302100585
4.28 6.408217303834307 11.212946879583445
4.3 6.629824311701295 11.04380998484303
4.32 6.847961326867409 10.838409866204607
4.34 7.061904762286587 10.594279037294784
4.36 7.270881990236388 10.309142429248146
4.38 7.474075208737856 9.980976602709115
4.4 7.67062649019596 9.608071908401444
4.42 7.859644070735695 9.18909614680063
4.44 8.040209909707825 8.723157940078872
4.46 8.211388513104161 8.209867714563927
4.48 8.372236972606396 7.649393924851086
4.5 8.521816124653556 7.042511959353414
4.5200000000000005 8.659202682767518 6.390643079978098
4.54 8.783502143509676 5.695880792311647
4.5600000000000005 8.89386221447012 4.961002238720564
4.58 8.989486464658667 4.189462568437292
4.6000000000000005 9.069647856884657 3.3853707680293335
4.62 9.133701791528633 2.5534461209449524
4.64 9.181098274637677 1.6989552793234055
4.66 9.211392823094906 0.8276308332732049
4.68 9.224255737465384 -0.05442680228231811
4.7 9.219479409655875 -0.940861459134337
4.72 9.19698338711834 -1.825183716204798
4.74 9.156816985999845 -2.7008996573184945
4.76 9.099159329121827 -3.5616412914180033
4.78 9.024316776582834 -4.40129387498602
4.8 8.932717811958332 -5.214115487871729
4.82 8.824905539984904 -5.994844597287229
4.84 8.701528036824431 -6.738791972572705
4.86 8.563326866626326 -7.441914137070645
4.88 8.41112413423122 -8.100866499401706
4.9 8.245808480850846 -8.713035322333873
4.92 8.068320446225679 -9.276548689925793
4.94 7.879637617380211 -9.790267555413768
4.96 7.680759962287021 -10.253758738119865
4.98 7.472695709240926 -10.667252348298371
5.0 7.256448083041742 -11.031586533115497
5.0200000000000005 7.033003151118985 -11.348142651420213
5.04 6.80331897050238 -11.618774011817841
5.0600000000000005 6.568316163800503 -11.845731172276324
5.08 6.328869992358033 -12.031586532736123
5.1000000000000005 6.08580394012726 -12.179160591617315
5.12 5.839884774385159 -12.291451819588424
5.14 5.5918190103892345 -12.371571663493894
5.16 5.342250676848617 -12.42268575893508
5.18 5.0917602575668885 -12.44796202443326
5.2 4.840864671209739 -12.450525949521747
5.22 4.590018144961319 -12.433423083331688
5.24 4.339613837759582 -12.399588483490417
5.26 4.0899860736603095 -12.351822697139601
5.28 3.8414130544957192 -12.292773712898063
5.3 3.594119932264852 -12.224924238645475
5.32 3.348282134638087 -12.150583617788207
5.34 3.104028850742026 -12.071883688400964
5.36 2.861446598334673 -11.990777907666654
5.38 2.620582807058516 -11.909043101284995
5.4 2.3814493652867017 -11.828283247755191
5.42 2.1440260898956396 -11.749934765407072
5.44 1.9082640889519171 -11.675272831516185
5.46 1.6740889967229042 -11.605418324516611
5.48 1.44140406860513 -11.541345039822037
5.5 1.2100931305564666 -11.483886885442427
5.5200000000000005 0.9800233834942202 -11.433744814436519
5.54 0.7510480679782985 -11.39149329677665
5.5600000000000005 0.5230089984433881 -11.357586173327977
5.58 0.29573897938507454 -11.332361769564228
5.6000000000000005 0.06906411834645138 -11.316047176778492
5.62 -0.15719394761007444 -11.30876163447189
5.64 -0.38321589394193545 -11.31051896998419
5.66 -0.6091826003014986 -11.321229070996086
5.68 -0.8352730724924233 -11.340698384040337
5.7 -1.0616623752231409 -11.368629448375533
5.72 -1.2885195563681042 -11.404619490265315
5.74 -1.5160055439558282 -11.448158118616684
5.76 -1.7442709979289903 -11.498624179790314
5.78 -1.9734540998848427 -11.555281847899309
5.8 -2.203678265536614 -11.617276047714013
5.82 -2.4350497665842608 -11.68362733098554
5.84 -2.6676552511035303 -11.753226354107397
5.86 -2.9015591545238686 -11.824828135963324
5.88 -3.1368009968443973 -11.89704630982546
5.9 -3.3733925660145005 -11.968347622340804
5.92 -3.611314992464331 -12.037046975792595
5.94 -3.8505157256899065 -12.101303356437473
5.96 -4.090905430646398 -12.159117040881931
5.98 -4.332354829532223 -12.208328522757991
https://ideone.com/XC0UOQ
language:
Python 3 (python  3.12)
created:
14 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

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