Use Matlab to find the PA=LU of the vector A and forward substitution? to find the first 100 values of q for an evenly distributed pressure value and plot the q-values (on the same graph) as q is a funtion of .

What I have so far is:

A = [1 -1 0 0 -1 0 0; 0 1 -1 -1 0 0 0; 0 0 1 0 1 -1 0; 0 0 0 1 0 1 -1; 2 1 0 1 0 0 2; 0 -1 0 -1 1 1 0; 0 -1 -sqrt(2) 0 1 0 0];

for p0 = 0.035:0.0008585:0.12; L = 1; d = 0.2; f = 0.02; p = 998; mu = 1.002*10^(-3); QA = 48*10^(-5);

b5 = p0*(pi*d^4./(128*mu*L)); b = [0; 0; QA; 0; b5; 0; 0]; [L,U,P]=lu(A); Pb = P*b; c = L\Pb; Pb = P*b; c = L\Pb; q = U\c %first 100 values of q end

That gives the first 100 values of q but I have problem of plotting the graph

Po 0.035 Pa, 0.12Pa