function [rhoX rhoY] = SensVSpower(power, alphaX, alphaY)
%Linear regression of sensitivity on power.
%First, enter or import power and alpha vectors.
rhoX = ((2.0116e13.*alphaX).^0.5).*1e-6;
rhoY = ((2.0116e13.*alphaY).^0.5).*1e-6;

px = polyfit(power, rhoX, 1);
plot(power, rhoX, 'ob', 'MarkerSize',8)
hold on

xfit = 0:0.01:1;
yfitX = polyval(px, xfit);
plot(xfit,yfitX,'-b', 'LineWidth',2)

py = polyfit(power, rhoY, 1);
plot(power, rhoY, 'sr', 'MarkerSize',8)
hold on


yfitY = polyval(py, xfit);
plot(xfit,yfitY,'-r', 'LineWidth',2)

xlabel('Laser Power (proportion max)')
ylabel('sensitivity (V/micrometer')
title('Relationship between sensitivity (rho) and laser power')
legend('rho-X', ['rhoX = ', num2str(px(2), 3), '+ ', num2str(px(1), 3), 'power'],...
    'rho-Y', ['rhoY = ', num2str(py(2), 3), '+ ', num2str(py(1), 3), 'power'])
hold off
end