Описание синтеза бутадиенового каучука дифференциальными уравнениями

[Mirabella]

Добрый день! 
Необходимо математически описать синтез бутадиенового каучука, т.е. массовый баланс и тепловой для лабораторного или промышленного реактора. 
Уравнения получились такие:
изменение концентрации мономера  

d[btd]/dt = -kp*[NBL]^0.24*[btd] 

где [btd] - концентрация мономера, kp - константа скорости роста, [NBL] - концентрация катализатора [n-BuLi] (возведена в степень 0,24)

тепловой баланс
m*Cp*dT/dt = (-deltaH)*Rp*V - UA(T-T0) , где
m - масса реакционной смеси,
Сp - теплоемкость реакционной смеси
deltaH - энтальпия
Rp - скорость роста
V - объем
U - телпопроводность стенки
A - площадь
Т - температура реакционной смеси
T0 - температура окружающей среды.
Т.е. первое слагаемое, это то тепло, которое выделилось при реакции, а второе, то что ушло в окружающую среду.
Не могу найти ошибку в уравнениях, но графики изменения конверсии мономера и температуры со временем отчаянно не хотят строится.
За любую информацию (где ошибка, где посмотреть уравнения, где посмотреть числовые значения констант и пр.) буду очень благодарна. Модель нужная самая грубая и простая для оценки изменения концентрации и температуры реакции.