Написать программу на С++

Нужен код на С++ для подбора коэффициентов Ug, D1 и D2 в уравнение U = U g ? D1(T ? Tg) + D2(T ? Tg)^2 с применением полинома Ньютона для расчёта температурной зависимости энергии активации стеклования для листового силикатного стекла, используя температуру стеклования температура стеклования(1578,0), логарифм вязкости(2,5) и ?0 – предэкспоненциальный множитель(-2,05), чтобы расчёты совпадали с расчётами по формуле U = 2.3RT[lg ?(T) - lg ?0]