我用GAMS软件解决StoNED模型，运行中出现错误报告（第91行），修改多次仍然没有解决。求高手帮忙！正确解决者奖励100论坛币。

GAMS Rev 240  WIN-VS8 24.0.2 x86/MS Windows             06/06/14 13:36:47 Page 1

G e n e r a l   A l g e b r a i c   M o d e l i n g   S y s t e m

C o m p i l a t i o n

   1   

   2   

   3   

   4  SETS

   5  i index of industries/i1*i14/

   6  t time periods/t1*t7/;

   7  alias(i,h)

   8  alias(t,s);

   9   

  10  PARAMETERS

  11  table Y(i,t) output

  12           t1             t2                t3                 t4           

            t5                     t6              t7

  13  i1     9304.4        9743.994043        9962.495575        10679.50789   

          11151.10979        10942.60471        11494.88075

  14  i2     4236.3        4262.297394        4586.986742        4861.670498   

          5134.966965        4948.967829        4983.9177

  15  i3     12453.8       11974.63566        12242.55032        12898.87819   

          13922.95563        13749.02765        14835.90626

  16  i4     3664.8        3735.446552        3792.517279        4005.1948      

          4357.807158        4445.062131        4580.137524

  17  i5     5393          5543.670246        5654.483261        6751.019039   

          7176.186508        7256.202689        7786.88816

  18  i6     7174.1        7345.07495         8001.209463        8568.017503   

          9057.58567         10301.87981        11896.25508

  19  i7     2627.5        2603.537009        3059.444973        3403.856345   

          3664.021704        3753.263513        4368.412084

  20  i8     1759.5        1833.072592        2166.83057         2495.035244   

          2608.979732        2862.360053        3163.044583

  21  i9     768.6         759.7831855        763.3340772        805.2842337   

          826.7910629        897.4583089        983.1592567

  22  i10    2481.5        2797.519779        2858.426685        2897.524326   

          3023.646916        3195.624876        3423.858705

  23  i11    4892.6        5056.400363        5170.929585        5577.721968   

          5806.457256        6354.162595        6757.206829

  24  i12    2620.7        2623.254048        2684.536514        2910.055593   

          3024.103686        3081.095616        3356.017854

  25  i13    1043.2        1062.205018        1099.784214        1182.712415   

          1255.962068        1352.458483        1400.473074

  26  i14    6141.4        6104.486107        7131.861596        7852.264797   

          9005.328133        9191.195141        9096.116944

  27   

  28  table K(i,t) capital

        。。。。

  45  table L(i,t) labor

       。。。

62  Parameter Tr(t) /t1 1,t2 2,t3 3,t4 4,t5 5,t6 6,t7 7/;

  63   

  64  Parameter Tr2(t)/t1 1,t2 4,t3 9,t4 16,t5 25,t6 36,t7 49/;

  65   

  66  VARIABLES

  67  SS  sum of squares of errors

  68  V(i,t)  residual

  69  ThetaK  capital trend coefficient

  70  PsiK   capital trend coefficient 2

  71  Thetal  labor trend coefficient

  72  PsiL   labor trend coefficient 2

  73  A(i) inefficiency constant

  74  B(i) inefficiency  slope

  75  C(i) inefficiency  quadratic;

  76   

  77  POSITIVE VARIABLES

  78  BK(i,t)   marginal product of capital

  79  BL(i,t)   marginal product of labor;

  80   

  81  EQUATIONS

  82  QSSE  objective=sum of squares of errors

  83  QREGRESSION(i,t)  regression equation

  84  QCONCAVITY(i,h,t,s)  concavity constraint

  85  QKCON(t)  nonnegative capital embodied technical progress

  86  QLCON(t)  nonnegative labor embodied technical progress

  87  QEFF(i,t) nonnegative inefficiency;

  88   

  89  QSSE..     SS=e=sum(t,sum(i,V(i,t)*V(i,t)));

  90  QREGRESSION(i,t)..   Y(i,t) =e= BK(i,t)*K(i,t)+BL(i,t)*L(i,t)+(ThetaK*Tr(t

      )+PsiK*Tr2(t))*K(i,t)+(Thetal*Tr(t)+PsiL*Tr2(t))*L(i,t)-(A(i)+B(i)*Tr(t)+C

      (i)*Tr2(t))*Y(i,t)+Y(i,t)*V(i,t);

  91  QCONCAVITY(i,h,t,s).. BK(i,t)*K(i,t)+BL(i,t)*L(i,t) =1= BK(h,s)*K(i,t)+BL(

****                                                        $37,409

      h,s)*L(i,t);

  92  QKCON(t)..              ThetaK*Tr(t)+PsiK*Tr2(t) =g= 0;

  93  QLCON(t)..              ThetaL*Tr(t)+PsiL*Tr2(t) =g= 0;

  94   

  95  MODEL    StoNED/all/

  96   

  97  SOLVE StoNED using NLP Minimizing SS;

****                                      $257

GAMS Rev 240  WIN-VS8 24.0.2 x86/MS Windows             06/06/14 13:36:47 Page 2

G e n e r a l   A l g e b r a i c   M o d e l i n g   S y s t e m

Error Messages

37  '=l=' or '=e=' or '=g=' operator expected

257  Solve statement not checked because of previous errors

409  Unrecognizable item - skip to find a new statement

       looking for a ';' or a key word to get started again

**** 3 ERROR(S)   0 WARNING(S)

COMPILATION TIME     =        0.000 SECONDS      3 Mb  WIN240-240 Feb 14, 2013

USER: GAMS Development Corporation, Washington, DC   G871201/0000CA-ANY

      Free Demo,  202-342-0180,  sales@gams.com,  www.gams.com   DC0000

**** FILE SUMMARY

Input      C:\Users\wmt\Desktop\Untitled_1234.gms

Output     C:\Users\wmt\Documents\gamsdir\projdir\Untitled_1234.lst

**** USER ERROR(S) ENCOUNTERED