function [] =quickcheck(A2,B2,C2,D2, A1,B1,C1,D1);

clf;

cA2=conv(conv(A2.^2,[1 1]/2),[1 1 1 1 1 1]/6);
cB2=conv(conv(A2.^2,[1 1]/2),[1 1 1 1 1 1]/6);
cC2=conv(conv(C2.^2,[1 1]/2),[1 1 1 1 1 1]/6);
cD2=conv(conv(D2.^2,[1 1]/2),[1 1 1 1 1 1]/6);
S2=(cA2+cB2+cC2+cD2).^0.5;


cA1=conv(conv(A1.^2,[1 1]/2),[1 1 1 1 1 1]/6);
cB1=conv(conv(B1.^2,[1 1]/2),[1 1 1 1 1 1]/6);
cC1=conv(conv(C1.^2,[1 1]/2),[1 1 1 1 1 1]/6);
cD1=conv(conv(D1.^2,[1 1]/2),[1 1 1 1 1 1]/6);
S1=(cA1+cB1+cC1+cD1).^0.5;


plot(S1,'.-b'); hold on; plot(S2,'.-r'); axis([1200 1800 0 6000]);