x = [1, 2.07; 1, 2.37; 1, 2.54];
y = [0.779; 0.916; 0.905];
%x = load('ex2x.dat');
%y = load('ex2y.dat');
m = length(y);
%x = [ones(m, 1), x];
theta = [0.0;0.0];
alpha = 0.0002;
for k=1:2000,
   acc = 0.0;
   for j = 1:m,
      acc = ((theta'*x(j,:)'-y(j,1))*x(j,:))'+acc;
   end
   theta = theta - (alpha / m)*acc;
end