Ugyanaz kevesebb számítással

    y=1/sqrt(1-v1*v1/(c*c))
    p1x=m*v1*y
    E1=m*c*c*y
    p1y=E1/c
    
    y=1/sqrt(1-v2*v2/(c*c))
    p2x=m*v2*y
    E2=m*c*c*y
    p2y=E2/c
    

    fi1=atan(v1/c)
    fi2=atan(v2/c)
    dfi=fi2-fi1

   py= cos(fi2)*p2y + sin(fi2)*p2x
    px=-sin(fi1)*p2y + cos(fi1)*p2x

 
    le=h/py
    d=h/px/2
    fi3=asin(le/(2*d))