%------------------------------------------------------------------------------
%
%  Syfte:    Löser, med cG(1)-metoden, differentialekvationen 
%                                                                              
%          -u''(x) + u(x) = f(x),    a < x < b,                    
%
% med Robinrandvillkor                                               
%           
%           u' = gamma*(u - g_D) + g_N, om x = a,   
%          -u' = gamma*(u - g_D) + g_N, om x = b.   
%
%-----------------------------------------------------------------------------

clear all        % Tar bort alla variabler från Matlab:s workspace.

h = 0.05;
p = 0:h:1;

[S, M, R, v, r] = PoissonAssembler1D(p);

A = S + M + R;
l = v + r;

U = A\l;
 
plot(p,U)
xlabel('x')
title('U(x)')

%------------------------------------------------------------------------------
%
% Pröva ditt program genom att lösa 
% 
%           -u''+ u = x-x.^2+2, u(0) = u(1) = 0. 
%
% Jämför med den exakta lösningen u(x) = x(1-x).