content_image/371637.jpg content_image/371643.jpg content_image/371642.jpg content_image/371641.jpg content_image/371640.jpg content_image/371639.jpg content_image/371638.jpg
where α = A, B, C . One possible reversible measure-preserving integration algorithm for the BKNH thermostat can then be derived from the following Trotter factorization: 

<img src='content_image/371637.jpg'>

The direct translation technique gives the following pseudo-code: 

<img src='content_image/371643.jpg'>

<img src='content_image/371642.jpg'>

<img src='content_image/371641.jpg'>

<img src='content_image/371640.jpg'>

<img src='content_image/371639.jpg'>

<img src='content_image/371638.jpg'>

