
[x, y] = m2samplenoisydata();

plot(x, y, 'b.')
hold on

% Here goes the code similar to the one from HW07
%a1 = m2svdfit(x,y)

% Next line only as a demo, not accepted on the test
a1 = [0.19 0.21 0.22];

yy1 = a1(1) ./ x + a1(2)*cos(x) + a1(3)*exp(-x.^2/2);
plot(x, yy1, 'r--', 'LineWidth', 2)

hold off
