%RUN_SR_ENU
%
f0='%8.0f %12.4e\n';
N=40; 
disp('       m      err')
for m=1:5
  ab1=sr_Enu(60,32,N,m); 
  ab2=sr_Enu(68,32,N,m);
  err=max(max(subs(abs((ab1-ab2)./ab2))));
  fprintf(f0,m,err)
end
%ab1
%xw=sgauss(60,10,ab1(1:10,:))
%xw=sgauss(60,20,ab1(1:20,:))
