clc; clear all; close all;
data=load('A04_VAN0_99-02.dat');
subplot(1,3,1);
errorbar(data(:,1), data(:,2),data(:,3));
pn=polyfit(data(:,1),data(:,2), 1)
hold on;

pnx=@(f)(pn(1)*f+pn(2))
pnxx=pnx(data(:,1));
plot(data(:,1),pnxx)
title('norr');
ylabel('Change in mm');
xlabel('year');

subplot(1,3,2);
errorbar(data(:,1), data(:,4),data(:,5));
po=polyfit(data(:,1),data(:,4), 1)
hold on;

pox=@(f)(pn(1)*f+pn(2))
poxx=pox(data(:,1));
plot(data(:,1),poxx)

title('öst');
ylabel('Change in mm');
xlabel('year');

subplot(1,3,3);
errorbar(data(:,1), data(:,6),data(:,7));
ph=polyfit(data(:,1),data(:,6), 1)
hold on;

phx=@(f)(pn(1)*f+pn(2))
phxx=phx(data(:,1));
plot(data(:,1),phxx)

title('höjd');
ylabel('Change in mm');
xlabel('year');