float scpval(float[],long,float);
void spfit(long,float[],float[],float[],long,LOGICAL32,LOGICAL32,LOGICAL32,float[],
      long*,float*,float*);