'peaklock': determines peacklocking errors from velocity histograms. ------------------------------------------------------- first smooth the input histogram 'histu' in such a way that the integral over n-n+1 is preserved, then deduce the peaklocking 'error' function of the pixcel displacement 'x'. [histinter,x,error]=peaklock(nbb,minim,maxim,histu) OUTPUT: histinter: smoothed interpolated histogram x: vector of displacement values. error: vector of estimated errors corresponding to x INPUT: histu=vector representing the values of histogram of measured velocity ; minim, maxim: extremal values of the measured velocity (absica for histu) nbb: number of bins inside each integer interval for the histograms SUBROUTINES INCLUDED: spline4.m% spline interpolation at 4th order splinhist.m: give spline coeff cc for a smooth histo (call spline4) histsmooth.m(x,cc): calculate the smooth histo for any value x histder.m(x,cc): calculate the derivative of the smooth histo