#polynomials.maple interface(quiet=true): # # This file contains eliminant polynomials for each of the 8 variables # in the 4 X 2 matrix K representing the feedback law in the example # of a real linear system that cannot be controlled by real output # feedback in "Some remarks on real and complex output feedback". # # It also contains the eliminant of various sums of the coordinates. # It solves these, and uses them to match up the values of the # different coordinates to assemble the solutions. # # It first checks that these eliminants have no real roots, thus # none of the entries in any of the matrices Ki are real. # # Next, it solves each of the eliminants, and then matches up the # solutions to assemble the matrices Ki, only half of which are # output, the others are their complex conjugates. # Digits:= 20: Poly11:=371620181654204560521275971849420243774149894588726336736578305621579445723745129144613359751169912163917499745958842036719288702230689803523657660147971531540167731283457807751203903152436389631857433379185920*X11^14-835669327812986155811399651078284646254382379774596970028329067394168278964152580597730101791385474189211441556673564203095712869977110640979061981940816049039969148671665003268313121797877691057872285919382814016*X11^13+1621402472485906772639804814367742248235851171051322215373612499221127411551840643933185354876639005000482826517792638526611555167346260604838061268619713796791703743767503734967877368786787518139111063642462486136656*X11^12-7776104598136691762613375710428980161472594500156499113669366169406052328457899619112437281317840693235618978416211955913451551380629317360070978893763025482378215062996433642784488510930915273082115434779180570010815324*X11^11+4861675254743905151620587073915508634260822428163863201715360440966102547087694982593677199253374258334202667839235249070770275391947804772717522922164662343792163882503604411576478818596392806681828031939548421132440450566*X11^10+11683431781649064168733182206010695434799143888318589530664586558075874589391768719848657472567688725589221721241102802460319841781311349074759858136188058146051034844547522822771674955172142951613642666173263478912849018921163*X11^9+53480774475234315028537707008438319295106743642264408656865894669263700199506387442347721430141195624304315830177525195989673990077510733638989590183207241660874443223040599317557489731906940393120458306863564605253487503893888518*X11^8+84994339745668535044466718472264721978517018553861943503166546469527783400423157628030616073270032752603202099585903716481529388277852645755434812186903029858953818294758923634635657349359576075684026515586568162785616578006614304091*X11^7+181468418641108934857392298552334506137263132634457498698326839856002277593755113842724731279353748731740593205735093889986828729189316061406896382732083791996793802600520023601112113973827126822622006770017078376960258100230883433506143*X11^6+70484304162143408959854037631030157425168718746440427111543081321132502329500829438378823192177363755248990704723583078353838025154665766144055474708248061190924576130475868797558895543188225553420192887151726737557229339477226502959533382*X11^5+136652190598663833006149977636575217674659186174681048574559779859347704992142351156235848334782701289226648213113042542895746654747699630561315986295374172826619934117416782171929532674461549363994695599428576031907510675453564283069465224301*X11^4-27511444135126642954601192805469248242970437815047519893977695017044351621689663402273202563769961428536751688852715464426771784210523456438835404463554171677408460175449132822739567794258970205189911391638203671552390147359383472241270543180492*X11^3+8601761082873237418326934913830925390987933961457742280801485156192592958853695730306266414191280358888473640725309523458448426741469974271282744084345654582382010582413277153312677177432856635135886318120456908193091921967774037454725489162384484*X11^2+1357446086401725343182051207368460161735965006282574245770507428775681922569061145027287774152930836067252657092496975129967639059318553317367843462906736666265704853642177807580940191744750127392540475480972633153563840496737731491907677816059108608*X11+48550689289372764989455064011294994963575206837869211644859672292904958051207787054705352810325682817791983844637606196009869685559069336422029473629004907269140161462624256049691075202911947095092842103415285329304778174788941319486530083006637981126592: ############################################################### Poly12:=58065653383469462581449370601471913089710921029488490115090360253371788394335176428845837461120298775612109335306069068237388859723545281800571509398120551803151208013040282461125609867568185879977723965497800*X12^14+59322279974063387969702265447308158034464379404939090054847464378533416979839689749486677989193723561561145544391430897784591962576611683155728389232790112631219718561489775477357392290629496837925266359342760580*X12^13-1408034663755441026449709567295680051946738365497385037247337243997261501210400656529995843851645254011511643467079707164613509694883207225326724981284662820748117892432506350685290988537013867260368814653746335481038*X12^12-5368296300287791684592578773021012187731880685532657852495535800213651377566292341424015949456590467883024974914311334903821460244680679082686785983649043787336368529385632754411264248867513731762793312213495621784497033*X12^11+29527144794933605841678596640073094948616674953484187999774157488330892135032914785070760021139240209894490164323827919363523315940629782494100742428854020097336364831964256100454955831226906296050390305446684936529441077001*X12^10+44470565793125536713910528059744394106227790575287999197845862769744493891841986988761246769446307454022132148689405526511725147853644788856394028975359400089237974374948648904998776327595548168305578773326934435686028494615985*X12^9-215130536997581940535372020648604037820966626082856430406129977565788141736825393255268228655573241080665216202872432263817103682012335929656635205398730501117462305327397014354397328073872390844193888808303990818157145776836236835*X12^8-155731313678627347295713914564183136218457704215150902886574387206771549754076332317340144154937783540763460272032670284155306939546010517041070680199597791883645404522211361983243880280475428881107635929026330970463002779963059911859*X12^7+1035828934498592908835923965441508740556907061268039228576947227499575464736079622970572826568869796150983695525092134046341188603495862903624044435789460031189388148358497897398679620686748314660854735008979153849724421931083866056310919*X12^6-1918569550859523280216782456095660040331429329977534421060801838383557844655996043232308202595236225105743856142351394096767624480933272253184704921719259394745890939354129680026751173477702324487761788951761799928188512877380582230237891153*X12^5+4542684927198435713536878320284484976030167478512071623818303700053124947639375893419952176667454031633544690404729102974254897435363349333672803334382434816382442628613430765720545854979327229657857730763172138757838648230809653640647325611921*X12^4-7996757479571801867928682713053847476235039264713706775120184637274318305622852304739078388449462802531632216887049515175976776472647524780383929922755512970354394655739291253617581617349507970068611674300952059543475558031607011627214195028356328*X12^3+8160917265915953372869176432491306173992984298038900580555604402251026916223565595965618982545869263581969524586006430694924555992737406671178975191020476998368359707066203869651826189968132447803391286552568974506630371581321233158871988553889293352*X12^2-4707006381635795934256113467218427469920505769366895061196680416430589512030517880482222528464997943903007929080909131615058776299123370322513368989044551960463608469405075382955332632933783861966763843030074035788548993879943103088726573199110839083664*X12+1240526051289702614189864911330753720168521496989670853503761489308466509942521872997045957082380320186139403394667815570552998936079019086276000804236996364331506394832273707852790100207409399227552692069415680821720593974531674426168881191985806179854736: ############################################################### Poly21:=3344581634887841044691483746644782193967349051298537030629204750594215011513706162301520237760529209475257497713629578330473598320076208231712918941331743783861509581551120269760835128371927506686716900412673280*X21^14-2560449628871048673606256678034776883370393435323840152812413687124448708077390362251770593327927041596948421139767014659602388419986916289646215864519831764650605982401172960418825711638778796520847153020914340288*X21^13+1960206021254275775905977575090007852485807474089045199518354029849177252166272406550909563737882921005636925821025213137849934748463706339913695527954036639841372735940843014518128006798885809261094140577693441198560*X21^12-678314828549852450766458815607282805054752338054557232853606693863074027315422428732854154933331041651519925100528677633430621785104468328686806262642603537518801566310583593357947149904853041679237546625231995647089632*X21^11+150042629268301301548347639377105802866892059735527378603825818379790898964869161089167584232755720880790082355439468325868700955455232826627941069613830180584948031922065379865521920302024653904377437761275945194907340248*X21^10+45767746808149843011215624771607672525605799142688551755499457975488447792398447480578931757750928131932346978768438345154961214197958825759942167186347864040065147717288531442041027453720107693648692366678038113283366231644*X21^9-58787541910478119170276338326218004922087236756470941252027095226470978530972132381987249348623010945290702824652216374734503944033577868904554503108703025218648050553387634990222649773534892762456548953125466657893266485395192*X21^8+26510881684710593163083933040213918893560513557263026066086649469535312948644075040310537812524526686048381914890574089491669573855595974680819520542198480942781357973641506597426743494104843893952861618435745689794481320061776032*X21^7-3656525615782331336146264369955064568610514467584546763201398352833813464654877984595317953127323044239467023062370452812872891184641025961738477590856552193497007101187494162458114649563626115514491469515077633136350664581915252267*X21^6-1140493468792249828611306809631998733860051556264138254682709258002341707823187080558108526571004350196746819276107785395277377937990124950616809706008504853980503927217468342839399070420669189195572674911859074827027890169049115787943*X21^5+660299002914598958042358589477886685202752322737424882295704188859373587585560234375985719093580090405912374199120790046657210126578148355857207236542881538678455562187892786275656815394373173490105396373134994410769411527911609529370949*X21^4-120168467749291415236903067624306160929340777116212600332897361946944212870406808058771964308574045190500391486704900950919847926021352222644824255811255204903856456614856278478387709617642255398844838596347572520318250110653816573131148782*X21^3+8063879546423809341978140593403340037859563182680032587119172203845534291702667723447866314918413527108305912738397493581395787206933342207344320227133407388309465927285389318239953950389451736639925689822716645150094254364637152603246515775*X21^2+819147968966384810354869505267846759820609474144078367661093036491133441467105050107134947848179411774375375700230552175374633198554498306142962606662187099493840645528699605991285230690637092402995892604318367400255106565914001440762230418257*X21+15266621121975524042656070592859072915083988824421084298794256692969674919976000562089666154446714252463200486713944464032548457283245599933995640767658610346771269724832437693560323692706467738922873097479390828722768333980158809832303742608639: ############################################################################# Poly22:=2090363521804900652932177341652988871229593157061585644143252969121384382196066351438450148600330755922035936071018486456545998950047630144820574338332339864913443488469450168600521955232454691679198062757920800*X22^14-7026423347446501938064401372355586805896829260550100159333576671300085300403146783621968755016627645062676725294052482612802411220409785203982976070396844006410137872671071019798453540865427505543551070187578516480*X22^13+8557868070017458002560103150099697411565854361486496910309703907364392663629467655017359082493566219061599970541747197950887772658901906428027536350464911310758568528782925316402624728887636078370943430688406342717024*X22^12+966912297032218651317926300263415769335899452774099578126585528165225252978539907710881243642277995445037021922736003683240076916240794357661417140339960353460324778724834774536412883546128905419649972336894278109544876*X22^11-7917646683880391969132222652745692223013720802732904608474896603542164153848645564062920659946685545355544259086825582190824812721917378274819618078850243688821750930162729000388161009802360639382258545725049634467709917440*X22^10+735477642043926373552644124106912185581402146178905328860855320956704628971344759941435610434061800385183921357603018868263196614667633981069076105627289155906202849546249144997162862229729088884717847322091384546915707411010*X22^9+3414186341726801007020107445348610183935339815251650777968200923644568335154291314496277994752555451205609445064799500457359810375425646111618352371385037866296340196949162298253913221842876779387692193997606645686799220002904355*X22^8+1022252594189032334240351684165726841996663626169798764075606416205808616958328570618620563290706725787050552911929192810526167369231430857083142393491115557740977803849719404198192904490105869515635068994518598055806612642508384303*X22^7+262841396288746674917026409152987431833224621223564287299251192017694724244755219622229264236457185027722458961343470328755273826140707951249246509411473365464029192560220117584705884730630378220408003544011234305056841011185884755076*X22^6+163512267744236948973308065160241679222666035411627217967979582122722989310395282471874026065414503457806749768219613686020728706106351745530372362673998516783584549332047259605748345946573435334959721184746227662028940702386886562955653*X22^5+49025675530252393605600309657469886731383139638813636896135332714143455090670737474838967204727930387916021089168655425284194614451584155698690149377796678530076351173951571789873927455503119245670806489685911363364674767099180990533857828*X22^4+8518630642069260896707019491882278366102227314185207179101623380157642112714285135708004323491457424966837488543985718319788740031045827955001437377646924494052482349303885197713056616579922760617960848978265755409123677727294296406278213493*X22^3+1383706902901530413071680200895755739633439816510500784680004254738989877852515350987188871702183718912928547396747212798128966666883508379532559088350434035258041218845863833320216628407039428342046958693311355552560679809669163980063775592592*X22^2+70318778763960836322996037677323272754269154255230680365345704581245886128499052882061214440890583404758819637715391609584083164303323006357914180345435721557516477118445310754461227160620050825074390300260318719392660645811434707281205952658409*X22+1571912669753803207890900389398439891345203008799955238627382376414570600926634466000085576579611187020803761011011623010570634496893930187992067985145523978199516803575793706029529841163189128255540249887309518199169650186759945973619010921323813: ############################################################################# Poly31:=3344581634887841044691483746644782193967349051298537030629204750594215011513706162301520237760529209475257497713629578330473598320076208231712918941331743783861509581551120269760835128371927506686716900412673280*X31^14+81009594427888204232405265834825278689706674320464465904117142445569959850774684756281077578312793842785218201246943460508414879050825989814061425101277720053475663777507216097320482903344722996966072837658903424*X31^13+357024240416353758618926265784644450855605561298121543627361381157310778748907061949183276157516564415121926186772554876288949638601587816649331515894696040257313494253847649204921671248625249104672649506473137650512*X31^12+1069839286317027404759782099907890565650062908727019179066399846889534185825807602172769553151243447641703289374705680383903206971198865317652405182712654686665922177689861826852380315553756486881285621815831873561440160*X31^11+458092809310151006139189368714655711540741266029711589705695164117017631796827206766372844876048583809710825744926134730835179705775538060518098698509711638896904703915965422611115011700320954860377375328382843332544369648*X31^10-1092681452777960051074506324993584183943505465452824973242524337478487124363250013906567818662051007215389488713045278800584172779191918364658791837612908801835114831430883477584868280483986141780362469681606481217534809436110*X31^9+1102200752058282133751637492495911416486476397730877557807420059512816440957336651696087520490976014929923903820190214429833074085452473529935524927535244397023082298567621789309155224189330554818863217220424718443092246488227230*X31^8-635705666172506990890540435752526836955729508627994252704505906075533822090137294725491576003211971167967474113243266465823837471589917853849529536702612784829468706073299652286975997578884192540549655344211327913436020800296704718*X31^7+248710834083584747979791872084280795002604174678805380035983451222399854125111987471409074624972957387972251276110712727892064536562111832195223233373474466594769432283730924285565192717996967312677410250069757027173597076216289415327*X31^6-71424429036859417272423613234674604294350322286590456634642352304308881867378925649373425198338188700336465505017453630684908033769540867869237541727300485594264851418626472004825608523177708763522548106685210111411491718711758581284482*X31^5+16084656855292222585623438283007024123305386721273322160464854700908570134761689027963673068809449569119648259441809187018132105932892379664942140823106431015414185520618829929762338161844850761786096941685415334179525438788486203948207240*X31^4-2805159502338610991919480457991409770334679377080287897518869313367039209279265330783531097951907415058672284451172865997126811233514930579904808966716466583855131182832671645598212002648710468690435672035960303901863534527465885228014166406*X31^3+479692664993982233433471695264131868451088288305045243011604198071055194741020786473283889597209076179967350979466634208351432741408758894038136553670172861063284725515649572179698657462716501043506708039198931772545752003667796770346691055955*X31^2-84887112336548803827806851791067773383142751845425603287594026748066310765598982900011785263128310962992098228089028022909363805421503402219899134298286612103932114704678986802358715886799682655827218225824390296743401191850850076575430844357332*X31+9512575189302352898652562276101607232681904023885113218563151549806506957514834239197189470219284348746756178897044138201748051842987415406658001262169266039257327803311153171857864950307911206540226732905832305656957544592328821970533054019265872: ############################################################################# Poly32:=522590880451225163233044335413247217807398289265396411035813242280346095549016587859612537150082688980508984017754621614136499737511907536205143584583084966228360872117362542150130488808113672919799515689480200*X32^14-1329455244529399018841317085582283523919457112484430423671206564783513235854260114448112649921828580119054199949419896131188034983270325608192761842029432311900896637432132137198152393093663565222921529732364262100*X32^13-900831239439408927252783964342286997652468805739233726425394649128641264026482074507674814735691736749960276332665125572122977443761580463728137845492124823808035113167068151036870685499377698961691810445235184154778*X32^12+3372799756089362540404482540996274742657909000964776813983424215566125578331730267377048810196273013237158539657516084434950796104330939528555775134678400063192821229045039859055762137024122123984952592579100855751940996*X32^11+1412168143025388376508021945622185377224092384068251146438718577859370961961507037306843730817769901067191099307335697566910838217621279422448702249396450086745790269939480240807635347408022536244619181841495488463628433082*X32^10-3781758277314052804174699037911369484155143724684839266491522042417648762867714302224544795645787736357111356915556177534422458616518665257011195486234445723620730266138787188685854432432996620178198741194431034025494377890161*X32^9-1056148569633411972709980247026590145462878544719430390021183664034055906126099122435863356434074549275843811977180267467268606797697714294243745400764077922743008396208288683874896201078080848868457778991993049007797501942480749*X32^8+909584113647756954032613712148869246685357541757152254263772257799424422549161265317355525990006247522224249996857287202604153159770771656974144316197495875538455801184998889383587209636385661664903689047877510499235234938427218488*X32^7+501041530980124256193145267914353414014514191433432294508400810685072899162072790772827695697566043847417858672640219286062365093536309929081030714251169053811443250790076840038425537936069315772144034665018613450781601408072788074251*X32^6+256956205179818628016391360995319917563665835559732458886933706988750474373043916005986405317303948471740143644870125142180427681583808042184727922548917251231834494974798251328750193731163936663732969159908675410314631835127694654009900*X32^5+177252803099410731474948082864662059516167249904761784544627223085404619760245565776993169473221200876309367672016250230218419375511222081352405757827968551003647955617749877942952474203073171324837382794617031083098696187033319411229096309*X32^4+79062290140730793943992830748402117603274220562778233868690262399629416585528815718809017375843442775652215155396125456052180380008110123414376198516587542974172683727105020800516183862196986380282738002717431866408741027751425524910597405064*X32^3+20910823485093143640790809700258560764293167490992523631813636018426656084056115139943127495975331974925850699553663123438732498166977796233142100361209970697370321107383627238265168153758500602786302572671790071144433924326246370316150089615995*X32^2+3163123021762920395152017640311665085057709329309347713156078133422688640193351941978880073146304188112584146020291348851286250883675833087956471901492621096829321312572754679682340124321419745043652214805250418737215440815740187317502779163959861*X32+244588588530390409846734292148163222869124089881437614895478440426608103848547443352184158493633245017703781169635390092547916326797459470530251354570777765297634015477270782907559964317223373197974320159598875225556947119225745302810105130068845994: ############################################################################# Poly41:=836145408721960261172870936661195548491837262824634257657301187648553752878426540575380059440132302368814374428407394582618399580019052057928229735332935945965377395387780067440208782092981876671679225103168320*X41^14+3781178006591860094832137861040362386270078781683620761418742793984023583436462888232322041703098228144025706398222778623848486722440658627232298083260977397504754983420158971140377648255782535837902039777102734144*X41^13+8171397946819686274024238108823394287376540876332805434044761796132974859492720851488883724435480932805770699608564716383372079018092820316245998018751960166611584583618887560435458056526985411114987121326515862689856*X41^12+14119786744937619626564418097932884214867935318927477764515185995030361678784381641387691629630021596207040548876040393669557403203425662469603266299589299187715278100414788386589638149941082916385627951729653780954792848*X41^11+21165678879400731725566326699169076994265993769105062010543711580748269467661019415962204983264909918300968238059242277141520259784259012567046920173765146549007617313190325325787500188775701600728743829697838618210364924468*X41^10+19493563189013925299330736636063825919662042914517454169195529990787803247935448968663566303133307489458398369937258194615805361079938787132734932774885897404274619848011681960689835361881650099368360806121680382218347349827916*X41^9+5988614355292132880267278484654644490979221834422388151987515624703480552309163046091313407095402867977544717298286559670461261858100774149523680847684440444962999783201718463225016139299892508841815904350526632822668360568363254*X41^8-3444889641983424619929907372066041315078702084119848155493101998050526925887348588680146785996179436157800122563079065330194790382154673918122261756430664496515171805102246493652028077129792572460795695095379234155129280036932061765*X41^7+4754603835776305127520080300985597800384687238418249003680371680980175472540128566218704572645821297658083434960718978248675578271778586100801888849831459769230353267259974268126875387537888858509357444250648148815270657412306791918075*X41^6+9173683723395026908470674361220147772449263534143709730770107840971777620996492573826163513091016933499126430076516061301505970814669787106344258184040680912134210296681906211845131984545733073233093711552986367544096991548106200623178487*X41^5+10270966161083139004867681601766471212629683117798827777207115407077354247087813431175937116363213465487320818856097239943813077713113857182386012793576350034973498991726568530271068751211933417627956054226754265240185516359437411231337845254*X41^4+9385615576853458231179946055256277796079581491987835780015012504308280351367799977054107269832568425491340422837528455601134634053088372424358608601776024129004994832222307433561105428431698237998356786501909356404795600399349704009814287770300*X41^3+6063849785255590436335574771067712393381577316966366092006486098563393200748597527994296388658463277801880390705959026323776477847804161457694324381750385953642320150995002089112381455137676926284651842576104239184145135814043675243532466745268888*X41^2+2233061091054271247886614164063702576217522853144448711175159627282141993304069279801434014330463724374362172748368631442971687396787302061143761021013680305455343796683050182586001869439687472630245009287876368435200544176377986022760361823851969056*X41+361575789658201589139429003038086222868990684829213611138671379711884874793905986133957529718800391368130925893089546875076041546190729041026852392808985064941707520372402966971105966516602553964029890114620371436078954468635969245116082662090939207424: ############################################################################# Poly42:=209036352180490065293217734165298887122959315706158564414325296912138438219606635143845014860033075592203593607101848645654599895004763014482057433833233986491344348846945016860052195523245469167919806275792080*X42^14+1668510246034704361191260505732168627252401597860532686682096507833716266936057706169627746835123391758708166299793263549150265199515887650491436133994049282833728727242295767005630949427409112457878121093761660584*X42^13+10419492150773719927897075312223981264087779775180648559406285891980252577466476348839370067389016014970827831050129136294241165941212778453161394414341721195554324448730853455859817019200042113701115487538495201974692*X42^12+54311744347110087704811052968911693488796168374870969581226094981211889883320153434667364019092342083349680281821328968175320933178182601790658493220424352665516553636323622108343448750722424460611961174080839905130057056*X42^11+203986940822378888168342771720714743727147857673641410615249825836587655478410376855701636830186870981079161814835319956511837319202304119187764525733171839467118124117836993058330243525539203646520624206051503096177682246384*X42^10+553619860650592832976501829888905594284225319399944011184339718734567591989334658992922943721045054363956303629810783378433102874298831152447369193363654822325461070343777306778558859331330173669859149662989153360334985521579005*X42^9+1098947093393135874616612676945212439425829874709983460518051845218258860353955742612166510875809983485418057107347909813546058054788390678876535274155251452651224809393154950514119097142657073469961251088045388815045485137365627665*X42^8+1609855693214382728830359738702933582047695122680801379947390289279152116701388482903552112870814964900774420038959463761664813630033178937227239338013225199315022593478918882820671202287229015198660657099184389208686368668144413841384*X42^7+1850855121962115455005229164339235018227157498172777267404189403519697665952015876399052688365256933163891127060023971635725634030067943067368482944317057458164584899299948386437978848467600767478118128439237187388399714859746666493814469*X42^6+1823923764017257161236880990325964475084805205507141654434339585173040808382978097896770487686841572010829936412849646847215533667576882579987126300824531612836810811358884191304371262635873920639760636898457960823933668587608940370435468497*X42^5+1516863028520555905611987318265399033648169167360749863293161261458065361197933588214259310655642195492077869109193195092822431819008325392709000211075026731650774564620606695177758354356528683660122506554262776996948925032071024302917092634336*X42^4+947395223909173645087840293811443930769935311084200564931223231908940171807507172550500950107180057434254824600241142541723846496734830978520209388719132536272882560241692865318848501561795727175527451010113607839201363040962861108971712691832060*X42^3+405870427629361268692091420094206673100413805204615397120544509274615878013583085428707470939391947199980658422928209414422969904982047853354026056700198286912643278263657575069223665211717811030972412897939321736227599428544766821742708628268880112*X42^2+110502408680475899988931591939570169931208864751465059856172814189560529460456538040889937634765599395752866316764448545577376300563330531116886922640485801560742126262765872360370870631266623869564227541823939014726043184726115084846385561376688263168*X42+14849330812354801591044155807197541403871712878868377378231058041436491363265948443249039450296332710928520940065532749337135554101115005091037430760979264348900679372524562717633502275342107827040177549730194800101261687513981439950710256825099885853696: ############################################################################# # # X11 + X12 # PolyA:=123873393884734853507091990616473414591383298196242112245526101873859815241248376381537786583723304054639166581986280678906429567410229934507885886715990510513389243761152602583734634384145463210619144459728640*A^14-152002245326326824268435050738504144944603450527662264559101765123851470097726188733671787553515214465073370024189468819091441436495598622927466763616986442733387649959376813404408603712616303765050193739863048768*A^13-2041643901091184903958139601943496426557203749407566671306353685599346430865888261753924130216979566795034134618579936888287982744974469634241886269170952078232291005220277705467042017189501486500909988310457508104208*A^12-31979238464216545524206886439255985831802211893811596886921692200621000888540944616576667770683610573402272296193650080338249707556265783611680068066887152337910479478082151687789757402566473557237123983878216288976832668*A^11+84405511602632884219465999344283268386456259913276991559148566342260518715946463607722710731790938734251761314110778585921467928044249799652490900096030564513254160050099475729228278716570381552456167254021138647023969254990*A^10+401513176429856561238807431248459261800275482662231905127511081094557637856594996612850851546868464754452600091754432842630556713962681770753790607333580238710344858754066706307213142986545561429115083644402204169752015868859345*A^9+2825486191675735666787934786930625564599726030100746001327395669108000990862480661424296600458301774050491311374120522782058213738096374883419007659364202404544625540341653860012430290423889894276347571699355348768237738854102257069*A^8+2383968694543157400894825636147212591602887985111825661849121156231774815475679418175649833630706462249597075227265518757142957980262241783100090537987614006049101635934204913830412467456541150709880508577377608728256309988688194074231*A^7-5267617932295789099619091507341298457914908350316991361271300799614012932434603945716190414712572567698353669340693297846811750108267441218706495551648274080543568497937929281317065845405661043486221588336936324789829150528769074224606614*A^6-59194850933148722992093095444025230660107031037050956042018211007342136074493216323226659828200306329673880691318773105058311245558648885679201472650103138748595119148691878036579158304170013843695571074286337476271933455508458464701483214772*A^5+289731056801123973758017065225697924819698560641241750826451064227637968001236236892520452403461032741339451577327348704883440283872282439314543107944736502074323608087768968923655725043857057867295503693276697260713514038064818499276741621245598*A^4-789732462344730411281163475355324280674748019389481893136393929897290608160415295842314855116362594137694925283212630197159413971981696162768187201173440795274610680207994676053846882238722672964627125018875962204113390757880690807551835628314521080*A^3+1308039248151002022561776846641426240971106705741494555595598245552405713455962615960788786851921343769961975914771869763448786661035613419260250816315809935016130056927299912882525152144474219929119627541220914423749646673414493527640540085016633125432*A^2-1079511683324890765273387837225599692130387564183464753578753850833809610441182061581545736107193722183310178830094071172257945206463619841125279509752090966786903918060197293600009791683153357199550971212023047080497707946603132091689940153132752992520384*A+368956087611388770004056910309466311374569463207268448004040608753781472337277789357335029803401477530336313448343676207546114660835037720391431629188440270482178033010301377400707645814073965208097802596277176269180277917057078430087675820252464825939959808: ############################################################################# # # X21 + X22 # PolyB:=69678784060163355097739244721766295707653105235386188138108432304046146073202211714615004953344358530734531202367282881884866631668254338160685811277744662163781449615648338953350731841081823055973268758597360*B^14-287556812183030245302277059870910745266777505254250008494711174191762191431717192000977512528219401535358982950213562225835130132763386929500395366190724628643892220389060137335340653687988808445636018027521665972*B^13+594651058545410061821173969201671245398536162544753717292220447018526836190545465018118976040979391510544670048199304102868362296220259939165029112326491075685113574977398322279555801028988281989536752707791898689200*B^12-461477632658984150099134073836744381068073031469857884378821245678688041646768822252627458836867562848026008369517728142730379609776717203398431432011726533051210962316835829904702226106132134567659269225780118642232213*B^11-223993913877866184587137033725169519725984381822274862092427434480675113749978432593459080971219089399632133971747335803751421169821592940818780175457511945512976996710609973718046293661579928935205082547419560874455726170*B^10+611816883319186731794170059059388698209942960228719118234653314543367434529565210355824994736027027635831107482628169884160125243112596112740837183944582842315704111838559324738812367048900765084731842432061454988634232950543*B^9-216094370625178583955000412593378658515183715203968339822671444947556989292579512456761898324062110606814702635970097097940182080223223032910369261418810341370205207302948283731214729933560220273446294546156797271759024297362369*B^8+21734739565382900021795058301050257909195777350763632624800992385773264797005529002455543696067018653896162957689995767489017217259385657983578223722187090517481985341804085307197265414121864926392372635939497150573527377111187554*B^7+269403589481568441396700054856676580700246876825291477583253889187531824114031286656702471653008466514104175308896698814165752429328091521741028406715598236348023870498935217401121709566650866755865895088091085035969158309663176230119*B^6+47446909897729455292261154547514084046854812990962999679511044002231108450935409977171022683556283406336064341243348824960407969126277297410608576900369387423473198060220612455989271440315433510343070374236065983943523528988055868471134*B^5+40179246187438254003038053628586868490611636194177685126592864482998066186282064644162094811757147339334610676521226683191813591649320175199591884920314152871151303366029164515116972450969105799223545026267193155559990792734247830252683640*B^4+16907791393359042740615486965940562761056701247785682469473413548860528884013722920807827887215248465319028767153336391220237872869934897221292940168098144458860835808332151313161730637860912842304318069234584140278689886400246013856777355960*B^3+3170925972022218582513660142403782147002137472326286827163733321605787211816694450534602490088307145016976430056038703519943958078942069470512022576952022119860868264229009180220965407154847679563317401324427330560441740264316717244829821072592*B^2+241399012931929866713418466373164074687370673841834907574468419645534146593561942760473642186363909616247809719187265064967711512020067405001455317502255607477552514324206663336523899290651166309919718739554101541160748926199443033377884107452544*B+7244655060589296993786378839453856235645933094880536326030440428245090661630692631370438475575104597714286940754424383986777080051958007880648093972921824633052913669833749144150332670076396029845855136515326786044239432774557138077525568835173120: ############################################################################# # # X31 + X32 # PolyC:=139357568120326710195478489443532591415306210470772376276216864608092292146404423429230009906688717061469062404734565763769733263336508676321371622555489324327562899231296677906701463682163646111946537517194720*C^14-351145998773344396514667670078824553099784118565828760232983536340371447900687085321318328413391254954965069228126682990795625375578302412609150598495295378171344283991172435915452284704004253934572154893728015584*C^13-452897129872874689549145651203391796510765226164894675783585298275155369084581764982127267478543971953094834054187752440684825676491281345496038086836545305910982923763365726765429797725012395584080791462098685625194*C^12+1628230035683557012432591275906785869751262343371346304033322240789426182710298069188981478165204475423263789340352599202004602147840516341670777131688789857532432764492174848266144181366499564680256564787385440611178521*C^11+869604512843820969396807254100418858216041566179494936403272322674815971713885949153589935433470649756451970533463791476710837512543955846207328651121911985430181801629572773607462219815260396465683344484092656809587520762*C^10-5262687883287253865623699586426779610520093450719506180708662946294091029770504887154301686682733072934021099592646186619130933813209727486260116803536693940786775490008226503998860872625688493727335384092207599250547483016893*C^9+6852226485828420651612526094734255963952672474281636111449757036293174946945387896798796995346011040150622907603673349753519583231444234321019492683877287930198341890984185787229078343857475233142558905697363726497443710464648722*C^8-5564779686057873127873559786476985704200596453960386109401420152450750284650235224891436516834881897572049334974127617154593449180892267165335042085383848870179399312548703703893726657664083143522434544979148886946259831511826860175*C^7+450249853058300422407745708467799225666346088263567596996525225134483491243802652978521156606461997746128871403641727605890876635326945838748410792578538793871309727164398714443355659469703410534767450728428229051691015729509864476922*C^6+303158117943071850919764602966980354799592251005409343994125947332617678697163552665108001011126730611823073772412928713846735440619735448739497046720659922000373989455887461799618670040542374222717435381013791065927749613376099229310728*C^5+1881814091001384478110963367272751516058883829168055527456823589366186475695237458618530268001298140673358874979348165559114796126490402669204623436924905850254140790298002531689226621806921045027247261571341552505013722058143254206802391252*C^4+1075938509638588053668499853319963519271342863350876471145560964389831810452457503446921042172743196987730817090820032121714930600276779083284668954662854213971297757388825183076381119301577169009004322179409814331143238125252341968276268040624*C^3+561538879804401966059848958300599737987430839691255402634125613657659279348887688245009322791461660033528485218143747960127018756569852070419839088999599727152185724946221041653341422323366634895139764174925125689860062467613672045906340658827456*C^2+113071124822754579143124134961272155364506187052173030441905900211979206342276625413948094284034625153057070629694138861704558434810958881120226878688565963702107548606660574433481610418739658129248895686660317596319800499175791259973621140263098752*C+9118443943785943927715189529917098860645655143069394747957151048284224043443308650432931524577223623874246783511464309468664618253317691349048732304306343813309651080477804474858755196014136687706034499138595656932509573474898996637951568745961385792: ############################################################################# # # X41 + X42 # PolyD:=13935756812032671019547848944353259141530621047077237627621686460809229214640442342923000990668871706146906240473456576376973326333650867632137162255548932432756289923129667790670146368216364611194653751719472*D^14+174253649845511292326619664732817281587994752885429191802452147088648144186344895215180550484059863252980972859956597213674159125341736820486634043653952908813994498206489033986048357432756983094490242069202489608*D^13+1460013395642643705650940116933332093329961825021423457414457022819585721761632982905427152568185916519005162063371151045136721468039697537539918211847056408403115659094191500096043301904047520504037593955737559546776*D^12+9369343486694616489032437170351428863693474644904636368897452085232970209385560884614149456358673999660430519778624135050804762994060564329558476993785043712190316054214657184322467335873283540992587015497034792440088558*D^11+43310992155759537901077665825204025160269924995073783769914504098352465306598866931860144639906900897391778702988708194488975341536701091368528742715089887060623919317015321924395654570469991067011031533179035469514336890081*D^10+155723623183316956264224298365127261959610991551698919772056551258066043488391640047400916498748522997055491430269672818062821361455343970253679506270047510638743518592179397348259073335051461917964835614694258100047650358820193*D^9+463304871731107566245182861561704529501144319373370247659886602913861204746806279647389839981762462522623858719558089407171092726736039268945746586637893290292906182260294019287151067130638800908672721686545047453532995748861065412*D^8+1129501274467970210045491560025312632412173792504123925747147429055565601841744024966740469864390925836717783963628001317509450842573517965784917270412069371542910281997307701341415556642657293741327519743341468342554408002685637101106*D^7+2294625726678156779482002448921497903592642205070997259959108149329518354708427527465930806870457021958625100543717965072463642217735896553521234237148425672161957956521487373970077213897406082937119225682913743646997559634298088844067089*D^6+3543254462112117998856075746522953205124152838721882137308790862957777875327530331927167609358789604873430910754658000356466388956046148287315886408414090260902772378792569066753409110443360078491702894458325664506338030603051083646032855539*D^5+4091993350785182036425927551684108975525571233627604578709088394717818316044947249994114886961038532337181399567744024761984847405381917140114329094026559794360381974809374664380014888646800683113740858369116646388598572044623501457779580762910*D^4+3608905857686648642384394358290018675490125538252295263928727839222407632653077795621072440941227193433738873216579574727396314164874751044625316935949323083955528611075516235919662951997710730784899200686870516680517824055952839071692398588663352*D^3+2243000810403665903772021457303266059674280943807604477835590601048386506128603528165431841720895396406768174682246587944521045145942846894042775520593000483701831780118683531898737450192205929517163899574345472277332482086959228951259252787685924032*D^2+824805934465597551116025726207704537059676723078506041015300067808264686945934608982603523011000524742977630883453782383941308929000084418178772712183248527723853848938426407880090394218881359795305520676805350028345801691967203223054114392578472462448*D+138073378409557870560796197736118751972173911132574707699078569897825934331002986805619991151434934982404312708506275272254231212726080417501947207356026564833106412374252985030844959377804613046190373197744156837096646072676306978066067361495193076128480: ############################################################################# # # X11 + X12 + X21 + X22 # PolyFSUM:=185810090827102280260637985924710121887074947294363168368289152810789722861872564572306679875584956081958749872979421018359644351115344901761828830073985765770083865641728903875601951576218194815928716689592960*FSUM^14-994821533810904223875391402430184871461645189802826752824549112197143048964501795103114381405524559125234009570187035830864175842112429746392254455267412007150460729309892253000854648736894611502604672016519015744*FSUM^13-120589847992694489022754747871869709516787966561952260876845071324977624387382476001075026119075271895888061088123769337033943999462414550511303574309357531086046918344086650546588941022603657315924305731285551026904*FSUM^12-11845129255706337244112802495102522778749879012738081300830816814796375375113074916387731002574144030346838015549815900386003291479788811035381098856697757599479271036208769876751050916590602111596646113831471644092112480*FSUM^11+100888198569770405982504051862803468004193706028923665540017958852140693887498484322219189982578794376288577288634430845020419827200107841440749204664233071520569033093431509281446309269064169563837432231458394541761627441812*FSUM^10+36479642733958812137534580951992445865132574974493857703773437941763292280543949779805185656107532782521170395676750782167486858363818572922233269941476427166321192649723799590143458049877088799550368481014789795553347131783823*FSUM^9+827723394262424412478370237061369834013770410768431761993637249035821762270141608414242676811975350001204251755326202089397575547358165290609605360681215990413931406786613836889412649128735237686076307872860794942995566224679178310*FSUM^8-1165158837055282005023765248911042563896523592432885859365841158195094518178277766115330795103981759359926661214832138437875467529302035900386403802868046339910760840484421518628419375946811553733018998929676673849687293672987374240259*FSUM^7-908761696248675973448699585243235587614350999610043027488781436040637379266505865130139611078632285797782669673026212696975464265183792802748559525753329587406248139492288999160972642717291988940290546828561224366903524926712661959379981*FSUM^6+267476952501223528934056856374156991139402930564364630186288173308649374845450869515562519385322309039027669095202272115627892300695143664995695103704720990644558953505699480504456741117060643937119872359954136183702604257754901442796260376*FSUM^5+15418217781397606702742949495470447027130987464850238720350040448906452252109365472349251449489722705584123316684054989704303393489457870625775551011609235726856777542981879341915391076629076015122430659709421853015983235313515132219906619201835*FSUM^4-49647681225468106524598371061837968027777619687126110791463537377954882258208797127445920671586491308904672178136821656430741985449034289022523699402895187521996642768340357067744971002169221373565879870736728684607659374654838896134092102281628160*FSUM^3+84003556207871638050349165553876364151554315168343601117565425110979177706838974756315129686502955370254953123640271419378058065991488846631265467594523643518285864792533157660478924061617886465845072849132403753020743091502186850440400562880557584544*FSUM^2-63566049832195688802165610378115615666943482625432106620061671278709233253389267257005553439130179245623200609524466708559340697626710006959822742405389865712474190772178384417311350984739988292654652902930264662831019484939787347134069394530093721174016*FSUM+19147984809382855996476126498136667410632609970271703555439704991618443760770797638035715299238738932481287462234272954409175865892940372672271437653729321807634063228485081317102082649152115973012134216962889265713166347255121946539243975753247146945198848: ############################################################################# # # X31 + X32 + X41 + X42 # PolyLSUM:=61936696942367426753545995308236707295691649098121056122763050936929907620624188190768893291861652027319583290993140339453214783705114967253942943357995255256694621880576301291867317192072731605309572229864320*LSUM^14+618395777636341567445123989888599227902294849017094736796239082020493330650116385257994300634314390033264293053973017398198207056817362574336528816908548315541600310254985735086680573166029145337924562577020835776*LSUM^13+4743817767266833887197018123179066495953071640770397787209428944203351611427156573369160810992775174981440661087065936244050110670755757277877897518153586226538013623692135032918560507786256306925024647291803507179960*LSUM^12+41811380110363275428822705579122342600906074301228906627895050665843876513598504646200246953233532042173484284245419369233595165361192634119645892947616764203413516594903191024710468099055293904514823077717387708083588752*LSUM^11+233689859714235360836485416319183399903633936156374305414216565477924903077053995835930777401037650568711582135709717287776256861217431236589666276865546259075585874106273253244342104748634336093313351511957353250142292862144*LSUM^10+969436489976271769551579760086395719765907876288726587436017609617955292740493888937796199363507275421138998912946286692285078502126994506474531707089149967664839419573756212102121345854956105971147573244995635033647920051136077*LSUM^9+4445189270561696328093541102237412450404590759367012437383771127957696849261528455562378325821440970132875525854817417442010274187028817535377322038364382572145663931205428997452891135606042350795136082546120800666628687656112200264*LSUM^8+13287048609469484675071419375940466582390248497733720798230017031094167638525283139571042532353196849848091687049283099817291527380920363205497339354082414711347948814871284931885247712563318113060890378825088919712034843145347615771895*LSUM^7+36897626489935570851482695514690266030034311991692394995703600688059616833556994159742411959046990462787689220942123154639487548548271113551634499541722364323567400101180734701689883489696687132402463789575910751582699370600002260576367961*LSUM^6+82993188267684199766077829638266621172448363488394950990287606302798384150825124231563048200572281971332730272251695703297167652153692141463070210790564683474258992957872044965444720004267590995163518615462554428773261720600082214948164473996*LSUM^5+140083793192089470089644874272047752058590769966923453623105781246289022920915527672428167608375804933177314384683143612583337767440747506817373491019438544355025725756958169724836394440006434175421751232715788247722017483622257780914651516012773*LSUM^4+185761426477937407733792385076959949630904085945011957595240748654900771042537788795571787622877304436120602253639296525572631451599417572732192555649848056685205022887615260708952711180660776132317150546174871239238870161924642227632194302141453392*LSUM^3+180567431426273058046780221181195309814902683875705704632484068758237395433834622654643999419403790426043499487482194423833609994206962530190801244400105084084941532852654822108395653085337918823654981522864287185016577019116518630076698362953836795408*LSUM^2+100675834858110397901863771050862112295034415588892182771533774502138974346025456910740364399473526015094787227622401772814663755477966682714707442427650553324648568277945475566770597094357911578777977716385628566656423714903340964869884877468277547155584*LSUM+25186730209653788100312620646398482737809896033336816996573102938599347893753722212097911981199477522925216467152666481209475599770325863894798215090298261468849749476457447841811205698946290268414663622164644990553008760840467447968439229981489591008960512: ############################################################################# # # X11+X12 + ... + X42 # PolySUM:=483880444862245521512078088345599275747591008579070750959086335444764903286126470240381978842669156463434244460883908901978240497696210681671429244984337931692926733442002353842713415563068215666481033045815*SUM^14+2240535935151355412656199393842741698555310826167857962406687848605044122358977501747053438795360966079580389561802125947214701375884525980774302071505814946547760941289231646508299663440606313831169311756623676*SUM^13+8264358627874934344570904185719393007315736978197526261760820481139763934952448409484185711465569245379385590792989228064680176660304605127362870497675827604450259083827837212339739683305709479264308354514363890586*SUM^12+12278803462892222214256270721787351437433477136280509604886527082037475912777285032412263267928986322708403564571931468607642801932201757424632605994544781713055184266418093987890333394373645273913521968976783820944842*SUM^11+5540314664408095686242680805559229126161347134272887259283188932946576842808836766190280002890483834738182926047876896011607279126370556128310692705031013418000961863913718141860389458380003207133443385312077229016032987*SUM^10-20521664063735143510123436552629895955719364797886184705012619439269698348288400144844013807191818125543115382692483678996521299362913421909962509882249532221638115499569725425799259410512417418214299053346030200626465877508*SUM^9+15999979998253784933879612924890700718733218372523674281527455210782913996066742737234728600751104739793773838535716539883174926322574278154856092193535084663117116383725403844832297369550173927613557529077456519686504315100888*SUM^8-979349275383126782238020826060842909233457818296526318929903280994388258060480660316088299189881933039178723969387390373195652226622678110475305546286120053509924032328293790335585658167881960658425641756102563878804277846737706*SUM^7-1576112368404431972730423073336027979493734626484131237483452159677640934052519032346903176348006144276799836615469883654474159938530926946956932711159735476319055035922227535487680069554093799760149194002652912456629537529731637436*SUM^6+178030324128221664931052098871389683428849870212084019298355291328632684353088116661959420232921561173045316480184891762202798723505375086571897951237407684314104132023457808684242731084353765304813383435126506414679015899572352343024*SUM^5+143702797994776225868199135517207699689804894346962153096186890255941187702402427930858406515984519341926877907706051451887454819914483935291797678162239192518727385164328582341534803271735664943257360996229369513308209685284847482441472*SUM^4-33343041510049184955203021161213120280168819198183338205452312496385450142577112912502327077963172311076553818146281339982786663042379878335815193044747024105160973792352394098503558956105463327642947709938610833823578540101904006421401856*SUM^3+7416180420082937747877155105356562008320989570775002173055121173986882889041474312012499376413246612869628063422649056800716616028412243102771307974826385520932724835823304803485977717924980406006485976059270237353102087341171710143322979328*SUM^2-536827942664978394488544443004066606101088491883776505153256468944209435201115339955590702439957027055293868432065287585546135961718154007961621999917294733031849983936382083506367113990805483759350527794544142946349167494488466174505733816320*SUM+15767719712200033437595611964331357415013628276855305846468375992399889792935540857431862781055733450308794381745802392923679637706258921391802801403963010101963816570524846620387204655619162470947642124068963471466404447958122471924512470335488: readlib( realroot ): # Checks to see that all of the coordinates of all of the matrices # are complex. realroot(Poly11, 1/10); realroot(Poly12, 1/10); realroot(Poly21, 1/10); realroot(Poly22, 1/10); realroot(Poly31, 1/10); realroot(Poly32, 1/10); realroot(Poly41, 1/10); realroot(Poly42, 1/10); # These are all the coordinates of the solutions; but they are not # matched up yet. R11:=[fsolve(Poly11,X11,complex)]: R12:=[fsolve(Poly12,X12,complex)]: R21:=[fsolve(Poly21,X21,complex)]: R22:=[fsolve(Poly22,X22,complex)]: R31:=[fsolve(Poly31,X31,complex)]: R32:=[fsolve(Poly32,X32,complex)]: R41:=[fsolve(Poly41,X41,complex)]: R42:=[fsolve(Poly42,X42,complex)]: # These are needed for efficient assembling of the solutions RA:=[fsolve(PolyA,A,complex)]: RB:=[fsolve(PolyB,B,complex)]: RC:=[fsolve(PolyC,C,complex)]: RD:=[fsolve(PolyD,D,complex)]: RFSUM:=[fsolve(PolyFSUM,FSUM,complex)]: RLSUM:=[fsolve(PolyLSUM,LSUM,complex)]: RSUM:=[fsolve(PolySUM,SUM,complex)]: S11:=vector(14,0): S12:=vector(14,0): S21:=vector(14,0): S22:=vector(14,0): S31:=vector(14,0): S32:=vector(14,0): S41:=vector(14,0): S42:=vector(14,0): SA:=vector(14,0): SB:=vector(14,0): SC:=vector(14,0): SD:=vector(14,0): SF:=vector(14,0): SL:=vector(14,0): # # Matches the sum of first half with sum of second half # for ii from 1 to 14 do for jj from 1 to 14 do for kk from 1 to 14 do top:=RFSUM[ii]+RLSUM[jj]-RSUM[kk]: if (abs(top)< 1/100000 ) then SF[kk]:= RFSUM[ii]: SL[kk]:= RLSUM[jj]: fi: od:od:od: print(`matched sum of first half with sum of second half`): # # Matches the sum of first two with sum of second two # for ii from 1 to 14 do for jj from 1 to 14 do for kk from 1 to 14 do top:=RA[ii]+RB[jj]-SF[kk]: if (abs(top)< 1/100000 ) then SA[kk]:= RA[ii]: SB[kk]:= RB[jj]: fi: od:od:od: print(`matched sum of first two with sum of second two`): # # Matches the sum of third two with sum of fourth two # for ii from 1 to 14 do for jj from 1 to 14 do for kk from 1 to 14 do top:=RC[ii]+RD[jj]-SL[kk]: if (abs(top)< 1/100000 ) then SC[kk]:= RC[ii]: SD[kk]:= RD[jj]: fi: od:od:od: print(`matched sum of third two with sum of fourth two`): # # Matches the first two coordinates # for ii from 1 to 14 do for jj from 1 to 14 do for kk from 1 to 14 do top:=R11[ii]+R12[jj]-SA[kk]: if (abs(top)< 1/100000 ) then S11[kk]:= R11[ii]: S12[kk]:= R12[jj]: fi: od:od:od: print(`matched first two`): # # Matches the second two coordinates # for ii from 1 to 14 do for jj from 1 to 14 do for kk from 1 to 14 do top:=R21[ii]+R22[jj]-SB[kk]: if (abs(top)< 1/100000 ) then S21[kk]:= R21[ii]: S22[kk]:= R22[jj]: fi: od:od:od: print(`matched second two`): # # Matches the third two coordinates # for ii from 1 to 14 do for jj from 1 to 14 do for kk from 1 to 14 do top:=R31[ii]+R32[jj]-SC[kk]: if (abs(top)< 1/100000 ) then S31[kk]:= R31[ii]: S32[kk]:= R32[jj]: fi: od:od:od: print(`matched third two`): # # Matches the last two coordinates # for ii from 1 to 14 do for jj from 1 to 14 do for kk from 1 to 14 do top:=R41[ii]+R42[jj]-SD[kk]: if (abs(top)< 1/100000 ) then S41[kk]:= R41[ii]: S42[kk]:= R42[jj]: fi: od:od:od: print(`matched fourth two`): for ii from 1 to 14 do K.ii:= matrix(4,2,0): od: for ii from 1 to 14 do K.ii[1,1]:= S11[ii]: K.ii[1,2]:= S12[ii]: K.ii[2,1]:= S21[ii]: K.ii[2,2]:= S22[ii]: K.ii[3,1]:= S31[ii]: K.ii[3,2]:= S32[ii]: K.ii[4,1]:= S41[ii]: K.ii[4,2]:= S42[ii]: od: print(K1); # print(K2); print(K3); # print(K4); print(K5); # print(K6); print(K7); # print(K8); print(K9); # print(K10); print(K11); # print(K12); print(K13); # print(K14); time(); quit