tempC = linspace(-40, 1670, 181);

te_voltage = [-0.194 -0.150 -0.103 -0.053 0 0.055 .113 .173 .235 .299 .365 .433 .502 .573 .646 .720 .795 .872 .950 1.029 1.110 1.191 1.273 1.357 1.441 1.526 1.612 1.698 1.786 1.874 1.962 2.052 2.141 2.232 2.323 2.415 2.507 2.599 2.692 2.786 2.880 2.974 3.069 3.164 3.259 3.355 3.451 3.548 3.645 3.742 3.840 3.938 4.036 4.134 4.233 4.332 4.432 4.532 4.632 4.732 4.833 4.934 5.035 5.137 5.239 5.341 5.443 5.546 5.649 5.753 5.857 5.961 6.065 6.170 6.275 6.381 6.486 6.593 6.699 6.806 6.913 7.020 7.128 7.236 7.345 7.454 7.563 7.673 7.783 7.893 8.003 8.114 8.226 8.337 8.449 8.562 8.674 8.787 8.900 9.014 9.128 9.242 9.357 9.472 9.587 9.703 9.819 9.935 10.051 10.168 10.285 10.403 10.520 10.638 10.757 10.875 10.994 11.113 11.232 11.351 11.471 11.590 11.710 11.830 11.951 12.071 12.191 12.312 12.433 12.554 12.675 12.796 12.917 13.038 13.159 13.280 13.402 13.523 13.644 13.786 13.887 14.009 14.130 14.251 14.373 14.494 14.615 14.736 14.857 14.978 15.099 15.220 15.341 15.461 15.582 15.702 15.822 15.942 16.062 16.182 16.301 16.420 16.539 16.658 16.777 16.895 17.013 17.131 17.249 17.366 17.483 17.600 17.717 17.8032 17.947 18.061 18.174 18.285 18.395 18.503 18.609];

tempC_ = linspace(0,1670,177);
te_voltage_ = te_voltage(:,5:181);

type_s_gain = 268;

type_k_gain = 92.5;

type_s_voltage = te_voltage_ * type_s_gain;

f1 = fit(tempC_.', type_s_voltage.', 'poly1')
f2 = fit(tempC_.', type_s_voltage.', 'poly2')
f3 = fit(tempC_.', type_s_voltage.', 'smoothingspline')



plot(f2,tempC_.',type_s_voltage.')