Вычисление двойного интеграла

from scipy.integrate import dblquad
import numpy as np
import math
from math import log
import pandas as pd
import decimal as dc
dc.getcontext().prec = 100
def f_x (x,l,h):
    numinator=(dc.Decimal(x)**(dc.Decimal(l-1)))*((dc.Decimal(math.e))**((dc.Decimal(-1)*dc.Decimal(x))/(dc.Decimal(2)*(dc.Decimal(1)+dc.Decimal(h)))))
    denuminator = ((dc.Decimal(2)**dc.Decimal(l))*(dc.Decimal(math.factorial(l-1)))*((dc.Decimal(1)+dc.Decimal(h))**dc.Decimal(l)))
    return numinator/denuminator
def f_y (y,l):
    numinator=(dc.Decimal(y)**(dc.Decimal(l)-dc.Decimal(1)))*((dc.Decimal(math.e))**(dc.Decimal((-1*y))/dc.Decimal(2)))
    denuminator = (dc.Decimal(math.factorial(l-1)))*(2**l)
    return numinator/denuminator
f = lambda y,x,l,h,m: f_x(x,l,h)*((f_y(y,l))**(m-1))
d={}
m=5000
a=lambda x:0
b=lambda x:x
for l in range(2,88):
    for h in np.arange(0.0, 60, 0.1):
        print(h)
        P_resault,a= dblquad(f,0,np.inf,a,b,args=[l,h,m]);
        print(P_resault)
        if P_resault>=0.99:
            d[l]=10*math.log(h,10)
            print (10*math.log(h,10))
            print(l)
            break
df = pd.DataFrame([d])
df.to_excel('1.xlsx',index=False)
print("finish")