\documentclass{article}
\usepackage{hyperref}
\usepackage[T1]{fontenc}
\usepackage{pythontex}
\usepackage{amsmath}
\title{Title of Document}
\author{Name of Author}
\begin{document}
\section{Python}
\begin{pycode}
import math
print(r'2 * Pi = {}'.format(2*math.pi))
\end{pycode}
The result of ... is \py{2+5}.
$$\frac{(1+a^2-b^2)}{1-a+a^2+b^2}$$
$\begin{aligned}
\nabla \cdot \mathbf{E} &= \frac{\rho}{\varepsilon_{0}} \\
\nabla \cdot \mathbf{B} &= 0 \\
\nabla \times \mathbf{E} &= -\frac{\partial \mathbf{B}}{\partial t} \\
\nabla \times \mathbf{B} &= \mu_{0}\mathbf{J} + \mu_{0}\varepsilon_{0}\frac{\partial \mathbf{E}}{\partial t}
\end{aligned}$
Note: The variables in Maxwell's equations are defined as follows:
\begin{itemize}
\item $\mathbf{E}$ is the electric field
\item $\mathbf{B}$ is the magnetic field
\item $\rho$ is the charge density
\item $\mathbf{J}$ is the current density
\item $\varepsilon_{0}$ is the permittivity of free space
\item $\mu_{0}$ is the permeability of free space
\item $\nabla$ is the del operator, representing spatial derivatives
\item $\frac{\partial}{\partial t}$ is the partial derivative with respect to time.
\end{itemize}
\end{document}
