gauss-jackson method