main.m

来自「for the parameter estimation and plot th」· M 代码 · 共 8 行

clear
clc

[time,sol]=ode45(@heatexco,[0:0.01:1],[0 0 0 0]);

sol
plot(time,sol)
legend('xgc','tgc','xgh','tgh');