function []=prob_9_3_10()
%
% Written by:
% -- 
% John L. Weatherwax                2006-12-20
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----
  
%clc; close all; 

nArray = [ 10 20 50 ]; 
for ni=1:length(nArray),
  % create the matrix of second derivatives:
  A = gen_der2_mat(nArray(ni));
  % create the requested right hand side:
  b = zeros(nArray(ni),1); b(1) = 1; 
  %tic; x_jac  = jacobi(A,b,[],[],10^(-3));, toc
  tic; x_gsed = gseidel(A,b,[],[],10^(-3));, toc
end