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

        J['y1','y2'] = -0.2
        J['y1','z'] = np.array([[2.0*params['z'][0], 1.0]])
        J['y1','x'] = 1.0

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