J = 1 upper= 0 ;

else

upper= t̲a(i, J‑1) end

lower= 0 for j = J:life

「lshi = 1 form=J:j

end

「lsh l二「lshi*(l+t̲r (i+m‑1)) end

if j

>

「lis l二「lshi/(l+t̲r (i+J‑1)) else

「lisi = 1 end

upper= upper+t w(i+j‑l)*e(j)/rishi lower= lower+(l+t tau(i+j‑1))/rishi* ..

( (1 +rho）《（J‑j)*riisi*(l+t tau(i+J‑1))/.. (l+t tau(i+j‑1））

γ

（1/ganma)

t̲c(i,J) =upper/lower;

明今期以降（J歳以降）の消費を求める． if JくIife

end

for j = (J+l) I ife rishi = 1 for m =(J+l) :j

rishi = rishi*(l+t r(i+m‑1)) end

t̲c(i, j) = t̲c(i, J)*((l+rho）《（J‑j) ＊「lshi* ... (l+t tau(i+J‑1))/(l+t tau(i+j‑1））

γ

（1/ganma) end

t C (x)

=

0 . for j = 1 Ii fe

t̲C (x) = t̲C (x) +t̲N (x

− 」＋

1)*t̲c(x‑j+1,j) end

明税収と政府支出が一致しているか確認する， differ (x‑1 i fe, I oop) = t̲tau (x) *t̲C (x) ‑t̲G (x)

if differ (x‑1 i fe, I oop)

>

eps i I on tau h = t tau (x)

else if differ (x‑1 ife, loop）く −epsi I on tau I = t tau (x)

社会資本の新設更新・維持管理政策の世代間厚生分析付録

end end

首ここまでで，今期（x期）のtau決定！！！！

明今期の各個人の貯蓄aの計算

t a(x, 1) = t w(x)*e(l)‑t c(x, 1)*(1+t tau(x)) for j = 2:life

t a (x‑j+l, j) = t a (x‑j+l, j‑1)*(1+t r (x)) .. +t̲w (x) *e （」）−t̲c(x−」＋1，」）* (1 +t̲tau (x)) end

明今期の家計貯蓄を計算する t A (x)

=

0 ;

for j = 1: I i fe

t̲A(x) = t̲A(x)+t̲N(x‑j+l)*t̲a(x−」＋1, j) end

明もし家計貯蓄が負になったら計算は続行できない．

if t̲A (x）くO

danger= input （貯蓄が負に！！） end

明期（

x

）が一つ次へ進む end

disp （移行期の計算．終了） di sp （program end')

首移行過程（一時的均衡）の計算終了判明結果まとめとグラフ化のセクション

判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明｛同協判明判明判

明判明分析結果計算削判

明判明百判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明判明｛同判明判明判明明グラフ化したい期間数

show= 100 if show>move

show= move end

︶ W

H n u LH p o n ud

︵＋

ρU

z−

− 一

ーρ

u n H n U

7 ι

﹈

− −

nL Au nu nu nL F E E

− − L

m

﹈

nu dq dh u

−

nv nk un 日ロロリ 4E EA HV AH V

Entntnt 軸

α o o

nv nk un

口

︑ ︐

︐ ︐

AH vn ud nu d nu n t

・

フ

一一一一

﹇

ユノ

二ze

− − nt

ぷ

a a n t

eem

V I

−v−E

w而 a a

−

明figure(l）を2*6に分け，グラフをまとめて表示する目 figure(l)

subp I ot (2, 6, 1) hold on

plot(areal, t G(zonel），．）

plot(areal, t I (zonel），ー）

p I ot (areal, t̲M (zonel，）：） p I ot (area 1, t̲chouse i (zone 1），一．）

title（公的支出［］・新設［一］・維持菅理［］・調整費［一．］）

xi im([I iml])

subp I ot (2, 6, 2) hold on

p I ot (areal, t Y (zonel），ー）

title（総生産 Y） xi im([I iml])

社会資本の新設更新・維持管理政策の世代間厚生分析付録

︶︶ー

ρu

n H n U

7 ι

︵ v l

＋L 1

2 u

ρu

V 1 2 u

︶ー

ρu

n H n V

7 ι

︵凶

＋L 1

2 u

ρいv

︑

﹁﹈

2 u

nn u

︑ ︑

AU Hy

﹃1

︑i

J︑J︑i

＋L

・

1﹂O﹁L一︑iJ︑iJnuEd

− −

一 ︑

︐ ︐

1J1JrDAU

− 旦

4141F﹁

m m

nvnvニヨ

J e e

− − −

l﹁L﹁L

﹈

阜市

y y﹇

﹇

︐ ︐

ふLふL

m m ρ u

n︑un︑u

tE IE

1 J e e m m

ili−

−E

−

・1 E

HY Hy lL r I L l

− − ElElE

Z−

− −

﹄

i︶︐︐︐︐

︶

h−

− −

1J4lnL e44唱l唱lqL︑j1jfkfkn

善三 m e e

lηL

X X

−

ィ一 l n n

fkfkAnAn

n H

1 1 ﹄

II EI

E・

1E VA VA

ιH rl

︑﹁LlLlLAHAH

︑ ︐

l︑ef︑ιHιHιHιH t t J u

x−

−

mrlrlrlrl

e e

− −

ふL l

ふLふLふLふL

s s nu ln HI E

−−Eρνρu

ρν

ρu

LH rl L

ふL X S S S S

M河川鴨川

subp I ot (2, 6, 4) hold on

plot(areal, t tau(zonel），ー）

title（税率ピ） xi im([I iml])

明ylim([O, 0. 15])

︶ −

E E J

＋︵

H U 2 u

＋L

＋L ＊︶

E E J

︵ p u

＋L ＋︶︵

n H JU LU ρ U

v

・

e

二nuZIE

m n v 1

1j ＋

l l l

ρ ν

二−−

︑ZIE4EEnH ・

︑1 E

I円u

l l l

二

1Jfk一﹁lL4lnHIJIE

二二

一 V E

n u n H

・1

EI E Z IE

l nu v t

一O

iE Zl

end end

subp I ot (2, 6, 5) hold on

plot(1990:2089, i bdn(I ife+(l .show）），） title（世代別個人の負担総額）

xi im([1990, 2090]) 首ylim([O, 300])

i esp

= ] [

for i

=

1 I i fe+move

i csp(i)

=

sum(t c(i, :))

︶ a u

Ub − ︵

／︶

a u

m m q a

− ︵

︽︶ −

E E J

︵ p u

ふL＊︶︶ 1 ・ E E J

︵︵

︽︶

n u ln v l

︵︶＋︵ー＋

川U

ρ ν

一

u v ρ u

・1 E

o z l

m川・・

l E

二＋

1l e A

︑JU−−

− − − −

l二

二川

U．

︐

．1J一1・l・J・l︵﹈﹁L

二二一O

u r

川

− −

EP十l

j nu H Hu v l n H

一n u

ρ u t E Z IE

end

︶︶︶

n u 凶

LH FO

1 ︵＋

ρu

︵

川U

n ud nu hu nt nu n ud n ud

− ︶︶

n u 凶

1 ︵＋

ρu

︵︶ nド

p o

− −

J G＋一﹁LJ用

＆﹄

nxu

d e ︐ ︐ ︐ AU

主主

nL 刊以

か生

n u J V

−

n ud

−−

︶

︵−

﹈

︶

﹁LnunU1﹂

y

事ミ

1J1J

nu dn ud AU

ふL豆Snvnv

o

商

＋ 2

；

2

ー

9

主主

nunu

江川

Mn同

J vn

同J

A HV

二ト

t e e nu dn ud

﹁L

f l l

1﹂︑JVJVJ﹁L﹁L

ふLふLm川

ρ M

L︑n︑un︑u

t E n H

1jJauaum川m川

1 1E

．J

岡国 n n

・l・ly︶L︷

i i

−

p o l n H

H1IL

− −

﹄

ロn︐︐︐︐︑

ρHU4EE

怜﹁

︑

︐︐

︐ ︐

4 EE

守

4lnL︑iJ︑iJ

︑

n L n H

e e

− −

v A

︑14レ

n n

f︑f︑

＋Lnl﹄・l・lvAVA

ιH rl

︑lLlLAHAH

︑

el nH ln HI

n−

n H

ふLnvJux

− −

f︑f︑f︑f︑e

lot

− −

AふLふLふLふLふLS

n H Hu nu ln HI Eρ

uρい

vρ

いvρいVeshHrlLふ

L S S S S

M鴨川

subp I ot (2, 6, 7) hold on

plot(areal, t KG(zonel），ー）

title（インフラ水準 KG）

xi im([I iml])

明ylim([0,40])

subp I ot (2, 6, 8) hold on

xi im([I iml]) 明ylim([O, 200])

subp I ot (2, 6, 9) hold on

plot(areal, t̲L(zonel)) title（労働 L） xi im([I iml])

subp I ot (2, 6, 10)

figure（position，［100 100 350 350]) f

= ] [

hold on for X

=

1 10

f(x)

=

plot(20:79, t̲c(I ife‑9+10*x，・）） set (f (x），color，（x‑1)*0. 1*[1 1 1]) end

首title （消費の流列（10世代ごと，黒→白） y I i m ( [O, 0. 8])

subp I ot (2, 6, 11) f

= ] [

hold on for X

=

1 10

f(x)

=

plot(20:79, t a(I ife‑9+10*x, :)) set (f (x，） color，（x‑1)*0. 1*[1 1 1]) end

首title （貯蓄の流列（10世代ごと，黒→白）首yI im（一［100,400])

︶ −

＋︵

＋L

＋L ＊︶

︵ p

＋L

一一

ρ u

uvρu︑

nuZIEIJ

m l

＋

I l l

e−−f︑ZIE4EEVA

・

1 E 2 u

二

＋

﹂

一﹇

1 J I

− − 一一

n U

− −E

a u

ふ

L r

一O

i E Z l

end end

subp I ot (2, 6, 12)

figure（position，［100 100 350 350]) f

= ] [

hold on for X

=

1 10

f(x)

=

plot(20:79, i̲tax(I ife‑9+10*x，・）） set (f (x），color，（x‑1)*0. 1*[1 1 1]) end

首title （税負担の流列（10世代ごと，黒→白） ylim([0,0.2])

百以上でfigure(l）の表示分を終了

明各期tauの収束の正しさのチェック figure (2)

elf

for X

=

1:25 Y

=

4*X ; subp I ot (5, 5, x) p I ot (tautau (y, : ) ) end

明移行第1期世代の賃金収入と税率の推移．

figure (3)

xxx

=

w (

I i fe+ (1 I i fe)). *e (1 : I i fe) vvv

=

t tau ( I i fe+ (1: I i fe))

社会資本の新設更新・維持管理政策の世代間厚生分析付録

p I ot (vvv)

明10世代毎の消費の流列 figure (4)

figure（position，［100 100 350 350]) f

= ] [

hold on for X

=

1 10

社会資本の新設更新・維持管理政策の世代間厚生分析付録

f(x)

=

plot((2010:2069)+10*(x‑1), t̲c(I ife‑9+10*x, :)) set (f (x），color，（x‑1)*0. 1*[1 1 1])

end

首title （消費の流列（10世代ごと，黒→白） y I i m ( [O, 0. 8])

xlabel （年）

明10世代毎の税負担の流列 figure (5)

figure（position，［100 100 350 350]) f

= ] [

hold on for X

=

1・10

end

f(x)

=

plot((2010:2069)+10*(x‑1), i̲tax(life‑9+10*x，・）） set (f (x），color，（x‑1)*0. 1*[1 1 1])

首title （税負担の流列（10世代ごと，黒→白） ylim([0,0.2])

xlabel （年） grid on

社会資本の新設更新・維持管理政策の世代間厚生分析付録

5 . MA'fLAB

ソース：最適支出配分の経済

問準備〜初期定常状態計算終了までのセクション

鴨1%%%%品輸品輸1%%掛町晴樹品1%%%1%%%%品輸品輸品晴樹品1%%%1%%%%品W品W品、輸品輸旬、

判明百準備間百

%%%%%%品輸品輸1%%掛町晴樹品W品1%%%%品輸品輸品輸1%%%W品1%%%%品W品W部旬、輸品輸旬、

clear al I clear al I for zxcv = 1 :10

figure (zxcv) clear Marker clear Linestyle elf

c I ose (zxcv)

︶ 1 ・

E EJ

︵＊

凋斗

ρ o

＋ n

︽︶

ハ日

ー＋

−E E

︵＊ nu

凋斗

− n

︶ n日ー＋

−E E

︵＊

﹃ ︐ I

ρ0

﹁D

− n

＋﹁D﹁D守I

﹁ ︐

︐

nu ．︐．︐

AH V

一

内J﹁D・

凋斗

二

−n t n U 4 l ハU ρ 0

・︐︑iJ一一一一一・lJ二

f︑二

a

﹁L・

J e

me

n u m川干

一一

VE nH ln Ha u

−−E

rσ61le

end

e (1 : 40) = e (1: 40) /max (e) e(41 :60) = 0.0 ;

eps i I on = 0. 00000001 move= 159

kappa= 0.00472 ; aa = 2. 610

b = 2.059

問問 2009年日本を定常状態と想定したときの定常状態計算目

w =

1. 0

r = 0.01 M

=

I = 1/0. 3704*M ;

fun = @ (KG) KG* (M/KG）《kappa+l‑KG; KG = fzero (fun, (I +M) *60)

chouse i = KG* (aa/2* (I /KG）《2+b*I/KG）一I G = l+chousei+M

Y = (I +M) /0. 0734 ;

=

alpha一I= 0. 0 ; a I pha̲h = 1. 0 ; for loop= 1 :1000

alpha= (alpha一l+alpha̲h)/2;

V = v+l

N = Y*(l‑alpha)/sum(e) ; L = N*sum(e)

K = 100*alpha/(1‑alpha)*L . bunnki = 0

if bunnki ^{二二} 1

zxcv = [] ; vvv = 0 . zzzz = 0. 01 : 0. 01: 10. 0 . for tau = zzzz

vvv = vvv+l

upper = 0 ; I ower = 0 ; for j = 1: I i fe

= ] [

c(l) = w*upper/lower , for j = 2:life

社会資本の新設更新・維持管理政策の世代間厚生分析付録

c(j) = c(1)*((1+r)/(1+rho））《（(j‑1)/gamma)

end

C = sum (N. *c) zxcv(vvv) = C*tau , end

figure (v) hold on p I ot (zzzz, G) p I ot (zzzz, zxcv) hold off

tau I = 0.0 ; tau h = 2.0 for loop̲b = 1: 100

end

tau= (tau一l+tau̲h)/2; upper = 0 ; I ower = 0 ; for j = 1・Iife

upper = upper+e(j)/(l+r）《（j‑1)

lower = lower+(1+tau)/(1+r）《（j‑1)* ( (l+r)

I

(l+rho））《（（j‑1)/gamma) end

= ] [

c(l) = w*upper/lower ; for

」＝

2：・Iife

c(j) = c(1)*((1＋「）

I

(l+rho））《（(j‑1) /gamma) end

C = sum (N. *c) if G/tau‑C

>

tau I = tau elseif G/tau‑CくO

tau h = tau else

break end

= ] [

a(l) = w*e(1)‑c(1)*(1+tau) for j = 2:life

a（」）= a(j‑1)*(1+r)+w*e（」）−c(j) * (l+tau) end

Knew = sum (N. *a) delta(v) = (Knew‑K)/M ;

if delta(v)

>

epsi Ion alpha一I= alpha ; elseif delta(v）く −epsi I on

alpha̲h = alpha ; else

break end end

figure(10) p I ot (de I ta)

B = 1. 0/((1‑alpha）《（1‑aI pha) *a I pha《alpha*100《alpha*KG《（1‑aI pha))

間最適支出配分に変更する figure(l)

ドキュメント内社会資本の新設更新・維持管理政策の世代間厚生分析 (ページ 71-77)

>

γ

γ

=

− 」 ＋

>

=

x

m

α o o

︑ ︐

︐ ︐

・

eem

︑

−

v

e

= ] [

=

=

︐

u r

y

o

2

9

f l l

= ] [

=

=

= ] [

=

=

一 一

﹂

− − 一 一

L r

= ] [

=

=

=

=

=

w (

=

= ] [

=

=

= ] [

=

=

5 . MA'fLAB

a

me

w =

=

=

= ] [

I

= ] [

」 ＝

I

>

= ] [

>

− 」＋

一一

− − 一一

」＝