How to find the mixed partial derivatives on interpolating points when performing local bicubic surface interporlation?

Question

I am trying to perform a local bicubic surface interpolation based on the method stated in the book The Nurbs Book (1995), pg.402.

I am having diffic