2015年12月27日日曜日

151227

Ruby


離散力学系における軌道の計算(1)

Rationalで計算すると大変なことになった。

def f(lambda, x)
  lambda * x * (1 - x)
end

N = 12
x = 1 / 3.to_r
p [x] + (1..N).map{x = f(4, x)}

出力結果
[(1/3), (8/9), (32/81), (6272/6561), (7250432/43046721), (1038154236987392/1853020188851841), (3383826162019367796397224108032/3433683820292512484657849089281), (674838593766753484487654913831820720085359667709963001167872/11790184577738583171520872861412518665678211592275841109096961), (30004257812217921448517555209678850165868237842249191697191933280958327232676392983676899191261907114558496045481444114432/139008452377144732764939786789661303114218850808529137991604824430036072629766435941001769154109609521811665540548899435521), (13082359825356910286492360061136736044180324498109817790128474991459314029228384187522308703184360391592495083267439167657929717062912601762329623472844083023729847910612245007924921595506191938616119697389548402100243733465436848681679275425792/19323349832288915105454068722019581055401465761603328550184537628902466746415537000017939429786029354390082329294586119505153509101332940884098040478728639542560550133727399482778062322407372338121043399668242276591791504658985882995272436541441), (326587507748564819526644454569376331474154420084818092745594332137561018974006085389560132992720304408663106445778503081330740875395946979612958407581856575179900521300573125323094574403291444363685870280626797364502957292617630209535691823652432585205978126151670505839933787574882669635101678944234697049051805139530490463071202426711136068343565226938926588080054391637396003154271711916707548475413061895252789685202028958639851477169623855847348893412975686499092668317042928917676032/373391848741020043532959754184866588225409776783734007750636931722079040617265251229993688938803977220468765065431475158108727054592160858581351336982809187314191748594262580938807019951956404285571818041046681288797402925517668012340617298396574731619152386723046235125934896058590588284654793540505936202376547807442730582144527058988756251452817793413352141920744623027518729185432862375737063985485319476416926263819972887006907013899256524297198527698749274196276811060702333710356481), (61142852306159761941719350858894288745543170044146323109689703291980924851796533997103855261908015575214297660329162668746580929102334010791256690499085996667702727608110600383881916661175401025841128350218665654536958173329093283646444406916989702358602504498990873573726862396234275982132488343451071715679594802189481207543737898095804664205350986023945673496564730371922126707688091178762709213338405622273864048733071827025714877491861045148462685583532205278217127011376328565887694222960373626180059928802236919349419068041418998434116106220266700986879354395786476811458677962796732587694611828952750818857523063965617113143694838317807352194745983507768142547294302047316329224316725717274586072428240970325222521933525573197993768182995402639191191324453820256624497823309547681905713129756569529235933403577032995411915451001642656423046719386092362811348998437241697365079123990708487380245174236936757338067180490209172293310254419238315917518121075773796329193472/139421472706236791468735287967015707232606232113998186759762284283203248267399327423485002476018294548027207658124580006956269122224761375615358997681237349433785401584729106693395255281491169142017082115307947771801954704567433211464560639140944673635364155016693389201613594038984686916178623702578352190636251776039742101384851474311083735102877993131212196962695464654469816719202970619252343026646014132485548307935847260900441555694943498191955672798596377459777392636575523873495086339994375985976543710570091438233257228420529594557457359609723639787519567621434999954843382100489400924995961904023259324301047965917783124275162880514542735528152228927483841372844810255470399471841348890379437250393401562476664019972812975660227449988177768002661436897984534184331079950399306756818061310948084534985374661203278436590349231431071569402389429768922180485838968764968523571261288038064310442585283933256404572748606640151307376800250725427496845607662821346676098703361), (19144712503407417567135059100458633701667309261608387794719701072085835139938108421290218174975016215714657693126232069759616256865782480902558184057624088768345367248914677728705523334349125924209805233089412155636967764580681507466229285675684099954946152999752510440892216943190218741400027174833668946677071863297897335955473049439573127667920418447654975692709588701834242287351039317582220066092409167756751240609167198807693453352489167398606261406447370695057064406265179059513400881523612801573817543697395645633573279331847873608996935144861669021453741792200011206868696965367089813764169742463266525029607727867386150092729846004573907465251484405890652377201376321459395896086671591945559497982887561609777822082124699779405715974263703118531857016577898505897381283351314353641134506307788914991746931629670739943951673167167541172574300926553165662576893798902217928419118126936416968892344654164860264317127874477993584853836862183965649667474310078569520451512245499722386379798234696618124800502349263039691768845705811932926646241786271093837316572273943077542327574039330422577103805612624251155221033086217586676029341508566076392648196659837797500695724040328070035578093492366392131689417025783967158724647867285956529352351500875068181118448882580408080438771351439926396969467728612146280442475235449540289675384105398167854906540445142977318423987318135885441075818173283389942560829552740288844103803489518606303098944536131089369039738367199369084865077552110767959726272052448395432901848203746166527324536275834701908222021024521049717000965869513618332488406972145618119726265214858358022562966536680948229538149778771163307213368496392381156029814746882202626326165357558439353593699685923189241967919428162720017105704489398074777319153750232775130376714213645787326097409226746826821996050111345129518976879954712719532335918095906765322748778286763898574001993813959499148858315486192477572653260400061077728656392978432/19438347051575930593026637277464327123266958586038369937823081254224453278630252497889372715554605593282381159872094158595095268067402968819922058060551822238741112165184374168564842547377441712304718992557393039038068650790181662446618501807237815659967187794033032803500035732383779283962952457019941126317411288547298348393098492554260561270997109425187315296683000918052454625638849801275961486965984453943973091802013130902022480064009582344873235465674303023321390779804986857735192246409512396859725643331241004565363356447570421115052411463010710895821609990116581002342162299225771869569017520145772647959223587680923278063396420933182948493511315641718579755044437441435175740859842860053607410165012361020931202956648767463066402467873454067221827562969068586208309387212456515799198207666199100272945024205012536374239885475968620882340375145207974046250462049126376560527290449130131440049788901462599304866873205740965010110139132541314345764876189581440737408106418855697237772150397259001202610130643763522631178062380060118593870971595168515254708733434522944375449004769094553723966227387818112951763132792281368178696361699983006097111869296586184614483820710005212223076649640648680453168817055219121408137639068898013300964241873653031196757067256269494984584288619038952494840890585035037799939467208769356259970350911013289050848552248531637772378630974861351496031252497009626563341895925393417896721361397607975167965296260228765252019053692679530118811766299692563834399397255178641804981662398606198713407736666476041324967934194436339007404569332651584643092199985299610486183872564728311389752139542321950489431449936370044181920096335116574849401654651260886859807973366610174486966928282281448812271499242148110214706385942214916820701400716027195159400919830048720354640368632751312692913386133718693853550810552709435439109936878872272169754302271563047123599935593148473062504803871685643230134991892294944498425472696321)]

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。