Starting diary file "diary-2016-03-02.txt" on Wed Mar 02 2016 12:27 PM.

A1 = laplacian1d(4)
A1 =
     2    -1     0     0
    -1     2    -1     0
     0    -1     2    -1
     0     0    -1     2

A1 = laplacian1d(8)
A1 =
     2    -1     0     0     0     0     0     0
    -1     2    -1     0     0     0     0     0
     0    -1     2    -1     0     0     0     0
     0     0    -1     2    -1     0     0     0
     0     0     0    -1     2    -1     0     0
     0     0     0     0    -1     2    -1     0
     0     0     0     0     0    -1     2    -1
     0     0     0     0     0     0    -1     2

Steady_State_Heat(0,10,0,11*exp(10),100);

edit Steady_State_Heat.m  % Uncommenting the nonzero source term & the analytic solution plot

Steady_State_Heat(0,10,0,11*exp(10),100);

Heat(-pi/2, pi/2, 100, .1, 5);

edit Heat % Changing the right boundary condition to 1

Heat(-pi/2, pi/2, 100, .1, 5);

A2 = laplacian2d_dense(3)
A2 =
  Columns 1 through 8
     4    -1     0    -1     0     0     0     0
    -1     4    -1     0    -1     0     0     0
     0    -1     4     0     0    -1     0     0
    -1     0     0     4    -1     0    -1     0
     0    -1     0    -1     4    -1     0    -1
     0     0    -1     0    -1     4     0     0
     0     0     0    -1     0     0     4    -1
     0     0     0     0    -1     0    -1     4
     0     0     0     0     0    -1     0    -1
  Column 9
     0
     0
     0
     0
     0
    -1
     0
    -1
     4
                   |
A2 = laplacian2d_sparse(10);
spy(A2)

BC = @(x,y,t) 0
BC = 
    @(x,y,t)0
Initial_Data = @(x,y,t) 0;
result = Heat_2D(.1, @Source_Term, Initial_Data, BC, -1, 1, -1, 1, 50, 50, .1, 4);