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

	if ((left > 1.0-fuzz && left < 1.0+fuzz) &&
	    (rite > 1.0-fuzz && rite < 1.0+fuzz))
		return(1.0);
	return(log(rite)/log(left));
}