@Bill: Спасибо, что указали это! Вот мой исправленный код:
dSes[T_, p_, alpha_] = D[Ses[T, p, alpha], T];
Ces[T_, p_, alpha_] := T*dSes[T, p, alpha];
теперь работает дифференциация. Но когда я запускаю подгонку, я получаю следующую ошибку: длинный список чисел, скобок, компиляция '$ ... не является списком реальных чисел с размерами {28} в {p, alpha} = {0.5.1.764 }.