global c alpha beta
N=20; eps0=.5e-4;
alpha=-1/2; 
beta=alpha;
%beta=-alpha;
cc=(linspace(0,.999))'; x=zeros(100,1); w=zeros(100,1);
for n=1:N
  for i=1:100
    c=cc(i);
    ab=r_subjacobi(N,eps0,c,alpha,beta); 
    xw=gauss(N,ab);
    x(i)=xw(n,1);
%    w(i)=xw(n,2);
  end
  plot(cc,x)
%  plot(cc,w)
  hold on
end
