clear; close all;

t = [-40 -30 -20 -10 0 10 20 25 30 40 50 60 70 80];
vo = [4757 4595 4361 4048 3657 3210 2736 2500 2270 1842 1469 1158 911 715];
vo =vo/1000;
plot(vo,t);

P = polyfit(vo,t,3)

ty = polyval(P,vo);

plot(vo,t,'.',vo,ty);