#include <ncbi.h>
#include "karlin.h"

/* stoprob -- calculate the probability of achieving a particular score */
double
stoprob(S, Lambda, K, H, dblen, qlen)
	Score_t	S;
	double	Lambda, K, H; /* parameters for Karlin statistics */
	unsigned long	dblen;	/* length of database */
	unsigned long	qlen;	/* length of query sequence */
{
	return etop(stoe(S, Lambda, K, H, dblen, qlen));
}