Some Matlab programs