* corrected S matrix, note that the minus sign is a bit weird
* Jacobian computation