(ステップ1)
\begin{equation}
px(n)=1
\end{equation}
(ステップ2)
\begin{equation}
p1=0.9*px(n)
\end{equation}
\begin{equation}
p2=1.1*px(n)
\end{equation}
(ステップ3)
\begin{equation}
k1=k(n)+k(n)^a-cx(n)
\end{equation}
\begin{equation}
n1=k1/h
\end{equation}
\begin{equation}
n2=int(n1)
\end{equation}
\begin{equation}
n3=n2+1
\end{equation}
\begin{equation}
c1=cx(n2)+(n1-n2)*(cx(n3)-cx(n2))
\end{equation}
\begin{equation}
pp=px(n2)+(n1-n2)*(px(n3)-px(n2))
\end{equation}
\begin{equation}
r1=a*k1^(a-1)*l1^(1-a)
\end{equation}
\begin{equation}
i1=(1+r1)*(pp/p1)-1
\end{equation}
\begin{equation}
um=\frac{\beta*i1}{c1}
\end{equation}
\begin{equation}
z1=um*\bar{m}-p1
\end{equation}
(ステップ4)
\begin{equation}
k1=k(n)+k(n)^a-cx(n)
\end{equation}
\begin{equation}
n1=k1/h
\end{equation}
\begin{equation}
n2=int(n1)
\end{equation}
\begin{equation}
n3=n2+1
\end{equation}
\begin{equation}
c1=cx(n2)+(n1-n2)*(cx(n3)-cx(n2))
\end{equation}
\begin{equation}
pp=px(n2)+(n1-n2)*(px(n3)-px(n2))
\end{equation}
\begin{equation}
r1=a*k1^(a-1)*l1^(1-a)
\end{equation}
\begin{equation}
i1=(1+r1)*(pp/p2)-1
\end{equation}
\begin{equation}
um=\frac{\beta*i1}{c1}
\end{equation}
\begin{equation}
z2=um*\bar{m}-p2
\end{equation}
(ステップ5)
\begin{equation}
p3=p2-z2*(p2-p1)/(z2-z1)
\end{equation}
\begin{equation}
p1=p2
\end{equation}
\begin{equation}
p2=p3
\end{equation}
\begin{equation}
p1=p2
\end{equation}
(ステップ6)
\begin{equation}
z1^2 < \bar{e}
\end{equation}
上記の式が成立する場合、ステップ3へ戻る。
(ステップ7)
\begin{equation}
ps(n)=p2
\end{equation}
(ステップ8)
\begin{equation}
e2=\sum (px(n)-ps(n))^2
\end{equation}
\begin{equation}
e2 < \bar{e}
\end{equation}
上記の式が成立する場合、プログラムを終了させる。
(ステップ9)
\begin{equation}
px(n)=ps(n)
\end{equation}
ステップ2へ戻る。
\end{document}
最終更新:2009年07月14日 22:38