Update psdestimative.m