这里有一个偏微分方程:C1*(dT/dt)-C2*C3^(1/T)=d2T/d2x T=T(t,x) 条件是:T(0,t)=20,dT(0,t)/dx=0

题目

这里有一个偏微分方程:C1*(dT/dt)-C2*C3^(1/T)=d2T/d2x T=T(t,x) 条件是:T(0,t)=20,dT(0,t)/dx=0

答案

这里有个例子,你看看吧
第一步,把你要积分的函数编写出来
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
第二步,进行积分
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);%设置步长
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);