J =
 
[ (r*((m*(b + h))/(k*(a*h - a*g*m*(b + h))) - 1))/((a*g*m*(b + h))/(a*h - a*g*m*(b + h)) + 1) - r*((2*m*(b + h))/(k*(a*h - a*g*m*(b + h))) - 1),       0, -(m*(b + h))/h]
[                                                                                                -r*((m*(b + h))/(k*(a*h - a*g*m*(b + h))) - 1), - b - h,  (m*(b + h))/h]
[                                                                                                                                             0,       h,             -m]