數學建模問題

%（3）微分表達 b*Δt=a*(t+x0)*ΔS

%（4）模型： ds/dt=k/(t+x0)

%k=a/b,S(0)=0，S(2)=2，S(4)=3

clear;clc

S=dsolve('DS=k/(t+x0)');%解模型的微分方程通解

equat=subs(S,'t',[0 2 4])-[0 2 3];%把S(0)=0，S(2)=2，S(4)=3代入

[C1,k,x0]=solve(equat(1),equat(2),equat(3));%解方程

x0=simplify(x0)%化簡

運行結果：

x0 =

[ 5^(1/2)-1]%不太壹樣啊

[ exp(1/2*(3+5^(1/2))^(1/3)*(5^(1/2)-1))-4]