I would do it as a sequence of substitutions (planck42's use of a single substitution is correct and faster but harder to see).
First, because that "[itex](x- z)^2[/itex]" is a nuisance, let u= x- z. Then du= -dz so the integral becomes
[tex]-\int \frac{du}{y^2+ u^2}[/tex].
Now, remembering that [itex]sin^2(\theta)+ cos^2(\theta)= 1[/itex] so that [itex]tan^1(\theta)+ 1= sec^2(\theta)[/itex] and that [itex]y^2+ u^2= y^2(1+ (u/y)^2[/itex], let [itex]u/y= tan(\theta)[/itex] so that [itex]du= y^2 sec^2(\theta)d\theta[/itex] and [itex]y^2+ u^2= y^2(1+ (u/y)^2)= y^1(sec^2(\theta))[/itex] (Again, Planck42's hyperbolic substitution works fine but I learned trig substitutions before hyperbolic substutions so I tend to think of them first!). That makes the integral
[tex]\int\frac{y^2 sec^2(\theta)d\theta}{y^2 sec^2(\theta)}[/tex]