McCullagh_second_order_pi_wrt_psi_c function

Second order derivative of pi[i, j] wrt psi[i1, j1] and c.