Correct option is (b) \(-\frac{∂^2Φ}{∂x∂z}-γx\)
The explanation is: It is convenient to reduce the compatibility and equilibrium equations into a single equation in terms of stress function Φ.
The shear stress function τxz is given by,
\(τ_{xz}=-\frac{∂^2Φ}{∂x∂z}-γx\)