function createFit(a,b)

f_ = clf;

figure(f_);

set(f_,'Units','Pixels','Position',[473 113 688 485]);

legh_ = [];

legt_ = {};

xlim_ = [Inf -Inf];

ax_ = axes;

set(ax_,'Units','normalized','OuterPosition',[0 0 1 1]);

set(ax_,'Box','on');

axes(ax_);

hold on;

a = a(:);

b = b(:);

h_ = line(a,b,'Parent',ax_,'Color',[0.333333 0 0.666667],...

'LineStyle','none', 'LineWidth',1,...

'Marker','o', 'MarkerSize',6);

xlim_(1) = min(xlim_(1),min(a));

xlim_(2) = max(xlim_(2),max(a));

legh_(end+1) = h_;

legt_{end+1} = 'b vs. a';

if all(isfinite(xlim_))

xlim_ = xlim_ + [-1 1] * 0.01 * diff(xlim_);

set(ax_,'XLim',xlim_)

else

set(ax_, 'XLim',[54.369999999999997, 118.63]);

end

if length(a)~=8

error( 'GenerateMFile:IncompatibleExclusionRule',...

'Exclusion rule ''%s'' is incompatible with ''%s''.',...

'first5b1', 'a' );

end

ex_ = false(length(a),1);

ex_([]) = 1;

ex_ = ex_ | (a >= 105);

ok_ = isfinite(a) & isfinite(b);

if ~all( ok_ )

warning( 'GenerateMFile:IgnoringNansAndInfs',...

'Ignoring NaNs and Infs in data.' );

end

st_ = [0 0 0 0.12822827157509359 ];

ft_ = fittype('fourier1');

if sum(~ex_(ok_))<2

error( 'GenerateMFile:NotEnoughDataAfterExclusionRule',...

'Not enough data left to fit ''%s'' after applying exclusion rule ''%s''.',...

'fit 1', 'first5b1' );

else

cf_ = fit(a(ok_),b(ok_),ft_,'Startpoint',st_,'Exclude',ex_(ok_));

end

h_ = plot(cf_,'predobs',0.95);

set(h_(1),'Color',[1 0 0],...

'LineStyle','-', 'LineWidth',2,...

'Marker','none', 'MarkerSize',6);

legend off;

legh_(end+1) = h_(1);

legt_{end+1} = 'fit 1';

if length(h_)>1

set(h_(2:end),'Color',[1 0 0],...

'LineStyle',':', 'LineWidth',1,'Marker','none');

legh_(end+1) = h_(2);

legt_{end+1} = 'Pred bnds (fit 1)';

end

hold off;

leginfo_ = {'Orientation', 'vertical', 'Location', 'SouthEast'};

h_ = legend(ax_,legh_,legt_,leginfo_{:});

set(h_,'Interpreter','none');

xlabel(ax_,'');

ylabel(ax_,'');