[PHI,PSI,XVAL] = wavefun('Db2',16);
 integralPHI=sum(PHI)/length(PHI)*max(XVAL)
 velkostPHI=sum(sqrt(PHI.*PHI))/length(PSI)*max(XVAL)
 energiaPHI=sum(PHI.*PHI)/length(PSI)*max(XVAL)
 integralPHI=sum(PHI)/length(PHI)*max(XVAL)
 integralPSI=sum(PSI)/length(PSI)*max(XVAL)
 velkostPSI=sum(sqrt(PSI.*PSI))/length(PSI)*max(XVAL)
 energiaPSI=sum(PSI.*PSI)/length(PSI)*max(XVAL)