[freq,freqS,amp,ampS,quali,qualiS,curr,currS] = WriteData('Elektro22.data');

figure(1)
plot(freq/1000,[amp,amp+ampS,amp-ampS])
xlabel('frequency [kHz]'); ylabel('amplitude')

figure(2)
plot(freq/1000,[quali,quali+qualiS,quali-qualiS])
xlabel('frequency [kHz]'); ylabel('Q-factor')