Энергетические характеристики лидара на python

print("Введите растояние между осями приемника и передатчика Po(Po):")
Po = input("\tPo = ")
Po_int = int(Po)
 
print("Введите угол между осями приемника и передатчика a(a):")
a= float(input("\ta = "))#угол между осями передатчика и приемника 
 
q=150
Rt=25 # Аппертура передатчика (мм)
Qt=150 #Расходимость передатчика
Rr=25 #Апертура приемника
Qr=150*10**-6 #Поле зрения приемника
dx=10000 # пространственное разрешение
Xmax=10000 #дистанция
N=1000 #число фотонов
Mi=0
import random
Nx = int(Xmax/dx)
N = int(N*1)
for i in range(1,Nx,10):
    x = i*dx
    for j in range (1,N,10): #Нахождение случайных точек по кол-ву N(числу фотонов)
        y1t_x= -Rt-x*Qt
        y2t_x= Rt+x*Qt
        y1r= Po_int-Rr
        y2r= Po_int+Rr
        z1t_x=-Rt-x*Qt
        z2t_x=Rt+x*Qt
        z1r=Po_int-Rr
        z2r=Po_int+Rr
        first=(Rt+x*Qt)-(-Rt-x*Qt)-Rt
        second=(Po_int+Rr)-(Po_int-Rr)-Rt+Po_int
        third=(Rt+x*Qt)-(-Rt-x*Qt)-Rt
        fourth=(Po_int+Rr)-(Po_int-Rr)-Rt+Po_int    
        yt= random.randint(0,first)
        yr= random.randint(0,second)
        zt= random.randint(0,third)
        zr= random.randint(0,fourth)
        Gauss= 2.72**(((yt**2)+(zt**2))/2*q*1) #Знаю что есть отдельная формула для гауса, но так мне привычней.
    if Qt==0*10**-6:
        wgt=Gauss
    else:
        wgt=((Rt+x*Qt)-(-Rt-x*Qt))/Qt**2
    if math.sqrt(((yr-yt)**2)+((zr-zt)**2)/x)-a <= Qr:
        Mi=Mi+wgt
    Mi=Mi/N
print=Mi