アルゴリズム２

(ステップ1）

$k(n)=n*h$

（ステップ2）

$lx(n)=(1-a)/(2-a)$

$cx(n)=k(n)^a*lx(n)^{1-a}$

（ステップ3）

$k1=k(n)+k(n)^a*lx(n)^{1-a}-cx(n)$

$n1=k1/h$

$n2=int(n1)$

$n3=n2+1$

$c1=cx(n2)+(n1-n2)*(cx(n3)-cx(n2))$

$l1=lx(n2)+(n1-n2)*(lx(n3)-lx(n2))$

$r1=a*k1^{a-1}*l1^{1-a}$

$cp(n)=\frac{c1}{\beta*(1+r1)}$

$w1=(1-a)*k(n)^{a}*lx(n)^{-a}$

$lp(n)=1-\frac{cx(n)}{w1}$

（ステップ4）

$e1=\sum (cp(n)-cx(n))^2+ (lp(n)-lx(n))^2$

（ステップ5）

$e1 &lt; e2$

この式が成立すれば、プログラムを終了する。

（ステップ6）

$cx(n)=cp(n)$

$lx(n)=lp(n)$

ステップ３へ戻る

最終更新：2009年02月13日 18:34

プログラム倉庫