ppppp - プログラム倉庫

(ステップ0）

任意の配分$c(m),y(m)$を適当に仮定する。

\begin{equation}
t=0
\end{equation}

(ステップ１）

\begin{equation}
c1=c(m)
\end{equation}

\begin{equation}
l1=y(m)/\theta(s)
\end{equation}

\begin{equation}
u1=\log(c1)-\frac{l1^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
u1>us　ならば、　ms=m　 　
\end{equation}

\begin{equation}
u1>us　ならば、　us=u1　 　
\end{equation}

\begin{equation}
gr(s)=m　 　
\end{equation}

(ステップ2）

\begin{equation}
c1=c(m-1)+n1*h
\end{equation}

\begin{equation}
c2=c(m)+n2*h
\end{equation}

\begin{equation}
c3=c(m)+n2*h
\end{equation}

\begin{equation}
l1=y(m-1)/\theta(s)
\end{equation}

\begin{equation}
u1=\log(c1)-\frac{l1^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
l2=y(m)/\theta(s)
\end{equation}

\begin{equation}
u2=\log(c2)-\frac{l2^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
l3=y(m+1)/\theta(s)
\end{equation}

\begin{equation}
u3=\log(c3)-\frac{l3^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
us=u1
\end{equation}

\begin{equation}
u2>us　ならば、　us=u2　 　
\end{equation}

\begin{equation}
u3>us　ならば、　us=u3　 　
\end{equation}

\begin{equation}
gr(s)=m　ならば、　fs=fs+us　 　
\end{equation}

\begin{equation}
f(m,n1,n2,n3)=fs
\end{equation}

(ステップ3）

\begin{equation}
c1=c(m-1)+n1*h
\end{equation}

\begin{equation}
c2=c(m)+n2*h
\end{equation}

\begin{equation}
c3=c(m)+n2*h
\end{equation}

\begin{equation}
l1=y(m-1)/\theta(s)
\end{equation}

\begin{equation}
u1=\log(c1)-\frac{l1^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
l2=y(m)/\theta(s)
\end{equation}

\begin{equation}
u2=\log(c2)-\frac{l2^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
l3=y(m+1)/\theta(s)
\end{equation}

\begin{equation}
u3=\log(c3)-\frac{l3^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
us=u1
\end{equation}

\begin{equation}
bx=y(m-1)-c(m-1)
\end{equation}

\begin{equation}
u2>us　ならば、bx=y(m)-c(m)
\end{equation}

\begin{equation}
u2>us　ならば、　us=u2　 　
\end{equation}

\begin{equation}
u3>us　ならば、bx=y(m+1)-c(m+1)
\end{equation}

\begin{equation}
gr(s)=m　ならば、　bs=bx+us　 　
\end{equation}

\begin{equation}
b(m,n1,n2,n3)=bs
\end{equation}

(ステップ4）

\begin{equation}
c1=c(1)+n1*h
\end{equation}

\begin{equation}
c2=c(2)+n2*h
\end{equation}

\begin{equation}
l1=y(m-1)/\theta(s)
\end{equation}

\begin{equation}
us=\log(c1)-\frac{l1^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
l2=y(m)/\theta(s)
\end{equation}

\begin{equation}
u2=\log(c2)-\frac{l2^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
u2>us　ならば、　us=u2　 　
\end{equation}

\begin{equation}
gr(s)=1　ならば、　fs=fs+us　 　
\end{equation}

\begin{equation}
fastu(n1,n2)=fs
\end{equation}

(ステップ5）

\begin{equation}
c1=c(1)+n1*h
\end{equation}

\begin{equation}
c2=c(2)+n2*h
\end{equation}

\begin{equation}
l1=y(1)/\theta(s)
\end{equation}

\begin{equation}
us=\log(c1)-\frac{l1^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
l2=y(2)/\theta(s)
\end{equation}

\begin{equation}
u2=\log(c2)-\frac{l2^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
bx=y(1)-c(1)　 　
\end{equation}

\begin{equation}
u2>us　ならば、　bx=y(2)-c(2)　 　
\end{equation}

\begin{equation}
gr(s)=1　ならば、　bs=bs+bx　 　
\end{equation}

\begin{equation}
fastb(n1,n2)=bs
\end{equation}

(ステップ6）

\begin{equation}
c1=c(M-1)+n1*h
\end{equation}

\begin{equation}
c2=c(M)+n2*h
\end{equation}

\begin{equation}
l1=y(M-1)/\theta(s)
\end{equation}

\begin{equation}
us=\log(c1)-\frac{l1^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
l2=y(M)/\theta(s)
\end{equation}

\begin{equation}
u2=\log(c2)-\frac{l2^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
u2>us　ならば、　us=u2　 　
\end{equation}

\begin{equation}
gr(s)=M　ならば、　fs=fs+us　 　
\end{equation}

\begin{equation}
endu(n1,n2)=fs
\end{equation}

(ステップ7）

\begin{equation}
c1=c(M-1)+n1*h
\end{equation}

\begin{equation}
c2=c(M)+n2*h
\end{equation}

\begin{equation}
l1=y(M-1)/\theta(s)
\end{equation}

\begin{equation}
us=\log(c1)-\frac{l1^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
l2=y(M)/\theta(s)
\end{equation}

\begin{equation}
u2=\log(c2)-\frac{l2^{1+\sigma(s)}}{1+\sigma(s)}
\end{equation}

\begin{equation}
bx=y(M-1)-c(M-1)　 　
\end{equation}

\begin{equation}
u2>us　ならば、　bx=y(M)-c(M)　 　
\end{equation}

\begin{equation}
gr(s)=M　ならば、　bs=bs+bx　 　
\end{equation}

\begin{equation}
endb(n1,n2)=bs
\end{equation}

(ステップ8)

\begin{equation}
q=fastb(n1,n2)
\end{equation}

\begin{equation}
v(1,n1,n2,q)=fastu(n1,n2)
\end{equation}

(ステップ9)

\begin{equation}
u1=u(m,nx,n1,n2)
\end{equation}

\begin{equation}
qs=q-b(m,nx,n1,n2)
\end{equation}

\begin{equation}
v1=u1+v(m-1,nx,n1,qs)
\end{equation}

\begin{equation}
v1>vs ならば、　nxs=nx
\end{equation}

\begin{equation}
v1>vs ならば、　qxs=qs
\end{equation}

\begin{equation}
v1>vs ならば、　vs=v1
\end{equation}

\begin{equation}
v(m,n1,n2,q)=vs
\end{equation}

\begin{equation}
goton(m,n1,n2,q)=nxs
\end{equation}

\begin{equation}
gotoq(m,n1,n2,q)=qxs
\end{equation}

(ステップ10)

\begin{equation}
u1=endu(n1,n2)
\end{equation}

\begin{equation}
qs=-endb(n1,n2)
\end{equation}

\begin{equation}
endv(n1,n2)=u1+v(M-1,n1,n2,qs)
\end{equation}

(ステップ11)

\begin{equation}
v1=endv(n1,n2)
\end{equation}

\begin{equation}
もし、v1>vs　ならば、nx1=n1
\end{equation}

\begin{equation}
もし、v1>vs　ならば、nx2=n2
\end{equation}

\begin{equation}
もし、v1>vs　ならば、qs=-endb(n1,n2)
\end{equation}

\begin{equation}
もし、v1>vs　ならば、vs=v1
\end{equation}

(ステップ12)

\begin{equation}
opc(M-1)=nx1
\end{equation}

\begin{equation}
opc(M)=nx2
\end{equation}

\begin{equation}
opq(M-1)=qs
\end{equation}

(ステップ13)

\begin{equation}
opc(m)=gotoc(m+1,opc(m),opc(m+1),opq(m+1))
\end{equation}

\begin{equation}
opq(m)=gotoq(m+1,opc(m),opc(m+1),opq(m+1))
\end{equation}

(ステップ14)

\begin{equation}
e=\sum opc(m)^2
\end{equation}

(ステップ15)

\begin{equation}
e=0 であれば、h=h/2
\end{equation}

(ステップ16)

\begin{equation}
h<\bar{h} であれば、t=1000
\end{equation}

(ステップ17)

\begin{equation}
c(m)=c(m)+opc(m)*h
\end{equation}

(ステップ18)

\begin{equation}
t>1000 であれば、終了です。
\end{equation}

(ステップ19)

\begin{equation}
t=t+1
\end{equation}

ステップ1へ戻る