Generate a fit, round to appropriate sig figs, print/savetxt with correct sig figs, plot rounded function, possibly retain original fit in save?