function [cA, cD1, cD2]=mwave2(X,LoD,HiD)
[C,L] = wavedec(X,2,LoD,HiD);
cA = appcoef(C,L,LoD,HiD);
cD1 = detcoef(C,L,1);
cD2 = detcoef(C,L,2);

