void midexp(double funk[],double aa,double bb,double s,int n)
{
	double b,a,s,sum,x,del,ddel;
	int it,tnm,j;
	func(x)=funk[-log(x)]/x;
	b=exp(-aa);
	a=0.0;
	if (n==1)
	{
		s=(b-a)*func(0.5*(a+b));
		it=1;
	}
	else
	{
		tnm=it;
		del=(b-a)/(3.0*tnm);
		ddel=del+del;
		x=a+0.5*del;
		sum=0.0;
		for (j=1; j<=it; j++)
		{
			sum=sum+func(x);
			x=x+ddel;
			sum=sum+func(x);
			x=x+del;
		}
		s=(s+(b-a)*sum/(double)tnm)/3.0;
		it=3*it;
	}
}