Implementation of the first heuristic

The first heuristic is a method to find the infinitesimals xi and eta, which
are solutions to the PDE eta.diff(x) + (eta.diff(y) - xi.diff(x))*h(x,y)
- xi.diff(y)*h**2 - xi*h.diff(x) - eta*h.diff(y) = 0.

There are four assumptions to eta and xi

1. eta = 0, xi = f(x)
2. eta = 0, xi = f(y)
3. eta = f(x), xi = 0
4. eta = f(y), xi = 0

After finding the infinitesimals, the next step would be to find the
necessary transformation for the simplified coordinates r,s by the following
PDE's

1.r.diff(x)*xi + r.diff(y)*eta = 0
2.s.diff(x)*xi + s.diff(y)*eta = 1