The correct option is (b) \(\frac{∂^2 Φ}{∂z^2}\)
To explain I would say: It is convenient to reduce the compatibility and equilibrium equations into a single equation in terms of stress function Φ.
The stress function for stress in the x-direction is given by,
\(\frac{∂^2 Φ}{∂z^2}.\)