Maths Help

Maths help Asap

function [L,U] = mylu(A);
% LU without pivoting
n = size(A,1); if (size(A,2) ~= n), error(‘A must be square.’); end;
L = eye(n); U = zeros(n,n);
for i=1:(n-1)
if (abs(A(i,i)) < 1e-10), error('Zero pivot. A cannot be factored as LU.'); end; U(i,i:n) = A(i,i:n); L((i+1):n,i) = A((i+1):n,i)/A(i,i); A((i+1):n,(i+1):n) = A((i+1):n, (i+1):n) - L((i+1):n,i)*U(i,(i+1):n); end; if (abs(A(n,n)) < 1e-10), error('Zero pivot. A cannot be factored as LU.'); end; U(n,n) = A(n,n);

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