function []=prob_9_1_4() % % Written by: % -- % John L. Weatherwax 2006-12-14 % % email: wax@alum.mit.edu % % Please send comments and especially bug reports to the % above email address. % %----- clc; close all; disp( 'the (8x8) Hilbert matrix' ); A = hilb(8) format long; disp( 'the eigenvalues of hilb(8)' ); [E]=eig(A) disp( 'the norm of hilb(8) is given by' ); max(E) disp( 'the eigenvalues of inv(hilb(8))' ); 1./E disp( 'the norm of invhilb(8) is given by' ); max(1./E)