約 2,579,091 件
https://w.atwiki.jp/ce00582/pages/32.html
Private Sub Command1_Click() Dim ks As Single Dim a As Single Dim beta As Single Dim r1 As Single Dim n As Single Dim cx(1 To 100) As Single Dim cp(1 To 100) As Single Dim lx(1 To 100) As Single Dim lp(1 To 100) As Single Dim k(1 To 100) As Single Dim h As Single Dim k1 As Single Dim n1 As Single Dim n2 As Single Dim n3 As Single Dim c1 As Single Dim e As Single Dim ls As Single ls = 0.5 beta = 0.95 a = 0.33 ks = ls * ((1 / beta - 1) / a) ^ (1 / (a - 1)) h = 2 * ks / 100 For n = 1 To 100 k(n) = n * h lx(n) = ls cx(n) = k(n) ^ a * lx(n) ^ (1 - a) Next t1 = 0 Do Until t1 100 For n = 10 To 90 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) = c1 / (beta * (1 + r1)) Next w1 = (1 - a) * k(n) ^ a * lx(n) ^ (-a) lp(n) = 1 - cx(n) / w1 e = 0 For n = 10 To 90 e = e + (cx(n) - cp(n)) ^ 2 + (lx(n) - lp(n)) ^ 2 Next For n = 10 To 90 cx(n) = cp(n) lx(n) = lp(n) Next If e 10 ^ (-5) Then t1 = 1000 t1 = t1 + 1 Loop Dim ms As Single Dim px(1 To 100) As Single Dim ps(1 To 100) As Single Dim p1 As Single Dim p2 As Single Dim p3 As Single Dim z1 As Single Dim z2 As Single Dim z3 As Single Dim pi As Single Dim pxs As Single ms = 50 For n = 1 To 100 px(n) = 1 Next t4 = 0 Do Until t4 100 For n = 10 To 90 p1 = 1.1 * px(n) p2 = 0.9 * px(n) k1 = k(n) + k(n) ^ a * lx(n) ^ (1 - a) - cx(n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 pxs = px(n2) + (n1 - n2) * (px(n3) - px(n2)) l1 = lx(n2) + (n1 - n2) * (lx(n3) - lx(n2)) r1 = a * k1 ^ (a - 1) * l1 ^ (1 - a) pi = pxs / p1 - 1 iss = (1 + pi) * (1 + r1) - 1 z1 = iss * ms / (cx(n) * (1 + iss)) - p1 t3 = 0 Do Until t3 100 pi = pxs / p2 - 1 iss = (1 + pi) * (1 + r1) - 1 z2 = iss * ms / (cx(n) * (1 + iss)) - p2 p3 = p2 - z2 * (p2 - p1) / (z2 - z1) If z2 ^ 2 10 ^ (-5) Then t3 = 1000 p1 = p2 p2 = p3 z1 = z2 t3 = t3 + 1 Loop ps(n) = p2 Next e = 0 For n = 10 To 90 e = e + (px(n) - ps(n)) ^ 2 Next For n = 10 To 90 px(n) = ps(n) Next If e 10 ^ (-5) Then t4 = 1000 t4 = t4 + 1 Loop For n = 10 To 90 Debug.Print px(n) Next End Sub
https://w.atwiki.jp/ce00582/pages/25.html
Private Sub Command1_Click() Dim ks As Single Dim a As Single Dim beta As Single Dim r1 As Single Dim n As Single Dim cx(1 To 10, 1 To 100) As Single Dim cp(1 To 10, 1 To 100) As Single Dim k(1 To 100) As Single Dim h As Single Dim k1 As Single Dim n1 As Single Dim n2 As Single Dim n3 As Single Dim c1 As Single Dim e As Single Dim m As Single Dim th(1 To 10) As Single For m = 1 To 10 th(m) = 0.95 + 0.01 * m Next beta = 0.95 a = 0.33 ks = ((1 / beta - 1) / a) ^ (1 / (a - 1)) Debug.Print ks h = 2 * ks / 100 For n = 1 To 100 k(n) = n * h For m = 1 To 10 cx(m, n) = th(m) * k(n) ^ a Next Next t1 = 0 Do Until t1 100 For m = 1 To 10 For n = 10 To 90 k1 = k(n) + th(m) * k(n) ^ a - cx(m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 uc = 0 For m1 = 1 To 10 c1 = cx(m1, n2) + (n1 - n2) * (cx(m1, n3) - cx(m1, n2)) r1 = a * k1 ^ (a - 1) uc = uc + (beta * (1 + r1)) / c1 Next uc = 0.1 * uc cp(m, n) = 1 / uc Next Next e = 0 For m = 1 To 10 For n = 10 To 90 e = e + (cx(m, n) - cp(m, n)) ^ 2 Next Next For m = 1 To 10 For n = 10 To 90 cx(m, n) = cp(m, n) Next Next If e 10 ^ (-5) Then t1 = 1000 t1 = t1 + 1 Debug.Print t1, e Loop End Sub
https://w.atwiki.jp/ce00582/pages/26.html
Private Sub Command1_Click() Dim ks As Single Dim a As Single Dim beta As Single Dim r1 As Single Dim n As Single Dim cx(1 To 10, 1 To 100) As Single Dim cp(1 To 10, 1 To 100) As Single Dim lx(1 To 10, 1 To 100) As Single Dim lp(1 To 10, 1 To 100) As Single Dim k(1 To 100) As Single Dim h As Single Dim k1 As Single Dim n1 As Single Dim n2 As Single Dim n3 As Single Dim c1 As Single Dim e As Single Dim m As Single Dim th(1 To 10) As Single For m = 1 To 10 th(m) = 0.95 + 0.01 * m Next beta = 0.95 a = 0.33 ls = 0.5 ks = ls * ((1 / beta - 1) / a) ^ (1 / (a - 1)) h = 2 * ks / 100 For n = 1 To 100 k(n) = n * h For m = 1 To 10 lx(m, n) = ls cx(m, n) = th(m) * k(n) ^ a * lx(m, n) ^ (1 - a) Next Next t1 = 0 Do Until t1 100 For m = 1 To 10 For n = 10 To 90 k1 = k(n) + th(m) * k(n) ^ a * lx(m, n) ^ (1 - a) - cx(m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 uc = 0 For m1 = 1 To 10 c1 = cx(m1, n2) + (n1 - n2) * (cx(m1, n3) - cx(m1, n2)) l1 = lx(m1, n2) + (n1 - n2) * (lx(m1, n3) - lx(m1, n2)) r1 = th(m1) * a * k1 ^ (a - 1) * l1 ^ (1 - a) uc = uc + (beta * (1 + r1)) / c1 Next uc = 0.1 * uc cp(m, n) = 1 / uc w1 = (1 - a) * th(m) + k(n) ^ a * lx(m, n) ^ (-a) lp(m, n) = 1 - cx(m, n) / w1 Next Next e = 0 For m = 1 To 10 For n = 10 To 90 e = e + (cx(m, n) - cp(m, n)) ^ 2 + (lx(m, n) - lp(m, n)) ^ 2 Next Next For m = 1 To 10 For n = 10 To 90 lx(m, n) = lp(m, n) cx(m, n) = cp(m, n) Next Next If e 10 ^ (-5) Then t1 = 1000 t1 = t1 + 1 Debug.Print t1, e Loop End Sub
https://w.atwiki.jp/ce00582/pages/789.html
Function countb(bx As Single, bs As Single, h As Single) As Single Dim zs As Single Dim z1 As Single Dim z2 As Single Dim ct As Single zs = (bx - bs) / h z1 = Int(zs) z2 = z1 + 1 ct = z1 If (zs - z2) ^ 2 (zs - z1) ^ 2 Then ct = z2 countb = ct End Function Private Sub Command1_Click() Dim m As Single Dim n1 As Single Dim n2 As Single Dim n3 As Single Dim q As Single Dim qx As Single Dim pp As Single Dim h As Single Dim u(1 To 9, -1 To 1, -1 To 1, -1 To 1) As Single Dim b(1 To 9, -1 To 1, -1 To 1, -1 To 1) As Single Dim fastu(-1 To 1, -1 To 1) As Single Dim fastb(-1 To 1, -1 To 1) As Single Dim fastv(-1 To 1, -1 To 1, -100 To 100) As Single Dim endu(-1 To 1, -1 To 1) As Single Dim endb(-1 To 1, -1 To 1) As Single Dim endv(-1 To 1, -1 To 1) As Single Dim v(1 To 9, -1 To 1, -1 To 1, -1 To 1, -100 To 100) As Single Dim bx As Single Dim bs As Single Open "c /m2.txt" For Input As #2 Do Until EOF(2) Input #2, a1, a2, a3, a4, a5, a6 m = a1 n1 = a2 n2 = a3 n3 = a4 u(m, n1, n2, n3) = a5 b(m, n1, n2, n3) = a6 Loop Close #2 Open "c /m3.txt" For Input As #3 Do Until EOF(3) Input #3, a1, a2, a3, a4 n1 = a1 n2 = a2 fastu(n1, n2) = a3 fastb(n1, n2) = a4 Loop Close #3 Open "c /m4.txt" For Input As #4 Do Until EOF(4) Input #4, a1, a2, a3, a4 n1 = a1 n2 = a2 endu(n1, n2) = a3 endb(n1, n2) = a4 Loop Close #3 For n1 = -1 To 1 For n2 = -1 To 1 For q = -100 To 100 fastv(n1, n2, q) = -999 Next Next Next h = 0.01 For n1 = -1 To 1 For n2 = -1 To 1 bx = fastb(n1, n2) bs = fastb(0, 0) q = countb(bx, bs, h) fastv(n1, n2, q) = fastu(n1, n2) Next Next For m = 1 To 9 For n1 = -1 To 1 For n2 = -1 To 1 For n3 = -1 To 1 For q = -100 To 100 v(m, n1, n2, n3, q) = -999 Next Next Next Next Next m = 1 For n1 = -1 To 1 For n2 = -1 To 1 For n3 = -1 To 1 For q = -100 To 100 u1 = u(m, n1, n2, n3) bx = b(m, n1, n2, n3) bs = b(m, 0, 0, 0) qx = q - countb(bx, bs, h) pp = 0 If qx 100 Then pp = 1 If qx -100 Then pp = 1 If pp = 1 Then qx = 0 vs = u1 + fastv(n1, n2, qx) If pp = 1 Then vs = -999 v(m, n1, n2, n3, q) = vs Next Next Next Next For m = 2 To 9 For n1 = -1 To 1 For n2 = -1 To 1 For n3 = -1 To 1 For q = -100 To 100 u1 = u(m, n1, n2, n3) bx = b(m, n1, n2, n3) bs = b(m, 0, 0, 0) qx = q - countb(bx, bs, h) pp = 0 If qx 100 Then pp = 1 If qx -100 Then pp = 1 If pp = 1 Then qx = 0 vs = -999 For nx = -1 To 1 v1 = u1 + v(m - 1, nx, n1, n2, qx) If v1 vs Then vs = v1 Next If pp = 1 Then vs = -999 v(m, n1, n2, n3, q) = vs Next Next Next Next Next For n1 = -1 To 1 For n2 = -1 To 1 endv(n1, n2) = -999 Next Next For n1 = -1 To 1 For n2 = -1 To 1 u1 = endu(n1, n2) bx = endb(n1, n2) bs = endb(0, 0) qx = bs - bx pp = 0 If qx 100 Then pp = 1 If qx -100 Then pp = 1 If pp = 1 Then qx = 0 vs = -999 For nx = -1 To 1 v1 = u1 + v(9, nx, n1, n2, qx) If v1 vs Then vs = v1 Next If pp = 1 Then vs = -999 endv(n1, n2) = vs Debug.Print n1, n2, endv(n1, n2) Next Next End Sub
https://w.atwiki.jp/ce00582/pages/38.html
Private Sub Command1_Click() Dim ks As Single Dim a As Single Dim beta As Single Dim r1 As Single Dim n As Single Dim cx(1 To 10, 1 To 100) As Single Dim cp(1 To 10, 1 To 100) As Single Dim k(1 To 100) As Single Dim h As Single Dim k1 As Single Dim n1 As Single Dim n2 As Single Dim n3 As Single Dim c1 As Single Dim e As Single Dim m As Single Dim th(1 To 10) As Single For m = 1 To 10 th(m) = 0.95 + 0.01 * m Next beta = 0.95 a = 0.33 ks = ((1 / beta - 1) / a) ^ (1 / (a - 1)) Debug.Print ks h = 2 * ks / 100 For n = 1 To 100 k(n) = n * h For m = 1 To 10 cx(m, n) = th(m) * k(n) ^ a Next Next t1 = 0 Do Until t1 100 For m = 1 To 10 For n = 10 To 90 k1 = k(n) + th(m) * k(n) ^ a - cx(m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 uc = 0 For m1 = 1 To 10 c1 = cx(m1, n2) + (n1 - n2) * (cx(m1, n3) - cx(m1, n2)) r1 = a * k1 ^ (a - 1) uc = uc + (beta * (1 + r1)) / c1 Next uc = 0.1 * uc cp(m, n) = 1 / uc Next Next e = 0 For m = 1 To 10 For n = 10 To 90 e = e + (cx(m, n) - cp(m, n)) ^ 2 Next Next For m = 1 To 10 For n = 10 To 90 cx(m, n) = cp(m, n) Next Next If e 10 ^ (-5) Then t1 = 1000 t1 = t1 + 1 Loop Dim ms As Single Dim px(1 To 10, 1 To 100) As Single Dim ps(1 To 10, 1 To 100) As Single Dim p1 As Single Dim p2 As Single Dim p3 As Single Dim z1 As Single Dim z2 As Single Dim z3 As Single Dim pi As Single Dim pxs As Single ms = 50 For m = 1 To 10 For n = 1 To 100 px(m, n) = 1 Next Next t4 = 0 Do Until t4 100 For m = 1 To 10 For n = 10 To 90 p1 = 1.1 * px(m, n) p2 = 0.9 * px(m, n) k1 = k(n) + th(m) * k(n) ^ a - cx(m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 z1 = 0 For m1 = 1 To 10 pxs = px(m1, n2) + (n1 - n2) * (px(m1, n3) - px(m1, n2)) r1 = a * k1 ^ (a - 1) pi = pxs / p1 - 1 iss = (1 + pi) * (1 + r1) - 1 z1 = z1 + iss * ms / (cx(m, n) * (1 + iss)) Next z1 = 0.1 * z1 z1 = z1 - p1 t3 = 0 Do Until t3 100 z2 = 0 For m1 = 1 To 10 pxs = px(m1, n2) + (n1 - n2) * (px(m1, n3) - px(m1, n2)) r1 = a * k1 ^ (a - 1) pi = pxs / p2 - 1 iss = (1 + pi) * (1 + r1) - 1 z2 = z2 + iss * ms / (cx(m, n) * (1 + iss)) Next z2 = 0.1 * z2 z2 = z2 - p2 p3 = p2 - z2 * (p2 - p1) / (z2 - z1) If z2 ^ 2 10 ^ (-5) Then t3 = 1000 p1 = p2 p2 = p3 z1 = z2 t3 = t3 + 1 Loop ps(m, n) = p2 Next Next e = 0 For m = 1 To 10 For n = 10 To 90 e = e + (px(m, n) - ps(m, n)) ^ 2 Next Next For m = 1 To 10 For n = 10 To 90 px(m, n) = ps(m, n) Next Next If e 10 ^ (-5) Then t4 = 1000 t4 = t4 + 1 Debug.Print t4, e Loop For n = 10 To 90 Debug.Print px(5, n) Next End Sub
https://w.atwiki.jp/ankundora/pages/21.html
「プログラムの流れ」のページです まずはこのプログラムを入力してみましょう #include stdio.h int main(){ printf("Hello World\n"); return 0; } そしたら"Ctrl"+"F5"を押してください imageプラグインエラー ご指定のURLはサポートしていません。png, jpg, gif などの画像URLを指定してください。 こう表示されたと思います まず、ソースコードから実行ファイルが作られる事を「ビルド」と言います また、ビルドはだいたいは「コンパイル」と「リンク」で出来ています コンパイル ソースコードを、コンパイラをつかって人間が読みやすいコンピューター言語から コンピュータ上で実行可能な形式（オブジェクトファイル）に変換すること オブジェクトファイル 機械語(コンピュータのCPUが直接実行できる言葉)で書かれたコード リンク オブジェクトファイルやライブラリをリンカしてまとめ、実行可能ファイルにすること ビルド コンパイルやリンクをしてソースコードを実行可能ファイルにすること また、Visual C++の上の方に「デバッグ」というところがあります 見てみると今回行った"Ctrl"+"F5"は「デバッグなしで開始」というのに対応していることが分かると思います デバッグ プログラムの誤り(バグ)を探すこと 「バグ」とは英語で虫という意味です 今、「デバッグ開始」を選ぶとウィンドウがすぐに閉じてしまうので「デバッグなしで開始」にしています では、ソースコードを見ていきます 注意：プログラムは「 " ～ " 」の中など以外は半角文字しか使いません なのでもし全角スペースを使ってしまったらエラーになり、 見つけ出すのに苦労してしまうので気を付けてください プログラムの基本の基本となることがあります それは基本的に、ソースコードの上から実行されてゆく ということです これを分かっていないとプログラムができないので、注意してください #include stdio.h この部分は「インクルード文」といいます インクルード文 「include」とは「含める」という意味で ソースファイルの先頭などで、別のソースファイルなどを読み込み利用できるようにする機能のこと この場合は「printf」というものを利用できるようにするために「stdio.h」と言うものを読み込んでいます 「studio」ではなく「stdio(スタンダード(標準) アイオー(Input Output 入出力)」なので間違えないようにしてください また、「#include」のように「#」で始まるものを「プリプロセッサ命令」と言います では次を見てみましょう int main(){ printf("Hello World\n"); return 0; } int main()はプログラムのメイン部分であり、このなかで色々な処理が行われます そのままですね 中括弧の間( {～} )が「int main()」の内容、つまりメイン部分となります printfは次の「printf」でreturnは今度説明します 今回は難しい言葉が沢山出て来ました まだ無理やり全部覚えようとしなくていいので、流れだけ覚えておいてください
https://w.atwiki.jp/ce00582/pages/46.html
Private Sub Command1_Click() Dim ks As Single Dim a As Single Dim beta As Single Dim r1 As Single Dim n As Single Dim cx(1 To 10, 1 To 100) As Single Dim cp(1 To 10, 1 To 100) As Single Dim k(1 To 100) As Single Dim h As Single Dim k1 As Single Dim n1 As Single Dim n2 As Single Dim n3 As Single Dim c1 As Single Dim e As Single Dim m As Single Dim de As Single Dim th(1 To 10) As Single de = 0.05 For m = 1 To 10 th(m) = 0.95 + 0.01 * m Next de = 0.05 beta = 0.95 a = 0.33 ks = ((1 / beta + de - 1) / a) ^ (1 / (a - 1)) Debug.Print ks h = 2 * ks / 100 For n = 1 To 100 k(n) = n * h For m = 1 To 10 cx(m, n) = th(m) * k(n) ^ a - de * k(n) Next Next t1 = 0 Do Until t1 100 For m = 1 To 10 For n = 10 To 90 k1 = k(n) + th(m) * k(n) ^ a - cx(m, n) - de * k(n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 uc = 0 For m1 = 1 To 10 c1 = cx(m1, n2) + (n1 - n2) * (cx(m1, n3) - cx(m1, n2)) r1 = a * k1 ^ (a - 1) - de uc = uc + (beta * (1 + r1)) / c1 Next uc = 0.1 * uc cp(m, n) = 1 / uc Next Next e = 0 For m = 1 To 10 For n = 10 To 90 e = e + (cx(m, n) - cp(m, n)) ^ 2 Next Next For m = 1 To 10 For n = 10 To 90 cx(m, n) = cp(m, n) Next Next If e 10 ^ (-5) Then t1 = 1000 t1 = t1 + 1 Debug.Print t1, e Loop End Sub
https://w.atwiki.jp/program/pages/14.html
概要 ヤフーチャット内に作られているチャットルームです。 作成者によって名前が適当に変わります。 「のんびりプロ㌘」が本来の名称ですが、管理人が作らなかったときは他の人が部屋を作るため、若干名称が変わります。 チャットの参加方法 ヤフーメッセンジャー等からチャットへ参加します。 ヤフーチャット→コンピューターカテゴリー内のユーザールームに部屋があります。 部屋がなくても泣かない。 部屋がないときは自分で作る。 IRCにチャンネルがあります irc.freenode.net 6667#prog_kansaiです。 Wikipedia チャットに詳しい説明があります。 チャンネル名を#prog_kansaiと読みかえると大体同じ手順で参加できます。 文字コードはiso-2022-jpにしてください。 Mixiにコミュニティーがあります。 ［mixi］のんびりプログラムです。 Link http //bahoosoft.net/ ページビュー -
https://w.atwiki.jp/ce00582/pages/4640.html
import java.io.*; import java.net.*; class game0322{ String urlx; String bun; int year1,month1,day1; int year2,month2,day2; String ycode; String[] xml=new String[10000]; String[] xmlxml=new String[10000]; int[] xmllong=new int[10000]; int[] x=new int[10000]; int[] y=new int[10000]; int xmlnumber; int s,sx,ex; int h; public static void main(String [] args) { game0322 f=new game0322(); } game0322(){ year1=2005; month1=10; day1=1; year2=2005; month2=11; day2=1; ycode="4755.q"; urlx=getyahoo(ycode,year1,month1,day1,year2,month2,day2); getdata(urlx); xmlmaker(bun); for(s=1;s xmlnumber;s++){ xmlxml[s]=xmlxmlmaker(s); } for(s=1;s xmlnumber;s++){ xmlxml[s]=xmlxmlmaker(s); xmllong[s]=xmlxml[s].length(); } System.out.println(doragon(year1,month1)); } int doragon(int year1,int month1){ int s,z1; z1=0; for(s=1;s xmlnumber;s++){ if(check(s,year1,month1) 50)z1=z1+1; } return z1; } int check(int s,int year1,int month1){ String syear,smonth; int z; z=0; syear=String.valueOf(year1); smonth=String.valueOf(month1); h=0; if(xmllong[s] 2)h=h+1; if(xmllong[s] 20)h=h+1; if(xmlxml[s].indexOf(syear) -1)h=h+1; if(xmlxml[s].indexOf(smonth) -1)h=h+1; if(h==4)z=100; return z; } String getyahoo(String ycode,int year1,int month1,int day1,int year2,int month2,int day2){ String url; String syear,smonth,sday; String eyear,emonth,eday; syear=String.valueOf(year1); smonth=String.valueOf(month1); sday=String.valueOf(day1); eyear=String.valueOf(year2); emonth=String.valueOf(month2); eday=String.valueOf(day2); url="http //table.yahoo.co.jp/t?c="; url=url+syear; url=url+" a="; url=url+smonth; url=url+" b="; url=url+sday; url=url+" f="; url=url+eyear; url=url+" d="; url=url+emonth; url=url+" e="; url=url+eday; url=url+" g=d s="; url=url+ycode; url=url+" y=0 z="; url=url+ycode; url=url+" x=sb"; return url; } String xmlxmlmaker(int s){ int p1,p2; String str; p1=y[s]; p2=x[s+1]; str=bun.substring(p1+1,p2); return str; } void xmlmaker(String bun){ int p1; int s; s=1; p1=bun.indexOf(" "); x[1]=p1; while(p1 -1){ p1=bun.indexOf(" ",p1+1); if(p1 -1)s=s+1; if(p1 -1)x[s]=p1; } xmlnumber=s; for(s=1;s xmlnumber+1;s++){ y[s]=bun.indexOf(" ",x[s]); } for(s=1;s xmlnumber+1;s++){ xml[s]=bun.substring(x[s]+1,y[s]); } } void getdata(String name){ int num; int n; String data[]=new String[1000]; num=0; try{ URL url = new URL(name); HttpURLConnection con = (HttpURLConnection)url.openConnection(); con.setRequestMethod("GET"); con.connect(); BufferedReader reader = new BufferedReader(new InputStreamReader(con.getInputStream(),"UTF-8")); n=0; while (true){ String line = reader.readLine(); if ( line == null ){ break; } data[n]=line; n=n+1; num=n; } reader.close(); con.disconnect(); } catch(Exception e) { System.err.println(e); } bun=""; for (n=0;n num;n++){ bun=bun+data[n]; } } }
https://w.atwiki.jp/ce00582/pages/48.html
Private Sub Command1_Click() Dim a As Single Dim beta As Single Dim phi As Single Dim sig As Single Dim barm(1 To 10) As Single Dim t2 As Single Dim t1 As Single Dim t3 As Single Dim t4 As Single Dim t5 As Single Dim k(1 To 100) As Single Dim cx(-5 To 5, 1 To 10, 1 To 100) As Single Dim cp(-5 To 5, 1 To 10, 1 To 100) As Single Dim lx(-5 To 5, 1 To 10, 1 To 100) As Single Dim lp(-5 To 5, 1 To 10, 1 To 100) As Single Dim px(-5 To 5, 1 To 10, 1 To 100) As Single Dim ps(-5 To 5, 1 To 10, 1 To 100) As Single Dim checkc(-5 To 5, 1 To 10, 1 To 100) As Single Dim checkl(-5 To 5, 1 To 10, 1 To 100) As Single Dim checkp(-5 To 5, 1 To 10, 1 To 100) As Single Dim gotos(-5 To 5, 1 To 10, 1 To 100, 1 To 10) As Single Dim maki(-5 To 5, 1 To 10, 1 To 100, 1 To 10) As Single Dim v(0 To 2, -5 To 5, 1 To 10, 1 To 100) As Single Dim vs(0 To 2, -5 To 5, 1 To 10, 1 To 100) As Single Dim wedge(-5 To 5) As Single Dim ls As Single Dim ks As Single Dim s As Single Dim m As Single Dim n As Single Dim q As Single Dim m1 As Single Dim h As Single Dim n1 As Single Dim n2 As Single Dim n3 As Single Dim s1 As Single Dim s2 As Single Dim s3 As Single Dim q1 As Single Dim q2 As Single Dim q3 As Single Dim r1 As Single Dim c1 As Single Dim l1 As Single Dim y1 As Single Dim i1 As Single Dim ds As Single Dim dn As Single Dim dq As Single Dim e As Single Dim uc As Single Dim um As Single Dim mc As Single Dim p1 As Single Dim p2 As Single Dim p3 As Single Dim z1 As Single Dim z2 As Single Dim pxs As Single Dim pp As Single Dim ys As Single Dim vv1 As Single Dim vv2 As Single Dim vv0 As Single Dim price(-5 To 5, 1 To 10, 1 To 100) As Single Dim maxp(-5 To 5, 1 To 10, 1 To 100) As Single Dim pi As Single Dim pipi As Single Dim pipipi As Single Dim r(0 To 2, -5 To 5, 1 To 10, 1 To 100) As Single Dim dv(0 To 2) As Single Dim dr(0 To 2) As Single Dim rank As Single Dim prob(1 To 10, 1 To 10) As Single Open "c /pro146.txt" For Input As #81 Do Until EOF(81) Input #81, a1, a2, a3 m = a1 m1 = a2 prob(m, m1) = a3 Loop Close #81 phi = 0.9 Open "c /pro1551.txt" For Input As #2 Do Until EOF(2) Input #2, a1, a2, a3, a4, a5, a6 s = a1 m = a2 n = a3 lx(s, m, n) = a4 cx(s, m, n) = a5 px(s, m, n) = a6 Loop Close #2 Open "c /pro1552.txt" For Input As #72 Do Until EOF(72) Input #72, a1, a2, a3, a4, a5 s = a1 m = a2 n = a3 m1 = a4 gotos(s, m, n, m1) = a5 Loop Close #72 phi = 0.9 For s = -5 To 5 sig = 1 / (1 - phi) px1 = 1 + 0.01 * s px2 = ((1 - 0.1 * px1 ^ (1 - sig)) / 0.9) ^ (1 / (1 - sig)) wedge(s) = phi * px2 Next beta = 0.95 a = 0.33 ls = ((1 - a) * phi) / ((1 - a) * phi + 1) ks = ls * ((1 / beta - 1) / (a * phi)) ^ (1 / (a - 1)) h = 2 * ks / 100 For n = 1 To 100 k(n) = n * h Next t5 = 0 Do Until t5 1 t1 = 0 Do Until t1 100 For s = -5 To 5 For n = 10 To 90 For m = 1 To 10 uc = 0 k1 = k(n) + k(n) ^ a * lx(s, m, n) ^ (1 - a) - cx(s, m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 For m1 = 1 To 10 s1 = gotos(s, m, n, m1) s2 = Int(s1) If s2 4 Then s2 = 4 s3 = s2 + 1 ds = (s1 - s2) * (cx(s3, m1, n2) - cx(s2, m1, n2)) dn = (n1 - n2) * (cx(s2, m1, n3) - cx(s2, m1, n2)) c1 = cx(s2, m1, n2) + ds + dn ds = (s1 - s2) * (lx(s3, m1, n2) - lx(s2, m1, n2)) dn = (n1 - n2) * (lx(s2, m1, n3) - lx(s2, m1, n2)) l1 = lx(s2, m1, n2) + ds + dn mc = wedge(s2) + (s1 - s2) * (wedge(s3) - wedge(s2)) r1 = mc * a * k1 ^ (a - 1) * l1 ^ (1 - a) uc = uc + prob(m, m1) * (beta * (1 + r1)) / c1 Next uc = uc cp(s, m, n) = 1 / uc w1 = wedge(s) * (1 - a) * k(n) ^ a * lx(s, m, n) ^ (-a) lp(s, m, n) = 1 - cx(s, m, n) / w1 Next Next Next e = 0 For s = -5 To 5 For n = 10 To 90 For m = 1 To 10 e = e + (cx(s, m, n) - cp(s, m, n)) ^ 2 + (lx(s, m, n) - lp(s, m, n)) ^ 2 Next Next Next For s = -5 To 5 For n = 10 To 90 For m = 1 To 10 cx(s, m, n) = cp(s, m, n) lx(s, m, n) = lp(s, m, n) Next Next Next If e 10 ^ (-5) Then t1 = 1000 Debug.Print 81, t1, e t1 = t1 + 1 Loop For m = 1 To 10 barm(m) = 20 * (0.95 + 0.01 * m) Next t2 = 0 Do Until t2 100 For s = -5 To 5 For n = 10 To 90 For m = 1 To 10 p1 = 0.9 * px(s, m, n) p2 = 1.1 * px(s, m, n) k1 = k(n) + k(n) ^ a * lx(s, m, n) ^ (1 - a) - cx(s, m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 um = 0 For m1 = 1 To 10 s1 = gotos(s, m, n, m1) s2 = Int(s1) s3 = s2 + 1 If s3 5 Then s3 = 5 ds = (s1 - s2) * (cx(s3, m1, n2) - cx(s2, m1, n2)) dn = (n1 - n2) * (cx(s2, m1, n3) - cx(s2, m1, n2)) c1 = cx(s2, m1, n2) + ds + dn ds = (s1 - s2) * (lx(s3, m1, n2) - lx(s2, m1, n2)) dn = (n1 - n2) * (lx(s2, m1, n3) - lx(s2, m1, n2)) l1 = lx(s2, m1, n2) + ds + dn ds = (s1 - s2) * (px(s3, m1, n2) - px(s2, m1, n2)) dn = (n1 - n2) * (px(s2, m1, n3) - px(s2, m1, n2)) pxs = px(s2, m1, n2) + ds + dn mc = wedge(s2) + (s1 - s2) * (wedge(s3) - wedge(s2)) r1 = mc * a * k1 ^ (a - 1) * l1 ^ (1 - a) pp = pxs / p1 i1 = (1 + r1) * pp - 1 um = um + prob(m, m1) * (beta * i1) / (c1 * pp) Next z1 = um * barm(m) - p1 t3 = 0 Do Until t3 100 k1 = k(n) + k(n) ^ a * lx(s, m, n) ^ (1 - a) - cx(s, m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 um = 0 For m1 = 1 To 10 s1 = gotos(s, m, n, m1) s2 = Int(s1) s3 = s2 + 1 If s3 5 Then s3 = 5 ds = (s1 - s2) * (cx(s3, m1, n2) - cx(s2, m1, n2)) dn = (n1 - n2) * (cx(s2, m1, n3) - cx(s2, m1, n2)) c1 = cx(s2, m1, n2) + ds + dn ds = (s1 - s2) * (lx(s3, m1, n2) - lx(s2, m1, n2)) dn = (n1 - n2) * (lx(s2, m1, n3) - lx(s2, m1, n2)) l1 = lx(s2, m1, n2) + ds + dn ds = (s1 - s2) * (px(s3, m1, n2) - px(s2, m1, n2)) dn = (n1 - n2) * (px(s2, m1, n3) - px(s2, m1, n2)) pxs = px(s2, m1, n2) + ds + dn mc = wedge(s2) + (s1 - s2) * (wedge(s3) - wedge(s2)) r1 = mc * a * k1 ^ (a - 1) * l1 ^ (1 - a) pp = pxs / p2 i1 = (1 + r1) * pp - 1 um = um + prob(m, m1) * (beta * i1) / (c1 * pp) Next z2 = um * barm(m) - p2 p3 = p2 - z2 * (p2 - p1) / (z2 - z1) p1 = p2 p2 = p3 z1 = z2 If (z2) ^ 2 10 ^ (-5) Then t3 = 10000 t3 = t3 + 1 Loop ps(s, m, n) = p2 Next Next Next e = 0 For s = -5 To 5 For n = 10 To 90 For m = 1 To 10 e = e + (px(s, m, n) - ps(s, m, n)) ^ 2 Next Next Next If e 10 ^ (-5) Then t2 = 1000 For s = -5 To 5 For n = 10 To 90 For m = 1 To 10 px(s, m, n) = ps(s, m, n) Next Next Next Debug.Print 82, t2, e t2 = t2 + 1 Loop For s = -5 To 5 For m = 1 To 10 For n = 10 To 90 ys = k(n) ^ a * lx(s, m, n) ^ (1 - a) pi = px(s, m, n) * ys - wedge(s) * px(s, m, n) * ys pipi = (phi - wedge(s)) * (ys / (phi - 1)) pipipi = (mr + wedge(s) * ys) / (px(s, m, n) * (phi - 1)) r(2, s, m, n) = pipipi / 2 r(1, s, m, n) = pipi - 2 * px(s, m, n) * r(2, s, m, n) r(0, s, m, n) = pi - r(1, s, m, n) * px(s, m, n) - r(2, s, m, n) * px(s, m, n) ^ 2 Next Next Next t4 = 0 Do Until t4 100 For s = -5 To 5 For m = 1 To 10 For n = 10 To 90 k1 = k(n) + k(n) ^ a * lx(s, m, n) ^ (1 - a) - cx(s, m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 vv0 = 0 vv1 = 0 vv2 = 0 For m1 = 1 To 10 s1 = gotos(s, m, n, m1) s2 = Int(s1) s3 = s2 + 1 If s3 5 Then s3 = 5 ds = (s1 - s2) * (lx(s3, m1, n2) - lx(s2, m1, n2)) dn = (n1 - n2) * (lx(s2, m1, n3) - lx(s2, m1, n2)) l1 = lx(s2, m1, n2) + ds + dn ds = (s1 - s2) * (px(s3, m1, n2) - px(s2, m1, n2)) dn = (n1 - n2) * (px(s2, m1, n3) - px(s2, m1, n2)) pxs = px(s2, m1, n2) + ds + dn mc = wedge(s2) + (s1 - s2) * (wedge(s3) - wedge(s2)) r1 = mc * a * k1 ^ (a - 1) * l1 ^ (1 - a) pp = pxs / px(s, m, n) i1 = (1 + r1) * pp - 1 For rank = 0 To 2 ds = (s1 - s2) * (v(rank, s3, m1, n2) - v(rank, s2, m1, n2)) dn = (n1 - n2) * (v(rank, s2, m1, n3) - v(rank, s2, m1, n2)) dv(rank) = v(rank, s2, m1, n2) + ds + dn Next For rank = 0 To 2 ds = (s1 - s2) * (r(rank, s3, m1, n2) - r(rank, s2, m1, n2)) dn = (n1 - n2) * (r(rank, s2, m1, n3) - r(rank, s2, m1, n2)) dr(rank) = r(rank, s2, m1, n2) + ds + dn Next vv0 = vv0 + prob(m, m1) * (dr(0) + 0.8 * dv(0)) / (1 + i1) vv1 = vv1 + prob(m, m1) * (dr(1) + 0.8 * dv(1)) / (1 + i1) vv2 = vv2 + prob(m, m1) * (dr(2) + 0.8 * dv(2)) / (1 + i1) Next vs(0, s, m, n) = vv0 vs(1, s, m, n) = vv1 vs(2, s, m, n) = vv2 Next Next Next e = 0 For s = -5 To 5 For m = 1 To 10 For n = 10 To 90 For q = 0 To 2 e = e + (v(q, s, m, n) - vs(q, s, m, n)) ^ 2 Next Next Next Next For s = -5 To 5 For m = 1 To 10 For n = 10 To 90 For q = 0 To 2 v(q, s, m, n) = vs(q, s, m, n) Next Next Next Next If e 10 ^ (-5) Then t4 = 1000 Debug.Print 83, t4, e t4 = t4 + 1 Loop For s = -5 To 5 For m = 1 To 10 For n = 10 To 90 price(s, m, n) = -v(1, s, m, n) / (2 * v(2, s, m, n)) If price(s, m, n) 1.05 * px(s, m, n) Then price(s, m, n) = 1.05 * px(s, m, n) If price(s, m, n) 0.95 * px(s, m, n) Then price(s, m, n) = 0.95 * px(s, m, n) Next Next Next For s = -5 To 5 For m = 1 To 10 For n = 10 To 90 p2 = price(s, m, n) p1 = (1 + 0.01 * s) * px(s, m, n) p3 = (0.8 * p1 ^ (1 - sig) + 0.2 * p2 ^ (1 - sig)) ^ (1 / (1 - sig)) k1 = k(n) + k(n) ^ a * lx(s, m, n) ^ (1 - a) - cx(s, m, n) n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 For m1 = 1 To 10 s1 = gotos(s, m, n, m1) s1 = gotos(s, m, n, m1) s2 = Int(s1) s3 = s2 + 1 If s3 5 Then s3 = 5 ds = (s1 - s2) * (px(s3, m1, n2) - px(s2, m1, n2)) dn = (n1 - n2) * (px(s2, m1, n3) - px(s2, m1, n2)) pxs = px(s2, m1, n2) + ds + dn gt = 100 * (p3 / pxs - 1) If gt 5 Then gt = 5 If gt -5 Then gt = -5 maki(s, m, n, m1) = gt Next Next Next Next For s = -5 To 5 For m = 1 To 10 For n = 10 To 90 For m1 = 1 To 10 gotos(s, m, n, m1) = maki(s, m, n, m1) Next Next Next Next e = 0 For s = -5 To 5 For m = 1 To 10 For n = 10 To 90 e = e + (checkc(s, m, n) - cx(s, m, n)) ^ 2 + (checkl(s, m, n) - lx(s, m, n)) ^ 2 + (checkp(s, m, n) - px(s, m, n)) ^ 2 e = e + (price(s, m, n) - maxp(s, m, n)) ^ 2 Next Next Next For s = -5 To 5 For n = 10 To 90 For m = 1 To 10 checkc(s, m, n) = cx(s, m, n) checkl(s, m, n) = lx(s, m, n) checkp(s, m, n) = px(s, m, n) maxp(s, m, n) = price(s, m, n) Next Next Next If e 10 ^ (-5) Then t5 = 1000 Debug.Print 84, t5, e t5 = t5 + 1 Loop Open "c /pro1551.txt" For Output As #1 For s = -5 To 5 For m = 1 To 10 For n = 1 To 100 Write #1, s, m, n, lx(s, m, n), cx(s, m, n), px(s, m, n) Next Next Next Close #1 Open "c /pro1552.txt" For Output As #12 For s = -5 To 5 For m = 1 To 10 For n = 1 To 100 For m1 = 1 To 10 Write #12, s, m, n, m1, gotos(s, m, n, m1) Next Next Next Next Close #12 Dim mm(1 To 100) As Single For t = 1 To 100 mm(t) = 5 Next mm(1) = 8 sx = 0 mx = 5 nx = 50 k1 = k(50) For t = 1 To 100 s1 = sx s2 = Int(s1) s3 = s2 + 1 n1 = nx n2 = Int(n1) n3 = n2 + 1 m1 = mm(t) ds = (s1 - s2) * (gotos(s3, mx, n2, m1) - gotos(s2, mx, n2, m1)) dn = (n1 - n2) * (gotos(s2, mx, n3, m1) - gotos(s2, mx, n2, m1)) s1 = gotos(s2, mx, n2, m1) + ds + dn s2 = Int(s1) If s2 4 Then s2 = 4 s3 = s2 + 1 ds = (s1 - s2) * (cx(s3, m1, n2) - cx(s2, m1, n2)) dn = (n1 - n2) * (cx(s2, m1, n3) - cx(s2, m1, n2)) c1 = cx(s2, m1, n2) + ds + dn ds = (s1 - s2) * (lx(s3, m1, n2) - lx(s2, m1, n2)) dn = (n1 - n2) * (lx(s2, m1, n3) - lx(s2, m1, n2)) l1 = lx(s2, m1, n2) + ds + dn y1 = k1 ^ a * l1 ^ (1 - a) k1 = k1 + k1 ^ a * l1 ^ (1 - a) - c1 sx = s1 mx = m1 nx = k1 / h Debug.Print y1 Next End Sub