Example
4. Solve
.
4 (a). Use the
boundary values
and
.
4 (b). Use the
boundary values
and
.
Solution 4 (b).
Form the linear operator, the differential equation and the boundary values.
![[Graphics:../Images/GalerkinMod_gr_348.gif]](../Images/GalerkinMod_gr_348.gif)
The first function
will be chosen to be the straight line between the boundary
points.
![[Graphics:../Images/GalerkinMod_gr_351.gif]](../Images/GalerkinMod_gr_351.gif)
The other n function
will be chosen to be zero at the endpoints.
![[Graphics:../Images/GalerkinMod_gr_354.gif]](../Images/GalerkinMod_gr_354.gif)
The residual is formed by substituting
into
.
![[Graphics:../Images/GalerkinMod_gr_358.gif]](../Images/GalerkinMod_gr_358.gif)
Galerkin's requirement is that the inner product of the residual
with
is
zero.
Form the equations
for
.
![[Graphics:../Images/GalerkinMod_gr_363.gif]](../Images/GalerkinMod_gr_363.gif)
We must solve for the coefficients
.
![[Graphics:../Images/GalerkinMod_gr_366.gif]](../Images/GalerkinMod_gr_366.gif)
Use
to
form the Galerkin solution.
![[Graphics:../Images/GalerkinMod_gr_369.gif]](../Images/GalerkinMod_gr_369.gif)
![[Graphics:../Images/GalerkinMod_gr_370.gif]](../Images/GalerkinMod_gr_370.gif)
We are done.
Aside. We can use
Mathematica to find the analytic solution. This is
just for fun !
![[Graphics:../Images/GalerkinMod_gr_372.gif]](../Images/GalerkinMod_gr_372.gif)
Plot the analytic solution.
![[Graphics:../Images/GalerkinMod_gr_374.gif]](../Images/GalerkinMod_gr_374.gif)
![[Graphics:../Images/GalerkinMod_gr_375.gif]](../Images/GalerkinMod_gr_375.gif)
Plot both the analytic and Galerkin solution.
![[Graphics:../Images/GalerkinMod_gr_377.gif]](../Images/GalerkinMod_gr_377.gif)
![[Graphics:../Images/GalerkinMod_gr_378.gif]](../Images/GalerkinMod_gr_378.gif)
Check out the difference between the analytic solution and the Galerkin solution.
![[Graphics:../Images/GalerkinMod_gr_379.gif]](../Images/GalerkinMod_gr_379.gif)
![[Graphics:../Images/GalerkinMod_gr_380.gif]](../Images/GalerkinMod_gr_380.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_383.gif]](../Images/GalerkinMod_gr_383.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_392.gif]](../Images/GalerkinMod_gr_392.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_398.gif]](../Images/GalerkinMod_gr_398.gif)
This is the same as we obtained above.
(c) John H. Mathews 2005