
function N1
syms  a b c d e f g h l m   x         

p=poly2sym([2,-1],x);
norma=int(p*p,x,0,1);
B=p/(sqrt(norma));
uno=int(B,x,0,1);
unon=int(B*B,x,0,1);

%------------------------
% continuity in  x=1/2 
c1=poly2sym([c,d],x);
c2=poly2sym([e,0.5*c+d-0.5*e],x);

normaC=simplify(int(c1*c1,x,0,1/2)+int(c2*c2,x,1/2,1))
ortC1=simplify(int(c1,x,0,1/2)+int(c2,x,1/2,1))
ortC2=simplify(int(B*c1,x,0,1/2)+int(B*c2,x,1/2,1))

%------------------------
d1=poly2sym([g,h],x);
d2=poly2sym([l,m],x);

normaD=simplify(int(d1*d1,x,0,1/2)+int(d2*d2,x,1/2,1))
ortD1=simplify(int(d1,x,0,1/2)+int(d2,x,1/2,1))
ortD2=simplify(int(B*d1,x,0,1/2)+int(B*d2,x,1/2,1))
ortD3=simplify(int(c1*d1,x,0,1/2)+int(c2*d2,x,1/2,1))

return