%tar in alla variabler
C=input('Lägg in en vektor för kapacitansen (F): ');
d=input('Lägg in en vektor för avstånden (m): ');
A=input('Lägg in ett värde på din area på plattorna (m^2): ');
%för om till ett linjärt samband och plottar

x=d.^-1;
plot(x,C,'b*')
hold on
%skapar och plottar en linjär anpassning

k=polyfit(x,C,1);
z=linspace(x(1),x(end));
y=polyval(k,z);
plot(z,y,'r')
%tar fram koefficienterna samt tar fram vilka konfidensintervall vi har

[xData, yData] = prepareCurveData(x,C);
linjfit=fit(xData,yData,'poly1')
coeff=coeffvalues(linjfit)
conf=confint(linjfit,0.68)
u_k=coeff(1)-conf(1);

%fixar till grafen
xlabel('inversen av avståndet (1/m)')
ylabel('kapacitans (F)')
title('linjär graf över plattornas avstånd och kapacitans')
legend('Datapunkter','Linjär anpassning')
%skriver ut permittiviteten
E_r=k(1)/(A*8.854*10^-12)
u_E_r=sqrt((u_k/(8.854*10^-12*A))^2-(coeff(1)*0.0005/(8.854*10^-12*A^2))^2)