Example 3.  Solve Laplace's equation with the Neuman boundary condition over a  9 by 9  grid with boundary conditions
    Top:            180
    Left:              80
    Bottom:    [Graphics:Images/EllipticPDEMod_gr_54.gif] = 0
    Right:             0

Solution 3.

First set up the boundary values.

 

[Graphics:../Images/EllipticPDEMod_gr_55.gif]

Next, solve it.

[Graphics:../Images/EllipticPDEMod_gr_56.gif]

[Graphics:../Images/EllipticPDEMod_gr_57.gif]

Plot the solution.

[Graphics:../Images/EllipticPDEMod_gr_58.gif]


[Graphics:../Images/EllipticPDEMod_gr_59.gif]

[Graphics:../Images/EllipticPDEMod_gr_60.gif]

[Graphics:../Images/EllipticPDEMod_gr_61.gif]


To see the numerical values enter the command:

[Graphics:../Images/EllipticPDEMod_gr_62.gif]



[Graphics:../Images/EllipticPDEMod_gr_63.gif]

[Graphics:../Images/EllipticPDEMod_gr_64.gif]

[Graphics:../Images/EllipticPDEMod_gr_65.gif]

 

 

 

We can make a contour plot of the solution.  However, to get the orientation like the table and the 3D figure requires reversing the rows in the matrix.

[Graphics:../Images/EllipticPDEMod_gr_66.gif]


[Graphics:../Images/EllipticPDEMod_gr_67.gif]

[Graphics:../Images/EllipticPDEMod_gr_68.gif]

[Graphics:../Images/EllipticPDEMod_gr_69.gif]

This problem has smoother contours.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004