Как вы просили о похожих шаблонах на любом языке, вот код Mathematica (действительно легко, поскольку Wolfram Alpha основан на Mathematica)
Редактировать
Это действительно эффект округления:
Если мы установим:
и составить сюжет
Plot3D[f[x, y], {x, 7, 9}, {y, -8, -9},WorkingPrecision -> MachinePrecision]
Результат:
Но если мы увеличим точность графика до 30 цифр:
Plot3D[f[x, y], {x, 7, 9}, {y, -8, -9},WorkingPrecision -> 30]
Получаем
и шероховатость исчезла (из-за чего вы набросали рисунок)
Кстати, ваш f [x, y] - очень хорошая функция:
Так что, если мне удалось скопировать ваши формулы без ошибок (что следует считать чудом), обе стороны вашего уравнения равны только в определенных периодических диапазонах по x, возможно, в форме [2 n Pi, (2 n + 1) Пи]