function y=es_bvp_bc(u,v)
% Definisce la funzione g(u,v) tale che le BC sono g(x0,x1)=0
% g(u,v)=B0*u+B1*v-b

l=10;
B0=zeros(3);
B0(1,1)=1;
B1=zeros(3);

B1=zeros(3); B1(2,1)=1; B1(3,2)=1;

% Le BC sono B0 x(0) + B1 x(1)= b
b1=(exp(-l)+exp(-l*2)+1)/(2+exp(-l))+1;
b2=0;
b3=(-l*exp(-l)+3*l)/(2+exp(-l));
b=[b1;b2;b3];
y=B0*u+B1*v-b;