86 6 r (6) y y d y = y 3 (64) y r y r y r ϕ(x, y, y,, y r ) n dy = f(x, y) (6) 6 Lipschitz 6 dy = y x c R y(x) y(x) = c exp(x) x x = x y(x ) = y (init
|
|
- ふさこ よせ
- 6 years ago
- Views:
Transcription
1 8 6 ( ) ( ) 6 ( ϕ x, y, dy ), d y,, dr y r = (x R, y R n ) (6) n r y(x) (explicit) d r ( y r = ϕ x, y, dy ), d y,, dr y r y y y r (6) dy = f (x, y) (63) = y dy/ d r y/ r
2 86 6 r (6) y y d y = y 3 (64) y r y r y r ϕ(x, y, y,, y r ) n dy = f(x, y) (6) 6 Lipschitz 6 dy = y x c R y(x) y(x) = c exp(x) x x = x y(x ) = y (initial value problem, IVP) dy = f(x, y) (66) y(x ) = y Lipschitz 6 (Lipschitz ) x [x, α] z, z R n f(x, z ) f(x, z ) L z z (67) L R (66) y(x)
3 63 : Euler,, Runge-Kutta 87 6 dy = y y() = L = Lipschitz (67) y(x) = exp(x) Lipschitz (66) 63 : Euler,, Runge-Kutta [x, α] l x i = x + i j= h j y(x i ) y i h i (i =,,, l ) (Step Size) (Stiff) y/h=(y(x + h) y(x))/h dy/ y i+ := y i + h i f(x i, y i ) (68) Euler ( Euler Euler ) δy/(h)=(y(x + h) y(x h))/(h) y i+ := y i + h i f(x i, y i ) (69) y = y(x h ) y Euler y 3 (Euler ) x i (i =,,, l ) y i+ := y i + h i f(x i, y i ) 36 ( ) y Euler x i (i =,, l ) y i+ := y i + h i f(x i, y i )
4 88 6 Euler ( ) Runge-Kutta 37 ( Runge-Kutta ) x i (i =,,, l ) (a) k, k, k 3, k 4 k = f(x i, y i ) k = f(x i + h i, y i + h ik ) k 3 = f(x i + h i, y i + h ik ) k 4 = f(x i + h i, y i + h i k 3 ) (b) y i+ := y i + 6 h i(k + k + k 3 + k 4 ) 4 63 dy = x + y y() = Euler Runge-Kutta 6 相対誤差 e-6 e-9 Euler, h=/ Euler, h=/ 4 Midpoint h=/ ERK44, h=/ MIdpoint, h=/ 4 e- e- ERK44, h=/ x 6: Euler,, Runge-Kutta
5 64 Runge-Kutta x = h = /,/4 Euler, Runge-Kutta 4 [Hint: f Taylor 4 ] 64 Runge-Kutta Runge-Kutta m Runge-Kutta (m Runge-Kutta ) [x, α] l x, x,, x l = α x i x i+ h i x i (i =,,, l ) (a) k, k,, k m k = f(x i + c h i, y i + h i m j= a j k j ) k = f(x i + c h i, y i + h i m j= a j k j ) k m = f(x i + c m h i, y i + h i m j= a m j k j ) (b) y i+ := y i + h i m j= w j k j c, c,, c m, a,, a m,m, w, w,, w m 6 O(h p+ ) m p Runge-Kutta ( 6) 3 Runge-Kutta c a c 3 a 3 a 3 c m a m a m a m,m w w w m w m Butcher Table
6 9 6 6: m Runge-Kutta c a a a m c a a a m c m a m a m a m,m w w w m c = a i j = (i j) k k m Runge-Kutta c a c a a c m a m a m a m,m w w w m a i j = (i < j) k j ( j =,, m) n 3 Runge-Kutta j k j = f(x i + c j h i, y i + h i a js k j + h i a j j k j ) (6) s= c a a a m c a a a m c m a m a m a m,m w w w m k,, k m mn k = f(x i + c h i, y i + h i m j= a j k j ) k = f(x i + c h i, y i + h i m j= a j k j ) k m = f(x i + c m h i, y i + h i m j= a m j k j ) (6)
7 6 Runge-Kutta 9 Runge-Kutta Butcher[] Lawson, Huta, Shanks, [3, 36, 37, 38] Runge-Kutta Runge-Kutta ( = ) A- Runge-Kutta Runge-Kutta Newton ( ) Runge-Kutta Stiff Runge-Kutta [9] Runge-Kutta 6 Runge-Kutta k k m (6, 6) 6 Newton Newton 39 ( Newton ) k (),, k() m for l =,,, (a) J( (b), k(l) k (l+) k (l+) k (l+) m,, k(l) m ) := m J (,, k(l) m ) f(x i + c h i, y i + h mj= i a j j ) f(x i + c h i, y i + h mj= i a j j ) m f(x i + c m h i, y i + h mj= i a m j j ) I n J J J m J(k (l, k(l),, k(l) m ) = J I n J J m J m J m I n J mm
8 9 6 J( J(, k(l) J pq = h i a pq y f(x i + c p h i, y i + h i, k(l) m j= a p j j ) M n (R),, k(l) m ) LU,, k(l) m ) z z z m = f(x i + c h i, y i + h mj= i a j j ) f(x i + c h i, y i + h mj= i a j j ) m f(x i + c m h i, y i + h mj= i a m j j ) mn [z z z n ] T Newton LU LU Jacobi [6] 66 Runge-Kutta Runge-Kutta Runge-Kutta Stiff Runge-Kutta (6) step Runge-Kutta Runge-Kutta Runge-Kutta Runge-Kutta 66 Runge-Kutta dy = Ay + g(x) = f(x, y) y(x ) = y (6) A M n (R) (6) (6) n mn r (I h i a rr A)k r = f(x i + c r h i, y i + h i a r j k j ) (63) j=
9 66 Runge-Kutta 93 I h i a A h i a A h i a m A h i a A h i a m,m A h i a m A h i a m,m I h i a mm A k k m = f(x i + c h i, y i ) f(x i + c m h i, y i ) (64) 66 Runge-Kutta Runge-Kutta 3 [] 4 : Butcher[] : Butcher[] : Optimal[] : (classical)runge-kutta ( 37)[] : Butcher[]
10 94 6 Butcher Jain[] Stiff 66 dy = y + y cos x dy = y 3y + 3 cos x sin x y () = y () = (6) 66 dy = y + y cos x dy = 998y 999y cos x sin x y () =, y () = (66) (6) (66) [] (67) y (x) = exp( x) y (x) = exp( x) + cos x (67) Stiff 6 6 [, ] (x l = ) 66( 6, 63) 66( 64, 6) 66 (6)(66) (67) 67 Rössler Rössler Model [33] d dt x y z = (y + z) x + αy β + z(x µ) (68) α = β = / µ 3, 4, [33]
11 67 Rössler 9 6: 66 : Runge-Kutta -:optimal 4-4:classical 7-6:Butcher y y y y y y / 8e+6 376e- 9948e+4 6e-3 644e+3 6e- / 4 639e+3 3e-3 4e+ 634e-6 98e- 9e-9 / 6 977e+ 63e- 836e- 978e-9 9e-4 9e- / 8 373e+ 3887e-6 368e e- 864e-8 68e-6 / 9e+ 47e-7 4e- 347e-3 976e-8 63e- / 986e- 9e-8 9e-7 7e- 9e-7 993e- / 4 68e-3 94e- 794e-7 34e- 697e-7 39e- / 6 386e-4 89e- 889e-7 74e-4 894e-7 74e-4 / 8 66e- 37e- 3e-6 7e-4 e-6 99e-4 63: 66 : Runge-Kutta - -4:Butcher 3-6:Butcher y y y y y y / 4e+4 894e-3 7e+3 86e- 49e+ 3e-8 / 4 49e+3 64e-4 434e+ 76e-7 66e-3 9e- / 6 737e+ 446e- 644e- 779e- 34e-7 33e- / 8 e+ 6e-6 4e-4 34e- 6e-8 e+ / 7849e- 6e e-7 99e-4 986e-8 493e- / 4939e- e-8 43e-7 363e- 37e-7 363e- / 4 393e-3 664e- 78e-7 364e- 773e-7 34e- / 6 938e-4 39e- 944e-7 74e-4 89e-7 687e-4 / 8 e- 47e- 3e-6 44e-4 6e-6 44e-4
12 : 66 : Runge-Kutta -:optimal 4-4:classical 7-6:Butcher y y y y y y / NaN NaN NaN NaN NaN NaN / 4 NaN NaN NaN NaN NaN NaN / 6 NaN NaN NaN NaN NaN NaN / 8 NaN NaN NaN NaN NaN NaN / 6776e- 796e-6 79e- 89e-7 89e-3 8e-8 / 89e-4 36e-8 86e- 873e- e-7 47e- / e- 66e- 63e-8 646e-3 63e-9 48e-6 / 6 33e e- 9e-9 3e- 349e-9 e+ / 8 48e-7 44e- 99e-9 68e-6 373e-9 68e-6 ( NaN ) 6: 66 : Runge-Kutta - -4:Butcher 3-6:Butcher y y y y y y / 7e+ 49e- 97e+ 98e e- 47e-6 / 4 6e+ 4889e-4 693e+ 69e- 68e-3 9e-8 / 6 74e+ 3e- 8e- 798e-7 933e-6 879e- / e- 99e-6 37e-4 38e-9 363e-8 4e-3 / 484e-3 93e-7 936e e- 797e- e- / 334e-4 748e e-9 367e-4 746e- 36e-6 / 4 898e- 466e- 6e-9 36e-6 773e-9 36e-6 / 6 86e-6 93e- 6e-9 e+ e-9 e+ / 8 76e-8 8e- e-9 68e-6 6e-9 44e-6
13 67 Ro ssler モデルの数値例 97 y y x [, ], h = /89 Explicit: 6th Order - x [, ], h = /89 Implicit: 6th Order 図 6: 6 次の陽的 R-K 法 (左) と陰的 R-K 法 (右) の数値解 (µ = 4) y y x [,], h = /89 Explicit: 6th Order - x [,], h = /89 Implicit: 6th Order 図 63: 6 次の陽的 R-K 法 (左) と陰的 R-K 法 (右) の数値解 (µ = )
14 98 6 [ ] T [, ] 6 Runge-Kutta Runge- Kutta (68) h h = /89 Fig6 µ = 4 Fig63 µ = Rössler Model µ 3, 4,, µ = Runge-Kutta t = x 7 6 Runge-Kutta 3 6 Runge-Kutta (68) IEEE74 Explicit : 7 step (6th order), µ = h = /496 h = /89 RN 48989e e + RZ 674e + 99e + RP 43638e + 384e + RM 63464e e + Implicit : 3 step (6th order), µ = h = /496 h = /89 RN e e + RZ 67438e + 678e + RP 3447e e + RM e e + µ = y RN, RM, RP, RZ ( ) [7] order order (68) [, ] Rössler Model IEEE74 µ = 7 z µ = 7 z Fig64 µ µ =, 7 µ = 4
15 67 Rössler 99 e- e- Roessler Model: gmp 8bit(38 digits) Extrapolation: 4 steps Relative Round-off Error: h=/4 \mu= \mu=7 Roessler Model: [, ], 8bit Extrapolation 4steps, h=/4 Relative Error e-3 e-3 e-4 \mu=4 z e-4 e- 3 4 t 3 4 t 64: µ = 4,, 7 z ( ) µ = 7 z( ) (a) dy = x y y() = (b) h = /4 Euler y() ỹ() (c) ỹ(), 994,
16 6, 4 Runge-Kutta A Hamilton EHairer, SPNørsett and GWanner Solving Ordinary Differential Equations, I Springer-Verlag, 993 EHairer and GWanner Solving Ordinary Differential Equations, II Springer-Verlag, 996 plain TEX TEX
d dt A B C = A B C d dt x = Ax, A 0 B 0 C 0 = mm 0 mm 0 mm AP = PΛ P AP = Λ P A = ΛP P d dt x = P Ax d dt (P x) = Λ(P x) d dt P x =
3 MATLAB Runge-Kutta Butcher 3. Taylor Taylor y(x 0 + h) = y(x 0 ) + h y (x 0 ) + h! y (x 0 ) + Taylor 3. Euler, Runge-Kutta Adams Implicit Euler, Implicit Runge-Kutta Gear y n+ y n (n+ ) y n+ y n+ y n+
More informationTrapezoidal Rule θ = 1/ x n x n 1 t = 1 [f(t n 1, x n 1 ) + f(t n, x n )] (6) 1. dx dt = f(t, x), x(t 0) = x 0 (7) t [t 0, t 1 ] f t [t 0, t 1 ], x x
University of Hyogo 8 8 1 d x(t) =f(t, x(t)), dt (1) x(t 0 ) =x 0 () t n = t 0 + n t x x n n x n x 0 x i i = 0,..., n 1 x n x(t) 1 1.1 1 1 1 0 θ 1 θ x n x n 1 t = θf(t n 1, x n 1 ) + (1 θ)f(t n, x n )
More informationA
A05-132 2010 2 11 1 1 3 1.1.......................................... 3 1.2..................................... 3 1.3..................................... 3 2 4 2.1............................... 4 2.2
More information数学Ⅱ演習(足助・09夏)
II I 9/4/4 9/4/2 z C z z z z, z 2 z, w C zw z w 3 z, w C z + w z + w 4 t R t C t t t t t z z z 2 z C re z z + z z z, im z 2 2 3 z C e z + z + 2 z2 + 3! z3 + z!, I 4 x R e x cos x + sin x 2 z, w C e z+w
More informationerror_g1.eps
Runge-Kutta Method Runge-Kutta Method x n+ = x n +hf(t n x n ) dx dt dx x(t+h) x(t) (t) = lim dt h h t = t n h > dx dt (t x(t n +h) x(t n ) n) = lim x(t n+) x(t n ) h h h x = f(tx) x(t n+ ) x(t n ) f(t
More informationp12.dvi
301 12 (2) : 1 (1) dx dt = f(x,t) ( (t 0,t 1,...,t N ) ) h k = t k+1 t k. h k k h. x(t k ) x k. : 2 (2) :1. step. 1 : explicit( ) : ξ k+1 = ξ k +h k Ψ(t k,ξ k,h k ) implicit( ) : ξ k+1 = ξ k +h k Ψ(t k,t
More informationn (1.6) i j=1 1 n a ij x j = b i (1.7) (1.7) (1.4) (1.5) (1.4) (1.7) u, v, w ε x, ε y, ε x, γ yz, γ zx, γ xy (1.8) ε x = u x ε y = v y ε z = w z γ yz
1 2 (a 1, a 2, a n ) (b 1, b 2, b n ) A (1.1) A = a 1 b 1 + a 2 b 2 + + a n b n (1.1) n A = a i b i (1.2) i=1 n i 1 n i=1 a i b i n i=1 A = a i b i (1.3) (1.3) (1.3) (1.1) (ummation convention) a 11 x
More informationuntitled
3 3. (stochastic differential equations) { dx(t) =f(t, X)dt + G(t, X)dW (t), t [,T], (3.) X( )=X X(t) : [,T] R d (d ) f(t, X) : [,T] R d R d (drift term) G(t, X) : [,T] R d R d m (diffusion term) W (t)
More informationI ( ) 1 de Broglie 1 (de Broglie) p λ k h Planck ( Js) p = h λ = k (1) h 2π : Dirac k B Boltzmann ( J/K) T U = 3 2 k BT
I (008 4 0 de Broglie (de Broglie p λ k h Planck ( 6.63 0 34 Js p = h λ = k ( h π : Dirac k B Boltzmann (.38 0 3 J/K T U = 3 k BT ( = λ m k B T h m = 0.067m 0 m 0 = 9. 0 3 kg GaAs( a T = 300 K 3 fg 07345
More informationf : R R f(x, y) = x + y axy f = 0, x + y axy = 0 y 直線 x+y+a=0 に漸近し 原点で交叉する美しい形をしている x +y axy=0 X+Y+a=0 o x t x = at 1 + t, y = at (a > 0) 1 + t f(x, y
017 8 10 f : R R f(x) = x n + x n 1 + 1, f(x) = sin 1, log x x n m :f : R n R m z = f(x, y) R R R R, R R R n R m R n R m R n R m f : R R f (x) = lim h 0 f(x + h) f(x) h f : R n R m m n M Jacobi( ) m n
More information: 1g99p038-8
16 17 : 1g99p038-8 1 3 1.1....................................... 4 1................................... 5 1.3.................................. 5 6.1..................................... 7....................................
More informationdy = sin cos y cos () y () 1 y = sin 1 + c 1 e sin (3) y() () y() y( 0 ) = y 0 y 1 1. (1) d (1) y = f(, y) (4) i y y i+1 y i+1 = y( i + ) = y i
007 8 8 4 1 1.1 ( ) (partial differential equation) (ordinary differential equation) 1 dy = f(, y) (1) 1 1 y() (1) y() (, y) 1 dy = sin cos y cos () y () 1 y = sin 1 + c 1 e sin (3) 1 1 5 y() () y() y(
More information研修コーナー
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 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 l l l l
More information1 8, : 8.1 1, 2 z = ax + by + c ax by + z c = a b +1 x y z c = 0, (0, 0, c), n = ( a, b, 1). f = n i=1 a ii x 2 i + i<j 2a ij x i x j = ( x, A x), f =
1 8, : 8.1 1, z = ax + by + c ax by + z c = a b +1 x y z c = 0, (0, 0, c), n = ( a, b, 1). f = a ii x i + i
More informationi
009 I 1 8 5 i 0 1 0.1..................................... 1 0.................................................. 1 0.3................................. 0.4........................................... 3
More informationW u = u(x, t) u tt = a 2 u xx, a > 0 (1) D := {(x, t) : 0 x l, t 0} u (0, t) = 0, u (l, t) = 0, t 0 (2)
3 215 4 27 1 1 u u(x, t) u tt a 2 u xx, a > (1) D : {(x, t) : x, t } u (, t), u (, t), t (2) u(x, ) f(x), u(x, ) t 2, x (3) u(x, t) X(x)T (t) u (1) 1 T (t) a 2 T (t) X (x) X(x) α (2) T (t) αa 2 T (t) (4)
More informationII 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 information1 n A a 11 a 1n A =.. a m1 a mn Ax = λx (1) x n λ (eigenvalue problem) x = 0 ( x 0 ) λ A ( ) λ Ax = λx x Ax = λx y T A = λy T x Ax = λx cx ( 1) 1.1 Th
1 n A a 11 a 1n A = a m1 a mn Ax = λx (1) x n λ (eigenvalue problem) x = ( x ) λ A ( ) λ Ax = λx x Ax = λx y T A = λy T x Ax = λx cx ( 1) 11 Th9-1 Ax = λx λe n A = λ a 11 a 12 a 1n a 21 λ a 22 a n1 a n2
More informationnote1.dvi
(1) 1996 11 7 1 (1) 1. 1 dx dy d x τ xx x x, stress x + dx x τ xx x+dx dyd x x τ xx x dyd y τ xx x τ xx x+dx d dx y x dy 1. dx dy d x τ xy x τ x ρdxdyd x dx dy d ρdxdyd u x t = τ xx x+dx dyd τ xx x dyd
More information1 4 1 ( ) ( ) ( ) ( ) () 1 4 2
7 1995, 2017 7 21 1 2 2 3 3 4 4 6 (1).................................... 6 (2)..................................... 6 (3) t................. 9 5 11 (1)......................................... 11 (2)
More information1 28 6 12 7 1 7.1...................................... 2 7.1.1............................... 2 7.1.2........................... 2 7.2...................................... 3 7.3...................................
More information211 kotaro@math.titech.ac.jp 1 R *1 n n R n *2 R n = {(x 1,..., x n ) x 1,..., x n R}. R R 2 R 3 R n R n R n D D R n *3 ) (x 1,..., x n ) f(x 1,..., x n ) f D *4 n 2 n = 1 ( ) 1 f D R n f : D R 1.1. (x,
More information() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi)
0. A A = 4 IC () det A () A () x + y + z = x y z X Y Z = A x y z ( 5) ( s5590) 0. a + b + c b c () a a + b + c c a b a + b + c 0 a b c () a 0 c b b c 0 a c b a 0 0. A A = 7 5 4 5 0 ( 5) ( s5590) () A ()
More information1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0
1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0 0 < t < τ I II 0 No.2 2 C x y x y > 0 x 0 x > b a dx
More informationII (10 4 ) 1. p (x, y) (a, b) ε(x, y; a, b) 0 f (x, y) f (a, b) A, B (6.5) y = b f (x, b) f (a, b) x a = A + ε(x, b; a, b) x a 2 x a 0 A = f x (
II (1 4 ) 1. p.13 1 (x, y) (a, b) ε(x, y; a, b) f (x, y) f (a, b) A, B (6.5) y = b f (x, b) f (a, b) x a = A + ε(x, b; a, b) x a x a A = f x (a, b) y x 3 3y 3 (x, y) (, ) f (x, y) = x + y (x, y) = (, )
More information17 1 25 http://grape.astron.s.u-tokyo.ac.jp/pub/people/makino/kougi/stellar_dynamics/index.html http://grape.astron.s.u-tokyo.ac.jp/pub/people/makino/talks/index-j.html d 2 x i dt 2 = j i Gm j x j x i
More information1. 1 A : l l : (1) l m (m 3) (2) m (3) n (n 3) (4) A α, β γ α β + γ = 2 m l lm n nα nα = lm. α = lm n. m lm 2β 2β = lm β = lm 2. γ l 2. 3
1. 1 A : l l : (1) l m (m 3) (2) m (3) n (n 3) (4) A 2 1 2 1 2 3 α, β γ α β + γ = 2 m l lm n nα nα = lm. α = lm n. m lm 2β 2β = lm β = lm 2. γ l 2. 3 4 P, Q R n = {(x 1, x 2,, x n ) ; x 1, x 2,, x n R}
More information#A A A F, F d F P + F P = d P F, F y P F F x A.1 ( α, 0), (α, 0) α > 0) (x, y) (x + α) 2 + y 2, (x α) 2 + y 2 d (x + α)2 + y 2 + (x α) 2 + y 2 =
#A A A. F, F d F P + F P = d P F, F P F F A. α, 0, α, 0 α > 0, + α +, α + d + α + + α + = d d F, F 0 < α < d + α + = d α + + α + = d d α + + α + d α + = d 4 4d α + = d 4 8d + 6 http://mth.cs.kitmi-it.c.jp/
More information17 ( ) II III A B C(100 ) 1, 2, 6, 7 II A B (100 ) 2, 5, 6 II A B (80 ) 8 10 I II III A B C(80 ) 1 a 1 = 1 2 a n+1 = a n + 2n + 1 (n = 1,
17 ( ) 17 5 1 4 II III A B C(1 ) 1,, 6, 7 II A B (1 ), 5, 6 II A B (8 ) 8 1 I II III A B C(8 ) 1 a 1 1 a n+1 a n + n + 1 (n 1,,, ) {a n+1 n } (1) a 4 () a n OA OB AOB 6 OAB AB : 1 P OB Q OP AQ R (1) PQ
More informationohpr.dvi
2003-08-04 1984 VP-1001 CPU, 250 MFLOPS, 128 MB 2004ASCI Purple (LLNL)64 CPU 197, 100 TFLOPS, 50 TB, 4.5 MW PC 2 CPU 16, 4 GFLOPS, 32 GB, 3.2 kw 20028 CPU 640, 40 TFLOPS, 10 TB, 10 MW (ASCI: Accelerated
More informationd ϕ i) t d )t0 d ϕi) ϕ i) t x j t d ) ϕ t0 t α dx j d ) ϕ i) t dx t0 j x j d ϕ i) ) t x j dx t0 j f i x j ξ j dx i + ξ i x j dx j f i ξ i x j dx j d )
23 M R M ϕ : R M M ϕt, x) ϕ t x) ϕ s ϕ t ϕ s+t, ϕ 0 id M M ϕ t M ξ ξ ϕ t d ϕ tx) ξϕ t x)) U, x 1,...,x n )) ϕ t x) ϕ 1) t x),...,ϕ n) t x)), ξx) ξ i x) d ϕi) t x) ξ i ϕ t x)) M f ϕ t f)x) f ϕ t )x) fϕ
More information(iii) 0 V, x V, x + 0 = x. 0. (iv) x V, y V, x + y = 0., y x, y = x. (v) 1x = x. (vii) (α + β)x = αx + βx. (viii) (αβ)x = α(βx)., V, C.,,., (1)
1. 1.1...,. 1.1.1 V, V x, y, x y x + y x + y V,, V x α, αx αx V,, (i) (viii) : x, y, z V, α, β C, (i) x + y = y + x. (ii) (x + y) + z = x + (y + z). 1 (iii) 0 V, x V, x + 0 = x. 0. (iv) x V, y V, x + y
More information85 4
85 4 86 Copright c 005 Kumanekosha 4.1 ( ) ( t ) t, t 4.1.1 t Step! (Step 1) (, 0) (Step ) ±V t (, t) I Check! P P V t π 54 t = 0 + V (, t) π θ : = θ : π ) θ = π ± sin ± cos t = 0 (, 0) = sin π V + t +V
More information20 9 19 1 3 11 1 3 111 3 112 1 4 12 6 121 6 122 7 13 7 131 8 132 10 133 10 134 12 14 13 141 13 142 13 143 15 144 16 145 17 15 19 151 1 19 152 20 2 21 21 21 211 21 212 1 23 213 1 23 214 25 215 31 22 33
More information基礎数学I
I & II ii ii........... 22................. 25 12............... 28.................. 28.................... 31............. 32.................. 34 3 1 9.................... 1....................... 1............
More informationma22-9 u ( v w) = u v w sin θê = v w sin θ u cos φ = = 2.3 ( a b) ( c d) = ( a c)( b d) ( a d)( b c) ( a b) ( c d) = (a 2 b 3 a 3 b 2 )(c 2 d 3 c 3 d
A 2. x F (t) =f sin ωt x(0) = ẋ(0) = 0 ω θ sin θ θ 3! θ3 v = f mω cos ωt x = f mω (t sin ωt) ω t 0 = f ( cos ωt) mω x ma2-2 t ω x f (t mω ω (ωt ) 6 (ωt)3 = f 6m ωt3 2.2 u ( v w) = v ( w u) = w ( u v) ma22-9
More informationChap11.dvi
. () x 3 + dx () (x )(x ) dx + sin x sin x( + cos x) dx () x 3 3 x + + 3 x + 3 x x + x 3 + dx 3 x + dx 6 x x x + dx + 3 log x + 6 log x x + + 3 rctn ( ) dx x + 3 4 ( x 3 ) + C x () t x t tn x dx x. t x
More information1 3 1.1.......................... 3 1............................... 3 1.3....................... 5 1.4.......................... 6 1.5........................ 7 8.1......................... 8..............................
More information1 (1) ( i ) 60 (ii) 75 (iii) 315 (2) π ( i ) (ii) π (iii) 7 12 π ( (3) r, AOB = θ 0 < θ < π ) OAB A 2 OB P ( AB ) < ( AP ) (4) 0 < θ < π 2 sin θ
1 (1) ( i ) 60 (ii) 75 (iii) 15 () ( i ) (ii) 4 (iii) 7 1 ( () r, AOB = θ 0 < θ < ) OAB A OB P ( AB ) < ( AP ) (4) 0 < θ < sin θ < θ < tan θ 0 x, 0 y (1) sin x = sin y (x, y) () cos x cos y (x, y) 1 c
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(
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 information1 1.1 ( ). z = a + bi, a, b R 0 a, b 0 a 2 + b 2 0 z = a + bi = ( ) a 2 + b 2 a a 2 + b + b 2 a 2 + b i 2 r = a 2 + b 2 θ cos θ = a a 2 + b 2, sin θ =
1 1.1 ( ). z = + bi,, b R 0, b 0 2 + b 2 0 z = + bi = ( ) 2 + b 2 2 + b + b 2 2 + b i 2 r = 2 + b 2 θ cos θ = 2 + b 2, sin θ = b 2 + b 2 2π z = r(cos θ + i sin θ) 1.2 (, ). 1. < 2. > 3. ±,, 1.3 ( ). A
More information1 (Berry,1975) 2-6 p (S πr 2 )p πr 2 p 2πRγ p p = 2γ R (2.5).1-1 : : : : ( ).2 α, β α, β () X S = X X α X β (.1) 1 2
2005 9/8-11 2 2.2 ( 2-5) γ ( ) γ cos θ 2πr πρhr 2 g h = 2γ cos θ ρgr (2.1) γ = ρgrh (2.2) 2 cos θ θ cos θ = 1 (2.2) γ = 1 ρgrh (2.) 2 2. p p ρgh p ( ) p p = p ρgh (2.) h p p = 2γ r 1 1 (Berry,1975) 2-6
More information1 1 x y = y(x) y, y,..., y (n) : n y F (x, y, y,..., y (n) ) = 0 n F (x, y, y ) = 0 1 y(x) y y = G(x, y) y, y y + p(x)y = q(x) 1 p(x) q(
1 1 y = y() y, y,..., y (n) : n y F (, y, y,..., y (n) ) = 0 n F (, y, y ) = 0 1 y() 1.1 1 y y = G(, y) 1.1.1 1 y, y y + p()y = q() 1 p() q() (q() = 0) y + p()y = 0 y y + py = 0 y y = p (log y) = p log
More informationf(x,y) (x,y) x (x,y), y (x,y) f(x,y) x y f x (x,y),f y (x,y) B p.1/14
B p.1/14 f(x,y) (x,y) x (x,y), y (x,y) f(x,y) x y f x (x,y),f y (x,y) B p.1/14 f(x,y) (x,y) x (x,y), y (x,y) f(x,y) x y f x (x,y),f y (x,y) f(x 1,...,x n ) (x 1 x 0,...,x n 0), (x 1,...,x n ) i x i f xi
More informationC による数値計算法入門 ( 第 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 informationNo δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i x j δx j (5) δs 2
No.2 1 2 2 δs δs = r + δr r = δr (3) δs δs = r r = δr + u(r + δr, t) u(r, t) (4) δr = (δx, δy, δz) u i (r + δr, t) u i (r, t) = u i δx j (5) δs 2 = δx i δx i + 2 u i δx i δx j = δs 2 + 2s ij δx i δx j
More information2 1 Introduction (1.1.2) Logistic ث Malthus (1.1.3) (( ) ث)( ) α = ( ) ( + ) [Verhulst 1845] 0 ( ) ( + ) lim ( ) = 0 t (1.1.4) (( ) ث)( ) α = ( ) Logi
1 1 Introduction 1.1 ( ) ( ) Malthus ( ) (1.1.1) ( ) α = ( ) ( + ) 0 ( ) ( + ) lim ( ) = 0 t ( ) α = ( ) αt م = ( ) ( ) = ( 0 ) 0 (1.1.2) αt م 0 = ( ) (1.1.1) ( ) ( α + (1 = ( + ) = 0 1 ( ( ) 2... ) geometric
More information最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 3 版 1 刷発行時のものです.
最新耐震構造解析 ( 第 3 版 ) サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/052093 このサンプルページの内容は, 第 3 版 1 刷発行時のものです. i 3 10 3 2000 2007 26 8 2 SI SI 20 1996 2000 SI 15 3 ii 1 56 6
More information第3章 非線形計画法の基礎
3 February 25, 2009 1 Armijo Wolfe Newton 2 Newton Lagrange Newton 2 SQP 2 1 2.1 ( ) S R n (n N) f (x) : R n x f R x S f (x ) = min x S R n f (x) (nonlinear programming) x 0 S k = 0, 1, 2, h k R n ɛ k
More informationi 18 2H 2 + O 2 2H 2 + ( ) 3K
i 18 2H 2 + O 2 2H 2 + ( ) 3K ii 1 1 1.1.................................. 1 1.2........................................ 3 1.3......................................... 3 1.4....................................
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 informationDecember 28, 2018
e-mail : kigami@i.kyoto-u.ac.jp December 28, 28 Contents 2............................. 3.2......................... 7.3..................... 9.4................ 4.5............. 2.6.... 22 2 36 2..........................
More informationII ( ) (7/31) II ( [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Re
II 29 7 29-7-27 ( ) (7/31) II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ (6/29)] Navier Stokes 3 [ (6/19)] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I Euler Navier
More information( ) sin 1 x, cos 1 x, tan 1 x sin x, cos x, tan x, arcsin x, arccos x, arctan x. π 2 sin 1 x π 2, 0 cos 1 x π, π 2 < tan 1 x < π 2 1 (1) (
6 20 ( ) sin, cos, tan sin, cos, tan, arcsin, arccos, arctan. π 2 sin π 2, 0 cos π, π 2 < tan < π 2 () ( 2 2 lim 2 ( 2 ) ) 2 = 3 sin (2) lim 5 0 = 2 2 0 0 2 2 3 3 4 5 5 2 5 6 3 5 7 4 5 8 4 9 3 4 a 3 b
More informationẍ = kx, (k > ) (.) x x(t) = A cos(ωt + α) (.). d/ = D. d dt x + k ( x = D + k ) ( ) ( ) k k x = D + i D i x =... ( ) k D + i x = or ( ) k D i x =.. k.
K E N Z OU 8 9 8. F = kx x 3 678 ẍ = kx, (k > ) (.) x x(t) = A cos(ωt + α) (.). d/ = D. d dt x + k ( x = D + k ) ( ) ( ) k k x = D + i D i x =... ( ) k D + i x = or ( ) k D i x =.. k. D = ±i dt = ±iωx,
More information,. Black-Scholes u t t, x c u 0 t, x x u t t, x c u t, x x u t t, x + σ x u t, x + rx ut, x rux, t 0 x x,,.,. Step 3, 7,,, Step 6., Step 4,. Step 5,,.
9 α ν β Ξ ξ Γ γ o δ Π π ε ρ ζ Σ σ η τ Θ θ Υ υ ι Φ φ κ χ Λ λ Ψ ψ µ Ω ω Def, Prop, Th, Lem, Note, Remark, Ex,, Proof, R, N, Q, C [a, b {x R : a x b} : a, b {x R : a < x < b} : [a, b {x R : a x < b} : a,
More informationA S hara/lectures/lectures-j.html ϵ-n 1 ϵ-n lim n a n = α n a n α 2 lim a n = 0 1 n a k n n k= ϵ
A S1-20 http://www2.mth.kyushu-u.c.jp/ hr/lectures/lectures-j.html 1 1 1.1 ϵ-n 1 ϵ-n lim n n = α n n α 2 lim n = 0 1 n k n n k=1 0 1.1.7 ϵ-n 1.1.1 n α n n α lim n n = α ϵ N(ϵ) n > N(ϵ) n α < ϵ (1.1.1)
More information: , 2.0, 3.0, 2.0, (%) ( 2.
2017 1 2 1.1...................................... 2 1.2......................................... 4 1.3........................................... 10 1.4................................. 14 1.5..........................................
More information1 (1) () (3) I 0 3 I I d θ = L () dt θ L L θ I d θ = L = κθ (3) dt κ T I T = π κ (4) T I κ κ κ L l a θ L r δr δl L θ ϕ ϕ = rθ (5) l
1 1 ϕ ϕ ϕ S F F = ϕ (1) S 1: F 1 1 (1) () (3) I 0 3 I I d θ = L () dt θ L L θ I d θ = L = κθ (3) dt κ T I T = π κ (4) T I κ κ κ L l a θ L r δr δl L θ ϕ ϕ = rθ (5) l : l r δr θ πrδr δf (1) (5) δf = ϕ πrδr
More informationwebkaitou.dvi
( c Akir KANEKO) ).. m. l s = lθ m d s dt = mg sin θ d θ dt = g l sinθ θ l θ mg. d s dt xy t ( d x dt, d y dt ) t ( mg sin θ cos θ, sin θ sin θ). (.) m t ( d x dt, d y dt ) = t ( mg sin θ cos θ, mg sin
More informationtakei.dvi
0 Newton Leibniz ( ) α1 ( ) αn (1) a α1,...,α n (x) u(x) = f(x) x 1 x n α 1 + +α n m 1957 Hans Lewy Lewy 1970 1 1.1 Example 1.1. (2) d 2 u dx 2 Q(x)u = f(x), u(0) = a, 1 du (0) = b. dx Q(x), f(x) x = 0
More informationNumerical Analysis II, Exam End Term Spring 2017
H. Ammari W. Wu S. Yu Spring Term 2017 Numerical Analysis II ETH Zürich D-MATH End Term Spring 2017 Problem 1 Consider dx = f(t, x), t [0, T ] dt x(0) = x 0 R [28 Marks] with f C subject to the Lipschitz
More informationIA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................
More information9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x
2009 9 6 16 7 1 7.1 1 1 1 9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x(cos y y sin y) y dy 1 sin
More informationA
A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................
More informationTOP URL 1
TOP URL http://amonphys.web.fc.com/ 1 19 3 19.1................... 3 19.............................. 4 19.3............................... 6 19.4.............................. 8 19.5.............................
More information20 6 4 1 4 1.1 1.................................... 4 1.1.1.................................... 4 1.1.2 1................................ 5 1.2................................... 7 1.2.1....................................
More information2014 S hara/lectures/lectures-j.html r 1 S phone: ,
14 S1-1+13 http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html r 1 S1-1+13 14.4.11. 19 phone: 9-8-4441, e-mail: hara@math.kyushu-u.ac.jp Office hours: 1 4/11 web download. I. 1. ϵ-δ 1. 3.1, 3..
More informationf(x) = f(x ) + α(x)(x x ) α(x) x = x. x = f (y), x = f (y ) y = f f (y) = f f (y ) + α(f (y))(f (y) f (y )) f (y) = f (y ) + α(f (y)) (y y ) ( (2) ) f
22 A 3,4 No.3 () (2) (3) (4), (5) (6) (7) (8) () n x = (x,, x n ), = (,, n ), x = ( (x i i ) 2 ) /2 f(x) R n f(x) = f() + i α i (x ) i + o( x ) α,, α n g(x) = o( x )) lim x g(x) x = y = f() + i α i(x )
More information30
3 ............................................2 2...........................................2....................................2.2...................................2.3..............................
More informationI A A441 : April 21, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka ) Google
I4 - : April, 4 Version :. Kwhir, Tomoki TA (Kondo, Hirotk) Google http://www.mth.ngoy-u.c.jp/~kwhir/courses/4s-biseki.html pdf 4 4 4 4 8 e 5 5 9 etc. 5 6 6 6 9 n etc. 6 6 6 3 6 3 7 7 etc 7 4 7 7 8 5 59
More informationI, II 1, A = A 4 : 6 = max{ A, } A A 10 10%
1 2006.4.17. A 3-312 tel: 092-726-4774, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office hours: B A I ɛ-δ ɛ-δ 1. 2. A 1. 1. 2. 3. 4. 5. 2. ɛ-δ 1. ɛ-n
More information2000年度『数学展望 I』講義録
2000 I I IV I II 2000 I I IV I-IV. i ii 3.10 (http://www.math.nagoya-u.ac.jp/ kanai/) 2000 A....1 B....4 C....10 D....13 E....17 Brouwer A....21 B....26 C....33 D....39 E. Sperner...45 F....48 A....53
More informationk m m d2 x i dt 2 = f i = kx i (i = 1, 2, 3 or x, y, z) f i σ ij x i e ij = 2.1 Hooke s law and elastic constants (a) x i (2.1) k m σ A σ σ σ σ f i x
k m m d2 x i dt 2 = f i = kx i (i = 1, 2, 3 or x, y, z) f i ij x i e ij = 2.1 Hooke s law and elastic constants (a) x i (2.1) k m A f i x i B e e e e 0 e* e e (2.1) e (b) A e = 0 B = 0 (c) (2.1) (d) e
More information2.2 ( y = y(x ( (x 0, y 0 y (x 0 (y 0 = y(x 0 y = y(x ( y (x 0 = F (x 0, y(x 0 = F (x 0, y 0 (x 0, y 0 ( (x 0, y 0 F (x 0, y 0 xy (x, y (, F (x, y ( (
(. x y y x f y = f(x y x y = y(x y x y dx = d dx y(x = y (x = f (x y = y(x x ( (differential equation ( + y 2 dx + xy = 0 dx = xy + y 2 2 2 x y 2 F (x, y = xy + y 2 y = y(x x x xy(x = F (x, y(x + y(x 2
More informationII A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka )
II 214-1 : October 2, 214 Version : 1.1 Kawahira, Tomoki TA (Kondo, Hirotaka ) http://www.math.nagoya-u.ac.jp/~kawahira/courses/14w-biseki.html pdf 1 2 1 9 1 16 1 23 1 3 11 6 11 13 11 2 11 27 12 4 12 11
More information1
1 1 7 1.1.................................. 11 2 13 2.1............................ 13 2.2............................ 17 2.3.................................. 19 3 21 3.1.............................
More informationTOSM kanie/tosm/ HP kanie/agora/ kanie/ h
1 1.1.......................................................... 3 3 : 4 4 4 4.1.............................. 4 4............................... 8 4.3....................... 10 5 13 5.1.....................................
More informationall.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 informationOnsager SOLUTION OF THE EIGENWERT PROBLEM (O-29) V = e H A e H B λ max Z 2 Onsager (O-77) (O-82) (O-83) Kramers-Wannier 1 1 Ons
Onsager 2 9 207.2.7 3 SOLUTION OF THE EIGENWERT PROBLEM O-29 V = e H A e H B λ max Z 2 OnsagerO-77O-82 O-83 2 Kramers-Wannier Onsager * * * * * V self-adjoint V = V /2 V V /2 = V /2 V 2 V /2 = 2 sinh 2H
More information(yx4) 1887-1945 741936 50 1995 1 31 http://kenboushoten.web.fc.com/ OCR TeX 50 yx4 e-mail: yx4.aydx5@gmail.com i Jacobi 1751 1 3 Euler Fagnano 187 9 0 Abel iii 1 1...................................
More informationJanuary 27, 2015
e-mail : kigami@i.kyoto-u.ac.jp January 27, 205 Contents 2........................ 2.2....................... 3.3....................... 6.4......................... 2 6 2........................... 6
More information1 911 9001030 9:00 A B C D E F G H I J K L M 1A0900 1B0900 1C0900 1D0900 1E0900 1F0900 1G0900 1H0900 1I0900 1J0900 1K0900 1L0900 1M0900 9:15 1A0915 1B0915 1C0915 1D0915 1E0915 1F0915 1G0915 1H0915 1I0915
More informationv v = v 1 v 2 v 3 (1) R = (R ij ) (2) R (R 1 ) ij = R ji (3) 3 R ij R ik = δ jk (4) i=1 δ ij Kronecker δ ij = { 1 (i = j) 0 (i
1. 1 1.1 1.1.1 1.1.1.1 v v = v 1 v 2 v 3 (1) R = (R ij ) (2) R (R 1 ) ij = R ji (3) R ij R ik = δ jk (4) δ ij Kronecker δ ij = { 1 (i = j) 0 (i j) (5) 1 1.1. v1.1 2011/04/10 1. 1 2 v i = R ij v j (6) [
More information2 (2016 3Q N) c = o (11) Ax = b A x = c A n I n n n 2n (A I n ) (I n X) A A X A n A A A (1) (2) c 0 c (3) c A A i j n 1 ( 1) i+j A (i, j) A (i, j) ã i
[ ] (2016 3Q N) a 11 a 1n m n A A = a m1 a mn A a 1 A A = a n (1) A (a i a j, i j ) (2) A (a i ca i, c 0, i ) (3) A (a i a i + ca j, j i, i ) A 1 A 11 0 A 12 0 0 A 1k 0 1 A 22 0 0 A 2k 0 1 0 A 3k 1 A rk
More informationSO(3) 49 u = Ru (6.9), i u iv i = i u iv i (C ) π π : G Hom(V, V ) : g D(g). π : R 3 V : i 1. : u u = u 1 u 2 u 3 (6.10) 6.2 i R α (1) = 0 cos α
SO(3) 48 6 SO(3) t 6.1 u, v u = u 1 1 + u 2 2 + u 3 3 = u 1 e 1 + u 2 e 2 + u 3 e 3, v = v 1 1 + v 2 2 + v 3 3 = v 1 e 1 + v 2 e 2 + v 3 e 3 (6.1) i (e i ) e i e j = i j = δ ij (6.2) ( u, v ) = u v = ij
More informationB ver B
B ver. 2017.02.24 B Contents 1 11 1.1....................... 11 1.1.1............. 11 1.1.2.......................... 12 1.2............................. 14 1.2.1................ 14 1.2.2.......................
More information2011de.dvi
211 ( 4 2 1. 3 1.1............................... 3 1.2 1- -......................... 13 1.3 2-1 -................... 19 1.4 3- -......................... 29 2. 37 2.1................................ 37
More informationDOPRI5.dvi
ODE DOPRI5 ( ) 16 3 31 Runge Kutta Dormand Prince 5(4) [1, pp. 178 179] DOPRI5 http://www.unige.ch/math/folks/hairer/software.html Fortran C C++ [3, pp.51 56] DOPRI5 C cprog.tar % tar xvf cprog.tar cprog/
More informationD 24 D D D
5 Paper I.R. 2001 5 Paper HP Paper 5 3 5.1................................................... 3 5.2.................................................... 4 5.3.......................................... 6
More information18 2 F 12 r 2 r 1 (3) Coulomb km Coulomb M = kg F G = ( ) ( ) ( ) 2 = [N]. Coulomb
r 1 r 2 r 1 r 2 2 Coulomb Gauss Coulomb 2.1 Coulomb 1 2 r 1 r 2 1 2 F 12 2 1 F 21 F 12 = F 21 = 1 4πε 0 1 2 r 1 r 2 2 r 1 r 2 r 1 r 2 (2.1) Coulomb ε 0 = 107 4πc 2 =8.854 187 817 10 12 C 2 N 1 m 2 (2.2)
More informationIntroduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math))
Introduction to Numerical Analysis of Differential Equations Naoya Enomoto (Kyoto.univ.Dept.Science(math)) 2001 1 e-mail:s00x0427@ip.media.kyoto-u.ac.jp 1 1 Van der Pol 1 1 2 2 Bergers 2 KdV 2 1 5 1.1........................................
More information1.2 y + P (x)y + Q(x)y = 0 (1) y 1 (x), y 2 (x) y 1 (x), y 2 (x) (1) y(x) c 1, c 2 y(x) = c 1 y 1 (x) + c 2 y 2 (x) 3 y 1 (x) y 1 (x) e R P (x)dx y 2
1 1.1 R(x) = 0 y + P (x)y + Q(x)y = R(x)...(1) y + P (x)y + Q(x)y = 0...(2) 1 2 u(x) v(x) c 1 u(x)+ c 2 v(x) = 0 c 1 = c 2 = 0 c 1 = c 2 = 0 2 0 2 u(x) v(x) u(x) u (x) W (u, v)(x) = v(x) v (x) 0 1 1.2
More information( ) ( 40 )+( 60 ) Schrödinger 3. (a) (b) (c) yoshioka/education-09.html pdf 1
2009 1 ( ) ( 40 )+( 60 ) 1 1. 2. Schrödinger 3. (a) (b) (c) http://goofy.phys.nara-wu.ac.jp/ yoshioka/education-09.html pdf 1 1. ( photon) ν λ = c ν (c = 3.0 108 /m : ) ɛ = hν (1) p = hν/c = h/λ (2) h
More information