Maths Help
function [A] = problem1data % sparse LU u=[7.1918848102658095e-01 1.1417728306073511e+00 1.5519362755905890e+00 1.3836301143044067e+00 -7.5809225438954519e-01 4.4266274156259322e-01 9.1109790133759794e-01 -1.0740856172401285e+00 2.0176190936920457e-01 7.6286322314263832e-01 -1.2881871974223793e+00 -9.5296182463183576e-01 7.7817450051825032e-01 -6.3310695599488910e-03 5.2448680956852800e-01 1.3642717954238290e+00 4.8203933049125486e-01 -7.8706584502001042e-01 7.5199926612499590e-01 -1.6688791003091180e-01 -8.1622811056064892e-01 2.0940651895061930e+00 8.0152532344038366e-02 -9.3729506929487394e-01 6.3573877984906246e-01 1.6820279944331387e+00 5.9363417229789772e-01 7.9015287591405581e-01 1.0525384669901849e-01 -1.5857856161489242e-01 8.7090745123482771e-01 -1.9475892828364041e-01 7.5474485412664938e-02 -5.2663484487992240e-01 -6.8548438052366756e-01 -2.6838764878490173e-01 -1.1883459645153542e+00 2.4857897968336215e-01 1.0245161722246231e-01 -4.1006876984961096e-02 -2.2475823134718933e+00 -5.1077647633279921e-01 2.4924264127232151e-01 3.6919671446379743e-01 1.7919662051852325e-01 -3.7283412276979996e-02 -1.6033098660032012e+00 3.3937203890256584e-01 -1.3113481457279527e-01 4.8518979244871324e-01 5.9875069775740231e-01 -8.6030886506518703e-02 3.2529204656648913e-01 -3.3514322960904835e-01 -3.2244920375023883e-01 -3.8237412519474984e-01 -9.5337093281709240e-01 2.3357559892250851e-01 1.2352446504741668e+00 -5.7853159974971058e-01 -5.0153736904257773e-01 7.2286351624396594e-01 3.9498485331621903e-02 1.5412787658547042e+00 -1.7010528404698730e+00 -1.0337413076645554e+00 -7.6370763669822572e-01 2.1764259627156815e+00 4.3161209992731958e-01 -4.4376518130617915e-01 2.9995868238615803e-02 -3.1567086967047631e-01 9.7784570477197974e-01 1.8294873913410118e-02 8.1796283374686962e-01 7.0234110036618524e-01 -2.3127056429901954e-01 -1.1368958008529322e-01 1.2794051860964067e-01 -7.9940969968013953e-01 -2.3861195623431747e-01 -8.9463237500257545e-02 -1.0232641710863453e+00 9.3753792337937292e-01 -1.1317192831521843e+00 -7.1070208602827178e-01 -1.1695008619456433e+00 1.0654371444339836e+00 -6.8039374385658702e-01 -1.7257726569231200e+00 8.1319959967304112e-01 1.4418666182923234e+00 6.7227220216041783e-01 1.3866495543565135e-01 -8.5953392675766072e-01 -7.5225055816553088e-01 1.2296150838290758e+00 1.1507543530981155e+00 -6.0802501127071640e-01]; v = [8.0615791615995780e-01 2.1713285248001582e-01 -3.7346107032556553e-01 -8.3203043468849469e-01 2.8686630098358784e-01 -1.8188916235640564e+00 -1.5730514019458928e+00 2.0156656809047844e+00 -7.1982023978393053e-02 2.6289094176857488e+00 -2.4331694903778381e-01 1.7327649426949135e-01 9.2320707695267612e-01 -1.7855299872084196e-01 -5.2170490682295867e-01 1.4319619733911875e+00 -8.7011744642428468e-01 8.0754179199900744e-01 -5.1063464437659145e-01 7.4351413343903050e-01 8.4789830498617613e-01 -8.2990071750100802e-01 5.3299443371653021e-01 1.0328484132990641e+00 -1.0520239654671704e+00 3.6211373842305622e-01 -3.6786616879762153e-02 -1.2276355758077513e+00 -2.7509887388406173e-01 -1.6043524146830923e-01 -1.0835748189281726e+00 -1.9542127496674984e+00 -9.0948725759204407e-01 -5.5787276532498651e-03 -1.7234898877785685e+00 1.2630772536949719e+00 -6.0043337396096219e-01 -2.0639248148829119e+00 1.1091094017167046e-01 1.4876142571023079e+00 5.3002141658683589e-02 1.6198074381301233e-01 -2.6877918764701345e-02 1.7357619145408637e-01 8.8216775480003351e-01 1.8229441850997191e-01 7.5529542234764091e-01 5.0803455271098141e-01 1.3188014285595623e-01 2.8010410501918309e-01 -9.8284755553118153e-01 -9.4408743135651363e-01 -1.3057515808979701e-02 3.5434527592419762e-01 -8.9470882629212700e-01 8.1211102159521542e-01 1.0953735238992339e-01 2.7316439695391792e+00 4.1107915462951367e-01 -1.3068619689101308e+00 3.8380649419224561e-01 4.9950367851879746e-01 -5.1078558701308308e-01 2.3492204825244120e-01 -5.9782511443613606e-01 2.0770818548438701e-02 4.1944275817818127e-01 1.1911037484264235e+00 7.7121424508625180e-01 -2.6442219216807445e+00 2.8542989794149343e-01 8.2609261897717423e-01 -8.1224925350071864e-03 8.5843816230125425e-01 7.7478824003235403e-01 1.3059453260206879e+00 1.2315025330368079e+00 9.5856366647100510e-01 -1.6545476908421763e+00 -9.9039601111699360e-01 6.8523602244719295e-01 -9.7487001383901362e-01 -6.0672563646930089e-01 6.8679374136400895e-01 2.0048856212854318e-02 1.0638005526661041e+00 -1.3410498594340896e+00 4.7951043863175041e-01 -1.6339739932203248e+00 -1.4426650201173274e+00 2.9378142164438520e-01 -1.4036408377460596e-01 -1.1303407455811170e+00 -2.9253805278518674e-01 -5.8253587427217057e-01 -8.9634827661834648e-01 2.4860108242350701e-01 -1.4896633913508939e+00 3.1350866229782715e-01]; d=[-1.4440040658779980e+00 6.1234029223737352e-01 -1.3235029618032093e+00 -6.6157667657724395e-01 -1.4611450503437784e-01 2.4808495822675208e-01 -7.6632557379239005e-02 1.7381699701538176e+00 1.6219718997877377e+00 6.2643589682312906e-01 9.1813590225913552e-02 -8.0760704650874726e-01 -4.6133709421872093e-01 -1.4059691942684509e+00 -3.7452963273531803e-01 -4.7091123634376558e-01 1.7512956264276767e+00 7.5322461750340575e-01 6.4989341476980364e-02 -2.9276425865359323e-01 8.2822847199424537e-02 7.6619148033876217e-01 2.2368498781009816e+00 3.2688672516205608e-01 8.6330394268156485e-01 6.7938660174195831e-01 5.5475762588310240e-01 1.0016304273089194e+00 1.2593651481189680e+00 4.4151065192012819e-02 -3.1413767474565191e-01 2.2670764087427511e-01 9.9669193290845393e-01 1.2159117190131969e+00 -5.4270241586089962e-01 9.1222834229924443e-01 -1.7214113232944259e-01 -3.3595463728095515e-01 5.4148677646411847e-01 9.3211121349110748e-01 -5.7025264752726357e-01 -1.4986053783994235e+00 -5.0345644019042325e-02 5.5302477378873527e-01 8.3497886682452727e-02 1.5775235124089682e+00 -3.3077420319113976e-01 7.9515483199829706e-01 -7.8480015048190344e-01 -1.2631214870993845e+00 6.6665498604277351e-01 -1.3926322050346449e+00 -1.3005624647030269e+00 -6.0502215908433565e-01 -1.4885645044331768e+00 5.5854295799593912e-01 -2.7735429146565144e-01 -1.2936847598522185e+00 -8.8843518137936861e-01 -9.8652007625449556e-01 -7.1617644436624828e-02 -2.4145909345801972e+00 -6.9434901111914316e-01 -1.3913890115665006e+00 3.2964779928237514e-01 5.9854447999166360e-01 1.4717539803566806e-01 -1.0143899070373294e-01 -2.6349808422200045e+00 2.8053414090014619e-02 -8.7631012622268989e-01 -2.6547744791284611e-01 -3.2757791031003242e-01 -1.1582474864169829e+00 5.8005321065211068e-01 2.3975572722305211e-01 -3.5088491838191171e-01 8.9209848362342525e-01 1.5782988358205232e+00 -1.1081739325708244e+00 -2.5931046975712430e-02 -1.1106277156089546e+00 7.5083417769790350e-01 5.0016703215394931e-01 -5.1726088476530196e-01 -5.5920948981830709e-01 -7.5337062086567663e-01 9.2581269305880642e-01 -2.4852031900480015e-01 -1.4983451724662422e-01 -1.2584153777730664e+00 3.1261987583866419e-01 2.6902769360063012e+00 2.8969629211422804e-01 -1.4228030065990676e+00 2.4678571548398118e-01 -1.4357728810904431e+00 1.4857312713838181e-01 -1.6930725688350658e+00]; A = [1 v'; u diag(d)]; end