Compare results/performance with tridiagonal Gaussian elimination solver for the problem arising from -y''=f on (0,1) with y(0)=0=y(1). You may also need to use sparse storage and MATLAB's built in '\' operator as a third solver.