Matlab code to play with the Pollutant Standards Index (PSI)