#NEWLINE# def linearize( self , params , unknowns , resids ) : #NEWLINE# #INDENT# ' Jacobian for Sellar discipline 1.' #NEWLINE# J = { #NEWLINE# #INDENT# #NEWLINE# #UNINDENT# } #NEWLINE# J [ ( 'y1' , 'y2' ) ] = ( - 0.2 ) #NEWLINE# J [ ( 'y1' , 'z' ) ] = np . array ( [ [ ( 2.0 * J [ 'z' ] [ 0 ] ) , 1.0 ] ] ) #NEWLINE# J [ ( 'y1' , 'x' ) ] = 1.0 #NEWLINE# return J