fork(1) download 
  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*q*dt
  5. 	return qt, pt
  6. dt= 0.01
  7. istep = 100
  8. q=0.0
  9. p=-0.1
  10. for i in range(1, istep): 
  11. 	t = i *dt
  12. 	qt,pt = inveuler(q,p,dt,t) 
  13. 	print(t,qt,pt)
  14. 	q=qt
  15. 	p=pt
Success #stdin #stdout 0.12s 14232KB
comments (?)
stdin 
Standard input is empty
stdout 
0.01 -0.0009999000099990003 -0.09999000099990002
0.02 -0.001999600059992001 -0.09997000499930009
0.03 -0.0029990002099640053 -0.09994001499720044
0.04 -0.0039980005598800215 -0.09990003499160163
0.05 -0.004996501259670071 -0.09985006997900493
0.06 -0.005994402519208199 -0.09979012595381284
0.07 -0.006991604618284499 -0.09972020990762999
0.08 -0.007988007916569142 -0.0996403298284643
0.09 -0.008983512863567428 -0.09955049469982863
0.1 -0.009978020008564858 -0.09945071449974298
0.11 -0.010971430010561231 -0.09934100019963736
0.12 -0.011963643648192787 -0.09922136376315543
0.13 -0.012954561829641376 -0.09909181814485903
0.14 -0.013944085602529714 -0.09895237728883373
0.15 -0.01493211616380167 -0.09880305612719571
0.16 -0.015918554869586668 -0.09864387057849984
0.17 -0.01690330324504716 -0.09847483754604937
0.18 -0.017886262994208233 -0.09829597491610728
0.19 -0.018867336009768327 -0.0981073015560096
0.2 -0.019846424382890134 -0.09790883731218068
0.21 -0.020823430412970644 -0.09770060300805097
0.22 -0.021798256617389412 -0.09748262044187708
0.23 -0.02277080574123406 -0.09725491238446475
0.24 -0.023740980767002005 -0.09701750257679473
0.25 -0.024708684924277528 -0.09677041572755196
0.26 -0.025673821699383106 -0.09651367751055813
0.27 -0.026636294845004186 -0.0962473145621081
0.28 -0.02759600838978629 -0.09597135447821023
0.29 -0.0285528666479036 -0.0956858258117312
0.3 -0.02950677422859805 -0.09539075806944522
0.31 -0.030457636045687932 -0.09508618170898833
0.32 -0.03140535732704511 -0.09477212813571788
0.33 -0.032349843624039885 -0.09444862969947747
0.34 -0.033291000820952564 -0.09411571969126795
0.35000000000000003 -0.03422873514435081 -0.09377343233982444
0.36 -0.035162953172431816 -0.09342180280810011
0.37 -0.03609356184432838 -0.09306086718965682
0.38 -0.03702046846937801 -0.09269066250496305
0.39 -0.037943580736354 -0.09231122669759952
0.4 -0.03886280672265773 -0.09192259863037294
0.41000000000000003 -0.03977805490347112 -0.09152481808133824
0.42 -0.040689234160868414 -0.09111792573972954
0.43 -0.041596253792886416 -0.09070196320180068
0.44 -0.042499023522552165 -0.09027697296657515
0.45 -0.04339745350686723 -0.08984299843150648
0.46 -0.044291454345747716 -0.089400083888049
0.47000000000000003 -0.04518093709091911 -0.08894827451713981
0.48 -0.04606581325476503 -0.08848761638459217
0.49 -0.04694599481912904 -0.08801815643640087
0.5 -0.047821394244068643 -0.08753994249396019
0.51 -0.048691924476560594 -0.08705302324919459
0.52 -0.04955749895915662 -0.08655744825960301
0.53 -0.050418031638588795 -0.08605326794321712
0.54 -0.05127343697432354 -0.08554053357347388
0.55 -0.052123629947063564 -0.08501929727400324
0.56 -0.052968526067196875 -0.08448961201333127
0.5700000000000001 -0.05380804138319187 -0.08395153159949935
0.58 -0.054642092489937864 -0.08340511067459996
0.59 -0.05547059653703016 -0.08285040470922966
0.6 -0.05629347123699876 -0.08228746999685967
0.61 -0.05711063487348001 -0.08171636364812486
0.62 -0.05792200630933033 -0.08113714358503156
0.63 -0.05872750499468118 -0.08054986853508475
0.64 -0.059527050974934534 -0.07995459802533542
0.65 -0.06032056489869802 -0.07935139237634843
0.66 -0.06110796802565893 -0.07874031269609184
0.67 -0.061889182234396406 -0.07812142087374788
0.68 -0.06266413003013087 -0.07749477957344657
0.6900000000000001 -0.0634327345524101 -0.07686045222792247
0.7000000000000001 -0.06419491958273105 -0.07621850303209517
0.71 -0.06495060955209679 -0.07556899693657419
0.72 -0.06569972954850768 -0.0749119996410891
0.73 -0.06644220532438613 -0.07424757758784524
0.74 -0.06717796330393419 -0.0735757979548059
0.75 -0.06790693059042321 -0.07289672864890166
0.76 -0.0686290349734149 -0.07221043829916751
0.77 -0.06934420493591298 -0.07151699624980838
0.78 -0.07005236966144492 -0.07081647255319393
0.79 -0.07075345904107275 -0.0701089379627832
0.8 -0.07144740368033255 -0.06939446392597987
0.81 -0.07213413490610174 -0.06867312257691885
0.8200000000000001 -0.07281358477339359 -0.0679449867291849
0.8300000000000001 -0.07348568607207823 -0.06721012986846413
0.84 -0.07415037233352952 -0.06646862614512883
0.85 -0.07480757783719709 -0.06572055036675685
0.86 -0.07545723761710293 -0.06496597799058582
0.87 -0.07609928746826196 -0.0642049851159032
0.88 -0.07673366395302568 -0.06343764847637294
0.89 -0.07736030440734867 -0.06266404543229945
0.9 -0.07797914694697697 -0.06188425396282968
0.91 -0.07859013047355791 -0.0610983526580941
0.92 -0.07919319468067078 -0.060306420711287397
0.93 -0.07978828005977767 -0.059508537910689614
0.9400000000000001 -0.08037532790609396 -0.058704784631628676
0.9500000000000001 -0.0809542803243778 -0.057895241828384895
0.96 -0.08152508023463818 -0.05707999102603851
0.97 -0.0820876713777608 -0.0562591143122609
0.98 -0.0826419983210513 -0.05543269432905038
0.99 -0.08318800646369544 -0.054600814264413426
https://ideone.com/kHA23g
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!