%用函數(shù)ode45()來(lái)解微分方程，并畫(huà)圖計(jì)算結(jié)果(注：ode45()采用龍格庫(kù)塔算法求解)
%example.m
%初始化x1，x2,...,x22。x1即xD,x22即xB
[t,x]=ode45(@ivpodefun,[0 500],[0.55;0.545;0.54;0.535;0.53;0.525;0.52;0.515;0.51;0.505;0.5;0.495;0.49;0.485;0.48;0.475;0.47;0.465;0.46;0.455;0.45;0.445]);
%畫(huà)圖
plot(t,x(:,1),'-',t,x(:,22),'--')
xlabel('t/min');
ylabel('x');
legend('xD','xB');