Example
4. Solve
.
4 (a). Use the
boundary values
and
.
Solution 4 (a).
Form the linear operator, the differential equation and the boundary values.
![[Graphics:../Images/GalerkinMod_gr_295.gif]](../Images/GalerkinMod_gr_295.gif)
The first function
will be chosen to be the straight line between the boundary
points.
![[Graphics:../Images/GalerkinMod_gr_298.gif]](../Images/GalerkinMod_gr_298.gif)
The other n function
will be chosen to be zero at the endpoints.
![[Graphics:../Images/GalerkinMod_gr_301.gif]](../Images/GalerkinMod_gr_301.gif)
The residual is formed by substituting
into
.
![[Graphics:../Images/GalerkinMod_gr_305.gif]](../Images/GalerkinMod_gr_305.gif)
Galerkin's requirement is that the inner product of the residual
with
is zero.
Form the equations
for
.
![[Graphics:../Images/GalerkinMod_gr_310.gif]](../Images/GalerkinMod_gr_310.gif)
We must solve for the coefficients
.
![[Graphics:../Images/GalerkinMod_gr_313.gif]](../Images/GalerkinMod_gr_313.gif)
Use
to
form the Galerkin solution.
![[Graphics:../Images/GalerkinMod_gr_316.gif]](../Images/GalerkinMod_gr_316.gif)
![[Graphics:../Images/GalerkinMod_gr_317.gif]](../Images/GalerkinMod_gr_317.gif)
We are done.
Aside. We can use
Mathematica to find the analytic solution. This is
just for fun !
![[Graphics:../Images/GalerkinMod_gr_319.gif]](../Images/GalerkinMod_gr_319.gif)
Plot the analytic solution.
![[Graphics:../Images/GalerkinMod_gr_321.gif]](../Images/GalerkinMod_gr_321.gif)
![[Graphics:../Images/GalerkinMod_gr_322.gif]](../Images/GalerkinMod_gr_322.gif)
Plot both the analytic and Galerkin solution.
![[Graphics:../Images/GalerkinMod_gr_324.gif]](../Images/GalerkinMod_gr_324.gif)
![[Graphics:../Images/GalerkinMod_gr_325.gif]](../Images/GalerkinMod_gr_325.gif)
Check out the difference between the analytic solution and the Galerkin solution.
![[Graphics:../Images/GalerkinMod_gr_327.gif]](../Images/GalerkinMod_gr_327.gif)
![]()
So the Galerkin solution appears to be good.
Caution. We proceed with caution when using Galerkin's method because the linear system might be ill conditioned.
![[Graphics:../Images/GalerkinMod_gr_330.gif]](../Images/GalerkinMod_gr_330.gif)
The condition number of the above system can be determined by Mathematica.
The condition number is small, we should expect a good solution.
Aside. We can check out the matrix form for setting up the Galerkin equations. This is just for fun !
The boundary values
are
used to form
and
there are
equations to solve
for
.
![[Graphics:../Images/GalerkinMod_gr_340.gif]](../Images/GalerkinMod_gr_340.gif)
Since the linear operator is
we have
which must be entered into the integrals on the right hand side.
for
.
The matrix form of the solution is
![[Graphics:../Images/GalerkinMod_gr_346.gif]](../Images/GalerkinMod_gr_346.gif)
This is the same as we obtained above.
(c) John H. Mathews 2005