diff options
author | Erich Eckner <git@eckner.net> | 2015-08-04 10:12:25 +0200 |
---|---|---|
committer | Erich Eckner <git@eckner.net> | 2015-08-04 10:12:25 +0200 |
commit | 16be431a71ff4d05d20bfcdf55627a2562506768 (patch) | |
tree | 659dac1518a9652ba6301e58a94170501d3b8efd | |
parent | 0484afec8922e5ac455b98bcd5516d53443a475c (diff) | |
download | Plasmapropagation-16be431a71ff4d05d20bfcdf55627a2562506768.tar.xz |
script zur Koeffizientenkonversion erstellt: rktopas
-rw-r--r-- | rk108.inc | 331 | ||||
-rw-r--r-- | rk1210.inc | 577 | ||||
-rw-r--r-- | rk1412.inc | 985 | ||||
-rwxr-xr-x | rktopas | 33 |
4 files changed, 1926 insertions, 0 deletions
diff --git a/rk108.inc b/rk108.inc new file mode 100644 index 0000000..de1c24e --- /dev/null +++ b/rk108.inc @@ -0,0 +1,331 @@ + felders[2].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + 0.1 * dT + ); + felders[2].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(2,mDT,dT) then + continue; + + felders[3].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + -0.91517656137529144052001501927534215431895138766436972056466 * dT, + 1.45453440217827322805250021715664459117622483736537873607016 * dT + ); + felders[3].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(3,mDT,dT) then + continue; + + felders[4].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + 0.202259190301118170324681949205488413821477543637878380814562 * dT, + 0.606777570903354510974045847616465241464432630913635142443687 * dT + ); + felders[4].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(4,mDT,dT) then + continue; + + felders[5].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[4], + 0.184024714708643575149100693471120664216774047979591417844635 * dT, + 0.197966831227192369068141770510388793370637287463360401555746 * dT, + -0.0729547847313632629185146671595558023015011608914382961421311 * dT + ); + felders[5].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(5,mDT,dT) then + continue; + + felders[6].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + 0.0879007340206681337319777094132125475918886824944548534041378 * dT, + 0.410459702520260645318174895920453426088035325902848695210406 * dT, + 0.482713753678866489204726942976896106809132737721421333413261 * dT + ); + felders[6].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(6,mDT,dT) then + continue; + + felders[7].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + felders[6], + 0.085970050490246030218848022594580840141113261563660022259388 * dT, + 0.330885963040722183948884057658753173648240154838402033448632 * dT, + 0.48966295730945019284450701113589820117801547843379009721079 * dT, + -0.0731856375070850736789057580558988816340355615025188195854775 * dT + ); + felders[7].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(7,mDT,dT) then + continue; + + felders[8].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[5], + felders[6], + felders[7], + 0.120930449125333720660378854927668953958938996999703678812621 * dT, + 0.260124675758295622809007617838335174368108756484693361887839 * dT, + 0.0325402621549091330158899334391231259332716675992700000776101 * dT, + -0.0595780211817361001560122202563305121444953672762930724538856 * dT + ); + felders[8].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(8,mDT,dT) then + continue; + + felders[9].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + 0.110854379580391483508936171010218441909425780168656559807038 * dT, + -0.0605761488255005587620924953655516875526344415354339234619466 * dT, + 0.32176370560177839010089879904987890408140436860307712925111 * dT, + 0.51048572560806303157775901228512341674467213703175235406759 * dT + ); + felders[9].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(9,mDT,dT) then + continue; + + felders[10].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + 0.112054414752879004829715002761802363003717611158172229329393 * dT, + -0.144942775902865915672349828340980777181668499748506838876185 * dT, + -0.333269719096256706589705211415746871709467423992115497968724 * dT, + 0.499269229556880061353316843969978567860276816592673201240332 * dT, + 0.50950460892968610423609869004538625398664323235298960218506 * dT + ); + felders[10].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(10,mDT,dT) then + continue; + + felders[11].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + 0.113976783964185986138004186736901163890724752541486831640341 * dT, + -0.0768813364203356938586214289120895270821349023390922987406384 * dT, + 0.239527360324390649107711455271882373019741311201004119339563 * dT, + 0.397774662368094639047830462488952104564716416343454639902613 * dT, + 0.0107558956873607455550609147441477450257136782823280838547024 * dT, + -0.327769124164018874147061087350233395378262992392394071906457 * dT + ); + felders[11].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(11,mDT,dT) then + continue; + + felders[12].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + 0.0798314528280196046351426864486400322758737630423413945356284 * dT, + -0.0520329686800603076514949887612959068721311443881683526937298 * dT, + -0.0576954146168548881732784355283433509066159287152968723021864 * dT, + 0.19478191571210416497630626214738287115614292135440936473809 * dT, + 0.145384923188325069727524825977071194859203467568236523866582 * dT, + -0.0782942710351670777553986729725692447252077047239160551335016 * dT, + -0.114503299361098912184303164290554670970133218405658122674674 * dT + ); + felders[12].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(12,mDT,dT) then + continue; + + felders[13].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + felders[12], + 0.985115610164857280120041500306517278413646677314195559520529 * dT, + 0.330885963040722183948884057658753173648240154838402033448632 * dT, + 0.48966295730945019284450701113589820117801547843379009721079 * dT, + -1.37896486574843567582112720930751902353904327148559471526397 * dT, + -0.861164195027635666673916999665534573351026060987427093314412 * dT, + 5.78428813637537220022999785486578436006872789689499172601856 * dT, + 3.28807761985103566890460615937314805477268252903342356581925 * dT, + -2.38633905093136384013422325215527866148401465975954104585807 * dT, + -3.25479342483643918654589367587788726747711504674780680269911 * dT, + -2.16343541686422982353954211300054820889678036420109999154887 * dT + ); + felders[13].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(13,mDT,dT) then + continue; + + felders[14].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[4], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + 0.89508029577163289104961313233658513814815627924156134599171 * dT, + 0.197966831227192369068141770510388793370637287463360401555746 * dT, + -0.0729547847313632629185146671595558023015011608914382961421311 * dT, + -0.851236239662007619739049371445966793289359722875702227166105 * dT, + 0.398320112318533301719718614174373643336480918103773904231856 * dT, + 3.63937263181035606029412920047090044132027387893977804176229 * dT, + 1.54822877039830322365301663075174564919981736348973496313065 * dT, + -2.12221714704053716026062427460427261025318461146260124401561 * dT, + -1.58350398545326172713384349625753212757269188934434237975291 * dT, + -1.71561608285936264922031819751349098912615880827551992973034 * dT, + -0.0244036405750127452135415444412216875465593598370910566069132 * dT + ); + felders[14].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(14,mDT,dT) then + continue; + + felders[15].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[5], + felders[7], + felders[13], + felders[14], + -0.91517656137529144052001501927534215431895138766436972056466 * dT, + 1.45453440217827322805250021715664459117622483736537873607016 * dT, + -0.777333643644968233538931228575302137803351053629547286334469 * dT, + -0.0910895662155176069593203555807484200111889091770101799647985 * dT, + 0.0910895662155176069593203555807484200111889091770101799647985 * dT, + 0.777333643644968233538931228575302137803351053629547286334469 * dT + ); + felders[15].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(15,mDT,dT) then + continue; + + felders[16].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[15], + 0.1 * dT, + -0.157178665799771163367058998273128921867183754126709419409654 * dT, + 0.157178665799771163367058998273128921867183754126709419409654 * dT + ); + felders[16].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(16,mDT,dT) then + continue; + + felders[17].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[3], + felders[5], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + 0.181781300700095283888472062582262379650443831463199521664945 * dT, + 0.675 * dT, + 0.34275815984718983994222055341385087174233873470395891993726 * dT, + 0.259111214548322744512977076191767379267783684543182428778156 * dT, + -0.358278966717952089048961276721979397739750634673268802484271 * dT, + -1.0459489594088330609505006875640990513158812317237848928608 * dT, + 0.930327845415626983292300564432428777137601651182965794680397 * dT, + 1.77950959431708102446142106794824453926275743243327790536 * dT, + 0.1 * dT, + -0.282547569539044081612477785222287276408489375976211189952877 * dT, + -0.159327350119972549169261984373485859278031542127551931461821 * dT, + -0.145515894647001510860991961081084111308650130578626404945571 * dT, + -0.259111214548322744512977076191767379267783684543182428778156 * dT, + -0.34275815984718983994222055341385087174233873470395891993726 * dT, + -0.675 * dT + ); + felders[17].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(17,mDT,dT) then + continue; + + felders[1-aktuelleFelder].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[3], + felders[5], + felders[7], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + 0.0333333333333333333333333333333333333333333333333333333333333 * dT, + 0.025 * dT, + 0.0333333333333333333333333333333333333333333333333333333333333 * dT, + 0.05 * dT, + 0.04 * dT, + 0.189237478148923490158306404106012326238162346948625830327194 * dT, + 0.277429188517743176508360262560654340428504319718040836339472 * dT, + 0.277429188517743176508360262560654340428504319718040836339472 * dT, + 0.189237478148923490158306404106012326238162346948625830327194 * dT, + -0.04 * dT, + -0.05 * dT, + -0.0333333333333333333333333333333333333333333333333333333333333 * dT, + -0.025 * dT, + 0.0333333333333333333333333333333333333333333333333333333333333 * dT + ); diff --git a/rk1210.inc b/rk1210.inc new file mode 100644 index 0000000..3410f78 --- /dev/null +++ b/rk1210.inc @@ -0,0 +1,577 @@ + felders[2].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + 0.2 * dT + ); + felders[2].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(2,mDT,dT) then + continue; + + felders[3].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + -0.216049382716049382716049382716049382716049382716049382716049 * dT, + 0.771604938271604938271604938271604938271604938271604938271605 * dT + ); + felders[3].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(3,mDT,dT) then + continue; + + felders[4].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + 0.208333333333333333333333333333333333333333333333333333333333 * dT, + 0.625 * dT + ); + felders[4].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(4,mDT,dT) then + continue; + + felders[5].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[4], + 0.193333333333333333333333333333333333333333333333333333333333 * dT, + 0.22 * dT, + -0.08 * dT + ); + felders[5].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(5,mDT,dT) then + continue; + + felders[6].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + 0.1 * dT, + 0.4 * dT, + 0.5 * dT + ); + felders[6].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(6,mDT,dT) then + continue; + + felders[7].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + felders[6], + 0.103364471650010477570395435690481791543342708330349879244197 * dT, + 0.124053094528946761061581889237115328211074784955180298044074 * dT, + 0.48317116756103289928883648045196250872410925751728917730238 * dT, + -0.038753024569476325208568144376762058039573330234136803880429 * dT + ); + felders[7].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(7,mDT,dT) then + continue; + + felders[8].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[5], + felders[6], + felders[7], + 0.12403826143183332408190458598017516814002467069863361229248 * dT, + 0.217050632197958486317846256953159942875916353757734167684657 * dT, + 0.0137455792075966759812907801835048190594443990939408530842918 * dT, + -0.0661095317267682844455831341498149531672668252085016565917546 * dT + ); + felders[8].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(8,mDT,dT) then + continue; + + felders[9].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + 0.0914774894856882983144991846980432197088832099976660100090486 * dT, + -0.00544348523717469689965754944144838611346156873847009178068318 * dT, + 0.0680716801688453518578515120895103863112751730758794372203952 * dT, + 0.408394315582641046727306852653894780093303185664924644551239 * dT + ); + felders[9].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(9,mDT,dT) then + continue; + + felders[10].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + 0.0890013652502551018954509355423841780143232697403434118692699 * dT, + 0.0049952822664553236019779340842069280040589114940681409195581 * dT, + 0.397918238819828997341739603001347156083435060931424970826304 * dT, + 0.427930210752576611068192608300897981558240730580396406312359 * dT, + -0.0865117637557827005740277475955029103267246394128995965941585 * dT + ); + felders[10].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(10,mDT,dT) then + continue; + + felders[11].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + 0.0695087624134907543112693906409809822706021061685544615255758 * dT, + 0.129146941900176461970759579482746551122871751501482634045487 * dT, + 1.53073638102311295076342566143214939031177504112433874313011 * dT, + 0.577874761129140052546751349454576715334892100418571882718036 * dT, + -0.951294772321088980532340837388859453930924498799228648050949 * dT, + -0.408276642965631951497484981519757463459627174520978426909934 * dT + ); + felders[11].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(11,mDT,dT) then + continue; + + felders[12].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + 0.0444861403295135866269453507092463581620165501018684152933313 * dT, + -0.00380476867056961731984232686574547203016331563626856065717964 * dT, + 0.0106955064029624200721262602809059154469206077644957399593972 * dT, + 0.0209616244499904333296674205928919920806734650660039898074652 * dT, + -0.0233146023259321786648561431551978077665337818756053603898847 * dT, + 0.00263265981064536974369934736325334761174975280887405725010964 * dT, + 0.00315472768977025060103545855572111407955208306374459723959783 * dT + ); + felders[12].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(12,mDT,dT) then + continue; + + felders[13].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + 0.0194588815119755475588801096525317761242073762016273186231215 * dT, + 0.0000678512949171812509306121653452367476194364781259165332321534 * dT, + -0.000042979585904927362327100533023016234356886338772488360367555 * dT, + 0.0000176358982260285155407485928953302139937553442829975734148981 * dT, + 0.0653866627415027051009595231385181033549511358787382098351924 * dT + ); + felders[13].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(13,mDT,dT) then + continue; + + felders[14].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + 0.206836835664277105916828174798272361078909196043446411598231 * dT, + 0.0166796067104156472828045866664696450306326505094792505215514 * dT, + -0.00879501563200710214457024178249986591130234990219959208704979 * dT, + 0.00346675455362463910824462315246379209427513654098596403637231 * dT, + -0.861264460105717678161432562258351242030270498966891201799225 * dT, + 0.908651882074050281096239478469262145034957129939256789178785 * dT + ); + felders[14].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(14,mDT,dT) then + continue; + + felders[15].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + 0.0203926084654484010091511314676925686038504449562413004562382 * dT, + 0.0869469392016685948675400555583947505833954460930940959577347 * dT, + -0.0191649630410149842286436611791405053287170076602337673587681 * dT, + 0.00655629159493663287364871573244244516034828755253746024098838 * dT, + 0.0987476128127434780903798528674033899738924968006632201445462 * dT, + 0.00535364695524996055083260173615567408717110247274021056118319 * dT, + 0.301167864010967916837091303817051676920059229784957479998077 * dT + ); + felders[15].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(15,mDT,dT) then + continue; + + felders[16].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + 0.228410433917778099547115412893004398779136994596948545722283 * dT, + -0.498707400793025250635016567442511512138603770959682292383042 * dT, + 0.134841168335724478552596703792570104791700727205981058201689 * dT, + -0.038745824405583415843990422692402923093516105914280680567436 * dT, + -1.27473257473474844240388430824908952380979292713250350199641 * dT, + 1.43916364462877165201184452437038081875299303577911839630524 * dT, + -0.214007467967990254219503540827349569639028092344812795499026 * dT, + 0.958202417754430239892724139109781371059908874605153648768037 * dT + ); + felders[16].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(16,mDT,dT) then + continue; + + felders[17].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + 2.00222477655974203614249646012506747121440306225711721209798 * dT, + 2.06701809961524912091954656438138595825411859673341600679555 * dT, + 0.62397813608613954195747127983149446615529231616702108066314 * dT, + -0.0462283685500311430283203554129062069391947101880112723185773 * dT, + -8.8497328836264961486007524672711894928660483545709270109463 * dT, + 7.74257707850855976227437225791835589560188590785037197433615 * dT, + -0.588358519250869210993353314127711745644125882130941202896436 * dT, + -1.10683733362380649395704708016953056176195769617014899442903 * dT, + -0.929529037579203999778397238291233214220788057511899747507074 * dT + ); + felders[17].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(17,mDT,dT) then + continue; + + felders[18].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + 3.1378953341207344293445160898988879680816125933032210026831 * dT, + 0.129146941900176461970759579482746551122871751501482634045487 * dT, + 1.53073638102311295076342566143214939031177504112433874313011 * dT, + 0.577874761129140052546751349454576715334892100418571882718036 * dT, + 5.42088263055126683050056840891857421941300558851862156403363 * dT, + 0.231546926034829304872663800877643660904880180835945693836936 * dT, + 0.0759292995578913560162301311785251873561801342333194895292058 * dT, + -12.3729973380186513287414553402595806591349822617535905976253 * dT, + 9.85455883464769543935957209317369202080367765721777101906955 * dT, + 0.085911143137043652957935770905236777288998049512232960115954 * dT, + -5.65242752862643921117182090081762761180392602644189218673969 * dT, + -1.94300935242819610883833776782364287728724899124166920477873 * dT, + -0.128352601849404542018428714319344620742146491335612353559923 * dT + ); + felders[18].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(18,mDT,dT) then + continue; + + felders[19].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + 1.38360054432196014878538118298167716825163268489922519995564 * dT, + 0.0049952822664553236019779340842069280040589114940681409195581 * dT, + 0.397918238819828997341739603001347156083435060931424970826304 * dT, + 0.427930210752576611068192608300897981558240730580396406312359 * dT, + -1.30299107424475770916551439123047573342071475998399645982146 * dT, + 0.661292278669377029097112528107513072734573412294008071500699 * dT, + -0.144559774306954349765969393688703463900585822441545655530145 * dT, + -6.96576034731798203467853867461083919356792248105919255460819 * dT, + 6.65808543235991748353408295542210450632193197576935120716437 * dT, + -1.66997375108841486404695805725510845049807969199236227575796 * dT, + 2.06413702318035263832289040301832647130604651223986452170089 * dT, + -0.674743962644306471862958129570837723192079875998405058648892 * dT, + -0.00115618834794939500490703608435907610059605754935305582045729 * dT, + -0.00544057908677007389319819914241631024660726585015012485938593 * dT + ); + felders[19].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(19,mDT,dT) then + continue; + + felders[20].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[5], + felders[6], + felders[7], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + 0.951236297048287669474637975894973552166903378983475425758226 * dT, + 0.217050632197958486317846256953159942875916353757734167684657 * dT, + 0.0137455792075966759812907801835048190594443990939408530842918 * dT, + -0.0661095317267682844455831341498149531672668252085016565917546 * dT, + 0.152281696736414447136604697040747131921486432699422112099617 * dT, + -0.33774101835759984080230079313399800435464342445753966767008 * dT, + -0.0192825981633995781534949199286824400469353110630787982121133 * dT, + -3.68259269696866809932409015535499603576312120746888880201882 * dT, + 3.16197870406982063541533528419683854018352080342887002331312 * dT, + -0.370462522106885290716991856022051125477943482284080569177386 * dT, + -0.0514974200365440434996434456698127984941168616474316871020314 * dT, + -0.000829625532120152946787043541792848416659382675202720677536554 * dT, + 0.00000279801041419278598986586589070027583961355402640879503213503 * dT, + 0.0418603916412360287969841020776788461794119440689356178942252 * dT, + 0.279084255090877355915660874555379649966282167560126269290222 * dT + ); + felders[20].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(20,mDT,dT) then + continue; + + felders[21].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + felders[6], + felders[8], + felders[10], + felders[11], + felders[18], + felders[19], + felders[20], + 0.103364471650010477570395435690481791543342708330349879244197 * dT, + 0.124053094528946761061581889237115328211074784955180298044074 * dT, + 0.48317116756103289928883648045196250872410925751728917730238 * dT, + -0.038753024569476325208568144376762058039573330234136803880429 * dT, + -0.438313820361122420391059788940960176420682836652600698580091 * dT, + -0.218636633721676647685111485017151199362509373698288330593486 * dT, + -0.0312334764394719229981634995206440349766174759626578122323015 * dT, + 0.0312334764394719229981634995206440349766174759626578122323015 * dT, + 0.218636633721676647685111485017151199362509373698288330593486 * dT, + 0.438313820361122420391059788940960176420682836652600698580091 * dT + ); + felders[21].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(21,mDT,dT) then + continue; + + felders[22].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[4], + felders[7], + felders[8], + felders[10], + felders[11], + felders[18], + felders[19], + felders[20], + felders[21], + 0.193333333333333333333333333333333333333333333333333333333333 * dT, + 0.22 * dT, + -0.08 * dT, + 0.0984256130499315928152900286856048243348202521491288575952143 * dT, + -0.196410889223054653446526504390100417677539095340135532418849 * dT, + 0.436457930493068729391826122587949137609670676712525034763317 * dT, + 0.065261372167572109856037093980555569835054381070841471673027 * dT, + -0.065261372167572109856037093980555569835054381070841471673027 * dT, + -0.436457930493068729391826122587949137609670676712525034763317 * dT, + 0.196410889223054653446526504390100417677539095340135532418849 * dT, + -0.0984256130499315928152900286856048243348202521491288575952143 * dT + ); + felders[22].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(22,mDT,dT) then + continue; + + felders[23].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[5], + felders[7], + felders[21], + felders[22], + -0.216049382716049382716049382716049382716049382716049382716049 * dT, + 0.771604938271604938271604938271604938271604938271604938271605 * dT, + -0.666666666666666666666666666666666666666666666666666666666667 * dT, + -0.390696469295978451446999802258495981249099665294395945559163 * dT, + 0.390696469295978451446999802258495981249099665294395945559163 * dT, + 0.666666666666666666666666666666666666666666666666666666666667 * dT + ); + felders[23].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(23,mDT,dT) then + continue; + + felders[24].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[23], + 0.2 * dT, + -0.164609053497942386831275720164609053497942386831275720164609 * dT, + 0.164609053497942386831275720164609053497942386831275720164609 * dT + ); + felders[24].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(24,mDT,dT) then + continue; + + felders[25].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[3], + felders[5], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + felders[24], + 1.47178724881110408452949550989023611293535315518571691939396 * dT, + 0.7875 * dT, + 0.421296296296296296296296296296296296296296296296296296296296 * dT, + 0.291666666666666666666666666666666666666666666666666666666667 * dT, + 0.348600717628329563206854421629657569274689947367847465753757 * dT, + 0.22949954476899484958289023371055544707382356966650670066251 * dT, + 5.79046485790481979159831978177003471098279506036722411333192 * dT, + 0.418587511856506868874073759426596207226461447604248151080016 * dT, + 0.307039880222474002649653817490106690389251482313213999386651 * dT, + -4.68700905350603332214256344683853248065574415794742040470287 * dT, + 3.13571665593802262152038152399873856554395436199962915429076 * dT, + 1.40134829710965720817510506275620441055845017313930508348898 * dT, + -5.5293110143949902362901030600576433642127605577765815640091 * dT, + -0.853138235508063349309546894974784906188927508039552519557498 * dT, + 0.103575780373610140411804607167772795518293914458500175573749 * dT, + -0.140474416950600941142546901202132534870665923700034957196546 * dT, + -0.418587511856506868874073759426596207226461447604248151080016 * dT, + -0.22949954476899484958289023371055544707382356966650670066251 * dT, + -0.348600717628329563206854421629657569274689947367847465753757 * dT, + -0.291666666666666666666666666666666666666666666666666666666667 * dT, + -0.421296296296296296296296296296296296296296296296296296296296 * dT, + -0.7875 * dT + ); + felders[25].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(25,mDT,dT) then + continue; + + felders[1-aktuelleFelder].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[3], + felders[5], + felders[7], + felders[8], + felders[10], + felders[11], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + felders[24], + felders[25], + 0.0238095238095238095238095238095238095238095238095238095238095 * dT, + 0.0234375 * dT, + 0.03125 * dT, + 0.0416666666666666666666666666666666666666666666666666666666667 * dT, + 0.05 * dT, + 0.05 * dT, + 0.1 * dT, + 0.0714285714285714285714285714285714285714285714285714285714286 * dT, + 0.138413023680782974005350203145033146748813640089941234591267 * dT, + 0.215872690604931311708935511140681138965472074195773051123019 * dT, + 0.24380952380952380952380952380952380952380952380952380952381 * dT, + 0.215872690604931311708935511140681138965472074195773051123019 * dT, + 0.138413023680782974005350203145033146748813640089941234591267 * dT, + -0.0714285714285714285714285714285714285714285714285714285714286 * dT, + -0.1 * dT, + -0.05 * dT, + -0.05 * dT, + -0.0416666666666666666666666666666666666666666666666666666666667 * dT, + -0.03125 * dT, + -0.0234375 * dT, + 0.0238095238095238095238095238095238095238095238095238095238095 * dT + ); diff --git a/rk1412.inc b/rk1412.inc new file mode 100644 index 0000000..c1334c1 --- /dev/null +++ b/rk1412.inc @@ -0,0 +1,985 @@ + felders[2].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + 0.111111111111111111111111111111111111111111111111111111111111 * dT + ); + felders[2].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(2,mDT,dT) then + continue; + + felders[3].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + -0.833333333333333333333333333333333333333333333333333333333333 * dT, + 1.38888888888888888888888888888888888888888888888888888888889 * dT + ); + felders[3].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(3,mDT,dT) then + continue; + + felders[4].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + 0.208333333333333333333333333333333333333333333333333333333333 * dT, + 0.625 * dT + ); + felders[4].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(4,mDT,dT) then + continue; + + felders[5].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[4], + 0.193333333333333333333333333333333333333333333333333333333333 * dT, + 0.22 * dT, + -0.08 * dT + ); + felders[5].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(5,mDT,dT) then + continue; + + felders[6].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + 0.1 * dT, + 0.4 * dT, + 0.5 * dT + ); + felders[6].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(6,mDT,dT) then + continue; + + felders[7].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + felders[6], + 0.103484561636679776672993546511910344499744798201971316606663 * dT, + 0.122068887306407222589644082868962077139592714834162134741275 * dT, + 0.482574490331246622475134780125688112865919023850168049679402 * dT, + -0.038140960001560699973088624000562020566411307247841147742197 * dT + ); + felders[7].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(7,mDT,dT) then + continue; + + felders[8].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[5], + felders[6], + felders[7], + 0.124380526654094412881516420868799316268491466359671423163289 * dT, + 0.226120282197584301422238662979202901196752320742633143965145 * dT, + 0.0137885887618080880607695837016477814530969417491493385363543 * dT, + -0.06722101339966844497493995074143058569500863415253821828562 * dT + ); + felders[8].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(8,mDT,dT) then + continue; + + felders[9].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + 0.0936919065659673815530885456083005933866349695217750085655603 * dT, + -0.00613406843450510987229498995641664735620914507128858871007099 * dT, + 0.216019825625503063708860097659866573490979433278117320188668 * dT, + 0.423695063515761937337619073960976753205867469544123532683116 * dT + ); + felders[9].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(9,mDT,dT) then + continue; + + felders[10].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + 0.0838479812409052664616968791372814085980533139224911131069335 * dT, + -0.0117949367100973814319755056031295775367961960590736150777613 * dT, + -0.247299020568812652339473838743194598325992840353340132697498 * dT, + 0.0978080858367729012259313014081291665503740655476733940756599 * dT, + 0.217590689243420631360008651767860318344168120024782176879989 * dT + ); + felders[10].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(10,mDT,dT) then + continue; + + felders[11].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + 0.061525535976942822795456238961431471433342396906482110745394 * dT, + 0.0059223278032450330804299000579804652473838956044425713683499 * dT, + 0.470326159963841112217224303205894113455362530746108825010848 * dT, + 0.299688863848679000853981837096192399136831121671781279184194 * dT, + -0.247656877593994914689992276329810825853958069263947095548189 * dT, + 0.110895029771437682893999851839061714522445173600678718208625 * dT + ); + felders[11].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(11,mDT,dT) then + continue; + + felders[12].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + 0.0419700073362782579861792864787277787213483656543104611245994 * dT, + -0.00317987696266205093901912847692712407988609169703103952205634 * dT, + 0.806397714906192077260821711520379506393543111567419750119748 * dT, + 0.0975983126412388979093522850684288851314672048003054550357187 * dT, + 0.778575578158398909027512446452927238999763460594181964958853 * dT, + 0.204890423831599428189499202098105603312029235081420653574829 * dT, + -1.56261579627468188307070943950527825211462892236424360892806 * dT + ); + felders[12].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(12,mDT,dT) then + continue; + + felders[13].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + 0.0437726782233730163574465242495339811688214967071614123256973 * dT, + 0.0062436502752019520879435862858093362528163121690309591720125 * dT, + 0.200043097109577314994435165469647856829066232218264969608768 * dT, + -0.00805328367804983036823857162048902911923392887337029314844206 * dT, + 0.0211517528067396521915711903523399601316877825157550573051221 * dT + ); + felders[13].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(13,mDT,dT) then + continue; + + felders[14].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + 0.0283499250363514563095023591920717312247137654896477097768495 * dT, + 0.00249163204855817407538949148805995149459884653585417680098222 * dT, + 0.0230138787854593149638399846373742768772087122638142234223658 * dT, + -0.00322155956692977098724476092467120878189463604760620461043308 * dT, + 0.00988442549447664668946335414487885256040819982786014648129297 * dT, + -0.0213010771328887351384307642875927384886634565429572466632092 * dT + ); + felders[14].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(14,mDT,dT) then + continue; + + felders[15].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + 0.343511894290243001049432234735147943083353174980701426268122 * dT, + 0.210451912023627385609097011999010655788807405225626700040882 * dT, + 1.03427452057230411936482926828825709938667999698324740166559 * dT, + 0.00600303645864422487051240448206640574939078092406156945568306 * dT, + 0.855938125099619537578012106002407728915062652616416005816477 * dT, + -0.977235005036766810872264852372525633013107656892839677696022 * dT, + -0.660026980479294694616225013856327693720573981219974874776419 * dT + ); + felders[15].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(15,mDT,dT) then + continue; + + felders[16].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + -0.0143574001672168069538206399935076366657755954378399880691949 * dT, + -0.0366253270049039970293685796848974791733119081733552207318285 * dT, + 0.0350254975636213681976849406979846524346789082471103574920148 * dT, + 0.0360946016362113508931786658758335239823689929864237671348749 * dT, + -0.0265219967553681106351595946834601923649627012457464284442911 * dT, + 0.0445699011305698119638911537508839908104336323082226770910408 * dT, + 0.124343093331358243286225595741786448038973408895106741855721 * dT, + 0.00413829693239480694403512496204335960426192908674476033832967 * dT + ); + felders[16].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(16,mDT,dT) then + continue; + + felders[17].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + 0.356032404425120290975609116398089176264106222379748802654822 * dT, + -0.450192758947562595966821779075956175110645100214763601190349 * dT, + 0.430527907083710898626656292808782917793030154094709462877146 * dT, + 0.511973029011022237668556960394071692077125787030651386389972 * dT, + 0.908303638886404260390159124638110213997496214819904630546596 * dT, + -1.2392109337193393175737246915153402885441388924860572618652 * dT, + -0.649048661671761465141672348879062553905402831967191097656668 * dT, + 0.251708904586819292210480529948970541404887852931447491219418 * dT, + 0.779906470345586398810756795282334476023540593411550187024263 * dT + ); + felders[17].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(17,mDT,dT) then + continue; + + felders[18].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + 0.0130935687406513066406881206418834980127470438213192487844956 * dT, + -0.0000932053067985113945908461962767108237858631509684667142124826 * dT, + 0.0505374334262299359640090443138590726770942344716122381702746 * dT, + 8.04470341944487979109579109610197797641311868930865361048975E-7 * dT, + 0.000591726029494171190528755742777717259844340971924321528178248 * dT, + -4.01614722154557337064691684906375587732264247950093804676867E-7 * dT + ); + felders[18].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(18,mDT,dT) then + continue; + + felders[19].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + 0.0207926484466053012541944544000765652167255206144373407979758 * dT, + 0.000582695918800085915101902697837284108951406103029871570103075 * dT, + -0.00801700732358815939083342186525852746640558465919633524655451 * dT, + 4.03847643847136940375170821743560570484117290330895506618968E-6 * dT, + 0.0854609998055506144225056114567535602510114622033622491802597 * dT, + -2.04486480935804242706707569691004307904442837552677456232848E-6 * dT, + 0.105328578824431893399799402979093997354240904235172843146582 * dT + ); + felders[19].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(19,mDT,dT) then + continue; + + felders[20].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + 1.40153449795736021415446247355771306718486452917597731683689 * dT, + -0.230252000984221261616272410367415621261130298274455611733277 * dT, + -7.21106840466912905659582237106874247165856493509961561958267 * dT, + 0.00372901560694836335236995327852132340217759566678662385552634 * dT, + -4.71415495727125020678778179392224757011323373221820091641216 * dT, + -0.00176367657545349242053841995032797673574903886695600132759652 * dT, + 7.64130548038698765563029310880237651185173367813936997648198 * dT, + 3.50602043659751834989896082949744710968212949893375368243588 * dT + ); + felders[20].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(20,mDT,dT) then + continue; + + felders[21].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + 11.9514650694120686799372385830716401674473610826553517297976 * dT, + 7.79480932108175968783516700231764388220284279598980948538579 * dT, + -56.4501393867325792523560991120904281440468100061340556540132 * dT, + 0.0912376306930644901344530449290276645709607450403673704844997 * dT, + -12.7336279925434886201945524309199275038162717529918963305155 * dT, + -0.0396895921904719712313542810939736674712383070433147873009352 * dT, + 54.4392141883570886996225765155307791861438378423305337073797 * dT, + -3.64411637921569236846406990361350645806721478409266709351203 * dT, + -0.804503249910509910899030787958579499315694913210787878260459 * dT + ); + felders[21].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(21,mDT,dT) then + continue; + + felders[22].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + -148.80942650710048842783886826864762556193061208214859707669 * dT, + -91.7295278291256484357935662402321623495228729036354276506427 * dT, + 707.656144971598359834575719286335716154821128966649565194286 * dT, + -1.1056361185748244090529696131159093080133830894263776955554 * dT, + 176.134591883811372587859898076055660406999516762301689616841 * dT, + 0.491384824214880662268898345164454557416884631402764792538746 * dT, + -684.278000449814944358237535610895081956077167893600278300805 * dT, + 27.9910604998398258984224332124380407446002518400668657974589 * dT, + 13.1939710030282333443670964371153238435064159623744975073252 * dT, + 1.25128781283980445450114974148056006317268830077396406361417 * dT + ); + felders[22].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(22,mDT,dT) then + continue; + + felders[23].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + -9.67307946948196763644126118433219395839951408571877262880482 * dT, + -4.46990150858505531443846227701960360497830681408751431146712 * dT, + 45.5127128690952681968241950400052751178905907817398483534845 * dT, + -0.0713085086183826912791492024438246129930559805352394367050813 * dT, + 11.227361406841274158259062447993938420782680077679448505154 * dT, + 0.126244376717622724516237912909138809361786889819105426371393 * dT, + -43.5439339549483313605810624907242107623814304467621407753424 * dT, + 0.787174307543058978398792994996550902064546091443233850464377 * dT, + 0.532264696744684215669300708603886690785395776821503851830821 * dT, + 0.422422733996325326010225127471388772575086538809603346825334 * dT, + 0.0859131249503067107308438031499859443441115056294154956487671 * dT + ); + felders[23].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(23,mDT,dT) then + continue; + + felders[24].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + -10.0664032447054702403396606900426891472202824757968765569183 * dT, + -0.0366253270049039970293685796848974791733119081733552207318285 * dT, + 0.0350254975636213681976849406979846524346789082471103574920148 * dT, + 0.0360946016362113508931786658758335239823689929864237671348749 * dT, + -0.0265219967553681106351595946834601923649627012457464284442911 * dT, + -6.27088972181464143590553149478871603839356122957396018530209 * dT, + 48.2079237442562989090702103008195063923492593141636117832993 * dT, + -0.0694471689136165640882395180583732834557754169149088630301342 * dT, + 12.6810690204850295698341370913609807066108483811412127009785 * dT, + 0.0119671168968323754838161435501011294100927813964199613229864 * dT, + -46.7249764992482408003358268242662695593201321659795608950429 * dT, + 1.33029613326626711314710039298216591399033511191227101321435 * dT, + 1.0076678750339829835343890361992665777116271779366171970837 * dT, + 0.0209512051933665091664122388475480702892770753864487241177616 * dT, + 0.0210134706331264177317735424331396407424412188443757490871603 * dT, + 0.00952196014417121794175101542454575907376360233658356240547761 * dT + ); + felders[24].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(24,mDT,dT) then + continue; + + felders[25].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + felders[24], + -409.478081677743708772589097409370357624424341606752069725341 * dT, + 0.210451912023627385609097011999010655788807405225626700040882 * dT, + 1.03427452057230411936482926828825709938667999698324740166559 * dT, + 0.00600303645864422487051240448206640574939078092406156945568306 * dT, + 0.855938125099619537578012106002407728915062652616416005816477 * dT, + -250.51699854744786049277765772931613038658405042078207596699 * dT, + 1946.42466652388427766053750328264758595829850895761428240231 * dT, + -3.04503882102310365506105809086860882786950544097602101685174 * dT, + 490.62637952828171352120826529916808384159854227406167157623 * dT, + 1.56647589531270907115484067013597445739595615245966775329993 * dT, + -1881.97428994011173362217267377035870619215906638453056643641 * dT, + 75.2592224724847175278837713643303149821620618914245864351135 * dT, + 34.5734356980331067622434344736554689696728644793551014989002 * dT, + 3.21147679440968961435417361847073755169022966748891627882572 * dT, + -0.460408041738414391307201404237058848867245095265382820823055 * dT, + -0.0870718339841810522431884137957986245724252047388936572215438 * dT, + -7.39351814158303067567016952195521063999185773249132944724553 * dT + ); + felders[25].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(25,mDT,dT) then + continue; + + felders[26].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[9], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + felders[24], + felders[25], + 3.43347475853550878921093496257596781120623891072008459930197 * dT, + 0.00249163204855817407538949148805995149459884653585417680098222 * dT, + 0.0230138787854593149638399846373742768772087122638142234223658 * dT, + -0.00322155956692977098724476092467120878189463604760620461043308 * dT, + 0.00988442549447664668946335414487885256040819982786014648129297 * dT, + 2.1625279937792250778830784190475735404575922533573270791653 * dT, + -16.2699864546457421328065640660139489006987552040228852402716 * dT, + -0.128534502120524552843583417470935010538029037542654506231743 * dT, + -8.98915042666504253089307820833379330486511746063552853023189 * dT, + -0.00348595363232025333387080201851013650192401767250513765000963 * dT, + 15.7936194113339807536235187388695574135853387025139738341334 * dT, + -0.574403330914095065628165482017335820148383663195675408024658 * dT, + -0.345602039021393296692722496608124982535237228827655306030152 * dT, + -0.00662241490206585091731619991383757781133067992707418687587487 * dT, + -0.0077778812924220416403254645860736430975934720962675911194615 * dT, + -0.00356084192402274913338827232697437364675240818791706587952939 * dT, + 4.79282506449930799649797749629840189457296934139359048988332 * dT, + 0.153725464873068577844576387402512082757034273069877432944621 * dT + ); + felders[26].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(26,mDT,dT) then + continue; + + felders[27].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[11], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + felders[24], + felders[25], + felders[26], + 32.3038520871985442326994734440031535091364975047784630088983 * dT, + -0.00317987696266205093901912847692712407988609169703103952205634 * dT, + 0.806397714906192077260821711520379506393543111567419750119748 * dT, + 0.0975983126412388979093522850684288851314672048003054550357187 * dT, + 0.778575578158398909027512446452927238999763460594181964958853 * dT, + 0.204890423831599428189499202098105603312029235081420653574829 * dT, + -1.56261579627468188307070943950527825211462892236424360892806 * dT, + 16.3429891882310570648504243973927174708753353504154550405647 * dT, + -154.544555293543621230730189631471036399316683669609116705323 * dT, + 1.56971088703334872692034283417621761466263593582497085955201 * dT, + 3.27685545087248131321429817269900731165522404974733504794135 * dT, + -0.0503489245193653176348040727199783626534081095691632396802451 * dT, + 153.321151858041665070593767885914694011224363102594556731397 * dT, + 7.17568186327720495846766484814784143567826308034865369443637 * dT, + -2.94036748675300481945917659896930989215320594380777597403592 * dT, + -0.0665845946076803144470749676022628870281920493197256887985612 * dT, + -0.0462346054990843661229248668562217261176966514016859284197145 * dT, + -0.0204198733585679401539388228617269778848579774821581777675337 * dT, + -53.3523106438735850515953441165998107974045090495791591218714 * dT, + -1.35548714715078654978732186705996404017554501614191325114947 * dT, + -1.57196275801232751882901735171459249177687219114442583461866 * dT + ); + felders[27].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(27,mDT,dT) then + continue; + + felders[28].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[10], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + felders[24], + felders[25], + felders[26], + felders[27], + -16.6451467486341512872031294403931758764560371130818978459405 * dT, + 0.0059223278032450330804299000579804652473838956044425713683499 * dT, + 0.470326159963841112217224303205894113455362530746108825010848 * dT, + 0.299688863848679000853981837096192399136831121671781279184194 * dT, + -0.247656877593994914689992276329810825853958069263947095548189 * dT, + 0.110895029771437682893999851839061714522445173600678718208625 * dT, + -0.491719043846229147070666628704194097678081907210673044988866 * dT, + -11.4743154427289496968389492564352536350842454130853175250727 * dT, + 80.2593166576230272541702485886484400152793366623589989106256 * dT, + -0.384132303980042847625312526759029103746926841342088219165648 * dT, + 7.28147667468107583471326950926136115767612581862877764249646 * dT, + -0.132699384612248379510571708176035274836827341616751884314074 * dT, + -81.0799832525730726674679289752255240006070716633632990308935 * dT, + -1.25037492835620639521768185656179119962253747492403205797494 * dT, + 2.59263594969543681023776379504377324994226447359296887778718 * dT, + -0.301440298346404539830163997260526875264431537275641495291993 * dT, + 0.22138446078983233745170645157277379169524683905731841430102 * dT, + 0.082757727477189293195598987097469315299627643542980989055121 * dT, + 18.9960662040611520464672450037243263998175161412237156872211 * dT, + 0.269231946409639685623468015128334167460051910348912845121977 * dT, + 1.6267482744706653746298936492962893398812502928418368027902 * dT, + 0.491719043846229147070666628704194097678081907210673044988866 * dT + ); + felders[28].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(28,mDT,dT) then + continue; + + felders[29].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[6], + felders[7], + felders[8], + felders[9], + felders[11], + felders[12], + felders[14], + felders[15], + felders[16], + felders[24], + felders[25], + felders[26], + felders[27], + felders[28], + 0.0838479812409052664616968791372814085980533139224911131069335 * dT, + -0.0117949367100973814319755056031295775367961960590736150777613 * dT, + -0.247299020568812652339473838743194598325992840353340132697498 * dT, + 0.0978080858367729012259313014081291665503740655476733940756599 * dT, + 0.217590689243420631360008651767860318344168120024782176879989 * dT, + 0.13758560676332522486565963219678774664744722297508486597544 * dT, + 0.0439870229715046685058790092341545026046103890294261359042581 * dT, + -0.513700813768193341957004456618630303738757363641964030086972 * dT, + 0.826355691151315508644211308399153458701423158616168576922372 * dT, + 25.7018139719811832625873882972519939511136556341960074626615 * dT, + -25.7018139719811832625873882972519939511136556341960074626615 * dT, + -0.826355691151315508644211308399153458701423158616168576922372 * dT, + 0.513700813768193341957004456618630303738757363641964030086972 * dT, + -0.0439870229715046685058790092341545026046103890294261359042581 * dT, + -0.13758560676332522486565963219678774664744722297508486597544 * dT + ); + felders[29].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(29,mDT,dT) then + continue; + + felders[30].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[5], + felders[6], + felders[7], + felders[10], + felders[11], + felders[12], + felders[14], + felders[15], + felders[16], + felders[24], + felders[25], + felders[26], + felders[27], + felders[28], + felders[29], + 0.124380526654094412881516420868799316268491466359671423163289 * dT, + 0.226120282197584301422238662979202901196752320742633143965145 * dT, + 0.0137885887618080880607695837016477814530969417491493385363543 * dT, + -0.06722101339966844497493995074143058569500863415253821828562 * dT, + -0.856238975085428354755349769879501772112121597411563802855067 * dT, + -1.96337522866858908928262850028093813988180440518267404553576 * dT, + -0.232332822724119401237246257308921847250108199230419994978218 * dT, + 4.30660719086453349461668936876562947772432562053478092626764 * dT, + -2.92722963249465482659787911202390446687687394950633612630592 * dT, + -82.3131666397858944454492334105458707735761966428138676971041 * dT, + 82.3131666397858944454492334105458707735761966428138676971041 * dT, + 2.92722963249465482659787911202390446687687394950633612630592 * dT, + -4.30660719086453349461668936876562947772432562053478092626764 * dT, + 0.232332822724119401237246257308921847250108199230419994978218 * dT, + 1.96337522866858908928262850028093813988180440518267404553576 * dT, + 0.856238975085428354755349769879501772112121597411563802855067 * dT + ); + felders[30].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(30,mDT,dT) then + continue; + + felders[31].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[4], + felders[5], + felders[6], + felders[8], + felders[10], + felders[11], + felders[14], + felders[15], + felders[16], + felders[24], + felders[25], + felders[26], + felders[28], + felders[29], + felders[30], + 0.103484561636679776672993546511910344499744798201971316606663 * dT, + 0.122068887306407222589644082868962077139592714834162134741275 * dT, + 0.482574490331246622475134780125688112865919023850168049679402 * dT, + -0.038140960001560699973088624000562020566411307247841147742197 * dT, + -0.550499525310802324138388507020508177411414311000037561712836 * dT, + -0.711915811585189227887648262043794387578291882406745570495765 * dT, + -0.584129605671551340432988730158480872095335329645227595707052 * dT, + 2.11046308125864932128717300046622750300375054278936987850718 * dT, + -0.0837494736739572135525742023001037992695260175335123517729291 * dT, + 5.10021499072320914075295969043344113107545060862804249161191 * dT, + -5.10021499072320914075295969043344113107545060862804249161191 * dT, + 0.0837494736739572135525742023001037992695260175335123517729291 * dT, + -2.11046308125864932128717300046622750300375054278936987850718 * dT, + 0.584129605671551340432988730158480872095335329645227595707052 * dT, + 0.711915811585189227887648262043794387578291882406745570495765 * dT, + 0.550499525310802324138388507020508177411414311000037561712836 * dT + ); + felders[31].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(31,mDT,dT) then + continue; + + felders[32].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[4], + felders[7], + felders[8], + felders[10], + felders[11], + felders[14], + felders[15], + felders[16], + felders[24], + felders[25], + felders[26], + felders[28], + felders[29], + felders[30], + felders[31], + 0.193333333333333333333333333333333333333333333333333333333333 * dT, + 0.22 * dT, + -0.08 * dT, + 0.109993425580724703919462404865068340845119058295846426463652 * dT, + -0.254297048076270161384068506997153122141835626976703920846242 * dT, + 0.865570777116694254343770343821098281832847401233011859346737 * dT, + 3.324164491140930831067995527865720183368600929369864071602 * dT, + -12.0102223315977933882352385148661841260301942633996815127277 * dT, + 0.476601466242493239430442776862061899602963782003580209476163 * dT, + -29.0243011221036390525802623213654099596251221332470910692353 * dT, + 29.0243011221036390525802623213654099596251221332470910692353 * dT, + -0.476601466242493239430442776862061899602963782003580209476163 * dT, + 12.0102223315977933882352385148661841260301942633996815127277 * dT, + -3.324164491140930831067995527865720183368600929369864071602 * dT, + -0.865570777116694254343770343821098281832847401233011859346737 * dT, + 0.254297048076270161384068506997153122141835626976703920846242 * dT, + -0.109993425580724703919462404865068340845119058295846426463652 * dT + ); + felders[32].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(32,mDT,dT) then + continue; + + felders[33].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[5], + felders[7], + felders[31], + felders[32], + -0.833333333333333333333333333333333333333333333333333333333333 * dT, + 1.38888888888888888888888888888888888888888888888888888888889 * dT, + -0.75 * dT, + -0.492529543718026304422682049114021320200214681580657784719074 * dT, + 0.492529543718026304422682049114021320200214681580657784719074 * dT, + 0.75 * dT + ); + felders[33].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(33,mDT,dT) then + continue; + + felders[34].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[3], + felders[33], + 0.111111111111111111111111111111111111111111111111111111111111 * dT, + -0.222222222222222222222222222222222222222222222222222222222222 * dT, + 0.222222222222222222222222222222222222222222222222222222222222 * dT + ); + felders[34].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(34,mDT,dT) then + continue; + + felders[35].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[3], + felders[5], + felders[7], + felders[8], + felders[10], + felders[11], + felders[12], + felders[13], + felders[14], + felders[15], + felders[16], + felders[17], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + felders[24], + felders[25], + felders[26], + felders[27], + felders[28], + felders[29], + felders[30], + felders[31], + felders[32], + felders[33], + felders[34], + 0.285835140388971558796088842163836414852927537894596466840753 * dT, + 0.291666666666666666666666666666666666666666666666666666666667 * dT, + 0.21875 * dT, + 0.1640625 * dT, + 0.218194354945556658327188241581352107093288824322187941141516 * dT, + 0.180392898478697766863635221946775437719620053641849228562435 * dT, + 0.205713839404845018859120755122929542277570094982808905393991 * dT, + 0.242715791581770239970282927959446515762745971386670541948576 * dT, + 0.246465780813629305833609291181891407799228103869305705137021 * dT, + -3.44991940790890824979834154601622662060370460614931644223924 * dT, + 0.228875562160036081760729060738458584294220372552740218459295 * dT, + 0.283290599702151415321527419056733335978436595493855789831434 * dT, + 3.21085125837766640960131490544236787005557320332238705967955 * dT, + -0.223538777364845699920233756214162507964125230083674032084065 * dT, + -0.707121157204419073518727286207487212130091231955206160635271 * dT, + 3.21123345150287080408174729202856500893260034443022374267639 * dT, + 1.40954348309669766030414474301123175769045945573548986335553 * dT, + -0.151362053443742613121602276742518111090963026203676055891793 * dT, + 0.372350574527014276454724080214619984397121028202148298716575 * dT, + 0.252978746406361336722199907762141285915775728129414319261111 * dT, + -3.21085125837766640960131490544236787005557320332238705967955 * dT, + -0.283290599702151415321527419056733335978436595493855789831434 * dT, + -0.228875562160036081760729060738458584294220372552740218459295 * dT, + -0.246465780813629305833609291181891407799228103869305705137021 * dT, + -0.242715791581770239970282927959446515762745971386670541948576 * dT, + -0.205713839404845018859120755122929542277570094982808905393991 * dT, + -0.180392898478697766863635221946775437719620053641849228562435 * dT, + -0.218194354945556658327188241581352107093288824322187941141516 * dT, + -0.1640625 * dT, + -0.21875 * dT, + -0.291666666666666666666666666666666666666666666666666666666667 * dT + ); + felders[35].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax); + + if pruefeMaxDT(35,mDT,dT) then + continue; + + felders[1-aktuelleFelder].liKo( + felders[aktuelleFelder], + felders[aktuelleFelder], + felders[2], + felders[3], + felders[5], + felders[7], + felders[8], + felders[10], + felders[11], + felders[12], + felders[14], + felders[15], + felders[16], + felders[18], + felders[19], + felders[20], + felders[21], + felders[22], + felders[23], + felders[24], + felders[25], + felders[26], + felders[27], + felders[28], + felders[29], + felders[30], + felders[31], + felders[32], + felders[33], + felders[34], + felders[35], + 0.0178571428571428571428571428571428571428571428571428571428571 * dT, + 0.005859375 * dT, + 0.01171875 * dT, + 0.017578125 * dT, + 0.0234375 * dT, + 0.029296875 * dT, + 0.03515625 * dT, + 0.041015625 * dT, + 0.046875 * dT, + 0.052734375 * dT, + 0.05859375 * dT, + 0.064453125 * dT, + 0.10535211357175301969149603288787816222767308308052388404167 * dT, + 0.170561346241752182382120338553874085887555487802790804737501 * dT, + 0.206229397329351940783526485701104894741914286259542454077972 * dT, + 0.206229397329351940783526485701104894741914286259542454077972 * dT, + 0.170561346241752182382120338553874085887555487802790804737501 * dT, + 0.10535211357175301969149603288787816222767308308052388404167 * dT, + -0.064453125 * dT, + -0.05859375 * dT, + -0.052734375 * dT, + -0.046875 * dT, + -0.041015625 * dT, + -0.03515625 * dT, + -0.029296875 * dT, + -0.0234375 * dT, + -0.017578125 * dT, + -0.01171875 * dT, + -0.005859375 * dT, + 0.0178571428571428571428571428571428571428571428571428571428571 * dT + ); @@ -0,0 +1,33 @@ +#!/bin/bash + +for inputFile in $* +do + outputFile="$(echo ${inputFile} | sed "s/\.txt$/.inc/")" + + lines=$(grep --text "^ *[0-9]\+ \+[0-9]\+ \+[-0-9.^*]\+.$" ${inputFile} | sed "s/^ \+//" | sed "s/[^0-9]$//" | grep -v "^[0-9]\+ \+[0-9]\+ \+0\(\.0\+\)\?$" | sed "s/0\+$//" | awk '{print $1+1" "$2+1" "$3}') + + [ -e ${outputFile} ] && rm ${outputFile} + touch ${outputFile} + + for i in $(echo -ne "${lines}" | cut -d " " -f 1 | sort -nu) + do + echo -ne " felders[${i}].liKo(\n felders[aktuelleFelder]" >> ${outputFile} + echo -ne "${lines}" | grep "^${i} \+" | sed "s/^[0-9]\+ \+\([0-9]\+\) \+[-0-9.*^]\+$/,\n felders[\1]/" | tr "\n" ";" | sed "s/;,/,/g" | sed "s/;$//" | tr ";" "\n" | sed "s/felders\[1\]/felders[aktuelleFelder]/g" >> ${outputFile} + echo -ne "${lines}" | grep "^${i} \+" | sed "s/^[0-9]\+ \+[0-9]\+ \+\([-0-9.*^]\+\)$/,\n \1 * dT/" | tr "\n" ";" | sed "s/,;/,\n/g" | sed "s/;,/,/g" | tr ";" " " | sed "s/ $//" | sed "s/\*10^/E/g" >> ${outputFile} + echo -e "\n );" >> ${outputFile} + echo " felders[${i}].berechneAbleitungen(dT/2,dX,iDT,iDX,pDNMax);" >> ${outputFile} + echo >> ${outputFile} + echo " if pruefeMaxDT(${i},mDT,dT) then" >> ${outputFile} + echo " continue;" >> ${outputFile} + echo >> ${outputFile} + done + + rest=$[$(wc -l ${inputFile} | cut -d " " -f 1)-$(grep --text -n "^ *k\s\+c\[k].$" ${inputFile} | cut -d ":" -f 1)] + lines=$(tail -n${rest} ${inputFile} | sed "s/^ \+//" | sed "s/[^0-9]$//" | grep --text "^[0-9]\+ \+[-0-9.^*]\+$" | grep -v "^[0-9]\+ \+0\(\.0\+\)\?$" | sed "s/0\+$//" | awk '{print $1+1" "$2}') + + echo -e " felders[1-aktuelleFelder].liKo(\n felders[aktuelleFelder]," >> ${outputFile} + echo -e "${lines}" | sed "s/^\([0-9]\+\) \+[-0-9.*^]\+$/ felders[\1],/" | sed "s/felders\[1\]/felders[aktuelleFelder]/g" >> ${outputFile} + echo -ne "${lines}" | sed "s/^[0-9]\+ \+\([-0-9.*^]\+\)$/ \1 * dT,/" | tr "\n" "," | sed "s/\([^,]\),$/\1\n );\n/" | sed "s/,,/,\n/g" | sed "s/\*10^/E/g" >> ${outputFile} + +done + |