R= 8.3143;
F=96487

tempDegC=25:-1:4
S= ((R * (tempDegC - 273.15)) / F) * log(10);

plot(tempDegC, S)