float oscil(amp,si,farray,len,phs)
float amp,si,*farray,*phs;
register len;
{
	register i =  *phs;   
	*phs += si;            
	while(*phs >= len)
	       *phs -= len;     
	return(*(farray+i) * amp);
}