model
sets
fac/1..21/:li,a;
endsets
data
n=108;
x0=116.3914;
t0=12.68333333 12.73333333 12.78333333 12.83333 12.88333 12.93333 12.98333 13.03333333 13.08333 13.13333 13.18333 13.23333 13.28333 13.33333 13.38333 13.43333 13.48333 13.53333 13.58333 13.63333 13.68333;
li=1.247256205 1.22279459 1.198921486 1.175428964 1.152439573 1.12991747 1.10783548 1.086254206 1.065081072 1.044446265 1.024264126 1.004640314 0.985490908 0.966790494 0.948584735 0.930927881 0.91375175 0.897109051 0.880973762 0.865492259 0.850504468;
enddata
min=@sum((fac(i):li(i+1)/li(i)-@tan(a(i))/@tan(a(i+1)))^2);
@tan(r(i))=@asin(@sin(z)@sin(b)+@cos(z)@cos(b)@cos(y));
b=23.45*@sin((n-80.25)*(1-(n/9500)));
y=15*(t-12);
t0=t+1/15*(x0-x);
@bnd(-90,z,90);
t>=12;
end
https://blog.csdn.net/sunshine_r_007/article/details/77971910