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 =

Size: px
Start display at page:

Download "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 ="

Transcription

1 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+ 3.3 d dt d dt d dt [A] = k[a] [B] = k[a] k[b] [C] = k[b] 4

2 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 = P x P x = = C C exp( t) C 3 exp( t) A B C A + B + C = C A + B = C exp( t) A = C 3 exp( t) 5

3 A 0 B 0 C 0 = C 3 = C = C = mm 0 mm 0 mm A + B + C = A + B = exp( t) A = exp( t) A B C = exp( t) exp( t) exp( t) exp( t) + exp( t) mm : Euler Euler Taylor y(x 0 + h) = y(x 0 ) + h y (x 0 ) + h! y (x 0 ) + Euler y(x 0 + h) = y(x 0 ) + h y (x 0 ) Euler Taylor t S S(t) t t + t S S(t + t) Euler S(t + t) = S(t) + ts (t) S (t) S ds dt 6

4 S(t + t) = S(t) + t ds dt ds dt ds = k S dt S(t + t) = S(t) t k S S(t) = S(t + t) S(t) = t k S Euler t S t k S d[glucose] = k[glucose] dt t =.0sec t = sec k = 0.sec t = 0 sec [Glucose]=.0mM. () t = 5.0 sec ε = n X analytical i Xi numerical n X analytical i i= 3. t.0 sec 4. 7

5 3.4.3 : Euler α 0 y n + α y n + + α k y n k = h(β 0 f n + β f n + + β m f n m ) y n t n y y n k = y(t n kh) f = y 3.5. Euler RK 3.5. : Adams-Bashforth Gear Adams y n = α y n + h(β f n + β f n ) α, β Taylor y(t n ) = α y(t n ) + h(β f(t n, y(t n )) + β f(t n, y(t n ))) t n = t h Taylor y(t n ) α { y(t n ) h d dt y(t n) + h + hβ { f(t n ) h d dt f(t n) + h + hβ { f(t n ) h d dt f(t n) + 4h } d dt y(t n) d dt f(t n) } d dt f(t n) } 8

6 h h 0 : y(t n ) = α y(t n ) h : 0 = α hy (t n ) + hβ y (t n ) + hβ y (t n ) h : 0 = h α y (t n ) h β y (t n ) h β y (t n ) α, β = α 0 = α + β + β 0 = α β β Adams-Bashforth α =, β = 3, β = Adams-Bashforth ds dt S = t ( 3 y n = y n + h f n + ) f n = k S ( 3 k S t=t t ) k S t=t t : Adams-Moulton (predictor-corrector) Adams-Bashforth Adams-Moulton MATLAB ode3 ( ) ( ) Backward Differential Formula(ode5s) 9

7 Runge-Kutta 3.7 Runge-Kutta Runge-Kutta 3.7. LM LM RK Adams 3.7. s Runge-Kutta k = f(x 0, y 0 ) k = f(x 0 + c h, y 0 + ha k ) k 3 = f(x 0 + c 3 h, y 0 + h(a 3 k + a 3 k )). k s = f(x 0 + c s h, y 0 + h(a s k + a s,s k s )) y = y 0 + h(b k + + b s k s ) 0

8 c = a c 3 = a 3 + a 3. c s = a s + + a s,s s Runge-Kutta Butcher s Runge-Kutta Butcher k 0 k c a k 3 c 3 a 3 a k s c s a s a s a s,s y b b b s b s Runge-Kutta (RK) Runge-Kutta k = f(x 0, y 0 ) k = f(x 0 + h, y 0 + hk ) y = y 0 + hk Runge-Kutta y(x 0 + h) y Kh p+ Runge-Kutta p Runge-Kutta Taylor p

9 3.7.6 Taylor Euler y = y 0 + hf ( x 0 + h, y 0 + h ) f 0 Taylor f ( ( ) x 0 + h, y 0 + hf 0) = h f(x0, y 0 ) + + hf x 0 f(x y 0, y 0 ) ( + h + hf! x 0 y) f(x0, y 0 ) + = f(x 0, y 0 ) + h (f x + ff y ) (x 0, y 0 ) + ( h 4 + h f x 0 x + h f y 4 0 y ) f(x 0, y 0 ) + = f(x 0, y 0 ) + h (f x + ff y ) (x 0, y 0 ) + h (f 8 xx + ff xy + f f yy ) (x 0, y 0 ) + y = y 0 + hf(x 0, y 0 ) + h (f x + f y f) (x 0, y 0 ) + h3 (f 8 xx + f xy f + f yy f ) (x 0, y 0 ) + Taylor y(x 0 + h) = y 0 + hy (x 0, y 0 ) + h y (x 0, y 0 ) + h3 6 y (x 0, y 0 ) + = y 0 + hf(x 0, y 0 ) + h x + f y f)(x 0, y 0 ) + h3 6 (f xx + f xy f + f yy f + f x f y + f y f)(x 0, y 0 ) + y(x 0 + h) y = h3 4 {(f xx + f xy f + f yy f + 4(f x f y + f y f)}(x 0, y 0 ) + y(x 0 + h) y Kh 3 Runge-Kutta

10 3.7.7 dy dx = f(x, y) df = f x df dx = f x + f y = y y = f x + f y f dx + f y dy dy dx dy = y x dy dx dx + y y dy = y x + y y = dy dx x (f x + f y f) + y (f x + f y f) dy dx = {f xx + (f xy f + f y f x )} + {f xy + (f yy f + fy )}f y = f xx + f xy f + f yy f + f x f y + f y f p(ˆp) p ˆp ( : Dormand-Prince 5(4), Bogacki-Shampine 3() ode45, ode3 ) Dormand-Prince Runge-Kutta(5,4) 3.8. RK3() RK

11 k = f(x n, y n ) () k = f(x n + c h, y n + ha k ) () k 3 = f(x n + c h, y n + h(a 3 k + a 3 k )) (3) y n+ = y n + h(b k + b k + b 3 k 3 ) (4) k 3 h 3 k (4) h h ( ) k = f(x n, y n ) + h c! x k f y ( + h c! x k y) f + O(h 3 ) = f(x n, y n ) + h(c f x + a ff y ) + h (c f xx + c a ff xy + a f f yy ) + O(h 3 ) k 3 k 3 = f(x n, y n ) + h! + h! ( c 3 ) f + (a x 3k + a 3 k ) y ( c 3 + (a x 3k + a 3 k ) y) f + O(h 3 ) = f(x n, y n ) + h(c 3 f x + a 3 ff y + a 3 k f y ) + h {c 3f xx + c 3 (a 3 f + a 3 k )f xy + (a 3 f + a 3 k ) f yy + O(h 3 ) k f a 3 k a 3 k = a 3 { f + h(c f x + a ff y ) + O(h ) } a 3 k = a 3 (f + O(h)) h k 3 k 3 = f(x n, y n ) + h[c 3 f x + a 3 ff y + a 3 {f + h(c f x + a ff y ) + O(h )} f y ] + h [c 3f xx + c 3 {a 3 f + a 3 (f + O(h))}f xy + {a 3 f + a 3 (f + O(h))} f yy ] RK h RK 4

12 y(x + h) = y(x) + h dy dx + h d y! dx + h3 d 3 y 3! dx + 3 O(h4 ) y(x + h) y(x) = hf + h (f x + f y f) + h3 6 (f xx + f xy f + f yy f + f y f x + fy f) + O(h 4 ) d y dx = f x + f y f df dx = f x + f df y dx RK h RK : hf = h(b + b + b 3 )f : 3 : h (f x + f y f) = b h (c f x + a ff y ) + b 3 h (c 3 f x + a 3 ff y + a 3 ff y ) h 3 6 (f xx + f xy f + f yy f + f y f x + fy f) = b h 3 (c f xx + c a ff xy + a f f yy ) + b 3 h 3 a 3 (c f x + a ff y )f y + b 3h 3 {c 3f xx h +c 3 a 3 ff xy + c 3 a 3 ff xy +(a 3 + a 3 ) f f yy } b + b + b 3 = h h f x : f x = b h c f x + b 3 h c 3 f x = b c + b 3 c 3 h ff y : f yf = b h a ff y + b 3 h (a 3 + a 3 )ff y = b a + b 3 (a 3 + a 3 ) h 3 5

13 h 3 f xx : 6 f xx = b h 3 c f xx + b 3h 3 c 3f xx = 6 b c + b 3c 3 h 3 ff xy : 6 (f xyf) = b h 3 (c a ff xy ) + b 3h 3 (c 3 a 3 ff xy + c 3 a 3 ff xy ) = b c a + b 3 c 3 (a 3 + a 3 ) 3 f yy f h 3 : 6 f yyf = b h 3 a f f yy + b 3h 3 (a 3 + a 3 ) f f yy = 6 b a + b 3(a 3 + a 3 ) h 3 f y f x : 6 f yf x = b 3 h 3 a 3 c f x f y = b 3 a 3 c 6 fy h 3 f : 6 f y f = b 3 h 3 a 3 a ffy = b 3 a 3 a 6 c = a c 3 = a 3 + a 3 b + b + b 3 = b c + b 3 c 3 = b c + b 3 c 3 = 3 b 3 a 3 c = 6 8 (a, a 3, a 3, b, b, b 3, c, c 3 ) 6 c = u, c 3 = v { b u + b 3 v = b u + b 3 v = 3 { b u + b 3 uv = u b u + b 3 v = 3 6

14 (uv v )b 3 = u 3 b 3 = = u 3 uv v 3u 6v(u v) { b uv + b 3 v = v b u + b 3 v = 3 (uv u )b = v 3 3v b = 6u(v u) b 3 a 3 c = 6 3u 6v(u v) a 3u = 6 a 3 = v(u v) u(3u ) b = b b 3 a = c = u a 3 = c 3 a 3 = v a 3 (RK) Runge-Kutta b, b, b 3, ŷ n+ = y n + h(ˆb k + + ˆb s k s ) Runge- Kutta 7

15 ε tol h ( ) RK3() ˆb 3 = 0 y i ŷ i ε tol b u + b 3 v = ˆb 3 = 0 ˆb = ˆb ˆb u = ˆb = = u c ˆb = c (RK) 0 0 c =, a =, ˆb = 0, ˆb = RK3 0 c 3 a 3 a 3 b b b 3 RK3 Ralston (REF.FIXME) Ralston a 3 = 0, a 3 = 3, b 4 =, b 9 =, b 3 3 = Ralston RK3 4 9, c 3 = 3 4 (REF.FIXME) 4 9 8

16 3.8.3 y i ŷ i ε tol err = n n ( yi ŷ i ε i= tol ) err err Ch q+ Ch q+ opt err ( h h opt ) q+ h opt h( q+ err h opt. err h. err h opt h Gear (BDF: Backward Differential Formula) α 0 y n + α y n + + α k y n k = hf n Adams Taylor Gear Numerical initial value problems pp.4 stiffly stable method 9

17 L.H.S. α 0 y(t n ) + α { y(t n ) h d dt y(t n) + h + α { y(t n ) h d dt y(t n) + 4h R.H.S. = h d dt y(t n) } d dt y(t n) d dt y(t n) } h h 0 : α 0 y(t n ) + α y(t n ) + α y(t n ) = 0 h : α hy (t n ) hα y (t n ) = hy (t n ) h : α h y (t n ) + α h y (t n ) = 0 α 0 + α + α = 0 α α = + α = 0 Gear α 0 = 3, α =, α = Gear ds dt 3 y n y n + y n = hf n = k S 3 S t S t t + S t t = t k S t 3.9. (stiff equation) (stiffness ratio) = > 0 4 (stiff) 30

18 A BDF A (stiffly stable) 6 BDF 3.0 Further reading Gear, Springer s yellow book by Gear 3

p12.dvi

p12.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 information

Trapezoidal 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

Trapezoidal 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 information

DVIOUT

DVIOUT A. A. A-- [ ] f(x) x = f 00 (x) f 0 () =0 f 00 () > 0= f(x) x = f 00 () < 0= f(x) x = A--2 [ ] f(x) D f 00 (x) > 0= y = f(x) f 00 (x) < 0= y = f(x) P (, f()) f 00 () =0 A--3 [ ] y = f(x) [, b] x = f (y)

More information

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

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 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 86 6 r (6) y y d y = y 3 (64) y

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

error_g1.eps

error_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 information

.1 z = e x +xy y z y 1 1 x 0 1 z x y α β γ z = αx + βy + γ (.1) ax + by + cz = d (.1') a, b, c, d x-y-z (a, b, c). x-y-z 3 (0,

.1 z = e x +xy y z y 1 1 x 0 1 z x y α β γ z = αx + βy + γ (.1) ax + by + cz = d (.1') a, b, c, d x-y-z (a, b, c). x-y-z 3 (0, .1.1 Y K L Y = K 1 3 L 3 L K K (K + ) 1 1 3 L 3 K 3 L 3 K 0 (K + K) 1 3 L 3 K 1 3 L 3 lim K 0 K = L (K + K) 1 3 K 1 3 3 lim K 0 K = 1 3 K 3 L 3 z = f(x, y) x y z x-y-z.1 z = e x +xy y 3 x-y ( ) z 0 f(x,

More information

I, II 1, A = A 4 : 6 = max{ A, } A A 10 10%

I, 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 information

A

A 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

211 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

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) 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 information

1 28 6 12 7 1 7.1...................................... 2 7.1.1............................... 2 7.1.2........................... 2 7.2...................................... 3 7.3...................................

More information

untitled

untitled 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 information

dy = 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

dy = 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

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

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 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 information

1 4 1 ( ) ( ) ( ) ( ) () 1 4 2

1 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 information

CALCULUS II (Hiroshi SUZUKI ) f(x, y) A(a, b) 1. P (x, y) A(a, b) A(a, b) f(x, y) c f(x, y) A(a, b) c f(x, y) c f(x, y) c (x a, y b)

CALCULUS II (Hiroshi SUZUKI ) f(x, y) A(a, b) 1. P (x, y) A(a, b) A(a, b) f(x, y) c f(x, y) A(a, b) c f(x, y) c f(x, y) c (x a, y b) CALCULUS II (Hiroshi SUZUKI ) 16 1 1 1.1 1.1 f(x, y) A(a, b) 1. P (x, y) A(a, b) A(a, b) f(x, y) c f(x, y) A(a, b) c f(x, y) c f(x, y) c (x a, y b) lim f(x, y) = lim f(x, y) = lim f(x, y) = c. x a, y b

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

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P

9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P 9 (Finite Element Method; FEM) 9. 9. P(0) P(x) u(x) (a) P(L) f P(0) P(x) (b) 9. P(L) 9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L)

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

2 1 κ c(t) = (x(t), y(t)) ( ) det(c (t), c x (t)) = det (t) x (t) y (t) y = x (t)y (t) x (t)y (t), (t) c (t) = (x (t)) 2 + (y (t)) 2. c (t) =

2 1 κ c(t) = (x(t), y(t)) ( ) det(c (t), c x (t)) = det (t) x (t) y (t) y = x (t)y (t) x (t)y (t), (t) c (t) = (x (t)) 2 + (y (t)) 2. c (t) = 1 1 1.1 I R 1.1.1 c : I R 2 (i) c C (ii) t I c (t) (0, 0) c (t) c(i) c c(t) 1.1.2 (1) (2) (3) (1) r > 0 c : R R 2 : t (r cos t, r sin t) (2) C f : I R c : I R 2 : t (t, f(t)) (3) y = x c : R R 2 : t (t,

More information

( ) Loewner SLE 13 February

( ) Loewner SLE 13 February ( ) Loewner SLE 3 February 00 G. F. Lawler, Conformally Invariant Processes in the Plane, (American Mathematical Society, 005)., Summer School 009 (009 8 7-9 ) . d- (BES d ) d B t = (Bt, B t,, Bd t ) (d

More information

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

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 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 information

17 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 information

( ) ( )

( ) ( ) 20 21 2 8 1 2 2 3 21 3 22 3 23 4 24 5 25 5 26 6 27 8 28 ( ) 9 3 10 31 10 32 ( ) 12 4 13 41 0 13 42 14 43 0 15 44 17 5 18 6 18 1 1 2 2 1 2 1 0 2 0 3 0 4 0 2 2 21 t (x(t) y(t)) 2 x(t) y(t) γ(t) (x(t) y(t))

More information

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,. 24(2012) (1 C106) 4 11 (2 C206) 4 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 (). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5... 6.. 7.,,. 8.,. 1. (75%)

More information

5.. z = f(x, y) y y = b f x x g(x) f(x, b) g x ( ) A = lim h 0 g(a + h) g(a) h g(x) a A = g (a) = f x (a, b)

5.. z = f(x, y) y y = b f x x g(x) f(x, b) g x ( ) A = lim h 0 g(a + h) g(a) h g(x) a A = g (a) = f x (a, b) 5 partial differentiation (total) differentiation 5. z = f(x, y) (a, b) A = lim h 0 f(a + h, b) f(a, b) h............................................................... ( ) f(x, y) (a, b) x A (a, b) x

More information

.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =,

.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =, [ ] IC. r, θ r, θ π, y y = 3 3 = r cos θ r sin θ D D = {, y ; y }, y D r, θ ep y yddy D D 9 s96. d y dt + 3dy + y = cos t dt t = y = e π + e π +. t = π y =.9 s6.3 d y d + dy d + y = y =, dy d = 3 a, b

More information

II (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 (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 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)

() 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 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

f : 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

f : 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

統計的データ解析

統計的データ解析 ds45 xspec qdp guplot oocalc (Error) gg (Radom Error)(Systematc Error) x, x,, x ( x, x,..., x x = s x x µ = lm = σ µ x x = lm ( x ) = σ ( ) = - x = js j ( ) = j= ( j) x x + xj x + xj j x + xj = ( x x

More information

,.,. 2, R 2, ( )., I R. c : I R 2, : (1) c C -, (2) t I, c (t) (0, 0). c(i). c (t)., c(t) = (x(t), y(t)) c (t) = (x (t), y (t)) : (1)

,.,. 2, R 2, ( )., I R. c : I R 2, : (1) c C -, (2) t I, c (t) (0, 0). c(i). c (t)., c(t) = (x(t), y(t)) c (t) = (x (t), y (t)) : (1) ( ) 1., : ;, ;, ; =. ( ).,.,,,., 2.,.,,.,.,,., y = f(x), f ( ).,,.,.,., U R m, F : U R n, M, f : M R p M, p,, R m,,, R m. 2009 A tamaru math.sci.hiroshima-u.ac.jp 1 ,.,. 2, R 2, ( ).,. 2.1 2.1. I R. c

More information

: 1g99p038-8

: 1g99p038-8 16 17 : 1g99p038-8 1 3 1.1....................................... 4 1................................... 5 1.3.................................. 5 6.1..................................... 7....................................

More information

I No. sin cos sine, cosine : trigonometric function π : π =.4 : n =, ±, ±, sin + nπ = sin cos + nπ = cos sin = sin : cos = cos :. sin. sin. sin + π si

I No. sin cos sine, cosine : trigonometric function π : π =.4 : n =, ±, ±, sin + nπ = sin cos + nπ = cos sin = sin : cos = cos :. sin. sin. sin + π si I 8 No. : No. : No. : No.4 : No.5 : No.6 : No.7 : No.8 : No.9 : No. : I No. sin cos sine, cosine : trigonometric function π : π =.4 : n =, ±, ±, sin + nπ = sin cos + nπ = cos sin = sin : cos = cos :. sin.

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

ax 2 + bx + c = n 8 (n ) a n x n + a n 1 x n a 1 x + a 0 = 0 ( a n, a n 1,, a 1, a 0 a n 0) n n ( ) ( ) ax 3 + bx 2 + cx + d = 0 4

ax 2 + bx + c = n 8 (n ) a n x n + a n 1 x n a 1 x + a 0 = 0 ( a n, a n 1,, a 1, a 0 a n 0) n n ( ) ( ) ax 3 + bx 2 + cx + d = 0 4 20 20.0 ( ) 8 y = ax 2 + bx + c 443 ax 2 + bx + c = 0 20.1 20.1.1 n 8 (n ) a n x n + a n 1 x n 1 + + a 1 x + a 0 = 0 ( a n, a n 1,, a 1, a 0 a n 0) n n ( ) ( ) ax 3 + bx 2 + cx + d = 0 444 ( a, b, c, d

More information

(1) (2) (3) (4) HB B ( ) (5) (6) (7) 40 (8) (9) (10)

(1) (2) (3) (4) HB B ( ) (5) (6) (7) 40 (8) (9) (10) 2017 12 9 4 1 30 4 10 3 1 30 3 30 2 1 30 2 50 1 1 30 2 10 (1) (2) (3) (4) HB B ( ) (5) (6) (7) 40 (8) (9) (10) (1) i 23 c 23 0 1 2 3 4 5 6 7 8 9 a b d e f g h i (2) 23 23 (3) 23 ( 23 ) 23 x 1 x 2 23 x

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

ma22-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

ma22-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 information

23 7 28 i i 1 1 1.1................................... 2 1.2............................... 3 1.2.1.................................... 3 1.2.2............................... 4 1.2.3 SI..............................

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

D xy D (x, y) z = f(x, y) f D (2 ) (x, y, z) f R z = 1 x 2 y 2 {(x, y); x 2 +y 2 1} x 2 +y 2 +z 2 = 1 1 z (x, y) R 2 z = x 2 y

D xy D (x, y) z = f(x, y) f D (2 ) (x, y, z) f R z = 1 x 2 y 2 {(x, y); x 2 +y 2 1} x 2 +y 2 +z 2 = 1 1 z (x, y) R 2 z = x 2 y 5 5. 2 D xy D (x, y z = f(x, y f D (2 (x, y, z f R 2 5.. z = x 2 y 2 {(x, y; x 2 +y 2 } x 2 +y 2 +z 2 = z 5.2. (x, y R 2 z = x 2 y + 3 (2,,, (, 3,, 3 (,, 5.3 (. (3 ( (a, b, c A : (x, y, z P : (x, y, x

More information

IA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................

More information

II 1 II 2012 II Gauss-Bonnet II

II 1 II 2012 II Gauss-Bonnet II II 1 II 212 II Gauss-Bonnet II 1 1 1.1......................................... 1 1.2............................................ 2 1.3.................................. 3 1.4.............................................

More information

Morse ( ) 2014

Morse ( ) 2014 Morse ( ) 2014 1 1 Morse 1 1.1 Morse................................ 1 1.2 Morse.............................. 7 2 12 2.1....................... 12 2.2.................. 13 2.3 Smale..............................

More information

W 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)

W 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 information

u = u(t, x 1,..., x d ) : R R d C λ i = 1 := x 2 1 x 2 d d Euclid Laplace Schrödinger N := {1, 2, 3,... } Z := {..., 3, 2, 1,, 1, 2, 3

u = u(t, x 1,..., x d ) : R R d C λ i = 1 := x 2 1 x 2 d d Euclid Laplace Schrödinger N := {1, 2, 3,... } Z := {..., 3, 2, 1,, 1, 2, 3 2 2 1 5 5 Schrödinger i u t + u = λ u 2 u. u = u(t, x 1,..., x d ) : R R d C λ i = 1 := 2 + + 2 x 2 1 x 2 d d Euclid Laplace Schrödinger 3 1 1.1 N := {1, 2, 3,... } Z := {..., 3, 2, 1,, 1, 2, 3,... } Q

More information

ii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,.

ii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,. (1 C205) 4 10 (2 C206) 4 11 (2 B202) 4 12 25(2013) http://www.math.is.tohoku.ac.jp/~obata,.,,,..,,. 1. 2. 3. 4. 5. 6. 7. 8. 1., 2007 ( ).,. 2. P. G., 1995. 3. J. C., 1988. 1... 2.,,. ii 3.,. 4. F. ( ),..

More information

7. y fx, z gy z gfx dz dx dz dy dy dx. g f a g bf a b fa 7., chain ule Ω, D R n, R m a Ω, f : Ω R m, g : D R l, fω D, b fa, f a g b g f a g f a g bf a

7. y fx, z gy z gfx dz dx dz dy dy dx. g f a g bf a b fa 7., chain ule Ω, D R n, R m a Ω, f : Ω R m, g : D R l, fω D, b fa, f a g b g f a g f a g bf a 9 203 6 7 WWW http://www.math.meiji.ac.jp/~mk/lectue/tahensuu-203/ 2 8 8 7. 7 7. y fx, z gy z gfx dz dx dz dy dy dx. g f a g bf a b fa 7., chain ule Ω, D R n, R m a Ω, f : Ω R m, g : D R l, fω D, b fa,

More information

No δ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 δ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 information

. p.1/14

. p.1/14 . p.1/14 F(x,y) = (F 1 (x,y),f 2 (x,y)) (x,y). p.2/14 F(x,y) = (F 1 (x,y),f 2 (x,y)) (x,y) (x,y) h. p.2/14 F(x,y) = (F 1 (x,y),f 2 (x,y)) (x,y) (x,y) h h { F 2 (x+ h,y) F 2 2(x h,y) F 2 1(x,y+ h)+f 2 1(x,y

More information

v er.1/ c /(21)

v er.1/ c /(21) 12 -- 1 1 2009 1 17 1-1 1-2 1-3 1-4 2 2 2 1-5 1 1-6 1 1-7 1-1 1-2 1-3 1-4 1-5 1-6 1-7 c 2011 1/(21) 12 -- 1 -- 1 1--1 1--1--1 1 2009 1 n n α { n } α α { n } lim n = α, n α n n ε n > N n α < ε N {1, 1,

More information

A B P (A B) = P (A)P (B) (3) A B A B P (B A) A B A B P (A B) = P (B A)P (A) (4) P (B A) = P (A B) P (A) (5) P (A B) P (B A) P (A B) A B P

A B P (A B) = P (A)P (B) (3) A B A B P (B A) A B A B P (A B) = P (B A)P (A) (4) P (B A) = P (A B) P (A) (5) P (A B) P (B A) P (A B) A B P 1 1.1 (population) (sample) (event) (trial) Ω () 1 1 Ω 1.2 P 1. A A P (A) 0 1 0 P (A) 1 (1) 2. P 1 P 0 1 6 1 1 6 0 3. A B P (A B) = P (A) + P (B) (2) A B A B A 1 B 2 A B 1 2 1 2 1 1 2 2 3 1.3 A B P (A

More information

x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x i ) A(x i ) = h 2 {f(x i) + f(x i+1 ) = h {f(a + i h) + f(a + (i + 1) h), (2) 2 a b n A(x i )

x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x i ) A(x i ) = h 2 {f(x i) + f(x i+1 ) = h {f(a + i h) + f(a + (i + 1) h), (2) 2 a b n A(x i ) 1 f(x) a b f(x)dx = n A(x i ) (1) ix [a, b] n i A(x i ) x i 1 f(x) [a, b] n h = (b a)/n y h = (b-a)/n y = f (x) h h a a+h a+2h a+(n-1)h b x 1: 1 x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x

More information

18 ( ) I II III A B C(100 ) 1, 2, 3, 5 I II A B (100 ) 1, 2, 3 I II A B (80 ) 6 8 I II III A B C(80 ) 1 n (1 + x) n (1) n C 1 + n C

18 ( ) I II III A B C(100 ) 1, 2, 3, 5 I II A B (100 ) 1, 2, 3 I II A B (80 ) 6 8 I II III A B C(80 ) 1 n (1 + x) n (1) n C 1 + n C 8 ( ) 8 5 4 I II III A B C( ),,, 5 I II A B ( ),, I II A B (8 ) 6 8 I II III A B C(8 ) n ( + x) n () n C + n C + + n C n = 7 n () 7 9 C : y = x x A(, 6) () A C () C P AP Q () () () 4 A(,, ) B(,, ) C(,,

More information

untitled

untitled 20010916 22;1017;23;20020108;15;20; 1 N = {1, 2, } Z + = {0, 1, 2, } Z = {0, ±1, ±2, } Q = { p p Z, q N} R = { lim a q n n a n Q, n N; sup a n < } R + = {x R x 0} n = {a + b 1 a, b R} u, v 1 R 2 2 R 3

More information

R R 16 ( 3 )

R R 16   ( 3 ) (017 ) 9 4 7 ( ) ( 3 ) ( 010 ) 1 (P3) 1 11 (P4) 1 1 (P4) 1 (P15) 1 (P16) (P0) 3 (P18) 3 4 (P3) 4 3 4 31 1 5 3 5 4 6 5 9 51 9 5 9 6 9 61 9 6 α β 9 63 û 11 64 R 1 65 13 66 14 7 14 71 15 7 R R 16 http://wwwecoosaka-uacjp/~tazak/class/017

More information

5. [1 ] 1 [], u(x, t) t c u(x, t) x (5.3) ξ x + ct, η x ct (5.4),u(x, t) ξ, η u(ξ, η), ξ t,, ( u(ξ,η) ξ η u(x, t) t ) u(x, t) { ( u(ξ, η) c t ξ ξ { (

5. [1 ] 1 [], u(x, t) t c u(x, t) x (5.3) ξ x + ct, η x ct (5.4),u(x, t) ξ, η u(ξ, η), ξ t,, ( u(ξ,η) ξ η u(x, t) t ) u(x, t) { ( u(ξ, η) c t ξ ξ { ( 5 5.1 [ ] ) d f(t) + a d f(t) + bf(t) : f(t) 1 dt dt ) u(x, t) c u(x, t) : u(x, t) t x : ( ) ) 1 : y + ay, : y + ay + by : ( ) 1 ) : y + ay, : yy + ay 3 ( ): ( ) ) : y + ay, : y + ay b [],,, [ ] au xx

More information

renshumondai-kaito.dvi

renshumondai-kaito.dvi 3 1 13 14 1.1 1 44.5 39.5 49.5 2 0.10 2 0.10 54.5 49.5 59.5 5 0.25 7 0.35 64.5 59.5 69.5 8 0.40 15 0.75 74.5 69.5 79.5 3 0.15 18 0.90 84.5 79.5 89.5 2 0.10 20 1.00 20 1.00 2 1.2 1 16.5 20.5 12.5 2 0.10

More information

数値計算:常微分方程式

数値計算:常微分方程式 ( ) 1 / 82 1 2 3 4 5 6 ( ) 2 / 82 ( ) 3 / 82 C θ l y m O x mg λ ( ) 4 / 82 θ t C J = ml 2 C mgl sin θ θ C J θ = mgl sin θ = θ ( ) 5 / 82 ω = θ J ω = mgl sin θ ω J = ml 2 θ = ω, ω = g l sin θ = θ ω ( )

More information

x (x, ) x y (, y) iy x y z = x + iy (x, y) (r, θ) r = x + y, θ = tan ( y ), π < θ π x r = z, θ = arg z z = x + iy = r cos θ + ir sin θ = r(cos θ + i s

x (x, ) x y (, y) iy x y z = x + iy (x, y) (r, θ) r = x + y, θ = tan ( y ), π < θ π x r = z, θ = arg z z = x + iy = r cos θ + ir sin θ = r(cos θ + i s ... x, y z = x + iy x z y z x = Rez, y = Imz z = x + iy x iy z z () z + z = (z + z )() z z = (z z )(3) z z = ( z z )(4)z z = z z = x + y z = x + iy ()Rez = (z + z), Imz = (z z) i () z z z + z z + z.. z

More information

I No. sin cos sine, cosine : trigonometric function π : π =.4 : n = 0, ±, ±, sin + nπ = sin cos + nπ = cos : parity sin = sin : odd cos = cos : even.

I No. sin cos sine, cosine : trigonometric function π : π =.4 : n = 0, ±, ±, sin + nπ = sin cos + nπ = cos : parity sin = sin : odd cos = cos : even. I 0 No. : No. : No. : No.4 : No.5 : No.6 : No.7 : No.8 : No.9 : No.0 : I No. sin cos sine, cosine : trigonometric function π : π =.4 : n = 0, ±, ±, sin + nπ = sin cos + nπ = cos : parity sin = sin : odd

More information

1 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 = 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 information

1 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω 1 ω α V T m T m 1 100Hz m 2 36km 500Hz. 36km 1

1 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω 1 ω α V T m T m 1 100Hz m 2 36km 500Hz. 36km 1 sin cos P (primary) S (secondly) 2 P S A sin(ω2πt + α) A ω ω α 3 3 2 2V 3 33+.6T m T 5 34m Hz. 34 3.4m 2 36km 5Hz. 36km m 34 m 5 34 + m 5 33 5 =.66m 34m 34 x =.66 55Hz, 35 5 =.7 485.7Hz 2 V 5Hz.5V.5V V

More information

5.. z = f(x, y) y y = b f x x g(x) f(x, b) g x ( ) A = lim h g(a + h) g(a) h g(x) a A = g (a) = f x (a, b)............................................

5.. z = f(x, y) y y = b f x x g(x) f(x, b) g x ( ) A = lim h g(a + h) g(a) h g(x) a A = g (a) = f x (a, b)............................................ 5 partial differentiation (total) differentiation 5. z = f(x, y) (a, b) A = lim h f(a + h, b) f(a, b) h........................................................... ( ) f(x, y) (a, b) x A (a, b) x (a, b)

More information

) Binary Cubic Forms / 25

) Binary Cubic Forms / 25 2016 5 2 ) Binary Cubic Forms 2016 5 2 1 / 25 1 2 2 2 3 2 3 ) Binary Cubic Forms 2016 5 2 2 / 25 1.1 ( ) 4 2 12 = 5+7, 16 = 5+11, 36 = 7+29, 1.2 ( ) p p+2 3 5 5 7 11 13 17 19, 29 31 41 43 ) Binary Cubic

More information

6.1 (P (P (P (P (P (P (, P (, P.101

6.1 (P (P (P (P (P (P (, P (, P.101 (008 0 3 7 ( ( ( 00 1 (P.3 1 1.1 (P.3.................. 1 1. (P.4............... 1 (P.15.1 (P.15................. (P.18............3 (P.17......... 3.4 (P................ 4 3 (P.7 4 3.1 ( P.7...........

More information

K E N Z OU

K E N Z OU K E N Z OU 11 1 1 1.1..................................... 1.1.1............................ 1.1..................................................................................... 4 1.........................................

More information

II 2 II

II 2 II II 2 II 2005 yugami@cc.utsunomiya-u.ac.jp 2005 4 1 1 2 5 2.1.................................... 5 2.2................................. 6 2.3............................. 6 2.4.................................

More information

Untitled

Untitled II 14 14-7-8 8/4 II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ 6/ ] Navier Stokes 3 [ ] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I 1 balance law t (ρv i )+ j

More information

2.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 ( (

2.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 information

6.1 (P (P (P (P (P (P (, P (, P.

6.1 (P (P (P (P (P (P (, P (, P. (011 30 7 0 ( ( 3 ( 010 1 (P.3 1 1.1 (P.4.................. 1 1. (P.4............... 1 (P.15.1 (P.16................. (P.0............3 (P.18 3.4 (P.3............... 4 3 (P.9 4 3.1 (P.30........... 4 3.

More information

..3. Ω, Ω F, P Ω, F, P ). ) F a) A, A,..., A i,... F A i F. b) A F A c F c) Ω F. ) A F A P A),. a) 0 P A) b) P Ω) c) [ ] A, A,..., A i,... F i j A i A

..3. Ω, Ω F, P Ω, F, P ). ) F a) A, A,..., A i,... F A i F. b) A F A c F c) Ω F. ) A F A P A),. a) 0 P A) b) P Ω) c) [ ] A, A,..., A i,... F i j A i A .. Laplace ). A... i),. ω i i ). {ω,..., ω } Ω,. ii) Ω. Ω. A ) r, A P A) P A) r... ).. Ω {,, 3, 4, 5, 6}. i i 6). A {, 4, 6} P A) P A) 3 6. ).. i, j i, j) ) Ω {i, j) i 6, j 6}., 36. A. A {i, j) i j }.

More information

1 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 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 information

DE-resume

DE-resume - 2011, http://c-faculty.chuo-u.ac.jp/ nishioka/ 2 11 21131 : 4 1 x y(x, y (x,y (x,,y (n, (1.1 F (x, y, y,y,,y (n =0. (1.1 n. (1.1 y(x. y(x (1.1. 1 1 1 1.1... 2 1.2... 9 1.3 1... 26 2 2 34 2.1,... 35 2.2

More information

[ ] 0.1 lim x 0 e 3x 1 x IC ( 11) ( s114901) 0.2 (1) y = e 2x (x 2 + 1) (2) y = x/(x 2 + 1) 0.3 dx (1) 1 4x 2 (2) e x sin 2xdx (3) sin 2 xdx ( 11) ( s

[ ] 0.1 lim x 0 e 3x 1 x IC ( 11) ( s114901) 0.2 (1) y = e 2x (x 2 + 1) (2) y = x/(x 2 + 1) 0.3 dx (1) 1 4x 2 (2) e x sin 2xdx (3) sin 2 xdx ( 11) ( s [ ]. lim e 3 IC ) s49). y = e + ) ) y = / + ).3 d 4 ) e sin d 3) sin d ) s49) s493).4 z = y z z y s494).5 + y = 4 =.6 s495) dy = 3e ) d dy d = y s496).7 lim ) lim e s49).8 y = e sin ) y = sin e 3) y =

More information

(u(x)v(x)) = u (x)v(x) + u(x)v (x) ( ) u(x) = u (x)v(x) u(x)v (x) v(x) v(x) 2 y = g(t), t = f(x) y = g(f(x)) dy dx dy dx = dy dt dt dx., y, f, g y = f (g(x))g (x). ( (f(g(x)). ). [ ] y = e ax+b (a, b )

More information

6. Euler x

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

More information

2011de.dvi

2011de.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 information

A 2 3. m S m = {x R m+1 x = 1} U + k = {x S m x k > 0}, U k = {x S m x k < 0}, ϕ ± k (x) = (x 0,..., ˆx k,... x m ) 1. {(U ± k, ϕ± k ) 0 k m} S m 1.2.

A 2 3. m S m = {x R m+1 x = 1} U + k = {x S m x k > 0}, U k = {x S m x k < 0}, ϕ ± k (x) = (x 0,..., ˆx k,... x m ) 1. {(U ± k, ϕ± k ) 0 k m} S m 1.2. A A 1 A 5 A 6 1 2 3 4 5 6 7 1 1.1 1.1 (). Hausdorff M R m M M {U α } U α R m E α ϕ α : U α E α U α U β = ϕ α (ϕ β ϕβ (U α U β )) 1 : ϕ β (U α U β ) ϕ α (U α U β ) C M a m dim M a U α ϕ α {x i, 1 i m} {U,

More information

sin cos No. sine, cosine : trigonometric function π : π = 3.4 : n = 0, ±, ±, sin + nπ = sin cos + nπ = cos : parity sin = sin : odd cos = cos : even.

sin cos No. sine, cosine : trigonometric function π : π = 3.4 : n = 0, ±, ±, sin + nπ = sin cos + nπ = cos : parity sin = sin : odd cos = cos : even. 08 No. : No. : No.3 : No.4 : No.5 : No.6 : No.7 : No.8 : No.9 : No.0 : No. : sin cos No. sine, cosine : trigonometric function π : π = 3.4 : n = 0, ±, ±, sin + nπ = sin cos + nπ = cos : parity sin = sin

More information

応力とひずみ.ppt

応力とひずみ.ppt in yukawa@numse.nagoya-u.ac.jp 2 3 4 5 x 2 6 Continuum) 7 8 9 F F 10 F L L F L 1 L F L F L F 11 F L F F L F L L L 1 L 2 12 F L F! A A! S! = F S 13 F L L F F n = F " cos# F t = F " sin# S $ = S cos# S S

More information

x ( ) x dx = ax

x ( ) x dx = ax x ( ) x dx = ax 1 dx = a x log x = at + c x(t) = e at C (C = e c ) a > 0 t a < 0 t 0 (at + b ) h dx = lim x(t + h) x(t) h 0 h x(t + h) x(t) h x(t) t x(t + h) x(t) ax(t) h x(t + h) x(t) + ahx(t) 0, h, 2h,

More information

VI VI.21 W 1,..., W r V W 1,..., W r W W r = {v v r v i W i (1 i r)} V = W W r V W 1,..., W r V W 1,..., W r V = W 1 W

VI VI.21 W 1,..., W r V W 1,..., W r W W r = {v v r v i W i (1 i r)} V = W W r V W 1,..., W r V W 1,..., W r V = W 1 W 3 30 5 VI VI. W,..., W r V W,..., W r W + + W r = {v + + v r v W ( r)} V = W + + W r V W,..., W r V W,..., W r V = W W r () V = W W r () W (W + + W + W + + W r ) = {0} () dm V = dm W + + dm W r VI. f n

More information

変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy,

変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, 変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, z + dz) Q! (x + d x + u + du, y + dy + v + dv, z +

More information

mugensho.dvi

mugensho.dvi 1 1 f (t) lim t a f (t) = 0 f (t) t a 1.1 (1) lim(t 1) 2 = 0 t 1 (t 1) 2 t 1 (2) lim(t 1) 3 = 0 t 1 (t 1) 3 t 1 2 f (t), g(t) t a lim t a f (t) g(t) g(t) f (t) = o(g(t)) (t a) = 0 f (t) (t 1) 3 1.2 lim

More information

II Karel Švadlenka * [1] 1.1* 5 23 m d2 x dt 2 = cdx kx + mg dt. c, g, k, m 1.2* u = au + bv v = cu + dv v u a, b, c, d R

II Karel Švadlenka * [1] 1.1* 5 23 m d2 x dt 2 = cdx kx + mg dt. c, g, k, m 1.2* u = au + bv v = cu + dv v u a, b, c, d R II Karel Švadlenka 2018 5 26 * [1] 1.1* 5 23 m d2 x dt 2 = cdx kx + mg dt. c, g, k, m 1.2* 5 23 1 u = au + bv v = cu + dv v u a, b, c, d R 1.3 14 14 60% 1.4 5 23 a, b R a 2 4b < 0 λ 2 + aλ + b = 0 λ =

More information

1 nakayama/print/ Def (Definition ) Thm (Theorem ) Prop (Proposition ) Lem (Lemma ) Cor (Corollary ) 1. (1) A, B (2) ABC

1   nakayama/print/ Def (Definition ) Thm (Theorem ) Prop (Proposition ) Lem (Lemma ) Cor (Corollary ) 1. (1) A, B (2) ABC 1 http://www.gem.aoyama.ac.jp/ nakayama/print/ Def (Definition ) Thm (Theorem ) Prop (Proposition ) Lem (Lemma ) Cor (Corollary ) 1. (1) A, B (2) ABC r 1 A B B C C A (1),(2),, (8) A, B, C A,B,C 2 1 ABC

More information

lecture

lecture 5 3 3. 9. 4. x, x. 4, f(x, ) :=x x + =4,x,.. 4 (, 3) (, 5) (3, 5), (4, 9) 95 9 (g) 4 6 8 (cm).9 3.8 6. 8. 9.9 Phsics 85 8 75 7 65 7 75 8 85 9 95 Mathematics = ax + b 6 3 (, 3) 3 ( a + b). f(a, b) ={3 (a

More information

1 : f(z = re iθ ) = u(r, θ) + iv(r, θ). (re iθ ) 2 = r 2 e 2iθ = r 2 cos 2θ + ir 2 sin 2θ r f(z = x + iy) = u(x, y) + iv(x, y). (x + iy) 2 = x 2 y 2 +

1 : f(z = re iθ ) = u(r, θ) + iv(r, θ). (re iθ ) 2 = r 2 e 2iθ = r 2 cos 2θ + ir 2 sin 2θ r f(z = x + iy) = u(x, y) + iv(x, y). (x + iy) 2 = x 2 y 2 + 1.3 1.4. (pp.14-27) 1 1 : f(z = re iθ ) = u(r, θ) + iv(r, θ). (re iθ ) 2 = r 2 e 2iθ = r 2 cos 2θ + ir 2 sin 2θ r f(z = x + iy) = u(x, y) + iv(x, y). (x + iy) 2 = x 2 y 2 + i2xy x = 1 y (1 + iy) 2 = 1

More information

1. x { e 1,..., e n } x = x1 e1 + + x n en = (x 1,..., x n ) X, Y [X, Y ] Intrinsic ( ) Intrinsic M m P M C P P M P M v 3 v : C P R 1

1. x { e 1,..., e n } x = x1 e1 + + x n en = (x 1,..., x n ) X, Y [X, Y ] Intrinsic ( ) Intrinsic M m P M C P P M P M v 3 v : C P R 1 1. x { e 1,..., e n } x = x1 e1 + + x n en = (x 1,..., x n ) X, Y [X, Y ] Intrinsic ( ) Intrinsic M m P M C P P M P M v 3 v : C P R 1 f, g C P, λ R (1) v(f + g) = v(f) + v(g) (2) v(λf) = λv(f) (3) v(fg)

More information

Chap9.dvi

Chap9.dvi .,. f(),, f(),,.,. () lim 2 +3 2 9 (2) lim 3 3 2 9 (4) lim ( ) 2 3 +3 (5) lim 2 9 (6) lim + (7) lim (8) lim (9) lim (0) lim 2 3 + 3 9 2 2 +3 () lim sin 2 sin 2 (2) lim +3 () lim 2 2 9 = 5 5 = 3 (2) lim

More information

2 A id A : A A A A id A def = {(a, a) A A a A} 1 { } 1 1 id 1 = α: A B β : B C α β αβ : A C αβ def = {(a, c) A C b B.((a, b) α (b, c) β)} 2.3 α

2 A id A : A A A A id A def = {(a, a) A A a A} 1 { } 1 1 id 1 = α: A B β : B C α β αβ : A C αβ def = {(a, c) A C b B.((a, b) α (b, c) β)} 2.3 α 20 6 18 1 2 2.1 A B α A B α: A B A B Rel(A, B) A B (A B) A B 0 AB A B AB α, β : A B α β α β def (a, b) A B.((a, b) α (a, b) β) 0 AB AB Rel(A, B) 1 2 A id A : A A A A id A def = {(a, a) A A a A} 1 { } 1

More information

untitled

untitled 1 n m (ICA = independent component analysis) BSS (= blind source separation) : s(t) =(s 1 (t),...,s n (t)) R n : x(t) =(x 1 (t),...,x n (t)) R m 1 i s i (t) a ji R j 2 (A =(a ji )) x(t) =As(t) (1) n =

More information

2 (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

2 (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 information

f (x) f (x) f (x) f (x) f (x) 2 f (x) f (x) f (x) f (x) 2 n f (x) n f (n) (x) dn f f (x) dx n dn dx n D n f (x) n C n C f (x) x = a 1 f (x) x = a x >

f (x) f (x) f (x) f (x) f (x) 2 f (x) f (x) f (x) f (x) 2 n f (x) n f (n) (x) dn f f (x) dx n dn dx n D n f (x) n C n C f (x) x = a 1 f (x) x = a x > 5.1 1. x = a f (x) a x h f (a + h) f (a) h (5.1) h 0 f (x) x = a f +(a) f (a + h) f (a) = lim h +0 h (5.2) x h h 0 f (a) f (a + h) f (a) f (a h) f (a) = lim = lim h 0 h h 0 h (5.3) f (x) x = a f (a) =

More information