Here is an old XLM UDF which accepts polynomial coefficients and temperature
from a worksheet and returns the predicted materials property based on a
simplistic model.  I wrote this in 1993 or 94.  I won't try to attach the next
project I worked on, which was an elaborate XLM macro that input about 18
fitted coefficients and 4 more environmental conditions, and iterated through
SOLVER to predict metal fatigue behavior in turbine engine parts.

MDD.POL5PH.AVG
Properties 01 (Young's Modulus), 02 (Shear Modulus), 03 (Poisson's Ratio), 04
(CTE), 05 (Specific Heat), 06 (Thermal Cond.), 09 (), 16 (Mono. Hardening
Exp.), 29 (Cycl. Hardening Exp.), 30 (0.2% Cycl. YS), 33 ()
=ARGUMENT("Azero")
=ARGUMENT("Aone")
=ARGUMENT("Atwo")
=ARGUMENT("Athree")
=ARGUMENT("Afour")
=ARGUMENT("Afive")
=ARGUMENT("Delta")
=ARGUMENT("temp")
=SET.NAME("X",temp/Delta)
=SET.NAME("prop",Azero+Aone*X+Atwo*X^2+Athree*X^3+Afour*X^4+Afive*X^5)
=RETURN(prop)

The VBA equivalent to this XLM macro, which I translated in 1997, is shown
below.

Function pol5ph(wksht, colnumb, Temp)
    a_0 = wksht.Cells(7, colnum)
    a_1 = wksht.Cells(8, colnum)
    a_2 = wksht.Cells(9, colnum)
    a_3 = wksht.Cells(10, colnum)
    a_4 = wksht.Cells(11, colnum)
    a_5 = wksht.Cells(12, colnum)
    delta = wksht.Cells(13, colnum)
    XX = Temp / delta
    pol5ph = a_0 + a_1 * XX + a_2 * XX ^ 2 + a_3 * XX ^ 3 _
    + a_4 * XX ^ 4 + a_5 * XX ^ 5
End Function