function vm_parabola=vm_curve_fitter(I_macro_vector,var_macro_vector)
%in reality, since there are only two parameters, it might be better to do a brute force search over the whole parameter space to get the best fit

% optionz=optimset('display','off');
% 
% lower_bounds=[1 1000];
% upper_bounds=[100 1];

[fitted_vals,resnorm,residual]=lsqcurvefit(@lsqfit_vm_parabola_generator,[10 3],I_macro_vector,var_macro_vector);%,lower_bounds,upper_bounds,optionz)

vm_parabola=lsqfit_vm_parabola_generator(fitted_vals,I_macro_vector);