%KEGEL
%
f0='%8.0f %5.1f %12.8f %12.8f %12.8f\n';
f1='%27.8f %12.8f %12.8f\n';
f2='%8.0f %5.1f %12.8f %12.8f %12.8f %12.8f\n';
f3='%27.8f %12.8f %12.8f %12.8f\n';
disp('       m    c                      nodes/weights')
for m=1:2
  for c=[.1:.1:1 1.2:.2:4 4.5:.5:8 9:15]
    ab=r_Emfin(4,m,c);
    xw=gauss(3,ab); x=xw(:,1); w=xw(:,2);
    fprintf(f0,m,c,x(1),x(2),x(3))
    fprintf(f1,w(1),w(2),w(3))
    xw=gauss(4,ab); x=xw(:,1); w=xw(:,2);
    fprintf(f2,m,c,x(1),x(2),x(3),x(4))
    fprintf(f3,w(1),w(2),w(3),w(4))
    fprintf('\n')
  end
end
