这个模型的代码这样写正确吗，其中x代表投入变量，一共有三种投入（资本 劳动 能源消耗），y代表gdp ，u代表非期望产出。一共有30个城市model:
sets:
capital/cap1..cap30/:x1;
labor/l1..l30/:x2;
energy/e1..e30/:x3;
gdp/g1..g30/:y;
carbon/car1..car30/:u;
para/p1..p30/:h;
endsets
data:
x1=
12766.51
6172.83
16116.01
6283.22
6928.55
13380.59
5906.82
7441.68
16689.37
28438.95
24749.09
8512.51
9589.43
6904.46
29257.69
14432.21
11319.03
9708.65
28934.31
5849.81
1729.20
6509.42
13425.16
3984.79
6268.66
6695.94
3622.88
1381.04
1593.81
5526.31
;
x2=
878.00
542.52
3568.97
1500.20
1041.10
2120.30
1238.90
1748.90
969.24
4578.75
3100.76
3669.70
1868.50
2276.70
5840.70
5662.00
3537.00
3801.48
5022.97
2703.00
379.55
1456.30
4702.00
2220.00
2461.30
1976.00
1391.40
291.04
299.60
744.49
;
x3=
5522
4085
19836
12750
9666
13611
5315
8050
8225
17167
12032
6506
6142
24162
4286
14625
10082
9709
17921
4869
822
4943
11816
5641
6024
5571
4368
1670
2536
5506
;
y=
4430.27
2413.68
6907.59
2667.59
2450.81
5512.54
2511.40
3835.69
5884.25
11994.38
8336.83
3706.61
4481.96
2755.17
11909.48
7350.10
4689.34
4651.41
14102.20
2829.86
621.83
2373.55
5082.11
1432.10
2487.66
2631.13
1323.32
355.76
405.45
1825.27
;
u=
13679.7
15050.06
71583.2
74967.08
41111.84
58017.91
22037.87
30884.87
26567.86
60185.48
37206.68
25613.07
16627.93
14769.28
89919.95
55331.19
30451.44
27881.43
45963.36
12466.28
1481.58
13011.79
27350.23
22164.29
22209.25
22098.92
15519.58
3179.696
9625.67
17234.33
;
enddata
max=r;
@sum(para:h*x1)<=12766.51;
@sum(para:h*x2)<=878;
@sum(para:h*x3)<=5522;
@sum(para:h*y)>=4430.27*(1+r);
@sum(para:h*u)=13679.6988*(1-r);
End