IDENTIFY THE ERROR SO THAT THE CODE WILL HAVE NO ERROR
USING SUCCESSIVE OVER RELAXATION METHOD
function [x,numIter,omega] =
gaussSeidel(func,x,maxIter,epsilon)
% Solves Ax = b by Gauss-Seidel method with
relaxation.
% USAGE: [x,numIter,omega] =
gaussSeidel(func,x,maxIter,epsilon)
% INPUT:
% func = handle of function that returns improved x
using
% x = starting solution vector
% maxIter = allowable number of iterations (default is
500)
% epsilon = error tolerance (default is 1.0e-9)
% OUTPUT:
% x = solution vector
%...