function run_fitting()

ij=[1,1;
    1,2;
    1,3;
    2,2;
    2,3;
    3,3;
    4,4;
    5,5;
    6,6];
paramFitting = [];
for k=1:9
    ii=ij(k,1); jj=ij(k,2);
    disp([ii,jj])
    b_reg = fitting(ii,jj,5);
    paramFitting = [paramFitting,b_reg];
end
save FittingParam.mat
end