float spsi(float);
void spsie(float*,long*);