%inmatning av värden
x = [10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280 290 300 310 320 330 340 350 360 370 380 390 400];

y1 = [295 279 314 286 258 246 240 253 220 205 196 193 189 188 160 170 142 140 133 135 149 120 118 101 123 120 101 109 89 95 112 69 72 58 78 71 59 55 53 39];
y2 = [251 97 253 103 224 100 218 101 223 108 214 92 236 98 199 80 177 67 166 66 153 63 155 71 151 62 161 68 139 47 121 41 135 52 111 48 125 43 97 49];

%subtrahera bakgrundsstrålning
y1 = y1 - (1/3);
y2 = y2 - (1/3);

%beräkna medelvärde
y = (y1+y2)/2;

%logaritmera
y = log(y);

%anpassa en linje
p = polyfit(x,y,1);
ypass = polyval(p,x);

%rita diagram
figure(1);
plot(x,y,'+',x,ypass,'-');  %plotta
xlabel('t/s'); %namnge x-axeln
ylabel('ln(A)');  %namnge y-axlen
title('Diagram 2');

%beräkna k samt skriv ut det
k = p(1);
display(k);