wrote code for F1 and F2 - next step is to check it is correct by making sure F1=F2 at the analytic solution for simga = x^2/2