
  close all % close all plots
  
  plot(Vs,VR2,'o') % Vs is the x-input, VR2 is the y-output, 'o' plots bubbles
  hold on % keep the current plot
  
  coef = polyfit(Vs,VR2,1); % find the slope and intercept of the best fit
  m = coef(1) % assign the slope
  b = coef(2) % assign the intercept
  
  plot(Vs,m*Vs+b,'--') % plot the best fit line as a dashed line