Nice SciPy 2021 talk using SymPy. "Conformal mapping with sympy Towards python driven analytical modeling in physics" youtu.be/P5ybpjv2uDA

import numpy as np import matplotlib.pyplot as plt x = linspace(-5, 5, 10000) y = linspace(-5, 5, 10000) X, Y = meshgrid(x, y) eps = 1e-2 plt.contour(X, Y, (abs(np.sin(2*np.pi*X) - np.cos(2*np.pi*Y)) < eps) & (abs(np.cos(2*np.pi*X) - np.sin(2*np.pi*Y)) < eps), [1])

In SymPy you can use >>> from sympy import * >>> x, y = symbols('x y') >>> plot_implicit(And(Eq(sin(2*pi*x), cos(2*pi*y)), Eq(cos(2*pi*x), sin(2*pi*y))))

Excellent #SciPy2017 talk by Michael Lange on Devito, which is built on SymPy. youtube.com/watch?v=KinmqFTE…

SymPy beginner tutorial from SciPy 2016 video youtube.com/watch?v=AqnpuGbM… sympy.org/scipy-2016-tutoria…

Use init_printing() to get pretty printing of SymPy expressions in the terminal or notebook. docs.sympy.org/latest/tutori…

Have a happy holidays! (Image from Anna Bujak)