%(3) Differential expression b * δ t = a * (t+x0) * δ s

%(4) model: ds/dt=k/(t+x0)

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

Clear; Chromaticity control center

s = dsolve(' DS = k/(t+x0)')； General solution of differential equation of% solution model

equate = subs(S，' t '，[0 2 4]-[0 2 3]； % into S(0)=0, S(2)=2, S(4)=3.

[C 1，k，x0]=solve(equat( 1)，equat(2)，equat(3))； % solution equation

X0=simplify(x0)% Simplify

Running results:

x0 =

[5 (1/2)- 1]% is very different.

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