%
% risolve il laplaciano misto con i p1p0 (RT0)
%
% [x,s,u]=p1p0(a,b,n)
%
function [x,s,u]=p1p0(a,b,n)
h=(b-a)/n;
x=(a:h:b)';
%
% costruzione di A
%
A=zeros(n+1,n+1);
A=1/6*h*diag(ones(1,n),-1)+2/3*h*diag(ones(1,n+1))+1/6*h*diag(ones(1,n),1);
A(1,1)=1/3*h; A(n+1,n+1)=1/3*h;
%
% costruzione di B
%
B=zeros(n,n+1);
for i=1:n
  B(i,i)=-1;
  B(i,i+1)=1;
end
%
% la matrice M
%
M=[A B';B zeros(n,n)];
%
% il termine noto
%
F=zeros(n+1+n,1);
F(n+1+1:n+1+n)=-g((x(1:n)+x(2:n+1))/2)*h;
%
% il calcolo della soluzione
%
X=M\F;
s=X(1:n+1);
u=X(n+2:2*n+1);
