29楼: Originally posted by wurongjun at 2013-10-05 13:52:58
function fcqiujie
AB=textread('AB.txt');
N=size(AB,1);
A=AB(:,1);B=AB(:,2);
tx=zeros(N,4);
for i=1:N
    tx(i,=roots()';
end
x=tx+1;
for i=1:N
y(i,=x(i,-A(i)*(1-x(i,).^2;
end
for  ...

2楼: Originally posted by wurongjun at 2013-10-02 16:58:55
由第一式得:y=x-A*(1-x)^2
代入第二式得:1+2*(x-1)+(-2*A+1)*(x-1)^2+(-2*A+2*B*A)*(x-1)^3+(A^2-2*B*A^2)*(x-1)^4=0
上式是关于(x-1)的四次方程!可由四次方程求根公式求解!