orig = [   0, 50, 31;
          59, 61, 39;
         120, 79, 57;
         181, 88, 69;
         243, 82, 63;
         304, 62, 40; ];

y_h = orig(:,2);
y_l = orig(:,3);
% Design matrix  y = B_0 + B_1 cos((2 * pi * t) / 365) +
%                          B_2 sin((2 * pi * t) /365)

for i=1:size(orig,1)
    t = orig(i,1);
    X(i,1) = 1;
    X(i,2) = cos((2 * pi * t) / 365);
    X(i,3) = sin(( 2 * pi * t) / 365);
end

X