class function StdDev(a: TFloatMatrix; meanVector: TFloatArray): TFloatArray;

class procedure Matrix.Normalize(var a: TFloatMatrix;out meanVector: TFloatArray); var stdD :TFloatArray; n,m,mm,nn: integer; begin meanVector:=Mean(a); stdD:=StdDev(a,meanVector); n:=Length(a); m:=Length(a[0]); for nn:=0 to n-1 do for mm:=0 to m-1 do a[nn,mm]:=(a[nn,mm]-meanVector[mm])/stdD[mm]; end;