#NEWLINE# def ode_rhs( self , t , y , p ) : #NEWLINE# #INDENT# ydot = self . ydot #NEWLINE# ydot [ 0 ] = ( ( ( - p [ 0 ] ) * y [ 0 ] ) + ( ( p [ 2 ] * y [ 1 ] ) * y [ 2 ] ) ) #NEWLINE# ydot [ 1 ] = ( ( ( p [ 0 ] * y [ 0 ] ) - ( p [ 1 ] * pow ( y [ 1 ] , 2 ) ) ) - ( ( p [ 2 ] * y [ 1 ] ) * y [ 2 ] ) ) #NEWLINE# ydot [ 2 ] = ( p [ 1 ] * pow ( y [ 1 ] , 2 ) ) #NEWLINE# return ydot