For four interior points (h 3/5), the resulting system of equations is
Setting the determinant equal to zero and expanding it gives
which can be solved for the first four eigenvalues
(b) The system can be normalized so that the coefficient of the eigenvalue is unity:
MATLAB solution:
(c) Power method:
The first 4 iterations are
0 0.2 0.4 0.6 0.8 1
0
20
40
60
80
100
120
y1
y2