I have $L(u,v,s',t')=\int_{-\infty}^{\infty}E(s,t,z)\delta(\frac{z}{z_0}s+(1-\frac{z}{z_0})u-Ms', \frac{z}{z_0}t+(1-\frac{z}{z_0})v-Mt')dz$
Where $E(s,t,z)=\delta(s)\delta(t)f(z)$, basically a line in 3D volume
Simplifying the integral knowing s=0 and t=0 I get
$L(u,v,s',t')=\int_{-\infty}^{\infty}f(z)\delta((1-\frac{z}{z_0})u-Ms', (1-\frac{z}{z_0})v-Mt')dz$
So my confusion is I wish to numerically evaluate this integral in python but given the delta(x,y) are a function of z, when I solve for z I get
$z=z_0(\frac{Ms'}{u}-1)=z_0(\frac{Mt'}{v}-1)$ leaving me unsure what to use to simplify like $\int f(x)\delta(x_0)dx=f(x_0) $ would.