我可能在论坛里没有把问题描述清楚。 可以加我qq讨论一下 287780095

parameter
        TauZ1                output tax on Z1,
        TauZ2                output tax on Z2,
        TauH                        income tax;

TauZ1 = 100/900;
TauZ2 = 100/900;
TauH = 1/9;

$ontext
$model:tax
$peps:1e-10

$sectors:        
        s1
        s2
        
$commodities:
        p1
        p2
        pf1
        pf2
        
$consumers:
        H        ! Household
        G        ! Government

$prod:S1        s:0  va:1
        o:P1        q:900                a:G t:TauZ1
        i:P1        q:50
        i:P2        q:250
        i:PF1        q:(350/(1+TauH))        a:G t:TauH        va:
        i:PF2        q:(150/(1+TauH))        a:G t:TauH        va:

$prod:S2        s:0  va:1
        o:P2        q:900                a:G t:TauZ1
        i:P1        q:250
        i:P2        q:50
        i:PF1        q:(350/(1+TauH))        a:G t:TauH        va:
        i:PF2        q:(150/(1+TauH))        a:G t:TauH        va:

$demand:H        s:0.5
        d:P1        q:450
        d:P2        q:450
        e:PF1        q:(700/(1+TauH))
        e:PF2        q:(300/(1+TauH))

$demand:G  
        d:P1        q:150
        d:P2        q:150
        
$offtext
$sysinclude mpsgeset tax

tax.iterlim=0;
$include tax.gen
solve tax using mcp;

哎，格式全乱了，凑合看吧，你要做的就是算出income tax rate.

看格式，用的是MPSGE吧。

没错，mpsge