Question

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

Answer 1

这里有个例子,你看看吧

第一步,把你要积分的函数编写出来

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);