Ti = (90:1:120)';
P = [ 0 0 0 0 0 0 0 0 0 20 40 60 80 100 100 100 100 100 100 100]';

Tu_90  = [139 136 133 130 127 124 121 118 115 112 109 106 103 100 97 94 91 88 85 82]'-10;
Tu_100 = [139 136 133 130 127 124 121 118 115 112 109 106 103 100 97 94 91 88 85 82]';
Tu_105 = [139 136 133 130 127 124 121 118 115 112 109 106 103 100 97 94 91 88 85 82]'+5;

P_TU_90  = interp1 (Tu_90, P,Ti);
P_TU_100 = interp1 (Tu_100,P,Ti);
P_TU_105 = interp1 (Tu_105,P,Ti);


deltaT = [-10:1:5]';

P_deltaT = interp1([-10 0 5], [P_TU_90'; P_TU_100';P_TU_105'], deltaT)';