using scipy.optimize.leastsq
B_V fixed