IntParameter a=[1,50];
Parameters Cus,Ku,Kr,Kf,So;
Variables x,y;
StartProgram [Pascal];
Procedure MainModel;
var i, j, p: integer;
    temd, temd1, temd2, temd3: double;
Begin
    p := Round(a);
    for j := 0 to DataLength - 1 do begin
        temd1 := Cus*Ku*(1+Kr);
        temd2 := 0;
        for i := 1 to p + 1 do
            temd2 := temd2 + power(x[j],i);
        temd3 := 0;
        for i := 1 to p + 1 do
            temd3 := temd3 + power(x[j],i-p);
        y[j] := temd1*temd2/(Ku*(1+Kr)*temd2+temd3)+So*Kf*x[j]/(1+Kf*x[j]);
    end;
End;
EndProgram;
Data;
0.56509 0
2.30222 0.444696164
4.03935 0.841178994
5.96484 1.135857962
8.07869 1.532340792
10.98786 1.923468197
15.21557 2.368164361
20.21766 2.909308963
32.16827 4.195191377
39.68188 5.127456128
49.89535 6.365125209
60.10883 8.138568953
70.32231 10.50674057
80.13813 13.46427738
90.16325 16.22356915
94.60025 17.01117213

3楼: Originally posted by dingd at 2014-06-10 11:53:40
1stOpt可以很容易实现的。