untitled

Size: px
Start display at page:

Download "untitled"

Transcription

1 1 2 3 WRITE (*,*) WRITE (*,*) depth(x,y) = + + WRITE (*,*) WRITE (*,*) ( ) d1 WRITE (*,*) d1 = u * x ^r(x) WRITE (*,*) u = 0.1 r=2/3 WRITE (*,*) WRITE (*,*) d2 WRITE (*,*) d2 = b exp -(x-xb)^2 / (xb/2) WRITE (*,*), b = ; xb = WRITE (*,*) WRITE (*,*) WRITE (*,*) d3 = a(1-x/lb) * sin (2pi*(y-delta)/ramda) WRITE (*,*), a= ; ib = breaker length WRITE (*,*) delta = ( ) WRITE (*,*) DELTA = DELTA1 + DELTA2 WRITE (*,*) delta1 = for skewness WRITE (*,*) delta2 = for the oblique downstream orientation of the crest of unduration WRITE (*,*) delta1 = deltamax * sin( 2pi*(y-delta1)/ramda ) WRITE (*,*) = deltamax sin(2pi*y/ramda)*cos(2pi*delta1/ramda) WRITE (*,*) - cos(2pi*y/ramda)*sin(2pi*delta1/ramda) WRITE (*,*)

2

3

4

5 C C C ///////////////////////////////////////////////////////////// C ******************************************************* C C C from C * nearshore circulations under sea breeze conditions and C * wave current interactios in the surf zone C C ( C ******************************************************* C ///////////////////////////////////////////////////////////// C (, ),( i, j) DIMENSION TOP1(0:202,-1:203),TOP2(0:202,-1:203),TOP3(0:202,-1:203) DIMENSION DEP1(0:202,-1:203),DEP2(0:202,-1:203),DEP3(0:202,-1:203) DIMENSION DELTA1(-1:203), DELTA2(0:202), R1(0:202) DIMENSION DEP1A(0:202,-1:203) C PI = G = WRITE (*,*) * WRITE (*,*)+ * WRITE (*,*)+ * WRITE (*,*) * * WRITE (*,*)*Nearshore circulations under sea breeze conditions* WRITE (*,*) * and WAVE CURRENT INTERACTIONS IN THE SURF ZONE * WRITE (*,*)* * WRITE (*,*)* July * WRITE (*,*) + July * WRITE (*,*) * WRITE (*,*) WRITE (*,*) C READ (*,*) " (LG1) = " ;LGI ILG1 = 200 C READ (*,*) " (OG1) = " ;OGI IOG1=200 C READ (*,*) " = " ;XINT XINT = 10.0 C READ (*,*) " Y = " ;YINT YINT = 10.0 C ( ) DELTAX = XINT DELTAY = YINT ILG1 = ILG1 + 3 IOG1 = IOG1 + 1

6 WRITE (*,*) WRITE (*,*) depth(x,y) = WRITE (*,*) + WRITE (*,*) + WRITE (*,*) WRITE (*,*) ( ) d1 WRITE (*,*) d1 = u * x ^r(x) WRITE (*,*) u = 0.1 r=2/3 at onshore WRITE (*,*) r=1/2 at offshore WRITE (*,*) r(x) = 2/3 -(2/3-1/2)*x / xca1 C WRITE (*,*) r(x) = 4/6 -(4/6-3/6)*x / xca1 WRITE (*,*) WRITE (*,*) d2 WRITE (*,*) d2 = b exp -(x-xb)^2 / (xb/2) WRITE (*,*), b = ; xb = C READ (*,*) (m) = ;B B = 4.0 C READ (*,*) (m) = ;XB XB = WRITE (*,*) WRITE (*,*) WRITE (*,*) WRITE (*,*) WRITE (*,*) d3 = a(1-x/lb) * sin (2pi*(y-delta)/ramda) WRITE (*,*), a= ; ib = breaker length WRITE (*,*) delta = ( ) WRITE (*,*) DELTA = DELTA1 + DELTA2 WRITE (*,*) delta1 = for skewness WRITE (*,*) delta2 = for the oblique downstream WRITE (*,*) orientation of the crest of unduration WRITE (*,*) WRITE (*,*) delta1 = deltamax * sin( 2pi*(y-delta1)/ramda ) WRITE (*,*) = deltamax sin(2pi*y/ramda)*cos(2pi*delta1/ramda) WRITE (*,*) - cos(2pi*y/ramda)*sin(2pi*delta1/ramda) WRITE (*,*) C READ (*,*) " = " ;RAMDA RAMDA = C READ (*,*) " input deltamax "; DMAX DMAX = 2.0 C C C DO 100 I=1, 202 DO 110 J=-1, 202 TOP1(I,J) = 0.0 TOP2(I,J) = 0.0 TOP3(I,J) = 0.0 DEP1(I,J) = 0.0 DEP2(I,J) = 0.0 DEP3(I,J) = 0.0

7 110 CONTINUE 100 CONTINUE C WRITE (*,*) " " WRITE (*,*)" " WRITE (*,*) " " C " C " WRITE (*,*)" dep1 " XCAL = XINT * FLOAT(IOG1-1) DO 200 I=1, IOG1 C WRITE (*,*) I X = XINT * FLOAT(I) C WRITE (*,*) X C R = (4.0/6.0) - (1.0/6.0)*(X/XCAL) R = 2.0/3.0 R1A = 1.0/2.0 C WRITE (*,*) R DEPTH1 =.1 * X **R DEPTH2 = 0.1*X**R1A C WRITE (*,*) DEPTH1 DO 210 J = -1, ILG1 C WRITE (*,*) I, J DEP1(I,J) = DEPTH1 DEP1A(I,J) =DEPTH2 C WRITE (*,*) DEP1(I,J) 210 CONTINUE 200 CONTINUE C ********************************************** C C ********************************************** WRITE (*,*)" dep2 " DO 220 I=1, IOG1 DEP2(I,1) = -B* EXP( -(XINT *FLOAT(I)-XB)**2 / (XB/2)**2 ) DO 230 J=-1, ILG1 DEP2(I,J) = DEP2(I,1) 230 CONTINUE 220 CONTINUE C *********************************************** C C *********************************************** WRITE (*,*)" dep3 " DDELTA1 = DMAX DO 300 J =-1, ILG1 DELTA1(J) = 0.0 Y = DELTAY * FLOAT(J-1) 50 DOLD = DDELTA1 DNEW = DMAX * SIN( 2.0*PI*(Y-DOLD)/RAMDA ) DELTA1(J) = DNEW 300 CONTINUE

8 WRITE (*,*) dep3a C delta2 = x * tan(arpha) C READ (*,*) " " ;ALPHA ALPHA = 20.0 C READ (*,*) " " ;A A = 1.0 C READ (*,*) " " ;XLB XLB = ALPHA = *ALPHA / DO 320 I=1, IOG1 XX = XINT*FLOAT(I) DELTA2(I) = XX * TAN(ALPHA) DO 330 J=-1, ILG1 DELTA = DELTA1(J) + DELTA2(I) IF (XX.LE. XLB) THEN DEP3(I,J) = A*(1.0- XINT*FLOAT(I)/XLB)* & SIN(2.0*PI*(DELTAY*FLOAT(J-1)-DELTA) /RAMDA) ELSE DEP3(I,J) = 0.0 ENDIF 330 CONTINUE 320 CONTINUE C C C --> file dep1 C + --> file dep12 C > dep123 C DEP1 + DEP2 DO 400 I=1, IOG1 DO 410 J=-1, ILG1 TOP1(I,J) = DEP1(I,J)+DEP2(I,J) 410 CONTINUE 400 CONTINUE C DEP1 + DEP3 DO 430 I=1, IOG1 DO 440 J=-1, ILG1 TOP2(I,J) = DEP1(I,J) + DEP3(I,J) 440 CONTINUE 430 CONTINUE C DEP1 + DEP2 + DEP3 DO 450 I=1, IOG1 DO 460 J=-1, ILG1 TOP3(I,J) = TOP1(I,J) + DEP3(I,J) 460 CONTINUE 450 CONTINUE C C C

9 C C WRITE (*,*)" " WRITE (*,*) " date stock " OPEN (3, FILE=DEP1.DAT) DO 500 I=1, IOG1 DO 505 J=-1, ILG1 WRITE (3,*) FLOAT(I)*XINT, FLOAT(J)*YINT, -DEP1(I,J) 505 CONTINUE 500 CONTINUE CLOSE (3, STATUS=KEEP) OPEN (4, FILE=TOP1.DAT) DO 510 I=1, IOG1 DO 520 J=-1, ILG1 WRITE (4,*) FLOAT(I)*XINT, FLOAT(J)*YINT, -TOP1(I,J) 520 CONTINUE 510 CONTINUE CLOSE (4, STATUS=KEEP) OPEN (5, FILE=TOP2.DAT) DO 540 I=1, IOG1 DO 550 J=-1, ILG1 WRITE (5,*) FLOAT(I)*XINT, FLOAT(J)*YINT, -TOP2(I,J) 550 CONTINUE 540 CONTINUE CLOSE (5, STATUS=KEEP) OPEN (6, FILE=TOP3.DAT) DO 560 I=1, IOG1 DO 570 J=1, ILG1 WRITE (6,*) FLOAT(I)*XINT, FLOAT(J)*YINT, -TOP3(I,J) 570 CONTINUE 560 CONTINUE CLOSE (6, STATUS=KEEP) OPEN (7, FILE=R1.DAT) DO 580 I=1, IOG1 WRITE (7,*) FLOAT(I)*XINT, DEP1(I,1), DEP2(I,1) 580 CONTINUE CLOSE (7, STATUS=KEEP) OPEN (8, FILE=R2.DAT) DO 590 I=1, IOG1 DO 600 J=1, ILG1 WRITE (8,*) FLOAT(I)*XINT, FLOAT(J-1)*YINT, -DEP3(I,J) 600 CONTINUE 590 CONTINUE CLOSE (8, STATUS=KEEP) STOP END

10

11

12 C ********************************************** C C ********************************************** WRITE (*,*)" dep2 " AMPLB = DO 220 I=1, IOG1 DO 230 J=-1, ILG1 SHI = AMPLB*SIN(2.0*PI*(DELTAY*FLOAT(J-1)) /RAMDA) C SHI = 0.0 for straight longshore bar DEP2(I,J)=-B * EXP( -((XINT *FLOAT(I)-SHI)-XB)**2/(XB/2)**2) 230 CONTINUE 220 CONTINUE

13

橡3章波浪取扱.PDF

橡3章波浪取扱.PDF 2 2-2 * -* *EM EM -1- (2.1) Dally x κ ( F cos θ ) + ( F sin θ ) = ( F F s ) 2.1 y d 2-7 F F Stable wave energy flux d F (2.2) F s = E s gs = 1 8 ρgh 2 s gh = 1 8 ρg(γ h) 2 gh 2.2 η (2.3) ds xx = ρgd d

More information

Microsoft Word - 資料 docx

Microsoft Word - 資料 docx y = Asin 2πt T t t = t i i 1 n+1 i i+1 Δt t t i = Δt i 1 ( ) y i = Asin 2πt i T 29 (x, y) t ( ) x = Asin 2πmt y = Asin( 2πnt + δ ) m, n δ (x, y) m, n 30 L A x y A L x 31 plot sine curve program sine implicit

More information

情報活用資料

情報活用資料 y = Asin 2πt T t t = t i i 1 n+1 i i+1 Δt t t i = Δt i 1 ( ) y i = Asin 2πt i T 21 (x, y) t ( ) x = Asin 2πmt y = Asin( 2πnt + δ ) m, n δ (x, y) m, n 22 L A x y A L x 23 ls -l gnuplot gnuplot> plot "sine.dat"

More information

コンピュータ概論

コンピュータ概論 4.1 For Check Point 1. For 2. 4.1.1 For (For) For = To Step (Next) 4.1.1 Next 4.1.1 4.1.2 1 i 10 For Next Cells(i,1) Cells(1, 1) Cells(2, 1) Cells(10, 1) 4.1.2 50 1. 2 1 10 3. 0 360 10 sin() 4.1.2 For

More information

P-12 P-13 3 4 28 16 00 17 30 P-14 P-15 P-16 4 14 29 17 00 18 30 P-17 P-18 P-19 P-20 P-21 P-22

P-12 P-13 3 4 28 16 00 17 30 P-14 P-15 P-16 4 14 29 17 00 18 30 P-17 P-18 P-19 P-20 P-21 P-22 1 14 28 16 00 17 30 P-1 P-2 P-3 P-4 P-5 2 24 29 17 00 18 30 P-6 P-7 P-8 P-9 P-10 P-11 P-12 P-13 3 4 28 16 00 17 30 P-14 P-15 P-16 4 14 29 17 00 18 30 P-17 P-18 P-19 P-20 P-21 P-22 5 24 28 16 00 17 30 P-23

More information

1 u t = au (finite difference) u t = au Von Neumann

1 u t = au (finite difference) u t = au Von Neumann 1 u t = au 3 1.1 (finite difference)............................. 3 1.2 u t = au.................................. 3 1.3 Von Neumann............... 5 1.4 Von Neumann............... 6 1.5............................

More information

. (.8.). t + t m ü(t + t) + c u(t + t) + k u(t + t) = f(t + t) () m ü f. () c u k u t + t u Taylor t 3 u(t + t) = u(t) + t! u(t) + ( t)! = u(t) + t u(

. (.8.). t + t m ü(t + t) + c u(t + t) + k u(t + t) = f(t + t) () m ü f. () c u k u t + t u Taylor t 3 u(t + t) = u(t) + t! u(t) + ( t)! = u(t) + t u( 3 8. (.8.)............................................................................................3.............................................4 Nermark β..........................................

More information

OHP.dvi

OHP.dvi 0 7 4 0000 5.. 3. 4. 5. 0 0 00 Gauss PC 0 Gauss 3 Gauss Gauss 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 u [] u [3] u [4] u [4] P 0 = P 0 (),3,4 (,), (3,), (4,) 0,,,3,4 3 3 3 3 4 4 4 4 0 3 6 6 0 6 3 6 0 6

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 0 1 2 3 4 5 6 1964 1978 7 0.0015+0.013 8 1 π 2 2 2 1 2 2 ( r 1 + r3 ) + π ( r2 + r3 ) 2 = +1,2100 9 10 11 1.9m 3 0.64m 3 12 13 14 15 16 17 () 0.095% 0.019% 1.29% (0.348%) 0.024% 0.0048% 0.32% (0.0864%)

More information

砂浜砕波帯における流れと地形変化

砂浜砕波帯における流れと地形変化 47 * Currents and Morphological Changes in the Surf Zone on a Sandy Beach Yoshiaki KURIYAMA, Littoral Drift Division, Port and Airport Research Institute 1 1 1) ρd M M x y 2 + + = (1) t x y 2.1 M x UM

More information

PowerPoint プレゼンテーション

PowerPoint プレゼンテーション 計算機実習 Ⅰ FORTRAN 担当 2018.05.29 本日の課題 プログラムの基本ルールを理解し 以下が含まれるプログラムを作成する (1) 文法の基礎 ( フローチャートなど ) (2) 変数宣言 (3) 入出力 (4) 四則演算 (5) 組込関数 (6) 判定文 (7) リダイレクション PROGRAM MAIN INTEGER I, J, K REAL A, B, C CHARACTER

More information

n 第1章 章立ての部分は、書式(PC入門大見出し)を使います

n 第1章 章立ての部分は、書式(PC入門大見出し)を使います FORTRAN FORTRAN FORTRAN ) DO DO IF IF FORTRAN FORTRAN(FORmula TRANslator)1956 IBM FORTRAN IV FORTRAN77 Fortran90 FORTRAN77 FORTRAN FORTARN IF, DO C UNIX FORTRAN PASCAL COBOL PL/I BASIC Lisp PROLOG Lisp

More information

y = x 4 y = x 8 3 y = x 4 y = x 3. 4 f(x) = x y = f(x) 4 x =,, 3, 4, 5 5 f(x) f() = f() = 3 f(3) = 3 4 f(4) = 4 *3 S S = f() + f() + f(3) + f(4) () *4

y = x 4 y = x 8 3 y = x 4 y = x 3. 4 f(x) = x y = f(x) 4 x =,, 3, 4, 5 5 f(x) f() = f() = 3 f(3) = 3 4 f(4) = 4 *3 S S = f() + f() + f(3) + f(4) () *4 Simpson H4 BioS. Simpson 3 3 0 x. β α (β α)3 (x α)(x β)dx = () * * x * * ɛ δ y = x 4 y = x 8 3 y = x 4 y = x 3. 4 f(x) = x y = f(x) 4 x =,, 3, 4, 5 5 f(x) f() = f() = 3 f(3) = 3 4 f(4) = 4 *3 S S = f()

More information

最 新 測 量 学 ( 第 3 版 ) サンプルページ この 本 の 定 価 判 型 などは, 以 下 の URL からご 覧 いただけます. このサンプルページの 内 容 は, 第 3 版 1 刷 発 行 時 の

最 新 測 量 学 ( 第 3 版 ) サンプルページ この 本 の 定 価 判 型 などは, 以 下 の URL からご 覧 いただけます.  このサンプルページの 内 容 は, 第 3 版 1 刷 発 行 時 の 最 新 測 量 学 ( 第 3 版 ) サンプルページ この 本 の 定 価 判 型 などは, 以 下 の URL からご 覧 いただけます. http://www.morikita.co.jp/books/mid/047143 このサンプルページの 内 容 は, 第 3 版 1 刷 発 行 時 のものです. 3 10 GIS 3 1 2 GPS GPS GNSS GNSS 23 3 3 2015

More information

!!! 10 1 110 88 7 9 91 79 81 82 87 6 5 90 83 75 77 12 80 8 11 89 84 76 78 85 86 4 2 32 64 10 44 13 17 94 34 33 107 96 14 105 16 97 99 100 106 103 98 63 at 29, 66 at 58 12 16 17 25 56

More information

(2-1) x, m, 2 N(m, 2 ) x REAL*8 FUNCTION NRMDST (X, M, V) X,M,V REAL*8 x, m, 2 X X N(0,1) f(x) standard-norm.txt normdist1.f x=0, 0.31, 0.5

(2-1) x, m, 2 N(m, 2 ) x REAL*8 FUNCTION NRMDST (X, M, V) X,M,V REAL*8 x, m, 2 X X N(0,1) f(x) standard-norm.txt normdist1.f x=0, 0.31, 0.5 2007/5/14 II II agata@k.u-tokyo.a.jp 0. 1. x i x i 1 x i x i x i x x+dx f(x)dx f(x) f(x) + 0 f ( x) dx = 1 (Probability Density Funtion 2 ) (normal distribution) 3 1 2 2 ( x m) / 2σ f ( x) = e 2πσ x m

More information

untitled

untitled 3,,, 2 3.1 3.1.1,, A4 1mm 10 1, 21.06cm, 21.06cm?, 10 1,,,, i),, ),, ),, x best ± δx 1) ii), x best ), δx, e,, e =1.602176462 ± 0.000000063) 10 19 [C] 2) i) ii), 1) x best δx

More information

9 8 7 (x-1.0)*(x-1.0) *(x-1.0) (a) f(a) (b) f(a) Figure 1: f(a) a =1.0 (1) a 1.0 f(1.0)

9 8 7 (x-1.0)*(x-1.0) *(x-1.0) (a) f(a) (b) f(a) Figure 1: f(a) a =1.0 (1) a 1.0 f(1.0) E-mail: takio-kurita@aist.go.jp 1 ( ) CPU ( ) 2 1. a f(a) =(a 1.0) 2 (1) a ( ) 1(a) f(a) a (1) a f(a) a =2(a 1.0) (2) 2 0 a f(a) a =2(a 1.0) = 0 (3) 1 9 8 7 (x-1.0)*(x-1.0) 6 4 2.0*(x-1.0) 6 2 5 4 0 3-2

More information

plotwsx PLOT-WSX X11 X spp[2]./a.out PLOT-WSX PLOT IPENS -13 IPENS 0 IPENS 0 2 s p PLOT MRI 3.2 PLOT-WSX PLOT IPE- NS -13 IPENS 0 1 Continue PLO

plotwsx PLOT-WSX X11 X spp[2]./a.out PLOT-WSX PLOT IPENS -13 IPENS 0 IPENS 0 2 s p PLOT MRI 3.2 PLOT-WSX PLOT IPE- NS -13 IPENS 0 1 Continue PLO PLOT-WSX 3 1 EWS FORTRAN PLOT-WSX FORTRAN spp 2 2.1 1 1 1 1 Y (Xmax,Ymax) (0,Ymax) 2.2 PLOT-WSX NEWPEN 1 1. NEWPEN 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 2.3 PLOT-WSX PLOTS 4.1.1 PLOT-WSX outle.ps 3 X (0,0)

More information

Microsoft Word - DF-Salford解説09.doc

Microsoft Word - DF-Salford解説09.doc Digital Fortran 解説 2009/April 1. プログラム形態とデ - タ構成 最小自乗法プログラム (testlsm.for) m 組の実験データ (x i,y i ) に最も近似する直線式 (y=ax+b) を最小自乗法で決定する 入力データは組数 mと m 組の (x i,y i ) 値 出力データは直線式の係数 a,bとなる 入力データ m=4 (x i,y i ) X=1.50

More information

診療ガイドライン外来編2014(A4)/FUJGG2014‐01(大扉)

診療ガイドライン外来編2014(A4)/FUJGG2014‐01(大扉) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

More information

caim03

caim03 ImageToolBox.swift fillrect fillcolor x1,y1,x2,y2 static func fillrect(_ img:caimimage, x1:int, y1:int, x2:int, y2:int, color:caimcolor) { // let mat = img.matrix // let wid = img.width // let hgt = img.height

More information

このcsvファイルを GraphR で 表 示 する あるいはエクセルで 読 み 込 んで 処 理 できる BMP 形 式 のファイルは Windows のソフトで 表 示 できる Mercury CCD では1ピクセルが2バイトで 記 述 されているが BMP でコンパクトに 表 すため 1 ピク

このcsvファイルを GraphR で 表 示 する あるいはエクセルで 読 み 込 んで 処 理 できる BMP 形 式 のファイルは Windows のソフトで 表 示 できる Mercury CCD では1ピクセルが2バイトで 記 述 されているが BMP でコンパクトに 表 すため 1 ピク Mercuryの img ファイルの 画 像 データを ASCI 変 換 し excel や GraphR にて 表 示 する 実 行 例 (アンダーライン 部 が 入 力 ) C:\Documents and Settings\My Documents >img-bmp_csv.exe > Input file

More information

joho09.ppt

joho09.ppt s M B e E s: (+ or -) M: B: (=2) e: E: ax 2 + bx + c = 0 y = ax 2 + bx + c x a, b y +/- [a, b] a, b y (a+b) / 2 1-2 1-3 x 1 A a, b y 1. 2. a, b 3. for Loop (b-a)/ 4. y=a*x*x + b*x + c 5. y==0.0 y (y2)

More information

(Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1

(Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1 (Basic Theory of Information Processing) Fortran Fortan Fortan Fortan 1 17 Fortran Formular Tranlator Lapack Fortran FORTRAN, FORTRAN66, FORTRAN77, FORTRAN90, FORTRAN95 17.1 A Z ( ) 0 9, _, =, +, -, *,

More information

N88 BASIC 0.3 C: My Documents 0.6: 0.3: (R) (G) : enterreturn : (F) BA- SIC.bas 0.8: (V) 0.9: 0.5:

N88 BASIC 0.3 C: My Documents 0.6: 0.3: (R) (G) : enterreturn : (F) BA- SIC.bas 0.8: (V) 0.9: 0.5: BASIC 20 4 10 0 N88 Basic 1 0.0 N88 Basic..................................... 1 0.1............................................... 3 1 4 2 5 3 6 4 7 5 10 6 13 7 14 0 N88 Basic 0.0 N88 Basic 0.1: N88Basic

More information

情報活用資料-03-20150604

情報活用資料-03-20150604 cp hello.f90 echo.f90 mv echo.f90 echofile.f90 cp echofile.f90 echo.f90 7 8 9 Echo key input program echo character(80):: A read (5,*) A write (6,*) A stop end program echo chracter read 10 Echo key input

More information

Autumn 2005 1 9 13 14 16 16 DATA _null_; SET sashelp.class END=eof; FILE 'C: MyFiles class.txt'; /* */ PUT name sex age; IF eof THEN DO; FILE LOG; /* */ PUT '*** ' _n_ ' ***'; END; DATA _null_;

More information

caim03

caim03 ImageToolBox.swift fillrect fillcolor x1,y1,x2,y2 static func fillrect(_ img:caimimage, x1:int, y1:int, x2:int, y2:int, color:caimcolor) { // let mat = img.matrix // let wid = img.width // let hgt = img.height

More information

untitled

untitled Y = Y () x i c C = i + c = ( x ) x π (x) π ( x ) = Y ( ){1 + ( x )}( 1 x ) Y ( )(1 + C ) ( 1 x) x π ( x) = 0 = ( x ) R R R R Y = (Y ) CS () CS ( ) = Y ( ) 0 ( Y ) dy Y ( ) A() * S( π ), S( CS) S( π ) =

More information

Evoltion of onentration by Eler method (Dirihlet) Evoltion of onentration by Eler method (Nemann).2 t n =.4n.2 t n =.4n : t n

Evoltion of onentration by Eler method (Dirihlet) Evoltion of onentration by Eler method (Nemann).2 t n =.4n.2 t n =.4n : t n 5 t = = (, y, z) t (, y, z, t) t = κ (68) κ [, ] (, ) = ( ) A ( /2)2 ep, A =., t =.. (69) 4πκt 4κt = /2 (, t) = for ( =, ) (Dirihlet ondition) (7) = for ( =, ) (Nemann ondition) (7) (68) (, t) = ( ) (

More information

2 H23 BioS (i) data d1; input group patno t sex censor; cards;

2 H23 BioS (i) data d1; input group patno t sex censor; cards; H BioS (i) data d1; input group patno t sex censor; cards; 0 1 0 0 0 0 1 0 1 1 0 4 4 0 1 0 5 5 1 1 0 6 5 1 1 0 7 10 1 0 0 8 15 0 1 0 9 15 0 1 0 10 4 1 0 0 11 4 1 0 1 1 5 1 0 1 1 7 0 1 1 14 8 1 0 1 15 8

More information

< 中略 > 24 0 NNE 次に 指定した日時の時間降水量と気温を 観測地点の一覧表に載っているすべての地点について出力するプログラムを作成してみます 観測地点の一覧表は index.txt というファイルで与えられています このファイルを読みこむためのサブルーチンが AMD

< 中略 > 24 0 NNE 次に 指定した日時の時間降水量と気温を 観測地点の一覧表に載っているすべての地点について出力するプログラムを作成してみます 観測地点の一覧表は index.txt というファイルで与えられています このファイルを読みこむためのサブルーチンが AMD 地上気象観測データの解析 1 AMeDAS データの解析 研究を進めるにあたって データ解析用のプログラムを自分で作成する必要が生じることがあります ここでは 自分で FORTRAN または C でプログラムを作成し CD-ROM に入った気象観測データ ( 気象庁による AMeDAS の観測データ ) を読みこんで解析します データを読みこむためのサブルーチンや関数はあらかじめ作成してあります それらのサブルーチンや関数を使って自分でプログラムを書いてデータを解析していきます

More information

Microsoft Word - 資料 (テイラー級数と数値積分).docx

Microsoft Word - 資料 (テイラー級数と数値積分).docx δx δx n x=0 sin x = x x3 3 + x5 5 x7 7 +... x ak = (-mod(k,2))**(k/2) / fact_k ( ) = a n δ x n f x 0 + δ x a n = f ( n) ( x 0 ) n f ( x) = sin x n=0 58 I = b a ( ) f x dx ΔS = f ( x)h I = f a h h I = h

More information

ex01.dvi

ex01.dvi ,. 0. 0.0. C () /******************************* * $Id: ex_0_0.c,v.2 2006-04-0 3:37:00+09 naito Exp $ * * 0. 0.0 *******************************/ #include int main(int argc, char **argv) double

More information

Contents P. P. 1

Contents P. P. 1 Contents P. P. 1 P. 2 TOP MESSAGE 3 4 P. P. 5 P. 6 7 8 9 P. P. P. P. P. 10 11 12 Economy P. P. 13 14 Economy P. 1,078 1,000 966 888 800 787 716 672 600 574 546 556 500 417 373 449 470 400 315 336 218 223

More information

6 1

6 1 (c) Masaya Kasuga Shaltics 2001 6 1 1 1 2 1 1.1 USO 2 3 4 EPR 5 6 27 2 (Unfinded Superconducting Object) 3 - - - 342K (Physica C 351 (2001) pp.78-81) 4 40K 5 Einstein-Podolsky-Rosen s paradox 1/2 ( 0)

More information

t 2 2 t 2 t F ( ) p- 2 2 F 2 G F ( ) 2 2 F 2 G F ( ) 2 2 2

t 2 2 t 2 t F ( ) p- 2 2 F 2 G F ( ) 2 2 F 2 G F ( ) 2 2 2 1 2 2 0 1 2 2 2 2 2 2 2 2.1 2 2 F={f ij }, G {g ij } t f ij t g ij = 1 f ij < t g ij = 0 t p- p S 0 S p = S 0 /S t p 2 t 1 t 2 2 t 2 t 2 2 3 3 1 2 F ( ) p- 2 2 F 2 G 3 2 2 F ( ) 2 2 F 2 G 3 3 2 F ( ) 2

More information

Appendix A BASIC BASIC Beginner s All-purpose Symbolic Instruction Code FORTRAN COBOL C JAVA PASCAL (NEC N88-BASIC Windows BASIC (1) (2) ( ) BASIC BAS

Appendix A BASIC BASIC Beginner s All-purpose Symbolic Instruction Code FORTRAN COBOL C JAVA PASCAL (NEC N88-BASIC Windows BASIC (1) (2) ( ) BASIC BAS Appendix A BASIC BASIC Beginner s All-purpose Symbolic Instruction Code FORTRAN COBOL C JAVA PASCAL (NEC N88-BASIC Windows BASIC (1 (2 ( BASIC BASIC download TUTORIAL.PDF http://hp.vector.co.jp/authors/va008683/

More information

.......p...{..P01-48(TF)

.......p...{..P01-48(TF) 1 2 3 5 6 7 8 9 10 Act Plan Check Act Do Plan Check Do 11 12 13 14 INPUT OUTPUT 16 17 18 19 20 21 22 23 24 25 26 27 30 33 32 33 34 35 36 37 36 37 38 33 40 41 42 43 44 45 46 47 48 49 50 51 1. 2. 3.

More information

11042 計算機言語7回目 サポートページ:

11042 計算機言語7回目  サポートページ: 11042 7 :https://goo.gl/678wgm November 27, 2017 10/2 1(print, ) 10/16 2(2, ) 10/23 (3 ) 10/31( ),11/6 (4 ) 11/13,, 1 (5 6 ) 11/20,, 2 (5 6 ) 11/27 (7 12/4 (9 ) 12/11 1 (10 ) 12/18 2 (10 ) 12/25 3 (11

More information

untitled

untitled CONTENTS 17 8 911 1213 1415 16 17 18 19 20 21 2012.7.20JULY 1 2012.7.20JULY 2 3 4 2012.7.20JULY 5 6 2012.7.20JULY 7 8 2012.7.20JULY 9 10 2012.7.20JULY 11 2012.7.20JULY 12 13 2012.7.20JULY 14 15 16 2012.7.20JULY

More information

main.dvi

main.dvi 1 F77 5 hmogi-2008f@kiban.civil.saitama-u.ac.jp 2013/5/13 1 2 f77... f77.exe f77.exe CDROM (CDROM D D: setupond E E: setupone 5 C:work\T66160\20130422>f77 menseki.f -o menseki f77(.exe) f77 f77(.exe) C:work\T66160\20130422>set

More information

+ + + + n S (n) = + + + + n S (n) S (n) S 0 (n) S (n) 6 4 S (n) S (n) 7 S (n) S 4 (n) 8 6 S k (n) 0 7 (k + )S k (n) 8 S 6 (n), S 7 (n), S 8 (n), S 9 (

+ + + + n S (n) = + + + + n S (n) S (n) S 0 (n) S (n) 6 4 S (n) S (n) 7 S (n) S 4 (n) 8 6 S k (n) 0 7 (k + )S k (n) 8 S 6 (n), S 7 (n), S 8 (n), S 9 ( k k + k + k + + n k 006.7. + + + + n S (n) = + + + + n S (n) S (n) S 0 (n) S (n) 6 4 S (n) S (n) 7 S (n) S 4 (n) 8 6 S k (n) 0 7 (k + )S k (n) 8 S 6 (n), S 7 (n), S 8 (n), S 9 (n), S 0 (n) 9 S (n) S 4

More information

Microsoft Word - 03-数値計算の基礎.docx

Microsoft Word - 03-数値計算の基礎.docx δx f x 0 + δ x n=0 a n = f ( n) ( x 0 ) n δx n f x x=0 sin x = x x3 3 + x5 5 x7 7 +... x ( ) = a n δ x n ( ) = sin x ak = (-mod(k,2))**(k/2) / fact_k 10 11 I = f x dx a ΔS = f ( x)h I = f a h I = h b (

More information

日本糖尿病学会誌第58巻第1号

日本糖尿病学会誌第58巻第1号 α β β β β β β α α β α β α l l α l μ l β l α β β Wfs1 β β l l l l μ l l μ μ l μ l Δ l μ μ l μ l l ll l l l l l l l l μ l l l l μ μ l l l l μ l l l l l l l l l l μ l l l μ l μ l l l l l l l l l μ l l l l

More information

A/B (2010/10/08) Ver kurino/2010/soft/soft.html A/B

A/B (2010/10/08) Ver kurino/2010/soft/soft.html A/B A/B (2010/10/08) Ver. 1.0 kurino@math.cst.nihon-u.ac.jp http://edu-gw2.math.cst.nihon-u.ac.jp/ kurino/2010/soft/soft.html 2010 10 8 A/B 1 2010 10 8 2 1 1 1.1 OHP.................................... 1 1.2.......................................

More information

0 1-4. 1-5. (1) + b = b +, (2) b = b, (3) + 0 =, (4) 1 =, (5) ( + b) + c = + (b + c), (6) ( b) c = (b c), (7) (b + c) = b + c, (8) ( + b)c = c + bc (9

0 1-4. 1-5. (1) + b = b +, (2) b = b, (3) + 0 =, (4) 1 =, (5) ( + b) + c = + (b + c), (6) ( b) c = (b c), (7) (b + c) = b + c, (8) ( + b)c = c + bc (9 1-1. 1, 2, 3, 4, 5, 6, 7,, 100,, 1000, n, m m m n n 0 n, m m n 1-2. 0 m n m n 0 2 = 1.41421356 π = 3.141516 1-3. 1 0 1-4. 1-5. (1) + b = b +, (2) b = b, (3) + 0 =, (4) 1 =, (5) ( + b) + c = + (b + c),

More information

情報処理概論(第二日目)

情報処理概論(第二日目) 情報処理概論 工学部物質科学工学科応用化学コース機能物質化学クラス 第 8 回 2005 年 6 月 9 日 前回の演習の解答例 多項式の計算 ( 前半 ): program poly implicit none integer, parameter :: number = 5 real(8), dimension(0:number) :: a real(8) :: x, total integer

More information

( ) 1 Windows HTML ( ) ( ) ( ) WWW 10 ( )

( ) 1 Windows HTML ( ) ( ) ( ) WWW 10 ( ) ( ) 1 Windows HTML ( ) ( ) ( ) 1. 2. 3. 4. WWW 10 ( ) 2 1. 2. 1 3. ( ) 4. 5. 3 Windows 2 7 8 MS Word MS Excel 1. MS Word 600 2. MS Excel 1 34 2 83 3 23 4 70 5 100 6 45 7 58 8 29 9 68 10 72 11 37 12 12

More information

平成19年度

平成19年度 1 2 3 4 H 3 H CC N + 3 O H 3 C O CO CH 3 CH O CO O CH2 CH 3 P O O 5 H H H CHOH H H H N + CHOH CHOH N + CH CH COO- CHOH CH CHOH 6 1) 7 2 ) 8 3 ) 4 ) 9 10 11 12 13 14 15 16 17 18 19 20 A A 0 21 ) exp( )

More information

c-all.dvi

c-all.dvi III(994) (994) from PSL (9947) & (9922) c (99,992,994,996) () () 2 3 4 (2) 2 Euler 22 23 Euler 24 (3) 3 32 33 34 35 Poisson (4) 4 (5) 5 52 ( ) 2 Turbo 2 d 2 y=dx 2 = y y = a sin x + b cos x x = y = Fortran

More information

untitled

untitled 1 4 4 6 8 10 30 13 14 16 16 17 18 19 19 96 21 23 24 3 27 27 4 27 128 24 4 1 50 by ( 30 30 200 30 30 24 4 TOP 10 2012 8 22 3 1 7 1,000 100 30 26 3 140 21 60 98 88,000 96 3 5 29 300 21 21 11 21

More information