added code to print out integral of f psisq