function [B0, B1, b]=bc(l)
% Calcola le matrici B0 e B1, e il vettore r tale che le BC del problema
% y(t)'=A(t,l)y(t)+q(t,l) 0< t< 1
% sono:
% B0 y(0)+ B1 y(1) = b

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

b1=(exp(-l)+2*exp(-2*l)+1)/(2+exp(-l))+1;
b2=(exp(-l)+3)/(2+exp(-l))-1;
b3=(-l*exp(-l)+5*l)/(2+exp(-l));
b= [b1; b2; b3];
