DEFINE_SOURCE(heat-gen,c,t,dS,eqn)
{
real A=;
real B=;
real C=;
real D=;
real T=C_T(c,t);
real source=A*T+ B*pow{T,2)+ C*pow{T,3)+ D*pow{T,4);
dS[eqn]=0;//關於dS系源項對溫度的壹階偏導,來源於傅裏葉展開,可以設置該值,簡化處理默認為零
return source;
}