>> help gauss_elim x = gauss_elim(A,b) takes a square matrix A and a vector b and returns the solution of Ax=b. % define a random 5x5 matrix >> M = rand(5,5) M = 1.5274e-02 4.1865e-01 8.3812e-01 5.0281e-01 1.9343e-01 7.4679e-01 8.4622e-01 1.9640e-02 7.0947e-01 6.8222e-01 4.4510e-01 5.2515e-01 6.8128e-01 4.2889e-01 3.0276e-01 9.3181e-01 2.0265e-01 3.7948e-01 3.0462e-01 5.4167e-01 4.6599e-01 6.7214e-01 8.3180e-01 1.8965e-01 1.5087e-01 % define a random right hand side >> b = rand(5,1) b = 6.9790e-01 3.7837e-01 8.6001e-01 8.5366e-01 5.9356e-01 % find the solution, I've asked to see the matrices along the way >> x = gauss_elim(M,b) A = 1.5274e-02 4.1865e-01 8.3812e-01 5.0281e-01 1.9343e-01 7.4679e-01 -1.9623e+01 -4.0958e+01 -2.3874e+01 -8.7752e+00 4.4510e-01 -1.1675e+01 -2.3742e+01 -1.4224e+01 -5.3340e+00 9.3181e-01 -2.5338e+01 -5.0752e+01 -3.0370e+01 -1.1259e+01 4.6599e-01 -1.2100e+01 -2.4738e+01 -1.5151e+01 -5.7505e+00 A = 1.5274e-02 4.1865e-01 8.3812e-01 5.0281e-01 1.9343e-01 7.4679e-01 -1.9623e+01 -4.0958e+01 -2.3874e+01 -8.7752e+00 4.4510e-01 -1.1675e+01 6.2618e-01 -1.9260e-02 -1.1315e-01 9.3181e-01 -2.5338e+01 2.1359e+00 4.5748e-01 7.1982e-02 4.6599e-01 -1.2100e+01 5.1877e-01 -4.2837e-01 -3.3927e-01 A = 1.5274e-02 4.1865e-01 8.3812e-01 5.0281e-01 1.9343e-01 7.4679e-01 -1.9623e+01 -4.0958e+01 -2.3874e+01 -8.7752e+00 4.4510e-01 -1.1675e+01 6.2618e-01 -1.9260e-02 -1.1315e-01 9.3181e-01 -2.5338e+01 2.1359e+00 5.2317e-01 4.5794e-01 4.6599e-01 -1.2100e+01 5.1877e-01 -4.1242e-01 -2.4553e-01 A = 1.5274e-02 4.1865e-01 8.3812e-01 5.0281e-01 1.9343e-01 7.4679e-01 -1.9623e+01 -4.0958e+01 -2.3874e+01 -8.7752e+00 4.4510e-01 -1.1675e+01 6.2618e-01 -1.9260e-02 -1.1315e-01 9.3181e-01 -2.5338e+01 2.1359e+00 5.2317e-01 4.5794e-01 4.6599e-01 -1.2100e+01 5.1877e-01 -4.1242e-01 1.1546e-01 x = 2.5807e+00 -1.1877e+00 2.2701e-01 3.7158e+00 -4.6678e+00 >> format short e % check it's a solution >> A*x-b ans = -1.1102e-16 -7.7716e-16 -4.4409e-16 -1.1102e-16 -3.3307e-16 >> quit 1180 flops.