I have attached the LU decomposition. I don't know how to construct a for loop to compute Ax=b

(c) In a for loop, use the LU decomposition found in part (d) to solve the system Ax=b for five different, values of b : b1=10000,b2=01000,b3=00100,b1=00010,b5=00001. Store cach solution x as a column in a square matrix called X_e, such that X4=(x1x2x3x4x5) LuA= 3.00000.33330001.00002.66670.37500001.00002.62500.38100001.00002.61900.38180001.00002.6182