x0 没有给出！！假x0=1953;function [xm,t]=nihe(p)
q=zeros(3,1);
xm=p(1);r=p(2);
x=[6.0 7.2 10.3];
t=[1953 1964 1982];
x0=x(1);
q(1)=xm/(1+(xm/x0-1)*exp(-r*t(1)))-x(1);
q(2)=xm/(1+(xm/x0-1)*exp(-r*t(2)))-x(2);
q(3)=xm/(1+(xm/x0-1)*exp(-r*t(3)))-x(3);
%命令窗口
[xm,t]=fsolve('nihe',xx)
