rewrite find_sigma function in postscripts to a utility function