smooth — smoothing by spline functions
[pt]=smooth(ptd [,step])
this function computes interpolation by spline functions for a given
set of points in the plane. The coordinates are
(ptd(1,i),ptd(2,i)). The components
ptd(1,:) must be in ascending order. The default value
for the step is
abs(maxi(ptd(1,:))-mini(ptd(1,:)))/100