maincount.m

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

clear
clc

solinit=bvpinit([0:0.01:1],[0.5 1.8 0.9 2.5])
sol=bvp4c(@heatexcount,@heatexcountbc,solinit);

xint=[0:0.01:1];
yint=deval(sol,xint);

plot(xint,yint)
legend('xgc','tgc','xgh','tgh');