n18=filter2(q, n18);n18=filter2(q, n18);n18=filter2(q, n18);
n19=filter2(q, n19);n19=filter2(q, n19);n19=filter2(q, n19);
n20=filter2(q, n20);n20=filter2(q, n20);n20=filter2(q, n20);
n21=filter2(q, n21);n21=filter2(q, n21);n21=filter2(q, n21);	% Apply Filter
n22=filter2(q, n22);n22=filter2(q, n22);n22=filter2(q, n22);
n23=filter2(q, n23);n23=filter2(q, n23);n23=filter2(q, n23);
n24=filter2(q, n24);n24=filter2(q, n24);n24=filter2(q, n24);
n25=filter2(q, n25);n25=filter2(q, n25);n25=filter2(q, n25);
n26=filter2(q, n26);n26=filter2(q, n26);n26=filter2(q, n26);
n27=filter2(q, n27);n27=filter2(q, n27);n27=filter2(q, n27);
n28=filter2(q, n28);n28=filter2(q, n28);n28=filter2(q, n28);
n29=filter2(q, n29);n29=filter2(q, n29);n29=filter2(q, n29);
n30=filter2(q, n30);n30=filter2(q, n30);n30=filter2(q, n30);
n31=filter2(q, n31);n31=filter2(q, n31);n31=filter2(q, n31);	% Apply Filter
n32=filter2(q, n32);n32=filter2(q, n32);n32=filter2(q, n32);
n33=filter2(q, n33);n33=filter2(q, n33);n33=filter2(q, n33);
n34=filter2(q, n34);n34=filter2(q, n34);n34=filter2(q, n34);
n35=filter2(q, n35);n35=filter2(q, n35);n35=filter2(q, n35);
n36=filter2(q, n36);n36=filter2(q, n36);n36=filter2(q, n36);
n37=filter2(q, n37);n37=filter2(q, n37);n37=filter2(q, n37);
n38=filter2(q, n38);n38=filter2(q, n38);n38=filter2(q, n38);
n39=filter2(q, n39);n39=filter2(q, n39);n39=filter2(q, n39);
n40=filter2(q, n40);n40=filter2(q, n40);n40=filter2(q, n40);
n41=filter2(q, n41);n41=filter2(q, n41);n41=filter2(q, n41);	% Apply Filter
n42=filter2(q, n42);n42=filter2(q, n42);n42=filter2(q, n42);
n43=filter2(q, n43);n43=filter2(q, n43);n43=filter2(q, n43);
n44=filter2(q, n44);n44=filter2(q, n44);n44=filter2(q, n44);
n45=filter2(q, n45);n45=filter2(q, n45);n45=filter2(q, n45);
n46=filter2(q, n46);n46=filter2(q, n46);n46=filter2(q, n46);
n47=filter2(q, n47);n47=filter2(q, n47);n47=filter2(q, n47);
n48=filter2(q, n48);n48=filter2(q, n48);n48=filter2(q, n48);
n49=filter2(q, n49);n49=filter2(q, n49);n49=filter2(q, n49);
n50=filter2(q, n50);n50=filter2(q, n50);n50=filter2(q, n50);
n51=filter2(q, n51);n51=filter2(q, n51);n51=filter2(q, n51);	% Apply Filter
n52=filter2(q, n52);n52=filter2(q, n52);n52=filter2(q, n52);
n53=filter2(q, n53);n53=filter2(q, n53);n53=filter2(q, n53);
n54=filter2(q, n54);n54=filter2(q, n54);n54=filter2(q, n54);
n55=filter2(q, n55);n55=filter2(q, n55);n55=filter2(q, n55);
n56=filter2(q, n56);n56=filter2(q, n56);n56=filter2(q, n56);
n57=filter2(q, n57);n57=filter2(q, n57);n57=filter2(q, n57);
n58=filter2(q, n58);n58=filter2(q, n58);n58=filter2(q, n58);
n59=filter2(q, n59);n59=filter2(q, n59);n59=filter2(q, n59);
n60=filter2(q, n60);n60=filter2(q, n60);n60=filter2(q, n60);
n61=filter2(q, n61);n61=filter2(q, n61);n61=filter2(q, n61);	% Apply Filter
n62=filter2(q, n62);n62=filter2(q, n62);n62=filter2(q, n62);
n63=filter2(q, n63);n63=filter2(q, n63);n63=filter2(q, n63);
n64=filter2(q, n64);n64=filter2(q, n64);n64=filter2(q, n64);
n65=filter2(q, n65);n65=filter2(q, n65);n65=filter2(q, n65);
n66=filter2(q, n66);n66=filter2(q, n66);n66=filter2(q, n66);
n67=filter2(q, n67);n67=filter2(q, n67);n67=filter2(q, n67);
n68=filter2(q, n68);n68=filter2(q, n68);n68=filter2(q, n68);
n69=filter2(q, n69);n69=filter2(q, n69);n69=filter2(q, n69);
n70=filter2(q, n70);n70=filter2(q, n70);n70=filter2(q, n70);
n71=filter2(q, n71);n71=filter2(q, n71);n71=filter2(q, n71);	% Apply Filter
n72=filter2(q, n72);n72=filter2(q, n72);n72=filter2(q, n72);
n73=filter2(q, n73);n73=filter2(q, n73);n73=filter2(q, n73);
n74=filter2(q, n74);n74=filter2(q, n74);n74=filter2(q, n74);
n75=filter2(q, n75);n75=filter2(q, n75);n75=filter2(q, n75);
n76=filter2(q, n76);n76=filter2(q, n76);n76=filter2(q, n76);
n77=filter2(q, n77);n77=filter2(q, n77);n77=filter2(q, n77);
n78=filter2(q, n78);n78=filter2(q, n78);n78=filter2(q, n78);
n79=filter2(q, n79);n79=filter2(q, n79);n79=filter2(q, n79);
n80=filter2(q, n80);n80=filter2(q, n80);n80=filter2(q, n80);
n81=filter2(q, n81);n81=filter2(q, n81);n81=filter2(q, n81);	% Apply Filter
n82=filter2(q, n82);n82=filter2(q, n82);n82=filter2(q, n82);
n83=filter2(q, n83);n83=filter2(q, n83);n83=filter2(q, n83);
n84=filter2(q, n84);n84=filter2(q, n84);n84=filter2(q, n84);
n85=filter2(q, n85);n85=filter2(q, n85);n85=filter2(q, n85);
n86=filter2(q, n86);n86=filter2(q, n86);n86=filter2(q, n86);
n87=filter2(q, n87);n87=filter2(q, n87);n87=filter2(q, n87);
n88=filter2(q, n88);n88=filter2(q, n88);n88=filter2(q, n88);
n89=filter2(q, n89);n89=filter2(q, n89);n89=filter2(q, n89);
n90=filter2(q, n90);n90=filter2(q, n90);n90=filter2(q, n90);
n91=filter2(q, n91);n91=filter2(q, n91);n91=filter2(q, n91);	% Apply Filter
n92=filter2(q, n92);n92=filter2(q, n92);n92=filter2(q, n92);
n93=filter2(q, n93);n93=filter2(q, n93);n93=filter2(q, n93);
n94=filter2(q, n94);n94=filter2(q, n94);n94=filter2(q, n94);
n95=filter2(q, n95);n95=filter2(q, n95);n95=filter2(q, n95);
n96=filter2(q, n96);n96=filter2(q, n96);n96=filter2(q, n96);
n97=filter2(q, n97);n97=filter2(q, n97);n97=filter2(q, n97);
n98=filter2(q, n98);n98=filter2(q, n98);n98=filter2(q, n98);
n99=filter2(q, n99);n99=filter2(q, n99);n99=filter2(q, n99);
n100=filter2(q, n100);n100=filter2(q, n100);n100=filter2(q, n100);

% ms1 = mspc(n1,k);								% Normalized Marginal Spectra
% ms2 = mspc(n2,k);
% ms3 = mspc(n3,k);
% ms4 = mspc(n4,k);
% ms5 = mspc(n5,k);
% ms6 = mspc(n6,k);
% ms7 = mspc(n7,k);
% ms8 = mspc(n8,k);
% ms9 = mspc(n9,k);
% ms10 = mspc(n10,k);

N=zeros(128,170, 100);			% Make & Assemble Data Volume

N(:, :, 1)=n1;
N(:, :, 2)=n2;
N(:, :, 3)=n3;
N(:, :, 4)=n4;
N(:, :, 5)=n5;
N(:, :, 6)=n6;
N(:, :, 7)=n7;
N(:, :, 8)=n8;
N(:, :, 9)=n9;
N(:, :, 10)=n10;
N(:, :, 11)=n11;
N(:, :, 12)=n12;
N(:, :, 13)=n13;
N(:, :, 14)=n14;
N(:, :, 15)=n15;
N(:, :, 16)=n16;
N(:, :, 17)=n17;
N(:, :, 18)=n18;
N(:, :, 19)=n19;
N(:, :, 20)=n20;
N(:, :, 21)=n21;
N(:, :, 22)=n22;
N(:, :, 23)=n23;
N(:, :, 24)=n24;
N(:, :, 25)=n25;
N(:, :, 26)=n26;
N(:, :, 27)=n27;
N(:, :, 28)=n28;
N(:, :, 29)=n29;
N(:, :, 30)=n30;
N(:, :, 31)=n31;
N(:, :, 32)=n32;
N(:, :, 33)=n33;
N(:, :, 34)=n34;
N(:, :, 35)=n35;
N(:, :, 36)=n36;
N(:, :, 37)=n37;
N(:, :, 38)=n38;
N(:, :, 39)=n39;
N(:, :, 40)=n40;
N(:, :, 41)=n41;
N(:, :, 42)=n42;
N(:, :, 43)=n43;
N(:, :, 44)=n44;
N(:, :, 45)=n45;
N(:, :, 46)=n46;
N(:, :, 47)=n47;
N(:, :, 48)=n48;
N(:, :, 49)=n49;
N(:, :, 50)=n50;
N(:, :, 51)=n51;
N(:, :, 52)=n52;
N(:, :, 53)=n53;
N(:, :, 54)=n54;
N(:, :, 55)=n55;
N(:, :, 56)=n56;
N(:, :, 57)=n57;
N(:, :, 58)=n58;
N(:, :, 59)=n59;
N(:, :, 60)=n60;
N(:, :, 61)=n61;
N(:, :, 62)=n62;
N(:, :, 63)=n63;
N(:, :, 64)=n64;
N(:, :, 65)=n65;
N(:, :, 66)=n66;
N(:, :, 67)=n67;
N(:, :, 68)=n68;
N(:, :, 69)=n69;
N(:, :, 70)=n70;
N(:, :, 71)=n71;
N(:, :, 72)=n72;
N(:, :, 73)=n73;
N(:, :, 74)=n74;
N(:, :, 75)=n75;
N(:, :, 76)=n76;
N(:, :, 77)=n77;
N(:, :, 78)=n78;
N(:, :, 79)=n79;
N(:, :, 80)=n80;
N(:, :, 81)=n81;
N(:, :, 82)=n82;
N(:, :, 83)=n83;
N(:, :, 84)=n84;
N(:, :, 85)=n85;
N(:, :, 86)=n86;
N(:, :, 87)=n87;
N(:, :, 88)=n88;
N(:, :, 89)=n89;
N(:, :, 90)=n90;
N(:, :, 91)=n91;
N(:, :, 92)=n92;
N(:, :, 93)=n93;
N(:, :, 94)=n94;
N(:, :, 95)=n95;
N(:, :, 96)=n96;
N(:, :, 97)=n97;
N(:, :, 98)=n98;
N(:, :, 99)=n99;
N(:, :, 100)=n100;

% --- Averaging Section --- %

N(:, :, 1) = (N(:, :, 1) + N(:, :, 2)) ./ 2.;
N(:, :, 100) = (N(:, :, 99) + N(:, :, 100)) ./ 2.;
N(:, :, 2) = (4. .* N(:, :, 2) + N(:, :, 1) + N(:, :, 3)) ./ 6.;
N(:, :, 99) = (4. .* N(:, :, 99) + N(:, :, 98) + N(:, :, 100)) ./ 6.;

for i = 3:98;
	N(:, :,i) = (6. .* N(:, :, i) + 2. .* (N(:, :, i-1) +  N(:, :, i+1)) ...
	             +  N(:, :, i-2) +  N(:, :, i+2)) ./ 12.;
end;

save Nvar N x k;

tt=linspace(15, 25, 100)';


% --- Change this Value -----> x.xxx <---------- Change this Value -- %
p = patch(isosurface(x,k,tt,N, 0.01));
% ------------------------------------------------------------------- %
isonormals(x,k, tt, N, p)
set(p, 'FaceColor', 'cyan', 'EdgeColor', 'none');
daspect([.5 .2 0.25])
campos([114 33.4 92]);
lighting phong;
grid on;
camlight left

% --- Change the Title ---> Case   Frames  Line No. Contour Level --- %
title( 'X Distance: 0-26 cm; Y Distance: 15-25 cm; Wavenumber: 0-4 1/cm')
toptitle('3D Hilbert Spectrum: NEH-07 F130 L193 to 492 C.01')

% --- Bigd3Hilbert.m Ends Normally --- %