Cr-Lohner algorithm