def linearize(self, params, unknowns, resids):
        J = {}

        J['y2', 'y1'] = .5*params['y1']**-.5
        J['y2', 'z'] = np.array([[1.0, 1.0]])

        return J
Jacobian for Sellar discipline 2.
dataset/ETHPy150Open OpenMDAO/OpenMDAO/openmdao/test/sellar.py/SellarDis2withDerivatives.linearize