Add a function for Pell equation x^2-Dy^2 = +-4.