Как получить промежуточное значение используя линейную интерполяцию?

short sArrayResistance [54] = {32554.0,  10000.0,  3582.96,  3428.57,  3298.32, 3175.51,  3059.52,  2949.81,  2831.46,  2733.58, 2640.57,  2539.79,  2456.08,  2365.13,  2278.85, 2206.90,  2128.44,  2053.73,  1982.51,  1914.53, 1849.58,  1787.47,  1728.00,  1671.02,  1616.37, 1563.91,  1511.06,  1459.13,  1412.74,  1368.06, 1325.00,  1283.48,  1243.42,  1204.74,  1167.37, 1131.25,  1096.31,  1058.35,  1025.74,  998.05, 967.31,   937.50,   908.58,   883.98,   856.62, 830.05,   807.42,   782.23,   760.76,   735.14, 716.44,   696.52,   677.05,   182.6};

short sArrayTemperature [54] = {0,       25,      52.81,   54.25,   55.44, 56.56,   57.69,   58.81,   59.94,   61.13, 62.25,   63.38,   64.50,   65.69,   66.81, 68.00,   69.12,   70.25,   71.37,   72.56, 73.69,   74.81,   76.00,   77.12,   78.31, 79.44,   80.56,   81.75,   82.87,   84.06, 85.25,   86.37,   87.56,   88.69,   89.87, 91.06,   92.25,   93.44,   94.56,   95.69, 96.87,   98.00,   99.19,   100.37,  101.50, 102.69,  103.87,  105.06,  106.19,  107.37, 108.50,  109.69,  110.87,  150};

// Если R строго внутри массива sArrayResistance, то как-то так
double getTemperature (double R) { 
    for (int i = 1; i < 54; i++)
        if (R>sArrayResistance[i])
            return sArrayTemperature[i-1] + (sArrayTemperature[i-1]-sArrayTemperature[i])/(sArrayResistance[i]-sArrayResistance[i-1])*(sArrayResistance[i-1]-R);
}