Я пытаюсь сделать это с симпатией:
from sympy.abc import *
from sympy import *
from sympy.plotting import *
plot3d(fourier_series(x**2*y,(x,-pi,pi),(y,-pi,pi)).truncate(2),(x,-pi,pi),(y,-pi,pi))
и то, что я получаю, кажется неверным.
документ не упоминает этот случай.