%Tidvattenståndet i Nordsjön
t=[0:2:10]';
H=[1 1.6 1.4 0.6 0.2 0.8]';
A=[ones(size(t)) sin(2*pi*t/12) cos(2*pi*t/12)];
c=A\H
tp=[0:0.1:10]';
Hp=c(1)+c(2)*sin(2*pi*tp/12)+c(3)*cos(2*pi*tp/12);
plot(t,H,'o',tp,Hp)
title('Tidvattenståndet i Nordsjön')
xlabel('t [h]'),ylabel('H [m]')
medel=c(1)
amplitud=sqrt(c(2)^2+c(3)^2)
