## Numerical solution of partial differential equations

### Quick description

Partial differential equations are considerably more complex to solve than ordinary differential equations, and there are a substantial number of special techniques developed to handle them. All involve some sort of "grid" that subdivides space into small packets. The different ways in which the packets are used gives the different methods.

### Prerequisites

Multivariate calculus, basic numerical analysis.

### Example 1: Galerkin

The most common (especially from a mathematical point of view) is the Galerkin method. The basic idea is that we set up the weak formulation of the partial differential equation. For a (linear) operator equation , where is an operator ( is the dual space to ), we take a basis of a finite dimensional subspace , with .

Then we choose where for . If is an elliptic partial differential operator, then the linear system generated by the Galerkin method is positive definite, and the linear system can be solved.

This method is also related to the Rayleigh-Ritz method, and the more general Petrov-Galerkin method.

The Galerkin method is commonly known as the finite element method (FEM), at least under common choices of the basis functions.