C:/大宮司先生原稿/数値積分数値微分.dvi

Size: px
Start display at page:

Download "C:/大宮司先生原稿/数値積分数値微分.dvi"

Transcription

1 6 u(x) 6. ::: x ; x x x ::: ::: u ; u u u ::: (dscrete) u = u(x ) ( = :::) x ;x = h =cnst: x <x<x x u(x) u(x) = x ;x h u x;x h u (x;x )(x;x ) x x u (x) = u ;u h Z x u(x)dx = x;x x (x;x ;x ) ; u u(x) ; (x;x ) x x = x x u = u = u ;u O(h) h Z x u(x)dx = x h(u u )O(h ) [ ]=[ ][ ] (trunctn errr) 6.:

2 x Tylr u(x) =u (x;x )u! (x;x )! (x;x ) u (x) =u (x;x )! (x;x ) Z x u(x)dx =(x;x )u x! (x;x ) u! (x;x ) 4! (x;x ) 4 u = u hu! (h)! (h) ( = :::) u u Tylr x;x x;x h O(h ), O(h ) u ; u u x ; <x<x x u(x) u(x) = h n (x;x )(x;x )u ; ;(x;x )(x;x ; )u (x;x ;)(x;x )u! (x;x ;)(x;x )(x;x ) u (x) = h (x;x ;x )u ; ;(x;x ; ;x )u (x;x ; ;x )u Z x! (x;x )(x;x )(x;x )(x;x ; )(x;x ; )(x;x ) u(x)dx = (x;x ;) (x;x ;x x; h )(u ; u );4(x;x ;h)u ; (x;x ;)(x;x ) h 4! (x;x ;) (x;x ) x u = h (;u ;u )O(h ) u = h (u ;;4u u )O(h ) Z ; u(x)dx = h(u ;4u u )O(h 4 ) 4 4 u

3 6. (numercl ntegrtn, numercl qudrture) u(x) dx u(x) u ( = ::: n) (dscrete) [ b] u(x) Newtn-Ctes Newtn-Ctes Guss-Legendre Guss Legendre u 6.. u(x) u u ::: u n u = u(x ) x = <x < <x n = b n n p n (x) n p n (x) = x x n x n 6. n u u ::: u n A A ::: A n p n (x) dx = A p n (x )A p n (x ) A n p n (x n ) (6.) ::: n n ::: n 6.: (dscntnutes) : : :

4 4 n dx =A A A n xdx=a x A x A n x n (6.) x n dx =A x n A x n A n x n n dx = b; xdx= b ; x n dx = bn ; n n A A A n (6.) x x x n :::::::::::::::::: x n x n x n n C A A A. A n C A = b; (b ; )=. (b n ; n )=(n) C A (6.) Vndermnde x x ::: x n A A ::: A n (6.) A A ::: A n u(x) dx = A u A u A n u n (6.4) n p n (x) b 6. u u ::: u n n (6.4) u = p n (x ) u = p n (x ) u n = p n (x n ) p n (x )= (x ) n (x )( = ::: n) (6.) u(x) dx = A p n (x )A p n (x )A n p n (x n ) = (A A A n ) (A x A x A n x n ) n (A x n A x n A n x n n ) = dx xdx n x n dx = p n (x)dx n = 4 A

5 5 SUBROUTINE NUMINT(x,u,f,s) DIMENSION x(:f),u(:f),a(:4),x(:4) REAL(8) v(5,6),w(6) w=x(f)-x() FORALL(=:f)x()=(x()-x())/w DO =,5 v(,6)=./float() IF(==)THEN FORALL(j=:5)v(,j)=. ELSE FORALL(j=:5)v(,j)=x(j-)*v(-,j) ENDIF ENDDO det=. e=. CALL GAUSS(v,w,5,6,det,e) FORALL(=:f)A()=w*v(,6) s=. DO =,f s=sa()*u() ENDDO ENDSUBROUTINE! Slutn f smultneus lner equtns by Gussn elmntn SUBROUTINE GAUSS(,w,n,n,det,e) REAL(8) (n,n),w(n) det=. e=max(e,.e-5) cycle_: DO k=,n l=k w=(k,k) IF(ABS(w)>e) GOTO l=k l=l IF(l>n)STOP 5555 w=(l,k) IF(ABS(w)<=e) GOTO FORALL(j=k:n) w(j)=(k,j) (k,j)=(l,j) (l,j)=w(j) ENDFORALL det=-det det=w*det FORALL(j=k:n)(k,j)=(k,j)/w IF(k==n)CYCLE cycle_ FORALL(=k:n,j=k:n)(,j)=(,j)-(,k)*(k,j) ENDDO cycle_ DO k=n,,- FORALL(=:k-)(,n)=(,n)-(,k)*(k,n) ENDDO END 6. = x! b = x 4! x s GAUSS Guss,v! n=n w det e

6 6 6.. Newtn-Ctes A A ::: A n A = h ( = ::: n) h n / / 4 Smpsn / /8 Smpsn /8 4 / / / (6/) ( ) (Weddle ) Newtn-Ctes [ b] n N h =(b ; )=N x = x h ( = ::: N) x = x N = b u = u(x ) n = ( trpezdl lw) Z xn ; u(x) dx = h x u u u u N; u N n = Smpsn / N =7 Smpsn /8 Z xn x u(x) dx = h(u 4u u 4u u 4 ) 8 h(u 4u 5 u 6 u 7 ) Newtn-Ctes (trunctn errr) Smpsn / u(x) Tylr u(x) =u (x;x )u! (x;x )! (x;x ) 4! (x;x ) 4 u (4) [x x ] Z x u(x) dx =hu h u 4 x h h4 4 5 h5 u (4) [x x ] Smpsn / Tylr h(u 4u u )=hu h u 4 h h4 5 8 h5 u (4) Z x u(x) dx = x h(u 4u u ) ; 9 h5 u (4) Smpsn / Smpsn / N =(b;)=h e t = ; b; 8 h4 u (4)

7 7 Smpsn / 4 O(h 4 ) Newtn-Ctes n 4 5 ;e t b; h b; 8 h4 u (4) b; 8 h4 u (4) (b;) 945 h6 u (6) 55(b;) 96 h6 u (6) Tylr Tylr n n Newtn-Ctes A n u(x) Newtn-Ctes N Smpsn / Z xn x u(x) dx = h(u 4u u 4u u 4 u N; 4u N; ) Smpsn / Z xn x u(x) dx = h(u 4u u 4u 4 u 5 u N; 4u N ) Z xn x u(x) dx = h(u u u u u N; ) n = n n 6.. Guss-Legendre [x x ] u u u(x = ) h [x x ] Smpsn / u u u Guss-Legendre Guss n Legendre n n (smth) (sucently smth)

8 8 Legendre [; ] [ b] ~x = b b; x [; ] ~u(~x) d~x = b; Z ; u(x) dx n u(x ) ( = ::: n) n n p n (x) p n (x) =p n (x)(x;x )(x;x ) (x;x n ) q n (x) (6.5) q n (x) n [; ] Z Z Z p n (x) dx = p n (x) dx (x;x )(x;x ) (x;x n ) q n (x) dx (6.6) ; ; ; x x ::: x n n Legendre P n (x) (x;x )(x;x ) (x;x n )=P n (x) (6.6) (6.4) (6.5) p n (x k )=p n (x k ) (k = ::: n) Z ; p n (x) dx = nx = A p n (x ) (6.6) Legendre Z ; x r P m (x) dx = (r = ::: m;) Guss-Legendre Z ; p n (x) dx = nx = A p n (x ) (6.7) x n Legendre A (6.) x A n x A Guss-Legendre [ b] e t = (b ; )n n (n )! (n )(n )(n ) u (n) (n )! (6.8)

9 9 I = Z dx x =:5 5 7 Smpsn / 4 Smpsn /8 n = Guss-Legendre trpezdl (n=) 7.55 :5 Smpsn / (n=).546 : : :9 Smpsn /8 (n=) 4.59 : ;:5 Newtn-Ctes (n=4) 5.54 :4 Newtn-Ctes (n=6) ;:444 Guss-Legendre (n=) ;: Guss-Legendre (n=) ;: Guss-Legendre I = Z ; dx (:5:5x) = :55555 (:;:7746) :88888 :55555 = :49988 : (::7746) x Smpsn /8 7 Newtn-Ctes (n=6) Guss-Legendre 6. u = u(x ) u(x) u(x) Lgrnge 6.. Lgrnge Lgrnge Lgrnge

10 Lgrnge u(x) = X k= C k() u ;k C () =; (h ;)(h h ;) h (h h )(h h h ) C () = (h )(h ;)(h h ;) h h (h h ) C () = (h )(h h ;) C (h )(h ;) () =; (h h )h h (h h h )(h h )h = x;x h = x ;x ; h = x ;x h = x ;x u(x) =(; )u u (; ; )( h ; ) h ( h )( h h ) u (; ; )(; h h ; ) h ( h ) (; )( h ; h ) (; )( h ( h u ; h ) ) h ( h )( h h ) u = =h h = h =h h = h =h [x x ] Z x x u(x) dx = h (u u ) h ; h h h ( h ) u h ; h h ( h ) u ; h ; h h ( h )( h h ) u ; h h ( h )( h h ) u u (6.9) (6.) [x x ] Z x u(x) dx = x h (u u ) h h 4h ; h ( h )( h h ) u em 4 h h h ( h h ) u ; h h h h ( h ) u h h ( h h )( h h ) u (6.b) = =h h = h =h h = h =h h = h = h = h Z x u(x) dx = h x 4 (;u ;u u ;u ) (6.) Z x u(x) dx = h x 4 (9u 9u ;5u u ) (6.b) Z xn n 8 u(x) dx = h x 4 u 4 u 4 u 5 4 u u 4 u 5 u N;4 5 4 u N; 4 u N; 4 u N; 8 4 u N (6.c)

11 [x x ] S(x) = (;) (;) (6.) S(x ) S(x ) S (x ) S (x ) S(x) = x (;) (;) u ; x 6 6 u u ; x 6 u (6.) x = x ;x =(x;x )=x u [x x ] Z x Z u(x)dx x S()d x = x 4 (u ) x nu ; x 6 u u ; x 6 u = x (u u ); x 4 (u ) (6.4) 4.. x 5 u(4) =5! S(x) = snh ; (;)x p snh(x p ) (;) u(x) = p ; snh snh (;)x x p u p (;) u ; p u ; p snh(x p ) (6.5) p x Z x u(x)dx csh(x x x p ); (u u ) p snh(x p ) ; x ( p ) (6.6) 4.. (6.) u(x) = 6 x u (;) 6 x u (;) (6.7)

12 [x x ] j 4.4. Z x u(x)dx x ) x ( ) (6.8) x 4 (u u(x) (;)u u R; R ; R; ; u (6.9) =(x;x )=h u = = u ;u u = u = ;u ;= R = u = u [x x ] Z x u(x) dx h ; ; u x u R = lg R ; ; u R; R; = h hu u R; lg R ; R u (6.) R; [x ; x ] [x x ] Z x dx h hu ;u x ;u(x) R; R lg R R; u R; Z x h u(x) dx h ; u u R x R; lg R ; R; u R; (6.b) (6.c) Rmberg u(x) T () b; = h T () = h fu(x )u(x h )g = I (b;)h u (x )O(h ) (6.) I O(h ) T () ( h =(b;)=) T () = h fu(x )u(x h )u(x h )g = I (b;)h u (x )O(h ) (6.) T () T () h h h! T () = () (4T ;T () )=I O(h ) (6.)

13 (6.) (6.) O(h ) T () =4 ( h =(b;)=4 =h =) T () = h fu(x )u(x h )u(x h )u(x h )u(x 4h )g (6.4) = I (b;)h u (x )O(h ) T () T () h! T () T () = () (4T ;T () )=I O(h ) (6.5) T () T () h! T () T () = T () ;T () ; = I O(h 4 ) (6.6) k k T (k) k; ( h k =(b;)= k; = h k; =) T (k) = h kfu(x )u(x h k )u(x h k )u(x k; h k )g (6.7) = I (b;)h k u (x )O(h ) T (k;) T (k) h! T (k) T (k) = T (k) ;T (k;) ; = I O(h ) (6.8) T (k;) T (k) h! T (k) T (k) k T (k) m = m T (k) (k;) m; ;T m; m = I O(h m ) ( <mk) (6.9) ; m= k; k T () T () T () T () T () T ()... k; T (k;) T (k;) T (k;) T (k;) k; k. T (k) T (k) T (k) T (k) k; T (k) k jt (k) k ;T (k) (k) j < T k; k

14 4 6. (nte-derence frmuls) Tylr n k n;k Lgrnge u(x) =(;)u u (;)h h = x ;x =(x;x )=h d=dx =(=h)d=d u (x) = h (;u u ) ; ; hu (6.) O(h) == O(h ) = x = x u = u = h u = O(h) u = h u ;= O(h) (6.) h = x ;x ; 6.. u(x) =u e u=! (;) e u! (m; )(;)h eu ;= = u ;u ; e u = m u= e = u ; ;u u= e = u ;u m m ; ( u= e ; u;= e ) e u = m ( u= e ; u= e ) h = x = m ; = x ;= =h m = x = =h u (x) = h n eu=! (;) e u! (m ; )(;)(;) h (6.) (x) = h e u (m; ;)h (6.b)

15 5 O(h ) O(h) = u = h eu=; e u O(h ) (6.) = h e u O(h) (6.b) e u e u u = h eu= ; e u O(h ) (6.c) = h e u O(h) (6.d) - ; m ; - - m u(x) =u e u= (;) n e u (m; ) e u = 4! (m; )(;)(;;m )h 4 u (4) e u = = m ; m ( e u ; e u u=@x = ;m ; (;)(m; ) e u = u(( m )h) =u u (m )m (m ; m ) m ; m m m = == Lgrnge 4 u ; u u u x ; x x x 4 (m ; )(;)(;;m ) Tylr h 4 u (4) 4 (=4!)h 4 u (4) 4

16 6 u (x) = hn eu=! (;) e u! C() e u = 4! (m ; )(;)C()(;;m ) h u (4) (6.4) (x) = h n e u (m; ;) e u = 4! C()(m ; ;)(;;m ) h u (4) (6.4b) (x) = h e u = 4 (4m; ;;m )hu (4) (6.4c) C() =(;)(m ; )(;) O(h ) O(h ) = u = hn eu= ;! e u ;! m; e u = O(h ) (6.5) = h n e u (m; ;) e u = O(h ) (6.5b) - ; m ; - - m - m - m - Lgrnge u(x) =u e u= (;) e u ;;m e u = O(h 4 ) eu 5= = m (u ;u ) e u = e u = = m ( e u ; e u ) m m ( e u5= ; e u= ) m =(x ;x )=h m =m m u (x) = hn eu=! (;) e u! C () e u = O(h ) (6.6) (x) = h n e u (;;m ) e u = O(h ) (6.6b) (x) = h e u = O(h) (6.6c) C () =(;)(;)(;;m ) = u = hn eu= ;! e u! (m ) e u = O(h ) (6.7) = h n e u ; (m ) e u = O(h ) (6.7b)

17 u(x) =u u= e h (;) e u n (m; ) e u = 4 (;;m ) e 4 u 5! (m;; m ; )(m ; )(;)(;;m )h 5 u (5) eu ;= = m ;; (u ;;u ; ) e u ; = m ;; m ; ( e u;= ; e u;= ) e u ;= = m ; ( e u ; e u ; ) e 4 u = 4 m ( e u = ; e u ;= ) m ;; =(x ; ;x ; )=h m ; = m ;; m ; m = m ; m u (x) = hn eu=! (;) e u! C() e u = ; (m ; )(;)C()(;;m ) e 4 u O(h 4 ) 4! (x) = n e u h (m; ;) e u = ; C()(m ; ;)(;;m ) e 4 u O(h ) (6.8) (6.8b) (x) = h n e u = 4 (4m; ;;m ) e 4 u O(h ) (6.8c) u (4) (x) = h 4 e 4 u O(h) (6.8d) = u = h n eu=;! e u ;! m; e u = 4! m; (m ) e 4 u O(h 4 ) = h n e u (m; ;) e u = ; ; m ; (m ; ;)(m ) e 4 u O(h ) (6.9) (6.9b) = h e u = 4 (m; ;;m ) e 4 u O(h ) (6.9c) 4 e 4 u e 4 u u (x) = hn eu=! (;) e u! C() e u = ; (m ; )(;)C()(;;m ) e 4 u O(h 4 ) 4! (x) = n e u h (m; ;) e u = ; C()(m ; ;)(;;m ) e 4 u O(h ) (6.4) (6.4b) (x) = h n e u = 4 (4m; ;;m ) e 4 u O(h ) (6.4c) u (4) (x) = h 4 e 4 u O(h) (6.4d)

18 8 = u = hn eu= ;! e u ;! m; e u = 4! m; (m ) e 4 u = h n e u (m; ;) e u = ; O(h 4 ) ; m ; (m ; ;)(m ) e 4 u O(h ) (6.4) (6.4b) = h e u = 4 (m; ;;m ) e 4 u O(h ) (6.4c) eu 5= = m (u ;u ) e u = m m ( e u5= ; e u= ) e u = = m ( e u ; e u ) e 4 u = 4^m ( e u = ; e u = ) m =(x ;x )=h ^m = m ; m - ; ; 4 m ;; - m ; - - m - m -m - m; - m; - m - m - m - ^m - m - 4 Lgrnge u(x) =u u= e n (;) e u ;;m e u ;m = 4 eu 7= = m (u 4;u ) e u = m m ( e u7= ; e u5= ) e u 5= = m ( e u ; e u ) e 4 u = 4 m ( e u 5= ; e u = ) e 4 u O(h 5 ) m =(x 4 ;x )=h m =m u (x) = hn eu=! (;) e u! C () e u = ; (;)(;;m )C ()(;m ) e 4 u O(h 4 ) 4! (x) = n e u h (;;m ) e u = ; C ()(;;m )(;m ) e 4 u O(h ) (6.4) (6.4b) (x) = h n e u = 4 (4;;m ;m ) e 4 u O(h ) (6.4c) u (4) (x) = h 4 e 4 u O(h) (6.4d)

19 9 = O(h 4 ) u = hn eu= ; e u!! (m ) e u = ; 4! (m )m e 4 u = n e u h ; (m ) e u = ; m (m )m e 4 u O(h ) (6.4) (6.4b) = h n e u = ; 4 (m m ) e 4 u O(h ) (6.4c) Lgrnge (6.) (6.b) u j u = = n h; h ; h h (u ;u ) h (u ;u ; ) h ; n (u ;u ); (u ;u ; ) h ; h h h ; O(h ) O(h) h = x ;x u = n h (h h ) (u ;u ); h h h h (u ;u ) h (h h )h h h h h h h h h h h (u ;u ) h = h h h O(h ) 6..5 m ;; = m ; = m = = u = u ;u ; h O(h) u = u ;u h O(h) u = h (u ;;4u ; u )O(h ) u = h (;u ;u )O(h ) u = h (;u 4u ;u )O(h ) u = 6h (;u ;;u 6u ;u )O(h ) u = 6h (;u 8u ;9u u )O(h ) u = h (u ;;8u ; 8u ;u )O(h 4 ) u = h (;u ;;u 8u ;6u u )O(h 4 ) u = h (;5u 48u ;6u 6u ;u 4 )O(h 4 )

20 = h (u ;;u ; u )O(h) = h (u ;;u u )O(h ) = h (u ;u u )O(h) = h (u ; 5u 4u ;u )O(h ) = h (;u ;6u ; ;u 6u ;u )O(h 4 ) = h (u ;;u 6u 4u ; u )O(h ) = h (5u ;4u 4u ;56u u 4 )O(h ) = h (;u ;u ;u u )O(h) = h (;u u ;u u )O(h) = h (;u ;u ; ;u u )O(h ) = h (;u ;u ;u 6u ;u )O(h ) = h (;5u 8u ;4u 4u ;u 4 )O(h ) u (4) = h 4 (u ;;4u ; 6u ;4u u )O(h ) u (4) u (4) = h 4 (u ;;4u 6u ;4u u )O(h) = h 4 (u ;4u 6u ;4u u 4 )O(h) j n O(h j;n ) (6.b) = 6.4 Newtn

21 6.4. Newtn Lgrnge Strlng 4 u(x) =u u =u ;=! u! ( ;) u = u ;= 4! ( ;) 4 u O(h 5 ) (;= =) u (x) = hn u= u ;= u 6 ( ;) u = u ;= ( ;) 4 u O(h 4 ) (6.44) (x) = n u h u = u ;= (6 ;) 4 u O(h ) (6.44b) (x) = n u = u ;= 4 u h O(h ) (6.44c) u (4) (x) = h 4 4 u O(h) = u = hn u= u ;= ; 6 u = u ;= O(h 4 ) (6.44d) (6.45) = u h ; 4 u O(h ) (6.45b) = u = u ;= O(h ) h x = 4 Bessel (6.45c) u(x) = u u ; ; u=! (;) u u (;); ;! u = 4! ( ;)(;) 4 u 4 u O(h 5 ) ( ) == u = = u =; 4 u = O(h 4 ) = u u = ; 5 4 = = u = O(h ) = 4 u 4 u O(h) = u (4) 4 u 4 u O(h ) (6.46) (6.46b) (6.46c) (6.46d)

22 (6.) u (x) = x ;(;) ;u x 6 u u ; x 6 u (6.47) (x) =(;) (6.47b) x = x ;x =(x ; x )=x j (6.5) u (x) = (x) = p snh x p ;csh ; (;)x p u ; u ; x p u ; x p csh(x p ) (6.48) snh x p snh ; (;)x p u snh(x p ) (6.48b) p 4.. j (6.7) u (x) = x ;(;) (; ) (6.49) x (x) =(;) (6.49b) [x x ] j FORTRAN (6.9) u (x) = hu = lg R x R; R; R ; u (6.5) (x) = (lg R) x (R;) R u (6.5b) R = u = u <! ; = ; =

1 Edward Waring Lagrange n {(x i, y i )} n i=1 x i p i p i (x j ) = δ ij P (x) = p i p i (x) = n y i p i (x) (1) i=1 n j=1 j i x x j x i x j (2) Runge

1 Edward Waring Lagrange n {(x i, y i )} n i=1 x i p i p i (x j ) = δ ij P (x) = p i p i (x) = n y i p i (x) (1) i=1 n j=1 j i x x j x i x j (2) Runge Edwrd Wring Lgrnge n {(x i, y i )} n x i p i p i (x j ) = δ ij P (x) = p i p i (x) = n y i p i (x) () n j= x x j x i x j (2) Runge [] [2] = ξ 0 < ξ < ξ n = b [, b] [ξ i, ξ i ] y = /( + 25x 2 ) 5 2,, 0,,

More information

1. 1 BASIC PC BASIC BASIC BASIC Fortran WS PC (1.3) 1 + x 1 x = x = (1.1) 1 + x = (1.2) 1 + x 1 = (1.

1. 1 BASIC PC BASIC BASIC BASIC Fortran WS PC (1.3) 1 + x 1 x = x = (1.1) 1 + x = (1.2) 1 + x 1 = (1. Section Title Pages Id 1 3 7239 2 4 7239 3 10 7239 4 8 7244 5 13 7276 6 14 7338 7 8 7338 8 7 7445 9 11 7580 10 10 7590 11 8 7580 12 6 7395 13 z 11 7746 14 13 7753 15 7 7859 16 8 7942 17 8 Id URL http://km.int.oyo.co.jp/showdocumentdetailspage.jsp?documentid=

More information

表1_表4

表1_表4 HN- 95 HN- 93 HN- 90 HN- 87 HN- 85 HN- 82 HN- 80 HN- 77 HN- 75 HN- 72 HN- 70 HN- 67 HN- 65 HN- 60 HN- 55 HN- 50 HN- 45 HN- 40 HN- 35 HN- 30 HN- 25 HN- 20 HN- 15 HN- 10 H02-80H H02-80L H02-70T H02-60H H05-60F

More information

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2 II No.1 [n/] [1]H n x) H n x) = 1) r n! r!n r)! x)n r r= []H n x) n,, H n x) = 1) n H n x) [3] H n x) = 1) n dn x e dx n e x [4] H n+1 x) = xh n x) nh n 1 x) ) d dx x H n x) = H n+1 x) d dx H nx) = nh

More information

3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t + u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., ν. t +

3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t + u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., ν. t + B: 2016 12 2, 9, 16, 2017 1 6 1,.,,,,.,.,,,., 1,. 1. :, ν. 2. : t = ν 2 u x 2, (1), c. t + c x = 0, (2). e-mail: iwayama@kobe-u.ac.jp,. 1 3. :, c, ν. 4. Burgers : t + c x = ν 2 u x 2, (3), ν. 5. : t +

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 7 (ordinary dierential equation) [a x b] u = f(x u) (7.) u(x )=u (7.2) u(x) x = a f(x u) 2 4 2 (Legendre ) (Bessel ) (7.) 7. (7.) u x u (starting values) Runge-Kutta 7.. x u(x) u(x) =u +(x ; x )u + 2!

More information

N N 1,, N 2 N N N N N 1,, N 2 N N N N N 1,, N 2 N N N 8 1 6 3 5 7 4 9 2 1 12 13 8 15 6 3 10 4 9 16 5 14 7 2 11 7 11 23 5 19 3 20 9 12 21 14 22 1 18 10 16 8 15 24 2 25 4 17 6 13 8 1 6 3 5 7 4 9 2 1 12 13

More information

3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,.,

3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2, (4), c. 2 u t 2 = c2 2 u x 2, (5) (1) (4), (1 Navier Stokes,., B:,, 2017 12 1, 8, 15, 22 1,.,,,,.,.,,,., 1,. 1. :, ν. 2. : u t = ν 2 u x 2, (1), c. u t + c u x = 0, (2), ( ). 1 3. :, c, ν. 4. Burgers : u t + c u x = ν 2 u x 2, (3), ν. 5. : u t + u u x = ν 2 u x 2,

More information

6. Euler x

6. Euler x ...............................................................................3......................................... 4.4................................... 5.5......................................

More information

第10章 アイソパラメトリック要素

第10章 アイソパラメトリック要素 June 5, 2019 1 / 26 10.1 ( ) 2 / 26 10.2 8 2 3 4 3 4 6 10.1 4 2 3 4 3 (a) 4 (b) 2 3 (c) 2 4 10.1: 3 / 26 8.3 3 5.1 4 10.4 Gauss 10.1 Ω i 2 3 4 Ξ 3 4 6 Ξ ( ) Ξ 5.1 Gauss ˆx : Ξ Ω i ˆx h u 4 / 26 10.2.1

More information

ii

ii ii iii 1 1 1.1..................................... 1 1.2................................... 3 1.3........................... 4 2 9 2.1.................................. 9 2.2...............................

More information

untitled

untitled 40 4 4.3 I (1) I f (x) C [x 0 x 1 ] f (x) f 1 (x) (3.18) f (x) ; f 1 (x) = 1! f 00 ((x))(x ; x 0 )(x ; x 1 ) (x 0 (x) x 1 ) (4.8) (3.18) x (x) x 0 x 1 Z x1 f (x) dx ; Z x1 f 1 (x) dx = Z x1 = 1 1 Z x1

More information

応用数学III-4.ppt

応用数学III-4.ppt III f x ( ) = 1 f x ( ) = P( X = x) = f ( x) = P( X = x) =! x ( ) b! a, X! U a,b f ( x) =! " e #!x, X! Ex (!) n! ( n! x)!x! " x 1! " x! e"!, X! Po! ( ) n! x, X! B( n;" ) ( ) ! xf ( x) = = n n!! ( n

More information

2009 IA I 22, 23, 24, 25, 26, a h f(x) x x a h

2009 IA I 22, 23, 24, 25, 26, a h f(x) x x a h 009 IA I, 3, 4, 5, 6, 7 7 7 4 5 h fx) x x h 4 5 4 5 1 3 1.1........................... 3 1........................... 4 1.3..................................... 6 1.4.............................. 8 1.4.1..............................

More information

x y x-y σ x + τ xy + X σ y B = + τ xy + Y B = S x = σ x l + τ xy m S y = σ y m + τ xy l σ x σ y τ xy X B Y B S x S y l m δu δv [ ( σx δu + τ )

x y x-y σ x + τ xy + X σ y B = + τ xy + Y B = S x = σ x l + τ xy m S y = σ y m + τ xy l σ x σ y τ xy X B Y B S x S y l m δu δv [ ( σx δu + τ ) 1 8 6 No-tension 1. 1 1.1................................ 1 1............................................ 5.1 - [B].................................. 5................................. 6.3..........................................

More information

ケミカルエンジニアのためのExcelを用いた化学工学計算法

ケミカルエンジニアのためのExcelを用いた化学工学計算法 VBA 7.1 f ()= ( f ( )) y = f ()(1) y = f ( )( ) + f ( ) (1) 1 1 1 = f ( ) f ( ) (2) 1 n = = y = f() y = () 1 n+1 = n (f( n )f( n )) 1.7 1 2 log()2 145 7.2 f ( ) 1-1 -1 f ( ) (3) (4) 34 5 1 1) ( ) ( ) '(

More information

Gauss

Gauss 15 1 LU LDL T 6 : 1g00p013-5 1 6 1.1....................................... 7 1.2.................................. 8 1.3.................................. 8 2 Gauss 9 2.1.....................................

More information

−g”U›ß™ö‡Æ…X…y…N…g…‰

−g”U›ß™ö‡Æ…X…y…N…g…‰ 1 / 74 ( ) 2019 3 8 URL: http://www.math.kyoto-u.ac.jp/ ichiro/ 2 / 74 Contents 1 Pearson 2 3 Doob h- 4 (I) 5 (II) 6 (III-1) - 7 (III-2-a) 8 (III-2-b) - 9 (III-3) Pearson 3 / 74 Pearson Definition 1 ρ

More information

sim98-8.dvi

sim98-8.dvi 8 12 12.1 12.2 @u @t = @2 u (1) @x 2 u(x; 0) = (x) u(0;t)=u(1;t)=0fort 0 1x, 1t N1x =1 x j = j1x, t n = n1t u(x j ;t n ) Uj n U n+1 j 1t 0 U n j =1t=(1x) 2 = U n j+1 0 2U n j + U n j01 (1x) 2 (2) U n+1

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

(1) 2000 ( ) ( ) 1000 2000 1000 0 http://www.spacepark.city.koriyama.fukushima.jp/ http://www.miraikan.jst.go.jp/ http://www.nasda.go.jp/ 3000 1 1 http://www.city.nara.nara.jp/citizen/jyugsidu/jgy/jsj/

More information

1/1 lim f(x, y) (x,y) (a,b) ( ) ( ) lim limf(x, y) lim lim f(x, y) x a y b y b x a ( ) ( ) xy x lim lim lim lim x y x y x + y y x x + y x x lim x x 1

1/1 lim f(x, y) (x,y) (a,b) ( ) ( ) lim limf(x, y) lim lim f(x, y) x a y b y b x a ( ) ( ) xy x lim lim lim lim x y x y x + y y x x + y x x lim x x 1 1/5 ( ) Taylor ( 7.1) (x, y) f(x, y) f(x, y) x + y, xy, e x y,... 1 R {(x, y) x, y R} f(x, y) x y,xy e y log x,... R {(x, y, z) (x, y),z f(x, y)} R 3 z 1 (x + y ) z ax + by + c x 1 z ax + by + c y x +

More information

6 6.1 B A: Γ d Q S(B) S(A) = S (6.1) T (e) Γ (6.2) : Γ B A R (reversible) 6-1

6 6.1 B A: Γ d Q S(B) S(A) = S (6.1) T (e) Γ (6.2) : Γ B A R (reversible) 6-1 6 6.1 B A: Γ d Q S(B) S(A) = S (6.1) (e) Γ (6.2) : Γ B A R (reversible) 6-1 (e) = Clausius 0 = B A: Γ B A: Γ d Q A + d Q (e) B: R d Q + S(A) S(B) (6.3) (e) // 6.2 B A: Γ d Q S(B) S(A) = S (6.4) (e) Γ (6.5)

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

2.

2. 2. 10 2. 2. 1995/12006/111995/42006/12 2. 10 1995120061119954200612 02505 025 05 025 02505 0303 02505 250100 250 200 100200 5010050 100200 100 100 50100 100200 50100 10 75100100 0250512 02505 1 025051205

More information

untitled

untitled 4 u(x) x 0 x x ::: x N u 0 u u ::: u N (discrete) u (interpolation) u (extrapolation) ( ) u(x) 3 n+ n (approximate functions) (polynomials) u(x) (discontinuities) u(x) 4. Lagrange u(x) n+ u 0 u ::: u n

More information

1 1.1 Excel Excel Excel log 1, log 2, log 3,, log 10 e = ln 10 log cm 1mm 1 10 =0.1mm = f(x) f(x) = n

1 1.1 Excel Excel Excel log 1, log 2, log 3,, log 10 e = ln 10 log cm 1mm 1 10 =0.1mm = f(x) f(x) = n 1 1.1 Excel Excel Excel log 1, log, log,, log e.7188188 ln log 1. 5cm 1mm 1 0.1mm 0.1 4 4 1 4.1 fx) fx) n0 f n) 0) x n n! n + 1 R n+1 x) fx) f0) + f 0) 1! x + f 0)! x + + f n) 0) x n + R n+1 x) n! 1 .

More information

2D-RCWA 1 two dimensional rigorous coupled wave analysis [1, 2] 1 ε(x, y) = 1 ε(x, y) = ϵ mn exp [+j(mk x x + nk y y)] (1) m,n= m,n= ξ mn exp [+j(mk x

2D-RCWA 1 two dimensional rigorous coupled wave analysis [1, 2] 1 ε(x, y) = 1 ε(x, y) = ϵ mn exp [+j(mk x x + nk y y)] (1) m,n= m,n= ξ mn exp [+j(mk x 2D-RCWA two dimensional rigoros copled wave analsis, 2] εx, εx, ϵ mn exp +jmk x x + nk ] m,n m,n ξ mn exp +jmk x x + nk ] 2 K x K x Λ x Λ ϵ mn ξ mn K x 2π Λ x K 2π Λ ϵ mn ξ mn Λ x Λ x Λ x Λ x Λx Λ Λx Λ

More information

B2 ( 19 ) Lebesgue ( ) ( ) 0 This note is c 2007 by Setsuo Taniguchi. It may be used for personal or classroom purposes, but not for commercia

B2 ( 19 ) Lebesgue ( ) ( ) 0 This note is c 2007 by Setsuo Taniguchi. It may be used for personal or classroom purposes, but not for commercia B2 ( 19) Lebesgue ( ) ( 19 7 12 ) 0 This note is c 2007 by Setsuo Taniguchi. It may be used for personal or classroom purposes, but not for commercial purposes. i Riemann f n : [0, 1] R 1, x = k (1 m

More information

a, b a bc c b a a b a a a a p > p p p 2, 3, 5, 7,, 3, 7, 9, 23, 29, 3, a > p a p [ ] a bp, b p p cq, c, q, < q < p a bp bcq q a <

a, b a bc c b a a b a a a a p > p p p 2, 3, 5, 7,, 3, 7, 9, 23, 29, 3, a > p a p [ ] a bp, b p p cq, c, q, < q < p a bp bcq q a < 22 9 8 5 22 9 29 0 2 2 5 2.............................. 5 2.2.................................. 6 2.3.............................. 8 3 8 4 9 4............................. 9 4.2 S(, a)..............................

More information

Fubini

Fubini 3............................... 3................................ 5.3 Fubini........................... 7.4.............................5..........................6.............................. 3.7..............................

More information

B000 B913 B913 S000 S500 L500 L913 B400 B913 B933 S320 L000 L913 492 498 P 38 5 P591 P595 P596 900 911 913 913 913 913 914 916 930 493 498 P 528 P594 P596 P597 910 913 913 913 913 913 914 918 700 723 746

More information

... 6

... 6 ... 6 1) 2) No. 01 02 03 04 05 06 07 08 09 10 11 12 No. 1 2 2 3 3cm 4

More information

( )

( ) ( ) () () 3 cm cm cm cm cm cm 1000 1500 50 500 1000 1000 1500 1000 10 50 300 1000 2000 1000 1500 50 10 1000 2000 300 50 1000 2000 1000 1500 50 10 1000 2000 300 30 10 300 1000 2000 1000 1500 1000

More information

-26-

-26- -25- -26- -27- -28- -29- -30- -31- -32- -33- -34- -35- -36- -37- -38- cm -39- -40- 1 2 3 4 4 3 2 1 5 5-41- -42- -43- -44- -45- -46- -47- -48- -49- -50- cm -51- -52- -53- -54- -55- -56- -57- -58- -59- -60-

More information

-1- 4 1 2 4-2- -- 2 2 cm 0 80cm 2 80cm 80cm 80cm 50cm 80cm 50cm 6 80cm 100 50 50 cm 10 6 4 50cm 4 4 50cm -4- -5- cm 50cm 4 4 4 50cm 50cm 4 80cm 50cm 80cm 50cm 6 cm -6- 20 250cm 1 2 1 4 0cm 60cm cm cm 1

More information

4 100g

4 100g 100g 10 20 30 40 50 60 70 80 4 5 7 9 12 15 19 24 60 100 10 80 100 20 10 5 20 195 20-1- 60 60 15 100 60 100 15 15 15 100 15 15 60 100 10 60 100 100 15 10 10 60 100 15 10 15 10 5-2- 80 80 24 100 80 100 24

More information

180 30 30 180 180 181 (3)(4) (3)(4)(2) 60 180 (1) (2) 20 (3)

180 30 30 180 180 181 (3)(4) (3)(4)(2) 60 180 (1) (2) 20 (3) 12 12 72 (1) (2) (3) 12 (1) (2) (3) (1) (2) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (1) (2) 180 30 30 180 180 181 (3)(4) (3)(4)(2) 60 180 (1) (2) 20 (3) 30 16 (1) 31 (2) 31 (3) (1) (2) (3) (4) 30

More information

untitled

untitled 1....1 2....2 2.1...2 2.2...2 3....14 3.1...14 3.2...14 4....15 4.1...15 4.2...18 4.3...21 4.4...23 4.5...26 5....27 5.1...27 5.2...35 5.3...54 5.4...64 5.5...75 6....79 6.1...79 6.2...85 6.3...94 6.4...

More information

113 120cm 1120cm 3 10cm 900 500+240 10 1 2 3 5 4 5 3 8 6 3 8 6 7 6 8 4 4 4 4 23 23 5 5 7

More information

untitled

untitled 21 14 487 2,322 2 7 48 4 15 ( 27) 14 3(1867) 3 () 1 2 3 ( 901923 ) 5 (1536) 3 4 5 6 7 8 ( ) () () 9 10 21 11 12 13 14 16 17 18 20 1 19 20 21 22 23 21 22 24 25 26 27 28 22 5 29 30cm 7.5m 1865 3 1820 5

More information

1948 1907 4024 1925 14 19281929 30 111931 4 3 15 4 16 3 15 4 161933 813 1935 12 17 11 17 1938 1945 2010 14 221 1945 10 1946 11 1947 1048 1947 1949 24

1948 1907 4024 1925 14 19281929 30 111931 4 3 15 4 16 3 15 4 161933 813 1935 12 17 11 17 1938 1945 2010 14 221 1945 10 1946 11 1947 1048 1947 1949 24 15 4 16 1988 63 28 19314 29 3 15 4 16 19283 15294 16 1930 113132 3 15 4 16 33 13 35 12 3 15 4 16 1945 10 10 10 10 40 1948 1907 4024 1925 14 19281929 30 111931 4 3 15 4 16 3 15 4 161933 813 1935 12 17 11

More information

裁定審議会における裁定の概要 (平成23年度)

裁定審議会における裁定の概要 (平成23年度) 23 23 23 4 24 3 10 11 12 13 14 () 1 23 7 21 23 12 14 (19 ) 30 1.876% 60 8 24 19 78 27 1 (10) 37 (3) 2 22 9 21 23 5 9 21 12 1 22 2 27 89 10 11 6 A B 3 21 12 1 12 10 10 12 5 1 9 1 2 61 ( 21 10 1 11 30 )

More information

Microsoft Word - 入居のしおり.doc

Microsoft Word - 入居のしおり.doc 1 1 2 2 2 3 2 4 3 5 3 6 3 7 3 8 4 1 7 2 7 3 7 4 8 5 9 6 9 7 10 8 10 9 11 10 11 11 11 12 12 13 13 1 14 2 17 3 18 4 19 5 20 6 22 (1) 24 (2) 24 (3) 24 (4) 24 (5) 24 (6) 25 (7) 25 (8) 25 (9) 25 1 29 (1) 29

More information

和県監査H15港湾.PDF

和県監査H15港湾.PDF ...1...1...1...1...1...1...1...1...2...2...2...3...3...3...5...5...10...11...12...13...13...13...14...14...14...14...14...14...15...15...15...15...15 ...16...17 14...17...18...18...19...21...23 2...25...27...27...28...28...28

More information

2002 (1) (2) (3) (4) (5) (1) (2) (3) (4) (5) (1) (2) (3) (4) (5) (6) (7) (8) (1) (2) (3) (4) (1) (2) (3) (4) (5) (6) (7) (8) No 2,500 3 200 200 200 200 200 50 200 No, 3 1 2 00 No 2,500 200 7 2,000 7

More information

untitled

untitled () () () () () ( ) () ( ) () ( ) () 2 () () 2 () () ( ) () () () 2 () () 2 3 ( ) () ( ) 2 3 4 () () 2 3 4 () () ( )( ) ( ) 2 ( ) 3 () () 2 3 () () 2 3 () () () () () () () () (( ) ( ) (( ))( )( ) ) 2 3

More information

河川砂防技術基準・基本計画編.PDF

河川砂防技術基準・基本計画編.PDF 4 1 1 1 1 1 2 1 2.1 1 2.2 2 2.3 2 2.4 2 3 2 4 3 2 4 1 4 1.1 4 1.2 4 2 4 2.1 4 2.2 4 2.3 5 2.4 5 2.5 5 2.5.1 5 2.5.2 5 2.6 5 2.6.1 5 2.6.2 5 2.6.3 5 2.6.4 5 2.6.5 6 2.7 6 2.7.1 6 2.7.2 6 2.7.3 6 2.7.4

More information

() n C + n C + n C + + n C n n (3) n C + n C + n C 4 + n C + n C 3 + n C 5 + (5) (6 ) n C + nc + 3 nc n nc n (7 ) n C + nc + 3 nc n nc n (

() n C + n C + n C + + n C n n (3) n C + n C + n C 4 + n C + n C 3 + n C 5 + (5) (6 ) n C + nc + 3 nc n nc n (7 ) n C + nc + 3 nc n nc n ( 3 n nc k+ k + 3 () n C r n C n r nc r C r + C r ( r n ) () n C + n C + n C + + n C n n (3) n C + n C + n C 4 + n C + n C 3 + n C 5 + (4) n C n n C + n C + n C + + n C n (5) k k n C k n C k (6) n C + nc

More information

1 yousuke.itoh/lecture-notes.html [0, π) f(x) = x π 2. [0, π) f(x) = x 2π 3. [0, π) f(x) = x 2π 1.2. Euler α

1   yousuke.itoh/lecture-notes.html [0, π) f(x) = x π 2. [0, π) f(x) = x 2π 3. [0, π) f(x) = x 2π 1.2. Euler α 1 http://sasuke.hep.osaka-cu.ac.jp/ yousuke.itoh/lecture-notes.html 1.1. 1. [, π) f(x) = x π 2. [, π) f(x) = x 2π 3. [, π) f(x) = x 2π 1.2. Euler dx = 2π, cos mxdx =, sin mxdx =, cos nx cos mxdx = πδ mn,

More information

M3 x y f(x, y) (= x) (= y) x + y f(x, y) = x + y + *. f(x, y) π y f(x, y) x f(x + x, y) f(x, y) lim x x () f(x,y) x 3 -

M3 x y f(x, y) (= x) (= y) x + y f(x, y) = x + y + *. f(x, y) π y f(x, y) x f(x + x, y) f(x, y) lim x x () f(x,y) x 3 - M3............................................................................................ 3.3................................................... 3 6........................................... 6..........................................

More information

all.dvi

all.dvi 29 4 Green-Lagrange,,.,,,,,,.,,,,,,,,,, E, σ, ε σ = Eε,,.. 4.1? l, l 1 (l 1 l) ε ε = l 1 l l (4.1) F l l 1 F 30 4 Green-Lagrange Δz Δδ γ = Δδ (4.2) Δz π/2 φ γ = π 2 φ (4.3) γ tan γ γ,sin γ γ ( π ) γ tan

More information

2001/Sep/12

2001/Sep/12 yamasita@imit.chiba-u.ac.jp 2001/Sep/12 i 1 1 2 1 2 2.1... 2 2.2... 3 2.3... 4 2.4 Roe... 5 2.5... 9 3 16 3.1 MUSCL... 16 3.1.1... 16 3.1.2... 18 3.1.3... 21 3.2 Chakravarthy-Osher... 22 4 2 33 5 3 40

More information

統計学のポイント整理

統計学のポイント整理 .. September 17, 2012 1 / 55 n! = n (n 1) (n 2) 1 0! = 1 10! = 10 9 8 1 = 3628800 n k np k np k = n! (n k)! (1) 5 3 5 P 3 = 5! = 5 4 3 = 60 (5 3)! n k n C k nc k = npk k! = n! k!(n k)! (2) 5 3 5C 3 = 5!

More information

1 filename=mathformula tex 1 ax 2 + bx + c = 0, x = b ± b 2 4ac, (1.1) 2a x 1 + x 2 = b a, x 1x 2 = c a, (1.2) ax 2 + 2b x + c = 0, x = b ± b 2

1 filename=mathformula tex 1 ax 2 + bx + c = 0, x = b ± b 2 4ac, (1.1) 2a x 1 + x 2 = b a, x 1x 2 = c a, (1.2) ax 2 + 2b x + c = 0, x = b ± b 2 filename=mathformula58.tex ax + bx + c =, x = b ± b 4ac, (.) a x + x = b a, x x = c a, (.) ax + b x + c =, x = b ± b ac. a (.3). sin(a ± B) = sin A cos B ± cos A sin B, (.) cos(a ± B) = cos A cos B sin

More information

A

A A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................

More information

72 5 f (x) f Tylor f (x) f (x) = f (x) + 2 f (x) + 2 3! f (x) + (5.) = f (x) + O() = f (x) 2 f (x) + 2 3! f (x) (5.2) = f (x) + O() δ f 2 = ( f (x) +

72 5 f (x) f Tylor f (x) f (x) = f (x) + 2 f (x) + 2 3! f (x) + (5.) = f (x) + O() = f (x) 2 f (x) + 2 3! f (x) (5.2) = f (x) + O() δ f 2 = ( f (x) + 7 5 Derivtives nd integrls re lso very useful s topics for tecing bout numericl computtion nd nlysis. Tey re esily understood, one cn present intuitive or pictoril motivtions, te lgebr is usully not too

More information

(2018 2Q C) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = ( ) a c b d (a c, b d) P = (a, b) O P ( ) a p = b P = (a, b) p = ( ) a b R 2 {( ) } R 2 x = x, y

(2018 2Q C) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = ( ) a c b d (a c, b d) P = (a, b) O P ( ) a p = b P = (a, b) p = ( ) a b R 2 {( ) } R 2 x = x, y (2018 2Q C) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = a c b d (a c, b d) P = (a, b) O P a p = b P = (a, b) p = a b R 2 { } R 2 x = x, y R y 2 a p =, c q = b d p + a + c q = b + d q p P q a p = c R c b

More information

C による数値計算法入門 ( 第 2 版 ) 新装版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 新装版 1 刷発行時のものです.

C による数値計算法入門 ( 第 2 版 ) 新装版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.  このサンプルページの内容は, 新装版 1 刷発行時のものです. C による数値計算法入門 ( 第 2 版 ) 新装版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009383 このサンプルページの内容は, 新装版 1 刷発行時のものです. i 2 22 2 13 ( ) 2 (1) ANSI (2) 2 (3) Web http://www.morikita.co.jp/books/mid/009383

More information

: BV15005

: BV15005 29 5 26 : BV15005 1 1 1.1............................................. 1 1.2........................................ 1 1.3........................................ 1 2 3 2.1.............................................

More information

(2016 2Q H) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = ( ) a c b d (a c, b d) P = (a, b) O P ( ) a p = b P = (a, b) p = ( ) a b R 2 {( ) } R 2 x = x, y

(2016 2Q H) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = ( ) a c b d (a c, b d) P = (a, b) O P ( ) a p = b P = (a, b) p = ( ) a b R 2 {( ) } R 2 x = x, y (2016 2Q H) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = a c b d (a c, b d) P = (a, b) O P a p = b P = (a, b) p = a b R 2 { } R 2 x = x, y R y 2 a p =, c q = b d p + a + c q = b + d q p P q a p = c R c b

More information

OpenMP¤òÍѤ¤¤¿ÊÂÎó·×»»¡Ê£²¡Ë

OpenMP¤òÍѤ¤¤¿ÊÂÎó·×»»¡Ê£²¡Ë 2013 5 30 (schedule) (omp sections) (omp single, omp master) (barrier, critical, atomic) program pi i m p l i c i t none integer, parameter : : SP = kind ( 1. 0 ) integer, parameter : : DP = selected real

More information

i

i i 3 4 4 7 5 6 3 ( ).. () 3 () (3) (4) /. 3. 4/3 7. /e 8. a > a, a = /, > a >. () a >, a =, > a > () a > b, a = b, a < b. c c n a n + b n + c n 3c n..... () /3 () + (3) / (4) /4 (5) m > n, a b >, m > n,

More information

.3. (x, x = (, u = = 4 (, x x = 4 x, x 0 x = 0 x = 4 x.4. ( z + z = 8 z, z 0 (z, z = (0, 8, (,, (8, 0 3 (0, 8, (,, (8, 0 z = z 4 z (g f(x = g(

.3. (x, x = (, u = = 4 (, x x = 4 x, x 0 x = 0 x = 4 x.4. ( z + z = 8 z, z 0 (z, z = (0, 8, (,, (8, 0 3 (0, 8, (,, (8, 0 z = z 4 z (g f(x = g( 06 5.. ( y = x x y 5 y 5 = (x y = x + ( y = x + y = x y.. ( Y = C + I = 50 + 0.5Y + 50 r r = 00 0.5Y ( L = M Y r = 00 r = 0.5Y 50 (3 00 0.5Y = 0.5Y 50 Y = 50, r = 5 .3. (x, x = (, u = = 4 (, x x = 4 x,

More information

68 A mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1

68 A mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1 67 A Section A.1 0 1 0 1 Balmer 7 9 1 0.1 0.01 1 9 3 10:09 6 A.1: A.1 1 10 9 68 A 10 9 10 9 1 10 9 10 1 mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1 A.1. 69 5 1 10 15 3 40 0 0 ¾ ¾ É f Á ½ j 30 A.3: A.4: 1/10

More information

2S III IV K A4 12:00-13:30 Cafe David 1 2 TA 1 appointment Cafe David K2-2S04-00 : C

2S III IV K A4 12:00-13:30 Cafe David 1 2 TA 1  appointment Cafe David K2-2S04-00 : C 2S III IV K200 : April 16, 2004 Version : 1.1 TA M2 TA 1 10 2 n 1 ɛ-δ 5 15 20 20 45 K2-2S04-00 : C 2S III IV K200 60 60 74 75 89 90 1 email 3 4 30 A4 12:00-13:30 Cafe David 1 2 TA 1 email appointment Cafe

More information

1F90/kouhou_hf90.dvi

1F90/kouhou_hf90.dvi Fortran90 3 33 1 2 Fortran90 FORTRAN 1956 IBM IBM704 FORTRAN(FORmula TRANslation ) 1965 FORTRAN66 1978 FORTRAN77 1991 Fortran90 Fortran90 Fortran Fortran90 6 Fortran90 77 90 90 Fortran90 [ ] Fortran90

More information

Part () () Γ Part ,

Part () () Γ Part , Contents a 6 6 6 6 6 6 6 7 7. 8.. 8.. 8.3. 8 Part. 9. 9.. 9.. 3. 3.. 3.. 3 4. 5 4.. 5 4.. 9 4.3. 3 Part. 6 5. () 6 5.. () 7 5.. 9 5.3. Γ 3 6. 3 6.. 3 6.. 3 6.3. 33 Part 3. 34 7. 34 7.. 34 7.. 34 8. 35

More information

&A : A = k j 1: 4-way., A set x, way y, way y LRU y, way., A (x,y).,,, L1( 1) L2, L3 3. L1., L2,L3., TLB(Translation Lookaside Buffer). OS,. TLB, ( ),

&A : A = k j 1: 4-way., A set x, way y, way y LRU y, way., A (x,y).,,, L1( 1) L2, L3 3. L1., L2,L3., TLB(Translation Lookaside Buffer). OS,. TLB, ( ), 1?,. 1,.,,. n-way (n ). 1, 4-way, n-way n (way).,., 1., ( set x ) (x), n., 2, 2 s, 2 l (, s, l )., s + l s., s,., n s. n. s + l way, (set,way)., way,. way, LRU(Least Recently Used, ). way. way, (,...).

More information

- II

- II - II- - -.................................................................................................... 3.3.............................................. 4 6...........................................

More information

構造と連続体の力学基礎

構造と連続体の力学基礎 II 37 Wabash Avenue Bridge, Illinois 州 Winnipeg にある歩道橋 Esplanade Riel 橋6 6 斜張橋である必要は多分無いと思われる すぐ横に道路用桁橋有り しかも塔基部のレストランは 8 年には営業していなかった 9 9. 9.. () 97 [3] [5] k 9. m w(t) f (t) = f (t) + mg k w(t) Newton

More information