約 2,225,481 件
https://w.atwiki.jp/ce00582/pages/788.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(2 To 9, -1 To 1, -1 To 1, -1 To 1) As Single Dim b(2 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, -50 To 50) 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(2 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 /m3.txt" For Input As #1 Do Until EOF(1) Input #1, a1, a2, a3, a4 n1 = a1 n2 = a2 fastu(n1, n2) = a3 fastb(n1, n2) = a4 Loop Close #1 For n1 = -1 To 1 For n2 = -1 To 1 For q = -50 To 50 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) Debug.Print n1, n2, q Next Next End Sub
https://w.atwiki.jp/ce00582/pages/755.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(0 To 10, 1 To 100) As Single Dim tk(0 To 10, 1 To 100) As Single Dim cs(0 To 10, 1 To 100) As Single Dim k(1 To 100) As Single Dim b(0 To 10) 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 g As Single g = 0.05 beta = 0.95 a = 0.33 ks = ((1 / beta - 1) / a) ^ (1 / (a - 1)) h = 2 * ks / 100 For n = 1 To 100 k(n) = n * h For m = 0 To 10 cx(m, n) = k(n) ^ a - g tk(m, n) = 0.05 Next Next For m = 0 To 10 b(m) = 0.01 * m Next t1 = 0 Do Until t1 100 For m = 0 To 10 For n = 10 To 90 k1 = k(n) + k(n) ^ a - cx(m, n) - g r1 = (1 - tk(m, n)) * a * k(n) ^ (a - 1) tax = tk(m, n) * a * k(n) ^ a m1 = (1 + r1) * b(m) + tax - g m2 = Int(m1) If m2 9 Then m2 = 9 If m2 0 Then m2 = 0 m3 = m2 + 1 n1 = k1 / h n2 = Int(n1) n3 = n2 + 1 dm = (m1 - m2) * (cx(m3, n2) - cx(m2, n2)) dn = (n1 - n2) * (cx(m2, n3) - cx(m2, n2)) c1 = cx(m2, n2) + dm + dn dm = (m1 - m2) * (tk(m3, n2) - tk(m2, n2)) dn = (n1 - n2) * (tk(m2, n3) - tk(m2, n2)) tk1 = tk(m, n2) + dm + dn r2 = (1 - tk1) * a * k1 ^ (a - 1) cs(m, n) = c1 / (beta * (1 + r2)) Next Next e = 0 For m = 0 To 10 For n = 10 To 90 e = e + (cx(m, n) - cs(m, n)) ^ 2 Next Next For m = 0 To 10 For n = 10 To 90 cx(m, n) = cs(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/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/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/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/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/47.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 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 beta = 0.95 a = 0.33 ls = 0.5 ks = ls * ((1 / beta + de - 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) - 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 * lx(m, n) ^ (1 - 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)) l1 = lx(m1, n2) + (n1 - n2) * (lx(m1, n3) - lx(m1, n2)) r1 = th(m1) * a * k1 ^ (a - 1) * l1 ^ (1 - a) - de 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/45.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 Dim ds As Single de = 0.05 ls = 0.5 beta = 0.95 a = 0.33 ks = ls * ((1 / beta + de - 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) - de * k(n) Next t1 = 0 Do Until t1 100 For n = 10 To 90 k1 = k(n) + k(n) ^ a * lx(n) ^ (1 - a) - cx(n) - de * k(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) - de cp(n) = c1 / (beta * (1 + r1)) w1 = (1 - a) * k(n) ^ a * lx(n) ^ (-a) lp(n) = 1 - cx(n) / w1 Next 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 Debug.Print t1, e Loop End Sub
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