用matlab對函數求導及其極值 ，可以先用diff（）函數求出y'，然後令y'=0，用vpasolve（）解出x值。實現代碼如下：
syms x
z=exp(4.1399-0.0235*x);
y=0.000157*z*(1-z)*(1+z)^3/((1+z)^4+0.00004476*z^2)^(3/2);
dy=diff(y,1)
x=vpasolve(dy==0)
運行結果
x=120.125
用matlab對函數求導及其極值 ，可以先用diff（）函數求出y'，然後令y'=0，用vpasolve（）解出x值。實現代碼如下：
syms x
z=exp(4.1399-0.0235*x);
y=0.000157*z*(1-z)*(1+z)^3/((1+z)^4+0.00004476*z^2)^(3/2);
dy=diff(y,1)
x=vpasolve(dy==0)
運行結果
x=120.125