function superN0

clear all
clc

syms  a1 b1 c1 d1  a2 b2 c2 d2  a3 b3 c3 d3  x         

A=1;
B=(2*x-1)*3^(1/2);

%------------------------
C1=poly2sym([0,a1],x);
C2=poly2sym([0,b1],x);
C3=poly2sym([0,c1],x);
C4=poly2sym([0,d1],x);

normaC=simplify(int(C1*C1,x,0,1/4)+int(C2*C2,x,1/4,1/2)+int(C3*C3,x,1/2,3/4)+int(C4*C4,x,3/4,1)) 
ortC1=simplify(int(C1,x,0,1/4)+int(C2,x,1/4,1/2)+int(C3,x,1/2,3/4)+int(C4,x,3/4,1)) 
ortC2=simplify(int(B*C1,x,0,1/4)+int(B*C2,x,1/4,1/2)+int(B*C3,x,1/2,3/4)+int(B*C4,x,3/4,1)) 

%------------------------
D1=poly2sym([0,a2],x);
D2=poly2sym([0,b2],x);
D3=poly2sym([0,c2],x);
D4=poly2sym([0,d2],x);

normaD=simplify(int(D1*D1,x,0,1/4)+int(D2*D2,x,1/4,1/2)+int(D3*D3,x,1/2,3/4)+int(D4*D4,x,3/4,1)) 
ortD1=simplify(int(D1,x,0,1/4)+int(D2,x,1/4,1/2)+int(D3,x,1/2,3/4)+int(D4,x,3/4,1)) 
ortD2=simplify(int(B*D1,x,0,1/4)+int(B*D2,x,1/4,1/2)+int(B*D3,x,1/2,3/4)+int(B*D4,x,3/4,1)) 
ortD3=simplify(int(C1*D1,x,0,1/4)+int(C2*D2,x,1/4,1/2)+int(C3*D3,x,1/2,3/4)+int(C4*D4,x,3/4,1)) 

%------------------------
E1=poly2sym([0,a3],x);
E2=poly2sym([0,b3],x);
E3=poly2sym([0,c3],x);
E4=poly2sym([0,d3],x);

normaE=simplify(int(E1*E1,x,0,1/4)+int(E2*E2,x,1/4,1/2)+int(E3*E3,x,1/2,3/4)+int(E4*E4,x,3/4,1)) 
ortE1=simplify(int(E1,x,0,1/4)+int(E2,x,1/4,1/2)+int(E3,x,1/2,3/4)+int(E4,x,3/4,1)) 
ortE3=simplify(int(C1*E1,x,0,1/4)+int(C2*E2,x,1/4,1/2)+int(C3*E3,x,1/2,3/4)+int(C4*E4,x,3/4,1)) 
ortE4=simplify(int(D1*E1,x,0,1/4)+int(D2*E2,x,1/4,1/2)+int(D3*E3,x,1/2,3/4)+int(D4*E4,x,3/4,1)) 

return