% zkpol(fctid,fctalp,xmin,xmax,ndeg,npts)
%
% fctid         = ID of interpolated function
% fctalp        = parameter alpha in interpolated function
% xmin,xmax     = interpolation interval (where we construct polynomials)
% aprmin,aprmax = evaluation interval (if > xmin,xmax, we extrapolate)
% ndeg(1:3)     = degrees of interpolation polynomials
% npts          = number of points for plotting

zkpol(1, 4.0, 0.0,6.0, -6.0,12.0, [2 4 10], 30);
pause

zkpol(2, 1.0, 0.0,1.0, -1.0,2.0, [3 6 12], 60);
pause

zkpol(3, 10.0, -2.0,2.0, -2.0,2.0, [5 8 11], 100);

