/* Copyright (C) 1994 by Thomas Glen Smith.  All Rights Reserved. */
/* logrithp APL2 V1.0.0 ************************************************
* Returns the base-n log of rite, left giving the base-n.              *
***********************************************************************/
#define INCLUDES MATH
#include "includes.h"
void logrithp(left,rite,ret)
double *left,*rite,*ret;
{
	Log;
	extern double fuzz;

	if ((*left > 1.0-fuzz && *left < 1.0+fuzz) &&
	    (*rite > 1.0-fuzz && *rite < 1.0+fuzz))
		*ret = 1.0;
	else *ret =(log(*rite)/log(*left));
}