solve xφ + y(-1/φ) - c - 3c = d and xφ^2 + y(-1/φ)^2 - 2c - 3c = d