之前从来没用过1stopt，下面是自己写的程序，用盗版程序根本运行不了。要拟合的参数是Cr1,Cr2,Cu1,Cu2,mu，自变量是t，因变量是xi1,xi2,lambda1,lambda2,s，现在实验数据是lambda1和s的。Title "lashennihe";Parameters Cr1,Cr2,Cu1,Cu2,mu;Variable t,xi1,xi2,lambda1,lambda2,s;ODEFunction s=Cr1*(2*lambda1-2*lambda1^(-3)*lambda2^(-2))+Cr2*(-2*lambda1^(-3)+2*lambda1*lambda2^2)+(Cu1-Cr1)*(2*lambda1*xi1^(-2)-2*lambda1^(-3)*lambda2^(-2)*xi1^2*xi2^2)+(Cu2-Cr2)*(-2*lambda1^(-3)*xi1^2+2*lambda1*lambda2^2*xi1^(-2)*xi2^(-2)); Cr1*(2*lambda2-2*lambda2^(-3)*lambda1^(-2))+Cr2*(-2*lambda2^(-3)+2*lambda2*lambda1^2)+(Cu1-Cr1)*(2*lambda2*xi2^(-2)-2*lambda2^(-3)*lambda1^(-2)*xi1^2*xi2^2)+(Cu2-Cr2)*(-2*lambda2^(-3)*xi2^2+2*lambda2*lambda1^2*xi1^(-2)*xi2^(-2))=0; lambda1=1.176471*t; xi1'=-1/mu*((Cu1-Cr1)*(-2*lambda1^2*xi1^(-3)+2*lambda1^(-2)*lambda2^(-2)*xi1*xi2^2))+(Cu2-Cr2)*(2*lambda1^(-2)*xi1-2*lambda1^2*lambda2^2*xi1^(-3)*xi2^(-2))); xi1'=-1/mu*((Cu1-Cr1)*(-2*lambda2^2*xi2^(-3)+2*lambda1^(-2)*lambda2^(-2)*xi2*xi1^2))+(Cu2-Cr2)*(2*lambda2^(-2)*xi2-2*lambda1^2*lambda2^2*xi2^(-3)*xi1^(-2)));Data;// s lambda1实验数据在excel表格里，有6638组，等级不够,发不上来啊

更正一下，代码有点乱，整理了一下，最后一个式子是xi2‘不是xi1’~谢谢！！
Title "lashennihe";
Parameters Cr1,Cr2,Cu1,Cu2,mu;
Variable t,xi1,xi2,lambda1,lambda2,s;
ODEFunction s=Cr1*(2*lambda1-2*lambda1^(-3)*lambda2^(-2))+Cr2(-2*lambda1^(-3)+2*lambda1*lambda2^2)(Cu1-Cr1)*(2*lambda1*xi1^(-2)-2*lambda1^(-3)*lambda2^(-2)*xi1^2*xi2^2)+(Cu2-Cr2)*(-2*lambda1^(-3)*xi1^2+2*lambda1*lambda2^2*xi1^(-2)*xi2^(-2));
Cr1*(2*lambda2-2*lambda2^(-3)*lambda1^(-2))+Cr2*(-2*lambda2^(-3)+2*lambda2*lambda1^2)+(Cu1-Cr1)*(2*lambda2*xi2^(-2)-2*lambda2^(-3)*lambda1^(-2)*xi1^2*xi2^2)+(Cu2-Cr2)*(-2*lambda2^(-3)*xi2^2+2*lambda2*lambda1^2*xi1^(-2)*xi2^(-2))=0;
lambda1=1.176471*t;
xi1'=-1/mu*((Cu1-Cr1)*(-2*lambda1^2*xi1^(-3)+2*lambda1^(-2)*lambda2^(-2)*xi1*xi2^2))+(Cu2-Cr2)*(2*lambda1^(-2)*xi1-2*lambda1^2*lambda2^2*xi1^(-3)*xi2^(-2)));
xi2'=-1/mu*((Cu1-Cr1)*(-2*lambda2^2*xi2^(-3)+2*lambda1^(-2)*lambda2^(-2)*xi2*xi1^2))+(Cu2-Cr2)*(2*lambda2^(-2)*xi2-2*lambda1^2*lambda2^2*xi2^(-3)*xi1^(-2)));
Data;
// s lambda1