function [nu, S] = innovation(xpred, Ppred, z, H, R)

     nu = z - H*xpred;                   %% innnvoation
     S = R + H*Ppred*H';                 %% innovation convariance



