error_g1.eps
|
|
- えりか いくのや
- 5 years ago
- Views:
Transcription
1 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 n x n ) h x(t n+ ) x(t n )+hf(t n x(t n )) x n+ = x n +hf(t n x n ) x n (explicit method) x n x n+ Jun 6 Version: naito@mathnagoya-uacjp
2 t = t n t = t n+ 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 n+ x n+ ) h x(t n+ ) x(t n )+hf(t n+ x(t n+ )) x n+ = x n +hf(t n+ x n+ ) (backward Euler method) x n x n+ x n (implicit method) x (t) = x(t) x n+ = x n +hx n+ x n x n+ x n+ x n+ = h x n x (t) = sin(x(t)) x n+ = x n +hsin(x n+ ) x n x n+ x n+ y = x n y k+ = x n +hsin(y k ) x n+ h Jun 6 Version: naito@mathnagoya-uacjp
3 x (t) = f(x(t)) x (t) = f(tx(t)) X(t) = X(t) x = f(x(t)) x x(t n+ ) = x(t n )+hx (t n )+ h x (t n )+O(h ) = x(t n )+hf(x(t n ))+ h f(x (t n ))f(x(t n ))+O(h ) x n+ = x n +hf(x n ) x(t n+ ) x n+ = x(t n ) x n +h(f(x(t n )) f(x n ))+ h f (x(t n ))f(x(t n ))+O(h ) f f(p) f(q) L p q f x(t n+ ) x n+ (+Lh) x(t n ) x n +Mh x(t n ) = x(t n+ ) hx (t n+ )+ h x (t n+ )+O(h ) = x(t n+ ) hf(x(t n+ ))+ h f(x (t n+ ))f(x(t n+ ))+O(h ) x n+ = x n +hf(x n+ ) x(t n+ ) x n+ = x(t n ) x n +h(f(x(t n+ )) f(x n+ ))+ h f (x(t n+ ))f(x(t n+ ))+O(h ) x(t n+ ) x n+ ( Lh) x(t n ) x n +Mh ǫ n+ = Aǫ n +Mh ǫ n = x(t n ) x n t = t n x n x(t n ) O(h ) Jun 6 Version: naito@mathnagoya-uacjp
4 4 t = T = Nh x(t) x N Che TL x(t) x N Ch h h h x(t) x N = O(h p ) p > x(t n+ ) = x(t n )+hx (t n )+ h x (t n )+O(h ) = x(t n )+hf(x(t n ))+ h f(x (t n ))f(x(t n ))+O(h ) h k k = x n+ = x n +hf(x n )+ h f (x n )f(x n ) O(h ) (h k O(h k+ ) k f x(t n+ ) = x(t n )+hx (t n )+ h x (t n )+ h 6 x() (t)+o(h 4 ) x(t n ) = x(t n ) hx (t n )+ h x (t n ) h 6 x() (t)+o(h 4 ) x(t n+ ) x(t n ) = hx (t n )+ h x() (t)+o(h 5 ) x n+ = x n +hf(x n ) O(h ) Jun 6 Version: naito@mathnagoya-uacjp
5 5 x (t) = x(t) x() = x n+ = x n hx n x = x = +h {x n } x n n x(t) = e t x(t) (t ) x(t) x n x n+ = x n hx n t +ht = α α α < < α < x n = A α n +A α n x n A = x = C x = Cα A x n a n+ = a n +a n Jun 6 Version: naito@mathnagoya-uacjp
6 6 h = x = x x = -x x (t) = x(t) x() = x (t) = x(t) x() = x = (-x^) x = sin(x) x (t) = ( x(t) ) x() = x (t) = sin(x(t)) x() = x = -x x = -x x (t) = x(t) x() = x () = Jun 6 Version: naito@mathnagoya-uacjp
7 7 x = -sin(x) x = -sin(x) - - x (t) = sin(x(t)) x() = x () = x = (-x^)x - x x = (-x^)x - x x (t) = ( x(t) )x (t) x(t) x() = 5 x () = X = -X/ X ^ X (t) = X(t)/ X(t) X() = () X () = () Double Pendulum Double Pendulum x () = x () = x () = x () = Jun 6 Version: naito@mathnagoya-uacjp
8 8 h = x = x x = -x x (t) = x(t) x() = x (t) = x(t) x() = x = (-x^) x = sin(x) x (t) = ( x(t) ) x() = x (t) = sin(x(t)) x() = x = -x x = -x x (t) = x(t) x() = x () = Jun 6 Version: naito@mathnagoya-uacjp
9 9 x = -sin(x) x = -sin(x) - - x (t) = sin(x(t)) x() = x () = x = (-x^)x - x x = (-x^)x - x x (t) = ( x(t) )x (t) x(t) x() = 5 x () = X = -X/ X ^ X (t) = X(t)/ X(t) X() = () X () = () Double Pendulum Double Pendulum x () = x () = x () = x () = Jun 6 Version: naito@mathnagoya-uacjp
10 h = x = x x = -x x (t) = x(t) x() = x (t) = x(t) x() = x = (-x^) x = sin(x) x (t) = ( x(t) ) x() = x (t) = sin(x(t)) x() = x = -x x = -x x (t) = x(t) x() = x () = Jun 6 Version: naito@mathnagoya-uacjp
11 x = -sin(x) x = -sin(x) - - x (t) = sin(x(t)) x() = x () = x = (-x^)x - x x = (-x^)x - x x (t) = ( x(t) )x (t) x(t) x() = 5 x () = X = -X/ X ^ X (t) = X(t)/ X(t) X() = () X () = () Double Pendulum Double Pendulum x () = x () = x () = x () = Jun 6 Version: naito@mathnagoya-uacjp
12 x = ±x x() = t = t = x = x relative error at t = x = -x relative error at t = Forward Euler Backward Euler Improved Euler Heun Heun Kutta Runge-Kutta Runge-Kutta-Gill e-6 Forward Euler Backward Euler Improved Euler Heun Heun Kutta Runge-Kutta Runge-Kutta-Gill e-6 e-8 e-8 e- e- e- e- e-4 e-4 e-6 e+6 e+7 e+8 e+9 e-6 e+6 e+7 e+8 e+9 x = x relative error at t = x = -x relative error at t = Forward Euler Backward Euler Improved Euler Heun Heun Kutta Runge-Kutta Runge-Kutta-Gill e-6 Forward Euler Backward Euler Improved Euler Heun Heun Kutta Runge-Kutta Runge-Kutta-Gill e-6 e-8 e-8 e- e- e- e- e-4 e-4 e-6 e-6 e+6 e+7 e+8 e+9 e-8 e+6 e+7 e+8 e+9 x (t) = x(t) x() = x (t) = x(t) x() = h 9 Jun 6 Version: naito@mathnagoya-uacjp
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 informationd 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 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.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 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 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 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 information1 28 6 12 7 1 7.1...................................... 2 7.1.1............................... 2 7.1.2........................... 2 7.2...................................... 3 7.3...................................
More information(, ) (, ) S = 2 = [, ] ( ) 2 ( ) 2 2 ( ) 3 2 ( ) 4 2 ( ) k 2,,, k =, 2, 3, 4 S 4 S 4 = ( ) 2 + ( ) ( ) (
B 4 4 4 52 4/ 9/ 3/3 6 9.. y = x 2 x x = (, ) (, ) S = 2 = 2 4 4 [, ] 4 4 4 ( ) 2 ( ) 2 2 ( ) 3 2 ( ) 4 2 ( ) k 2,,, 4 4 4 4 4 k =, 2, 3, 4 S 4 S 4 = ( ) 2 + ( ) 2 2 + ( ) 3 2 + ( 4 4 4 4 4 4 4 4 4 ( (
More information86 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 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 informationCALCULUS 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[] x < T f(x), x < T f(x), < x < f(x) f(x) f(x) f(x + nt ) = f(x) x < T, n =, 1,, 1, (1.3) f(x) T x 2 f(x) T 2T x 3 f(x), f() = f(t ), f(x), f() f(t )
1 1.1 [] f(x) f(x + T ) = f(x) (1.1), f(x), T f(x) x T 1 ) f(x) = sin x, T = 2 sin (x + 2) = sin x, sin x 2 [] n f(x + nt ) = f(x) (1.2) T [] 2 f(x) g(x) T, h 1 (x) = af(x)+ bg(x) 2 h 2 (x) = f(x)g(x)
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 informationt θ, τ, α, β S(, 0 P sin(θ P θ S x cos(θ SP = θ P (cos(θ, sin(θ sin(θ P t tan(θ θ 0 cos(θ tan(θ = sin(θ cos(θ ( 0t tan(θ
4 5 ( 5 3 9 4 0 5 ( 4 6 7 7 ( 0 8 3 9 ( 8 t θ, τ, α, β S(, 0 P sin(θ P θ S x cos(θ SP = θ P (cos(θ, sin(θ sin(θ P t tan(θ θ 0 cos(θ tan(θ = sin(θ cos(θ ( 0t tan(θ S θ > 0 θ < 0 ( P S(, 0 θ > 0 ( 60 θ
More informationf(x) = x (1) f (1) (2) f (2) f(x) x = a y y = f(x) f (a) y = f(x) A(a, f(a)) f(a + h) f(x) = A f(a) A x (3, 3) O a a + h x 1 f(x) x = a
3 3.1 3.1.1 A f(a + h) f(a) f(x) lim f(x) x = a h 0 h f(x) x = a f 0 (a) f 0 (a) = lim h!0 f(a + h) f(a) h = lim x!a f(x) f(a) x a a + h = x h = x a h 0 x a 3.1 f(x) = x x = 3 f 0 (3) f (3) = lim h 0 (
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 informationlecture
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) (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 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 information18 ( ) 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名古屋工業大の数学 2000 年 ~2015 年 大学入試数学動画解説サイト
名古屋工業大の数学 年 ~5 年 大学入試数学動画解説サイト http://mathroom.jugem.jp/ 68 i 4 3 III III 3 5 3 ii 5 6 45 99 5 4 3. () r \= S n = r + r + 3r 3 + + nr n () x > f n (x) = e x + e x + 3e 3x + + ne nx f(x) = lim f n(x) lim
More informationx 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 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 informationx A Aω ẋ ẋ 2 + ω 2 x 2 = ω 2 A 2. (ẋ, ωx) ζ ẋ + iωx ζ ζ dζ = ẍ + iωẋ = ẍ + iω(ζ iωx) dt dζ dt iωζ = ẍ + ω2 x (2.1) ζ ζ = Aωe iωt = Aω cos ωt + iaω sin
2 2.1 F (t) 2.1.1 mẍ + kx = F (t). m ẍ + ω 2 x = F (t)/m ω = k/m. 1 : (ẋ, x) x = A sin ωt, ẋ = Aω cos ωt 1 2-1 x A Aω ẋ ẋ 2 + ω 2 x 2 = ω 2 A 2. (ẋ, ωx) ζ ẋ + iωx ζ ζ dζ = ẍ + iωẋ = ẍ + iω(ζ iωx) dt dζ
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 information1 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 informationDVIOUT
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 informationx i [, b], (i 0, 1, 2,, n),, [, b], [, b] [x 0, x 1 ] [x 1, x 2 ] [x n 1, x n ] ( 2 ). x 0 x 1 x 2 x 3 x n 1 x n b 2: [, b].,, (1) x 0, x 1, x 2,, x n
1, R f : R R,.,, b R < b, f(x) [, b] f(x)dx,, [, b] f(x) x ( ) ( 1 ). y y f(x) f(x)dx b x 1: f(x)dx, [, b] f(x) x ( ).,,,,,., f(x)dx,,,, f(x)dx. 1.1 Riemnn,, [, b] f(x) x., x 0 < x 1 < x 2 < < x n 1
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 informationII 1 3 2 5 3 7 4 8 5 11 6 13 7 16 8 18 2 1 1. x 2 + xy x y (1 lim (x,y (1,1 x 1 x 3 + y 3 (2 lim (x,y (, x 2 + y 2 x 2 (3 lim (x,y (, x 2 + y 2 xy (4 lim (x,y (, x 2 + y 2 x y (5 lim (x,y (, x + y x 3y
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 informationr 1 m A r/m i) t ii) m i) t B(t; m) ( B(t; m) = A 1 + r ) mt m ii) B(t; m) ( B(t; m) = A 1 + r ) mt m { ( = A 1 + r ) m } rt r m n = m r m n B
1 1.1 1 r 1 m A r/m i) t ii) m i) t Bt; m) Bt; m) = A 1 + r ) mt m ii) Bt; m) Bt; m) = A 1 + r ) mt m { = A 1 + r ) m } rt r m n = m r m n Bt; m) Aert e lim 1 + 1 n 1.1) n!1 n) e a 1, a 2, a 3,... {a n
More informationslide1.dvi
1. 2/ 121 a x = a t 3/ 121 a x = a t 4/ 121 a > 0 t a t = a t t {}}{ a a a t 5/ 121 a t+s = = t+s {}}{ a a a t s {}}{{}}{ a a a a = a t a s (a t ) s = s {}}{ a t a t = a ts 6/ 121 a > 0 t a 0 t t = 0 +
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 information5.. 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春期講座 ~ 極限 1 1, 1 2, 1 3, 1 4,, 1 n, n n {a n } n a n α {a n } α {a n } α lim n an = α n a n α α {a n } {a n } {a n } 1. a n = 2 n {a n } 2, 4, 8, 16,
春期講座 ~ 極限 1 1, 1 2, 1 3, 1 4,, 1 n, n n {a n } n a n α {a n } α {a n } α lim an = α n a n α α {a n } {a n } {a n } 1. a n = 2 n {a n } 2, 4, 8, 16, 32, n a n {a n } {a n } 2. a n = 10n + 1 {a n } lim an
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 informationv_-3_+2_1.eps
I 9-9 (3) 9 9, x, x (t)+a(t)x (t)+b(t)x(t) = f(t) (9), a(t), b(t), f(t),,, f(t),, a(t), b(t),,, x (t)+ax (t)+bx(t) = (9),, x (t)+ax (t)+bx(t) = f(t) (93), b(t),, b(t) 9 x (t), x (t), x (t)+a(t)x (t)+b(t)x(t)
More information4 4 4 a b c d a b A c d A a da ad bce O E O n A n O ad bc a d n A n O 5 {a n } S n a k n a n + k S n a a n+ S n n S n n log x x {xy } x, y x + y 7 fx
4 4 5 4 I II III A B C, 5 7 I II A B,, 8, 9 I II A B O A,, Bb, b, Cc, c, c b c b b c c c OA BC P BC OP BC P AP BC n f n x xn e x! e n! n f n x f n x f n x f k x k 4 e > f n x dx k k! fx sin x cos x tan
More informationM3 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 informationFX ) 2
(FX) 1 1 2009 12 12 13 2009 1 FX ) 2 1 (FX) 2 1 2 1 2 3 2010 8 FX 1998 1 FX FX 4 1 1 (FX) () () 1998 4 1 100 120 1 100 120 120 100 20 FX 100 100 100 1 100 100 100 1 100 1 100 100 1 100 101 101 100 100
More informationA A = a 41 a 42 a 43 a 44 A (7) 1 (3) A = M 12 = = a 41 (8) a 41 a 43 a 44 (3) n n A, B a i AB = A B ii aa
1 2 21 2 2 [ ] a 11 a 12 A = a 21 a 22 (1) A = a 11 a 22 a 12 a 21 (2) 3 3 n n A A = n ( 1) i+j a ij M ij i =1 n (3) j=1 M ij A i j (n 1) (n 1) 2-1 3 3 A A = a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33
More information( ) 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 informationBlack-Scholes [1] Nelson [2] Schrödinger 1 Black Scholes [1] Black-Scholes Nelson [2][3][4] Schrödinger Nelson Parisi Wu [5] Nelson Parisi-W
003 7 14 Black-Scholes [1] Nelson [] Schrödinger 1 Black Scholes [1] Black-Scholes Nelson [][3][4] Schrödinger Nelson Parisi Wu [5] Nelson Parisi-Wu Nelson e-mail: takatoshi-tasaki@nifty.com kabutaro@mocha.freemail.ne.jp
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 informationMH MH 9.50 8.50 9.40 8.40 9.30 9.20 8.30 9.10 9.00 8.20 8.90 8.80 8.10 8.70 8.60 8.00 8.50 7.90 8.40 8.30 7.80 8.20 8.10 7.70 8.00 7.60 7.90 7.80 7.50 7.70 7.60 7.40 1 7.50 7.40 7.30 2 7.30 7.20 7.20
More informationTrv1_.xtv
第 0 回 PBF 選手権日本代表選抜試合ウィメン予選 ession 1 01 年 11 月 1 日 ( 土 ) 1 Board (1) - -W Contract L - -W IMP H D H C 0 Board () - -W Contract L - -W IMP H -1 D 0 H -1 D 0 Board () - -W Contract L - -W IMP T -1 D 0 D
More information2012 IA 8 I p.3, 2 p.19, 3 p.19, 4 p.22, 5 p.27, 6 p.27, 7 p
2012 IA 8 I 1 10 10 29 1. [0, 1] n x = 1 (n = 1, 2, 3,...) 2 f(x) = n 0 [0, 1] 2. 1 x = 1 (n = 1, 2, 3,...) 2 f(x) = n 0 [0, 1] 1 0 f(x)dx 3. < b < c [, c] b [, c] 4. [, b] f(x) 1 f(x) 1 f(x) [, b] 5.
More informationVersion C 1 2 3 4 5 1 2 3 4 5 6 7 8 9 0 A 1 2 1 3 4 5 1 1 2 1 1 1 2 4 5 6 7 8 3 1 2 C a b c d e f g A A B C B a b c d e f g 3 4 4 5 6 7 8 1 2 a b 1 2 a b 1 2 1 2 5 4 1 23 5 6 6 a b 1 2 e c d 3
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 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 informationx () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x
[ ] IC. f(x) = e x () f(x) f (x) () lim f(x) lim f(x) x + x (3) lim f(x) lim f(x) x + x (4) y = f(x) ( ) ( s46). < a < () a () lim a log xdx a log xdx ( ) n (3) lim log k log n n n k=.3 z = log(x + y ),
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 (
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平成20年5月 協会創立50年の歩み 海の安全と環境保全を目指して 友國八郎 海上保安庁 長官 岩崎貞二 日本船主協会 会長 前川弘幸 JF全国漁業協同組合連合会 代表理事会長 服部郁弘 日本船長協会 会長 森本靖之 日本船舶機関士協会 会長 大内博文 航海訓練所 練習船船長 竹本孝弘 第二管区海上保安本部長 梅田宜弘
More information
1 c Koichi Suga, ISBN
c Koichi Suga, 4 4 6 5 ISBN 978-4-64-6445- 4 ( ) x(t) t u(t) t {u(t)} {x(t)} () T, (), (3), (4) max J = {u(t)} V (x, u)dt ẋ = f(x, u) x() = x x(t ) = x T (), x, u, t ẋ x t u u ẋ = f(x, u) x(t ) = x T x(t
More informationNo.004 [1] J. ( ) ( ) (1968) [2] Morse (1997) [3] (1988) 1
No.004 [1] J. ( ) ( ) (1968) [2] Morse (1997) [3] (1988) 1 1 (1) 1.1 X Y f, g : X Y { F (x, 0) = f(x) F (x, 1) = g(x) F : X I Y f g f g F f g 1.2 X Y X Y gf id X, fg id Y f : X Y, g : Y X X Y X Y (2) 1.3
More informationx(t) + t f(t, x) = x(t) + x (t) t x t Tayler x(t + t) = x(t) + x (t) t + 1 2! x (t) t ! x (t) t 3 + (15) Eular x t Teyler 1 Eular 2 Runge-Kutta
6 Runge-KuttaEular Runge-Kutta Runge-Kutta A( ) f(t, x) dx dt = lim x(t + t) x(t) t 0 t = f(t, x) (14) t x x(t) t + dt x x(t + dt) Euler 7 t 1 f(t, x(t)) x(t) + f(t + dt, x(t + dt))dt t + dt x(t + dt)
More information1
1 1 7 1.1.................................. 11 2 13 2.1............................ 13 2.2............................ 17 2.3.................................. 19 3 21 3.1.............................
More information(3) (2),,. ( 20) ( s200103) 0.7 x C,, x 2 + y 2 + ax = 0 a.. D,. D, y C, C (x, y) (y 0) C m. (2) D y = y(x) (x ± y 0), (x, y) D, m, m = 1., D. (x 2 y
[ ] 7 0.1 2 2 + y = t sin t IC ( 9) ( s090101) 0.2 y = d2 y 2, y = x 3 y + y 2 = 0 (2) y + 2y 3y = e 2x 0.3 1 ( y ) = f x C u = y x ( 15) ( s150102) [ ] y/x du x = Cexp f(u) u (2) x y = xey/x ( 16) ( s160101)
More information5.. 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 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 information(1) 1 y = 2 = = b (2) 2 y = 2 = 2 = 2 + h B h h h< h 2 h
6 6.1 6.1.1 O y A y y = f() y = f() b f(b) B y f(b) f() = b f(b) f() f() = = b A f() b AB O b 6.1 2 y = 2 = 1 = 1 + h (1 + h) 2 1 2 (1 + h) 1 2h + h2 = h h(2 + h) = h = 2 + h y (1 + h) 2 1 2 O y = 2 1
More informationA
A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................
More informationhttp://know-star.com/ 3 1 7 1.1................................. 7 1.2................................ 8 1.3 x n.................................. 8 1.4 e x.................................. 10 1.5 sin
More information( ) ( ) ( ) i (i = 1, 2,, n) x( ) log(a i x + 1) a i > 0 t i (> 0) T i x i z n z = log(a i x i + 1) i=1 i t i ( ) x i t i (i = 1, 2, n) T n x i T i=1 z = n log(a i x i + 1) i=1 x i t i (i = 1, 2,, n) n
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 information< 1 > (1) f 0 (a) =6a ; g 0 (a) =6a 2 (2) y = f(x) x = 1 f( 1) = 3 ( 1) 2 =3 ; f 0 ( 1) = 6 ( 1) = 6 ; ( 1; 3) 6 x =1 f(1) = 3 ; f 0 (1) = 6 ; (1; 3)
< 1 > (1) f 0 (a) =6a ; g 0 (a) =6a 2 (2) y = f(x) x = 1 f( 1) = 3 ( 1) 2 =3 ; f 0 ( 1) = 6 ( 1) = 6 ; ( 1; 3) 6 x =1 f(1) = 3 ; f 0 (1) = 6 ; (1; 3) 6 y = g(x) x = 1 g( 1) = 2 ( 1) 3 = 2 ; g 0 ( 1) =
More informationuntitled
II(c) 1 October. 21, 2009 1 CS53 yamamoto@cs.kobe-u.ac.jp 3 1 7 1.1 : : : : : : : : : : : : : : : : : : : : : : 7 1.2 : : : : : : : : : : : : : : : : 8 1.2.1 : : : : : : : : : : : : : : : : : : : 8 1.2.2
More information2014計算機実験1_1
H26 1 1 1 seto@ics.nara-wu.ac.jp 数学モデリングのプロセス 問題点の抽出 定義 仮定 数式化 万有引力の法則 m すべての物体は引き合う r mm F =G 2 r M モデルの検証 モデルによる 説明 将来予測 解釈 F: 万有引力 (kg m s-2) G: 万有引力定数 (m s kg ) 解析 数値計算 M: 地球の質量 (kg) により解を得る m: 落下する物質の質量
More informationDE-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 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 informationD 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 informationB [ 0.1 ] x > 0 x 6= 1 f(x) µ 1 1 xn 1 + sin sin x 1 x 1 f(x) := lim. n x n (1) lim inf f(x) (2) lim sup f(x) x 1 0 x 1 0 (
. 28 4 14 [.1 ] x > x 6= 1 f(x) µ 1 1 xn 1 + sin + 2 + sin x 1 x 1 f(x) := lim. 1 + x n (1) lim inf f(x) (2) lim sup f(x) x 1 x 1 (3) lim inf x 1+ f(x) (4) lim sup f(x) x 1+ [.2 ] [, 1] Ω æ x (1) (2) nx(1
More information, x R, f (x),, df dx : R R,, f : R R, f(x) ( ).,, f (a) d f dx (a), f (a) d3 f dx 3 (a),, f (n) (a) dn f dx n (a), f d f dx, f d3 f dx 3,, f (n) dn f
,,,,.,,,. R f : R R R a R, f(a + ) f(a) lim 0 (), df dx (a) f (a), f(x) x a, f (a), f(x) x a ( ). y f(a + ) y f(x) f(a+) f(a) f(a + ) f(a) f(a) x a 0 a a + x 0 a a + x y y f(x) 0 : 0, f(a+) f(a)., f(x)
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 informationohp_06nov_tohoku.dvi
2006 11 28 1. (1) ẋ = ax = x(t) =Ce at C C>0 a0 x(t) 0(t )!! 1 0.8 0.6 0.4 0.2 2 4 6 8 10-0.2 (1) a =2 C =1 1. (1) τ>0 (2) ẋ(t) = ax(t τ) 4 2 2 4 6 8 10-2 -4 (2) a =2 τ =1!! 1. (2) A. (2)
More informationC10-075 26 2 12 1 1 4 1.1............................. 4 1.2............................ 5 1.3................................... 5 2 6 2.1............................ 6 2.2......................... 6
More information