function []=prob_9_3_12()
%
% 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; 

n = 10;
% create the matrix of second derivatives:
A = gen_der2_mat(n);
% create the requested right hand side:
b = zeros(n,1); b(1) = 1; 

omegaList = [1,1.4,1.8,2.2];

for oi=1:length(omegaList),
  fprintf('omega = %f...\n',omegaList(oi)); 
  tic; x_sor  = sor(A,b,[],[],10^(-3),omegaList(oi)), toc  
end