Consider a PDE of the form:
Consider the finite dimensional space of functions
We seek a solution of the form
In general, this does not exist,
because there is no reason for
So we only seek an approximation which satisfies:
Now we just need to pick our
and ideally an inner product where they are orthogonal.
We can use Fourier Series or Chebyshev Polynomials
For Chebyshev expansion, in general, boundary conditions may be a problem.
We might need to use
If there is time dependence on boundary conditions,
we may run into serious stability issues.