% vectorization procedure
robust = 1;
param = fitoptions('method','NonlinearLeastSquares');
param.Algorithm = 'Levenberg-Marquardt';
%param.DiffMaxChange = 0.5;
param.Display = 'Notify';
if robust == 1
    param.Robust = 'on';
end
get(param)
% Here get frames ;]
xdata=[-4;-3;-2;-1;0;1;2;3;4]
ydata=[16;9;4;1;0;1;4;9;16]
fresult = fit(xdata,ydata,'poly2',param)
for n=1:1000
fresult = fit(xdata,ydata,'poly2',param);
end
[p1, p2, p3] = fit(xdata,ydata,'poly2',param)