Choleski method