1 4 1 ( ) ( ) ( ) ( ) () 1 4 2
|
|
- あきお おおはし
- 5 years ago
- Views:
Transcription
1 7 1995, (1) (2) (3) t (1) (2) Euler (3) θ, (1) (Landau) O( ), o( ) (2) (3) (4) 1 (46) (5) (6) (Gauss) (7) (35) ( [1]) 7 ( ) 1
2 1 4 1 ( ) ( ) ( ) ( ) () 1 4 2
3 2 f x f (x) (1) f (x) = lim h 0 f(x + h) f(x) h h f (x) (2) f(x + h) f(x) h f (x) f Taylor f x (3) f(x + h) = f(x) + f (x)h + f (x) 2 h h (4) f(x h) = f(x) f (x)h + f (x) 2 (3) (5) f (x) = f(x + h) f(x) h h 2 + f (x) 3! h 2 f (x) 3! + O(h) (h 0) h 3 + f (4) (x) h 4 + 4! h 3 + f (4) (x) h ! O(h) O(h) Landau ( (1) (p. 16) ) (4) (6) f (x) = f(x) f(x h) h + O(h) (h 0) 1 f (x) O(h) (3) (4) (7) f (x) = f(x + h) f(x h) 2h + O(h 2 ) (h 0) 1 1 f (x) 1 O(h 2 ) (3) (4) : (8) f (x) = f(x + h) 2f(x) + f(x h) h 2 + O(h 2 ) (h 0). 1 2 Landau C 3
4 t O 1 x 1: [0, 1] [0, ) f x C 2 f(x + h) 2f(x) + f(x h) lim h 0 h 2 = f (x) 3 6 () (9) u t (x, t) = 2 u (x, t) x2 (t > 0, 0 < x < 1) (10) u(0, t) = u(1, t) = 0 (t > 0) (11) u(x, 0) = f(x) (0 x 1) u = u(x, t) [0, 1] [0, + ) 1 (x i, t n ) u (12) u n i := u(x i, t n ) Ui n x [0, 1] N x i (i = 0,, N) (stepsize) h = 1/N (13) x i = ih (i = 0, 1,, N). t τ (> 0) 1 (14) t n = nτ (n = 0, 1, 2, ) (x i, t n ) u u t 2 2 x2 (15) (16) u t (x i, t n ) = un+1 i u n i + O(τ) (τ +0), τ 2 u x (x i, t 2 n ) = un i+1 2u n i + u n i 1 + O(h 2 ) (h +0) h 2 4
5 u (9) h τ (17) U n i (18) u n+1 i U n+1 i u n i τ Ui n τ un i+1 2u n i + u n i 1 h 2 = U n i+1 2U n i + U n i 1 h 2 2 n 0 i U n i 1, U n i+1 i = 0, N (18) 1 i N 1 λ = τ/h 2 (19) U n+1 i U n i = λ(u n i+1 2U n i + U n i 1). (20) U n+1 i = (1 2λ)Ui n + λ ( ) Ui 1 n + Ui+1 n (10) (21) U n 0 = U n N = 0 (n = 1, 2, ), (11) (22) U 0 i = f(x i ) (0 i N) (1 i N 1; n = 0, 1, ) (9), (10), (11) {Ui n ; 0 i N, n 0} (23) (24) (25) U n+1 i = (1 2λ)Ui n + λ ( Ui 1 n + Ui+1 n U0 n = UN n = 0 (n = 1, 2, ) U 0 i = f(x i ) (0 i N) ) (1 i N 1; n = 0, 1, ) {U n i ; 0 i N, n 0} 1. n 0 {U n i ; 0 i N} (25) 2. n = 0, 1, {U n+1 i ; 0 i N} (a) 1 i N 1 U n+1 i (23) (b) i = 0, N (24) U n+1 0 = U n+1 N = 0. Euler U n+1 i ( ) 2 F n = F n 1 + F n 2 Fibonacci 5
6 N u(1/2, 0.1) (%) : (x, t) = (1/2, 0.1) (λ = 1/2). 4 (1) ( ) (9)-(11) u f f(x) = f 0 (x) sin πx u (26) u(x, t) = e π2t sin πx 6 (27) u(x, t) = b k e k2 π 2t sin kπx, b k = 2 1 k=1 0 f(x) sin kπx dx (k = 1, 2, ) b k b 1 = 1, b k = 0 (k = 2, 3, ) (26) u (9), (10), (11) 1 (x, t) = (1/2, 0.1) u(1/2, 0.1) = exp( 0.1π 2 ) sin(π/2) = N 2 ( (2) ) N N 2 (2) t t + u(, t) 6
7 Error N 2: (x, t) = (1/2, 0.1) (λ = 1/2). (Dirichlet) (9)-(11) (28) u(x, t) = b k e k2 π 2t sin kπx, b k = 2 k=1 t u 0 : 1 (29) lim t u(x, t) = 0 (0 x 1). ( (29) (3) ) f x (30) f 1 (x) = 1 x 0 f(x) sin kπx dx (k = 1, 2, ) ( 0 x 1 ) 2 ( ) 1 2 x 1 ( 3) (9)-(11) N = 50, λ = 1/2 4 7 n {U n i ; 0 i N} t = nτ = 0.01, 0.02, 0.03, 0.04 x u(x, t) 8 t = 0 t = (Neumann) (10) (31) u x (0, t) = u x (1, t) = 0 (t > 0) 6 (32) u(x, t) = 1 2 a 0 + a k e k2 π 2t cos kπx, a k = 2 1 k=1 0 7 f(x) cos kπx dx (k = 0, 1, )
8 theta=0, N=50, lambda=0.5 theta=0, N=50, lambda=0.5 t = 0.01 theta=0, N=50, lambda=0.5 t = : f = f 1 4: t = : t = 0.02 theta=0, N=50, lambda=0.5 t = 0.03 theta=0, N=50, lambda=0.5 t = 0.04 theta=0, N=50, lambda=0.5 Tmax = 0.3, interval = : t = : t = :, t = 0 0.3, t =
9 heat equation, Neumann B.C. theta=0, N=50, lambda=0.5 heat equation, Neumann B.C. theta=0, N=50, lambda=0.5 t = 0.01 heat equation, Neumann B.C. theta=0, N=50, lambda=0.5 t = : t = 0 10: t = : t = 0.02 heat equation, Neumann B.C. theta=0, N=50, lambda=0.5 t = 0.03 heat equation, Neumann B.C. theta=0, N=50, lambda=0.5 t = 0.04 heat equation, Neumann B.C. theta=0, N=50, lambda=0.5 Tmax = 0.2, interval = : t = : t = : t = t u : (33) lim t u(x, t) = a 0 2 (0 x 1) f 1 N = 50, λ = 1/2 u(, t) f = f 1 (34) a 0 2 = 1 0 f 1 (x) dx = 3 x = /4 = 0.25 ( ) (3) t 9
10 t \ x : (, θ = 0, N = 40, λ = 0.5) 6 f a (a k cos kπx + b k sin kπx), k=1 f ( ) f a k, b k a k, b k k : (35) f C l (l ) = a k, b k = o(k l ) (k + ). f f (= ) (36) k l ( a k + b k ) < = f C l. k=1 f 15 f 2 f k + a k, b k 3 f = f 2 (9)-(11) (27) k exp ( k 2 π 2 t) t > 0 k f = f 2 u () f = f t = () 3 a k, b k = o(1) (k ) k ( a k + b k ) = k=1 10
11 k sin kπx k exp ( k 2 π 2 t) t (27) t k = 2, 3, k = 1 b 1 exp( π 2 t) sin πx t (37) u(x, t) b 1 e π2t sin πx. 20 theta=0, N=100, lambda=0.25 theta=0, N=100, lambda=0.25 t = theta=0, N=100, lambda=0.25 t = : t = 0 theta=0, N=100, lambda=0.25 t = : t = theta=0, N=100, lambda=0.25 t = : t = theta=0, N=100, lambda=0.25 Tmax = 0.2, interval = : t = : t = : t = 0 0.2, t = (1) λ = τ/h 2 1/2 λ = 1/4 λ = 0.51 f = f 1 (9)-(11)
12 4 7 (N ) 3 λ 0 < λ 1/2 0 < λ 1/2 (38) Ui n max f(x) (0 i N; n = 0, 1, ) x [0,1] ( (5) 3 ) 0 < λ 1/2 T (39) lim max u(x i, t n ) Ui n = 0 N 0 i N 0 n T /τ ( (5) 2 ) 0 i N, 0 n T/τ (x i, t n ) [0, 1] [0, T ] theta=0, N=50, lambda=0.51 t = theta=0, N=50, lambda=0.51 t = theta=0, N=50, lambda=0.51 t = : t = : t = theta=0, N=50, lambda=0.51 t = : t = : t = (2) Euler 3 0 < λ 1/2 (= ) 12
13 3 u t (40) U n i U n 1 i τ = U n i+1 2U n i + U n i 1 h 2 (1 i N 1, n 1) n (41) U n+1 i Ui n τ = U n+1 i+1 2U n+1 i + U n+1 i 1 h 2 (1 i N 1, n 0) n n + 1 (42) (1 + 2λ)U n+1 i λ ( U n+1 i 1 + U ) n+1 i+1 = U n i (1 i N 1, n 0) Euler (42) 3 (23) n + 1 (23) U n+1 i (i = 1, 2,, N 1) (42) i (1 i N 1) 1 (23) ( ) (42) ( ) λ = (42) λ (> 0) λ = 5 (3) θ, - (41) (18) (41) 1 θ : θ (43) U n+1 i Ui n τ = (1 θ) U i+1 n 2Ui n + Ui 1 n + θ U n+1 i+1 h 2 2U n+1 i + U n+1 i 1 h 2 θ 0 θ 1 θ = 0 (18) θ = 1 (41) n + 1 (44) (1 + 2θλ)U n+1 i θλ ( U n+1 i 1 + U ) n+1 i+1 = {1 2(1 θ)λ} U n i + (1 θ)λ ( Ui 1 n + Ui+1) n (1 i N 1, n 0) ({Ui n ; 0 i N} ) N +1 {U n+1 i ; 0 i N} N 1 (45) U n+1 0 = U n+1 N = 0 13
14 theta=1, N=50, lambda=0.51 t = theta=1, N=50, lambda=0.51 t = theta=1, N=50, lambda=0.51 t = : t = : t = theta=1, N=50, lambda=0.51 t = : t = : t = theta=1, N=50, lambda=5 t = 0.01 theta=1, N=50, lambda=5 t = 0.02 theta=1, N=50, lambda=5 t = : t = : t = 0.02 theta=1, N=50, lambda=5 t = : t = : t =
15 2 U0 n+1, U n+1 N N 1 {U n+1 i ; 1 i N 1} N 1 1 (46) 1 + 2θλ θλ 0 U 1 n+1 θλ 1 + 2θλ θλ U n+1 θλ 1 + 2θλ θλ U i. n+1 θλ 1 + 2θλ θλ U n+1 N 2 0 θλ 1 + 2θλ U n+1 N 1 {1 2(1 θ)λ}u1 n + (1 θ)λ(u0 n + U2 n ) {1 2(1 θ)λ}u2 n + (1 θ)λ(u1 n + U3 n ) = {1 2(1 θ)λ}ui n + (1 θ)λ(ui 1 n + Ui+1) n {1 2(1 θ)λ}un 2 n + (1 θ)λ(un 3 n + U N 1 n ) {1 2(1 θ)λ}un 1 n + (1 θ)λ(un 2 n + U N n ) θ θ = 1/2 Crank-Nicolson (- ) ( Crank-Nicolson [2]). 6 6 (47) u t (x, t) = 2 u (x, t) + ku(x, t)(1 u(x, t)) x2 (t > 0, 0 < x < 1) u x (t, 0) = 1, u(t, 1) = 0 (t > 0) u(0, x) = f(x) (0 x 1) (48) U n+1 i Ui n τ = (1 θ) U n i+1 2U n i + U n i 1 + θ U n+1 i+1 h 2 2U n+1 i + U n+1 i 1 h 2 + ku n i (1 U n i ) 15
16 7 (1) (Landau) O( ), o( ) O(h) (h 0) o(n k ) (n ) O( ) C α x (49) f(x) C g(x) 4 (50) f(x) = O(g(x)) (x α) 5 α f g o( ) (51) lim x α f(x) g(x) = 0 (52) f(x) = o(g(x)) (x α) f a A (53) A = lim h 0 f(a + h) f(a) h (54) f(a + h) f(a) Ah = o(h) (h 0) f C k (55) f(a + h) = f(a) + f (a)h + f (a) h f (k 1) (a) 2! (k 1)! hk 1 + f (k) (a + θh) k! θ (0 < θ < 1) (56) f(a + h) = f(a) + f (a)h + f (a) h f (k 1) (a) 2! (k 1)! hk 1 + O(h k ) (h 0) 4 D f C > 0, ε > 0, x B(α; ε) D f(x) C g(x). 5 : f(x) O(g(x)) O(g(x)) 1 h k 16
17 α ± x 2 + x + 1 x = O(x) (x ). (2) 4 (1) x, y y = ax + b (a, b ) 1 a y = bx a (a, b, b > 0) y = bx a (57) log y = log bx a = a log x + log b. Y = log y, B = log b, X = log x Y = ax + B. X, Y 1 x, y X, Y a N E E = C/N 2 (C ) log N, log E 2 (3) u (58) lim t u(x, t) = 0, u (59) lim t u(x, t) = a f(x) dx (58) (0 ) (60) u(x, t) = b n e n2 π 2t sin nπx lim t e n2 π 2t = 0 (61) lim t n=1 b n e n2 π 2t sin nπx = n=1 17 lim b ne n2 π 2t sin nπx t n=1
18 n=1 lim t 0 1 t 0 > 0 (62) u(x, t) Ce π2 t (t t 0, 0 x 1) C f(t) C, τ (63) f(t) Ce t/τ ( t) Ce t/τ f(t) t u(x, t) 0 f M := max x [0,1] f(x) b n (64) b n = 2 (65) b n f(x) sin nπx dx 2 0 f(x) sin nπx dx 1 0 M 1 dx = 2M 1 0 dx = 2M. b n n 2M (66) u(x, t) = (67) u(x, t) = b n e n2 π 2t sin nπx n=1 b n e n2 π 2t sin nπx n=1 ( t t 0 2M n=1 e (n2 1)π 2 t (68) u(x, t) (69) ( ) 2M n=1 n=1 e π2t. ( ) 2M e n2 π 2t 1 = 2M ) e (n2 1)π 2 t 0 e π2t. n=1 e (n2 1)π 2 t 0 18 n=1 e n2 π 2 t
19 (70) e (n2 1)π 2 t 0 e (n 1)π2 t 0 = r n 1, r = e π2 t 0 < 1 r n 1, 0 < r < 1 n=1 (4) 1 (46) 5 (3) 1 (46) {U n+1 i } (, ) N 1 N 1 1 n A = (a ij ) (71) a ii > a ij (i = 1, 2,, n) 1 j n,j i A θλ 0 θλ, θλ (72) 1 + 2θλ θλ θλ = 1 + 2θλ 2(θλ) = 1 > A = (a ij ) Ax = 0 x ( 0) 6 (73) x = x k Ax = 0 k x 1. x n (74) n a kj x j = 0. j=1 6 n A A f : R n x Ax R n (i) A (A 1 ), (ii) f, (iii) f, (iv) f, (v) ker f = {0}, (vi) rank f = n, (vii) det A 0. n dim ker f = rank f ( ) 19
20 k (75) a kk x k = j k a kj x j. (76) a kk x k = a kj x j ( ) a kj x j a kj x k. j k j k j k x k (> 0) (77) a kk a jk. j k A (5) (9)-(11) 2 ( ) T (9)-(11) 0 x 1, 0 t T u θ λ 2 (i) 0 θ < 1, 0 < λ (ii) θ = 1, λ > (1 θ). θ {U n i } 0 i N,n 0 : (78) max Ui n u n i = O(τ + h 2 ) (N ). 0 i N 0 n J u n i = u(x i, t n ) J T/τ N (79) max Ui n u n i 0. 0 i N 0 n J (79)
21 3 ( ) {U n i ; 0 i N, 0 n J} (80) (1 + 2θλ)U n+1 i θλ(u n+1 i 1 + U n+1 i+1 ) = [1 2(1 θ)λ]ui n + (1 θ)λ(ui 1 n + Ui+1) n (1 i N 1, 0 n J 1) θ, λ 2 (i) 0 θ < 1, 0 < λ (ii) θ = 1, λ > (1 θ). : { = max (81) max Ui n 0 i N 0 n J max Uk 0, max U0, l max UN l 0 k N 0 l J 0 l J }. 2 (82) 2 u n+1 i u n i τ = u t (x i, t n ) + O(τ), (83) (84) u n i 1 2u n i + u n i+1 h 2 = u xx (x i, t n ) + O(h 2 ). u n+1 i u n i un i 1 2u n i + u n i+1 τ h 2 = u t (x i, t j ) u xx (x i, t j ) + O(τ + h 2 ) = O(τ + h 2 ). u u t = u xx (85) u n+1 i u n i τ un+1 i 1 2un+1 i + u n+1 i+1 = O(τ + h 2 ). h 2 {u n i } : (86) (1 + 2θλ)u n+1 i θλ(u n+1 i 1 + un+1 i+1 ) = [1 2(1 θ)λ]u n i + (1 θ)λ(u n i 1 + u n i+1) + O(τ 2 + τh 2 ). ( 1 i N 1, 0 n J 1 (i, n) ) Ui n ( ) (87) (1 + 2θλ)U n+1 i θλ(u n+1 i 1 + U n+1 i+1 ) = [1 2(1 θ)λ]u n i + (1 θ)λ(u n i 1 + U n i+1) 21
22 e n i := U n i u n i e n i (88) (1 + 2θλ)e n+1 i θλ(e n+1 i 1 + en+1 i+1 ) = [1 2(1 θ)λ]e n i + (1 θ)λ(e n i 1 + e n i+1) + O(τ 2 + τh 2 ) θ, λ 1 + 2θλ, θλ, 1 2(1 θ)λ, (1 θ)λ 0 (89) (1 + 2θλ) e n+1 i θλ( e n+1 i 1 + en+1 i+1 ) [1 2(1 θ)λ] e n i + (1 θ)λ( e n i 1 + e n i+1 ) + C(τ 2 + τh 2 ). C E n := max e n i 0 i N (1 + 2θλ) e n+1 i 2θλE n+1 (1 + 2θλ) e n+1 i θλ( e n+1 i 1 + en+1 i+1 ) (90) [1 2(1 θ)λ] e n i + (1 θ)λ( e n i 1 + e n i+1 ) + C(τ 2 + τh 2 ) [1 2(1 θ)λ]e n + 2(1 θ)λe n + C(τ 2 + τh 2 ) = E n + C(τ 2 + τh 2 ). e n+1 0 = e n+1 N = 0 i 1 i N 1 (91) (1 + 2θλ)E n+1 2θλE n+1 E n + C(τ 2 + τh 2 ). (92) E n+1 E n + C(τ 2 + τh 2 ). (93) E n E n 1 + C(τ 2 + τh 2 ) E n 2 + 2C(τ 2 + τh 2 ) E 0 + nc(τ 2 + τh 2 ) = nc(τ 2 + τh 2 ). i e 0 i = U 0 i u 0 i = f(x i ) f(x i ) = 0 E 0 = 0 nτ Jτ T (94) E n CT (τ + h 2 ). (95) max Ui n u n i CT (τ + h 2 ). 0 i N 0 n J 22
23 3 (81) U n i 0 n J 1 n (96) max U n+1 i max 0 i N { } max Uk n, U0 n+1, U n+1 N. 0 k N (97) max U n+1 i max Uk n 1 i N 1 0 k N i = i 0 (80) i = i 0 : max U n+1 i = U n+1 i 1 i N 1 0 (98) (1 + 2θλ)U n+1 i 0 θλ(u n+1 i U n+1 i 0 +1 ) = (1 2(1 θ)λ)u n i 0 + (1 θ)λ(u n i U n i 0 +1). θ, λ 1 + 2θλ, θλ, 1 2(1 θ)λ, (1 θ)λ 0 (99) (98) {1 2(1 θλ)} max Uk n + (1 θ)λ( max Uk n + max Uk n ) 0 k N 0 k N 0 k N = {[1 2(1 θλ)] + 2(1 θ)λ} max Uk n 0 k N = max Uk n. 0 k N (100) (98) (1 + 2θλ)U n+1 i 0 θλ{u n+1 i 0 + U n+1 i 0 } = U n+1 i 0. (99), (100) U n+1 i 0 max Uk n. (97) : 0 k N (101) max U n+1 i max Uk n. 1 i N 1 0 k N (96) (96) n = 0 { } (102) max Ui 1 max max Uk 0, U0 1, UN 1. 0 i N 0 k N (96) n = 1 { max (103) max Ui 2 0 i N (102), (103) { max (104) max Ui 2 0 i N max Uk 0, U0 1, U0 2, UN, 1 UN 2 0 k N max Uk 1, U0 2, UN 2 0 k N n = 1, 2,, J { } (105) max Ui n max max Uk 0, max U0, l max UN l 0 i N 0 k N 1 l n 1 l n (81) 23 } }. { } = max max Uk 0, max U0, l max UN l 0 k N 1 l 2 1 l 2
24 (6) (Gauss) 1 (Cramer) ( Jordan ) 8 2x 1 + 3x 2 x 3 = 5 (106) 4x 1 + 4x 2 3x 3 = 3 2x 1 + 3x 2 x 3 = , x 1 x 2 x 3 = x 1 + 3x 2 x 3 = 5, 2x 2 x 3 = 7, 5x 3 = 15 x 3 = 15 5 = 3, x 2 = 7 + x 3 2 = 2, x 1 = 5 3x 2 + x 3 2 = = 1 8 () ( ) 24
25 0 (forward elimination) (backward substitution) 3 n + 1 () n ( n ) 5 5 n 5, 6 θλ = 1 (7) (35) ( ) 2π f a k (f) = 1 π π π f(x) cos kx dx, b k (f) = 1 π f C l { a k (f) = ± 1 a k (f (l) ) (l ) k l b k (f (l) ) (l ), π π b k (f) = ± 1 k l { f (l) Riemann-Lebesgue lim a k(f (l) ) = lim b k (f (l) ) = 0 k k f(x) sin kx dx b k (f (l) ) a k (f (l) ) (l ) (l ) lim k kl a k (f) = lim k l b k (f) = 0. k 8 1. Neumann I [3] 25
26 2., von Neumann [4], I [3] 3. LU Gauss [4], 1 1 [5] 4. [6] 5., [7] 6., [8] 7. (?) Richtmyer-Morton [9], [10], Smith [11], [12] [1], (2000), 7 [2] Crank, J. and Nicolson, P.: A Practical Method for Numerical Evaluation of Solutions of Partial Differential Equations of the Heat-Conduction Type, Proc. Camb. Phil. Soc., Vol. 43, pp (1947), g j807/ reprint. [3] I, meiji.ac.jp/~mk/labo/text/heat-fdm-1.pdf (1998 ). [4], heat-fdm-0-add.pdf (1997 ). [5] 1 I, linear-eq-1.pdf (2002 ). [6], text/wave.pdf (1999? ). [7],, (1991). [8], (2006). [9] Richtmyer, R. D. and Morton, K. W.: Difference Methods for Initial-Value Problems, Krieger Pub Co, 2nd edition (1994). 26
27 [10], ( ), ( ), (1968), Finite-Difference Methods for Partial Differential Equations (1960) ( ). [11] Smith, G. D.: Numerical solution of partial differential equations third edition, Clarendon Press Oxford (1986), G. D.,,, (1996). [12],, (1969),
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 information1
1 1 7 1.1.................................. 11 2 13 2.1............................ 13 2.2............................ 17 2.3.................................. 19 3 21 3.1.............................
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微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.
微分積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. ttp://www.morikita.co.jp/books/mid/00571 このサンプルページの内容は, 初版 1 刷発行時のものです. i ii 014 10 iii [note] 1 3 iv 4 5 3 6 4 x 0 sin x x 1 5 6 z = f(x, y) 1 y = f(x)
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 informationz f(z) f(z) x, y, u, v, r, θ r > 0 z = x + iy, f = u + iv C γ D f(z) f(z) D f(z) f(z) z, Rm z, z 1.1 z = x + iy = re iθ = r (cos θ + i sin θ) z = x iy
f f x, y, u, v, r, θ r > = x + iy, f = u + iv C γ D f f D f f, Rm,. = x + iy = re iθ = r cos θ + i sin θ = x iy = re iθ = r cos θ i sin θ x = + = Re, y = = Im i r = = = x + y θ = arg = arctan y x e i =
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 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 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数値計算:有限要素法
( ) 1 / 61 1 2 3 4 ( ) 2 / 61 ( ) 3 / 61 P(0) P(x) u(x) P(L) f P(0) P(x) P(L) ( ) 4 / 61 L P(x) E(x) A(x) x P(x) P(x) u(x) P(x) u(x) (0 x L) ( ) 5 / 61 u(x) 0 L x ( ) 6 / 61 P(0) P(L) f d dx ( EA du dx
More information2 7 V 7 {fx fx 3 } 8 P 3 {fx fx 3 } 9 V 9 {fx fx f x 2fx } V {fx fx f x 2fx + } V {{a n } {a n } a n+2 a n+ + a n n } 2 V 2 {{a n } {a n } a n+2 a n+
R 3 R n C n V??,?? k, l K x, y, z K n, i x + y + z x + y + z iv x V, x + x o x V v kx + y kx + ky vi k + lx kx + lx vii klx klx viii x x ii x + y y + x, V iii o K n, x K n, x + o x iv x K n, x + x o x
More information第5章 偏微分方程式の境界値問題
October 5, 2018 1 / 113 4 ( ) 2 / 113 Poisson 5.1 Poisson ( A.7.1) Poisson Poisson 1 (A.6 ) Γ p p N u D Γ D b 5.1.1: = Γ D Γ N 3 / 113 Poisson 5.1.1 d {2, 3} Lipschitz (A.5 ) Γ D Γ N = \ Γ D Γ p Γ N Γ
More informationII 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K
II. () 7 F 7 = { 0,, 2, 3, 4, 5, 6 }., F 7 a, b F 7, a b, F 7,. (a) a, b,,. (b) 7., 4 5 = 20 = 2 7 + 6, 4 5 = 6 F 7., F 7,., 0 a F 7, ab = F 7 b F 7. (2) 7, 6 F 6 = { 0,, 2, 3, 4, 5 },,., F 6., 0 0 a F
More information1W II K =25 A (1) office(a439) (2) A4 etc. 12:00-13:30 Cafe David 1 2 TA appointment Cafe D
1W II K200 : October 6, 2004 Version : 1.2, kawahira@math.nagoa-u.ac.jp, http://www.math.nagoa-u.ac.jp/~kawahira/courses.htm TA M1, m0418c@math.nagoa-u.ac.jp TA Talor Jacobian 4 45 25 30 20 K2-1W04-00
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 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 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 information4 2016 3 8 2.,. 2. Arakawa Jacobin., 2 Adams-Bashforth. Re = 80, 90, 100.. h l, h/l, Kármán, h/l 0.28,, h/l.., (2010), 46.2., t = 100 t = 2000 46.2 < Re 46.5. 1 1 4 2 6 2.1............................
More informationA11 (1993,1994) 29 A12 (1994) 29 A13 Trefethen and Bau Numerical Linear Algebra (1997) 29 A14 (1999) 30 A15 (2003) 30 A16 (2004) 30 A17 (2007) 30 A18
2013 8 29y, 2016 10 29 1 2 2 Jordan 3 21 3 3 Jordan (1) 3 31 Jordan 4 32 Jordan 4 33 Jordan 6 34 Jordan 8 35 9 4 Jordan (2) 10 41 x 11 42 x 12 43 16 44 19 441 19 442 20 443 25 45 25 5 Jordan 26 A 26 A1
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 informationad bc A A A = ad bc ( d ) b c a n A n A n A A det A A ( ) a b A = c d det A = ad bc σ {,,,, n} {,,, } {,,, } {,,, } ( ) σ = σ() = σ() = n sign σ sign(
I n n A AX = I, YA = I () n XY A () X = IX = (YA)X = Y(AX) = YI = Y X Y () XY A A AB AB BA (AB)(B A ) = A(BB )A = AA = I (BA)(A B ) = B(AA )B = BB = I (AB) = B A (BA) = A B A B A = B = 5 5 A B AB BA A
More informationn ξ n,i, i = 1,, n S n ξ n,i n 0 R 1,.. σ 1 σ i .10.14.15 0 1 0 1 1 3.14 3.18 3.19 3.14 3.14,. ii 1 1 1.1..................................... 1 1............................... 3 1.3.........................
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() Remrk I = [0, ] [x i, x i ]. (x : ) f(x) = 0 (x : ) ξ i, (f) = f(ξ i )(x i x i ) = (x i x i ) = ξ i, (f) = f(ξ i )(x i x i ) = 0 (f) 0.
() 6 f(x) [, b] 6. Riemnn [, b] f(x) S f(x) [, b] (Riemnn) = x 0 < x < x < < x n = b. I = [, b] = {x,, x n } mx(x i x i ) =. i [x i, x i ] ξ i n (f) = f(ξ i )(x i x i ) i=. (ξ i ) (f) 0( ), ξ i, S, ε >
More informationmeiji_resume_1.PDF
β β β (q 1,q,..., q n ; p 1, p,..., p n ) H(q 1,q,..., q n ; p 1, p,..., p n ) Hψ = εψ ε k = k +1/ ε k = k(k 1) (x, y, z; p x, p y, p z ) (r; p r ), (θ; p θ ), (ϕ; p ϕ ) ε k = 1/ k p i dq i E total = E
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, 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 informationmugensho.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 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 information9. 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 informationi
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,. 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 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 informationall.dvi
5,, Euclid.,..,... Euclid,.,.,, e i (i =,, ). 6 x a x e e e x.:,,. a,,. a a = a e + a e + a e = {e, e, e } a (.) = a i e i = a i e i (.) i= {a,a,a } T ( T ),.,,,,. (.),.,...,,. a 0 0 a = a 0 + a + a 0
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 information: 1g99p038-8
16 17 : 1g99p038-8 1 3 1.1....................................... 4 1................................... 5 1.3.................................. 5 6.1..................................... 7....................................
More information2S III IV K A4 12:00-13:30 Cafe David 1 2 TA 1 appointment Cafe David K2-2S04-00 : C
2S III IV K200 : April 16, 2004 Version : 1.1 TA M2 TA 1 10 2 n 1 ɛ-δ 5 15 20 20 45 K2-2S04-00 : C 2S III IV K200 60 60 74 75 89 90 1 email 3 4 30 A4 12:00-13:30 Cafe David 1 2 TA 1 email appointment Cafe
More 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 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 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 informationA
A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................
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 informationA S- hara/lectures/lectures-j.html r A = A 5 : 5 = max{ A, } A A A A B A, B A A A %
A S- http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html r A S- 3.4.5. 9 phone: 9-8-444, e-mail: hara@math.kyushu-u.ac.jp, http://www.math.kyushu-u.ac.jp/ hara/lectures/lectures-j.html Office
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 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 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 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
[ ]. 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入試の軌跡
4 y O x 4 Typed by L A TEX ε ) ) ) 6 4 ) 4 75 ) http://kumamoto.s.xrea.com/plan/.. PDF) Ctrl +L) Ctrl +) Ctrl + Ctrl + ) ) Alt + ) Alt + ) ESC. http://kumamoto.s.xrea.com/nyusi/kumadai kiseki ri i.pdf
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 informationf(x) = e x2 25 d f(x) 0 x d2 dx f(x) 0 x dx2 f(x) (1 + ax 2 ) 2 lim x 0 x 4 a 3 2 a g(x) = 1 + ax 2 f(x) g(x) 1/2 f(x)dx n n A f(x) = Ax (x R
29 ( ) 90 1 2 2 2 1 3 4 1 5 1 4 3 3 4 2 1 4 5 6 3 7 8 9 f(x) = e x2 25 d f(x) 0 x d2 dx f(x) 0 x dx2 f(x) (1 + ax 2 ) 2 lim x 0 x 4 a 3 2 a g(x) = 1 + ax 2 f(x) g(x) 1/2 f(x)dx 11 0 24 n n A f(x) = Ax
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 information「産業上利用することができる発明」の審査の運用指針(案)
1 1.... 2 1.1... 2 2.... 4 2.1... 4 3.... 6 4.... 6 1 1 29 1 29 1 1 1. 2 1 1.1 (1) (2) (3) 1 (4) 2 4 1 2 2 3 4 31 12 5 7 2.2 (5) ( a ) ( b ) 1 3 2 ( c ) (6) 2. 2.1 2.1 (1) 4 ( i ) ( ii ) ( iii ) ( iv)
More information熊本県数学問題正解
00 y O x Typed by L A TEX ε ( ) (00 ) 5 4 4 ( ) http://www.ocn.ne.jp/ oboetene/plan/. ( ) (009 ) ( ).. http://www.ocn.ne.jp/ oboetene/plan/eng.html 8 i i..................................... ( )0... (
More informationall.dvi
72 9 Hooke,,,. Hooke. 9.1 Hooke 1 Hooke. 1, 1 Hooke. σ, ε, Young. σ ε (9.1), Young. τ γ G τ Gγ (9.2) X 1, X 2. Poisson, Poisson ν. ν ε 22 (9.) ε 11 F F X 2 X 1 9.1: Poisson 9.1. Hooke 7 Young Poisson G
More informationII K116 : January 14, ,. A = (a ij ) ij m n. ( ). B m n, C n l. A = max{ a ij }. ij A + B A + B, AC n A C (1) 1. m n (A k ) k=1,... m n A, A k k
: January 14, 28..,. A = (a ij ) ij m n. ( ). B m n, C n l. A = max{ a ij }. ij A + B A + B, AC n A C (1) 1. m n (A k ) k=1,... m n A, A k k, A. lim k A k = A. A k = (a (k) ij ) ij, A k = (a ij ) ij, i,
More informationi I II I II II IC IIC I II ii 5 8 5 3 7 8 iii I 3........................... 5......................... 7........................... 4........................ 8.3......................... 33.4...................
More information5. [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, 3, 6 = 3, 3,,,, 3,, 9, 3, 9, 3, 3, 4, 43, 4, 3, 9, 6, 6,, 0 p, p, p 3,..., p n N = p p p 3 p n + N p n N p p p, p 3,..., p n p, p,..., p n N, 3,,,,
6,,3,4,, 3 4 8 6 6................................. 6.................................. , 3, 6 = 3, 3,,,, 3,, 9, 3, 9, 3, 3, 4, 43, 4, 3, 9, 6, 6,, 0 p, p, p 3,..., p n N = p p p 3 p n + N p n N p p p,
More information(1) 3 A B E e AE = e AB OE = OA + e AB = (1 35 e ) e OE z 1 1 e E xy e = 0 e = 5 OE = ( 2 0 0) E ( 2 0 0) (2) 3 E P Q k EQ = k EP E y 0
(1) 3 A B E e AE = e AB OE = OA + e AB = (1 35 e 0 1 15 ) e OE z 1 1 e E xy 5 1 1 5 e = 0 e = 5 OE = ( 2 0 0) E ( 2 0 0) (2) 3 E P Q k EQ = k EP E y 0 Q y P y k 2 M N M( 1 0 0) N(1 0 0) 4 P Q M N C EP
More informationルベーグ積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.
ルベーグ積分 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/005431 このサンプルページの内容は, 初版 1 刷発行時のものです. Lebesgue 1 2 4 4 1 2 5 6 λ a
More information, 1 ( f n (x))dx d dx ( f n (x)) 1 f n (x)dx d dx f n(x) lim f n (x) = [, 1] x f n (x) = n x x 1 f n (x) = x f n (x) = x 1 x n n f n(x) = [, 1] f n (x
1 1.1 4n 2 x, x 1 2n f n (x) = 4n 2 ( 1 x), 1 x 1 n 2n n, 1 x n n 1 1 f n (x)dx = 1, n = 1, 2,.. 1 lim 1 lim 1 f n (x)dx = 1 lim f n(x) = ( lim f n (x))dx = f n (x)dx 1 ( lim f n (x))dx d dx ( lim f d
More informationy π π O π x 9 s94.5 y dy dx. y = x + 3 y = x logx + 9 s9.6 z z x, z y. z = xy + y 3 z = sinx y 9 s x dx π x cos xdx 9 s93.8 a, fx = e x ax,. a =
[ ] 9 IC. dx = 3x 4y dt dy dt = x y u xt = expλt u yt λ u u t = u u u + u = xt yt 6 3. u = x, y, z = x + y + z u u 9 s9 grad u ux, y, z = c c : grad u = u x i + u y j + u k i, j, k z x, y, z grad u v =
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 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 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 informationfunction2.pdf
2... 1 2009, http://c-faculty.chuo-u.ac.jp/ nishioka/ 2 11 38 : 5) i) [], : 84 85 86 87 88 89 1000 ) 13 22 33 56 92 147 140 120 100 80 60 40 20 1 2 3 4 5 7.1 7 7.1 1. *1 e = 2.7182 ) fx) e x, x R : 7.1)
More informationIA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................
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 information微分方程式の解を見る
( ) norikazu[at]ms.u-tokyo.ac.jp http://www.infsup.jp/saito/materials/829gairon.pdf 28 2 9 NS ( ) 28 2 9 / 5 Newton 2 3 Navier-Stokes 4 5 6 NS ( ) 28 2 9 2 / 5 2 x 2 + bx + c = x = b ± b 2 4c 2 3 x 3 +
More information(2018 2Q C) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = ( ) a c b d (a c, b d) P = (a, b) O P ( ) a p = b P = (a, b) p = ( ) a b R 2 {( ) } R 2 x = x, y
(2018 2Q C) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = a c b d (a c, b d) P = (a, b) O P a p = b P = (a, b) p = a b R 2 { } R 2 x = x, y R y 2 a p =, c q = b d p + a + c q = b + d q p P q a p = c R c b
More informationnumb.dvi
11 Poisson kanenko@mbkniftycom alexeikanenko@docomonejp http://wwwkanenkocom/ , u = f, ( u = u+f u t, u = f t ) 1 D R 2 L 2 (D) := {f(x,y) f(x,y) 2 dxdy < )} D D f,g L 2 (D) (f,g) := f(x,y)g(x,y)dxdy (L
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 information<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>
常微分方程式の局所漸近解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/007651 このサンプルページの内容は, 初版 1 刷発行当時のものです. i Leibniz ydy = y 2 /2 1675 11 11 [6] 100 Bernoulli Riccati 19 Fuchs
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 informationi ii iii iv v vi vii ( ー ー ) ( ) ( ) ( ) ( ) ー ( ) ( ) ー ー ( ) ( ) ( ) ( ) ( ) 13 202 24122783 3622316 (1) (2) (3) (4) 2483 (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 11 11 2483 13
More informationPart () () Γ 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 information1/1 lim f(x, y) (x,y) (a,b) ( ) ( ) lim limf(x, y) lim lim f(x, y) x a y b y b x a ( ) ( ) xy x lim lim lim lim x y x y x + y y x x + y x x lim x x 1
1/5 ( ) Taylor ( 7.1) (x, y) f(x, y) f(x, y) x + y, xy, e x y,... 1 R {(x, y) x, y R} f(x, y) x y,xy e y log x,... R {(x, y, z) (x, y),z f(x, y)} R 3 z 1 (x + y ) z ax + by + c x 1 z ax + by + c y x +
More informationkeisoku01.dvi
2.,, Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 5 Mon, 2006, 401, SAGA, JAPAN Dept.
More information(2 X Poisso P (λ ϕ X (t = E[e itx ] = k= itk λk e k! e λ = (e it λ k e λ = e eitλ e λ = e λ(eit 1. k! k= 6.7 X N(, 1 ϕ X (t = e 1 2 t2 : Cauchy ϕ X (t
6 6.1 6.1 (1 Z ( X = e Z, Y = Im Z ( Z = X + iy, i = 1 (2 Z E[ e Z ] < E[ Im Z ] < Z E[Z] = E[e Z] + ie[im Z] 6.2 Z E[Z] E[ Z ] : E[ Z ] < e Z Z, Im Z Z E[Z] α = E[Z], Z = Z Z 1 {Z } E[Z] = α = α [ α ]
More informationMicrosoft Word - 触ってみよう、Maximaに2.doc
i i e! ( x +1) 2 3 ( 2x + 3)! ( x + 1) 3 ( a + b) 5 2 2 2 2! 3! 5! 7 2 x! 3x! 1 = 0 ",! " >!!! # 2x + 4y = 30 "! x + y = 12 sin x lim x!0 x x n! # $ & 1 lim 1 + ('% " n 1 1 lim lim x!+0 x x"!0 x log 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 information1 1 u m (t) u m () exp [ (cπm + (πm κ)t (5). u m (), U(x, ) f(x) m,, (4) U(x, t) Re u k () u m () [ u k () exp(πkx), u k () exp(πkx). f(x) exp[ πmxdx
1 1 1 1 1. U(x, t) U(x, t) + c t x c, κ. (1). κ U(x, t) x. (1) 1, f(x).. U(x, t) U(x, t) + c κ U(x, t), t x x : U(, t) U(1, t) ( x 1), () : U(x, ) f(x). (3) U(x, t). [ U(x, t) Re u k (t) exp(πkx). (4)
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 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 informationchap1.dvi
1 1 007 1 e iθ = cos θ + isin θ 1) θ = π e iπ + 1 = 0 1 ) 3 11 f 0 r 1 1 ) k f k = 1 + r) k f 0 f k k = 01) f k+1 = 1 + r)f k ) f k+1 f k = rf k 3) 1 ) ) ) 1+r/)f 0 1 1 + r/) f 0 = 1 + r + r /4)f 0 1 f
More information1 Introduction 1 (1) (2) (3) () {f n (x)} n=1 [a, b] K > 0 n, x f n (x) K < ( ) x [a, b] lim f n (x) f(x) (1) f(x)? (2) () f(x)? b lim a f n (x)dx = b
1 Introduction 2 2.1 2.2 2.3 3 3.1 3.2 σ- 4 4.1 4.2 5 5.1 5.2 5.3 6 7 8. Fubini,,. 1 1 Introduction 1 (1) (2) (3) () {f n (x)} n=1 [a, b] K > 0 n, x f n (x) K < ( ) x [a, b] lim f n (x) f(x) (1) f(x)?
More information2010 II / y = e x y = log x = log e x 2. ( e x ) = e x 3. ( ) log x = 1 x 1.2 Warming Up 1 u = log a M a u = M a 0
2010 II 6 10.11.15/ 10.11.11 1 1 5.6 1.1 1. y = e x y = log x = log e x 2. e x ) = e x 3. ) log x = 1 x 1.2 Warming Up 1 u = log a M a u = M a 0 log a 1 a 1 log a a a r+s log a M + log a N 1 0 a 1 a r
More informationI , : ~/math/functional-analysis/functional-analysis-1.tex
I 1 2004 8 16, 2017 4 30 1 : ~/math/functional-analysis/functional-analysis-1.tex 1 3 1.1................................... 3 1.2................................... 3 1.3.....................................
More informationII R n k +1 v 0,, v k k v 1 v 0,, v k v v 0,, v k R n 1 a 0,, a k a 0 v 0 + a k v k v 0 v k k k v 0,, v k σ k σ dimσ = k 1.3. k
II 231017 1 1.1. R n k +1 v 0,, v k k v 1 v 0,, v k v 0 1.2. v 0,, v k R n 1 a 0,, a k a 0 v 0 + a k v k v 0 v k k k v 0,, v k σ kσ dimσ = k 1.3. k σ {v 0,...,v k } {v i0,...,v il } l σ τ < τ τ σ 1.4.
More information7. 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 information1 I
1 I 3 1 1.1 R x, y R x + y R x y R x, y, z, a, b R (1.1) (x + y) + z = x + (y + z) (1.2) x + y = y + x (1.3) 0 R : 0 + x = x x R (1.4) x R, 1 ( x) R : x + ( x) = 0 (1.5) (x y) z = x (y z) (1.6) x y =
More information4. ϵ(ν, T ) = c 4 u(ν, T ) ϵ(ν, T ) T ν π4 Planck dx = 0 e x 1 15 U(T ) x 3 U(T ) = σt 4 Stefan-Boltzmann σ 2π5 k 4 15c 2 h 3 = W m 2 K 4 5.
A 1. Boltzmann Planck u(ν, T )dν = 8πh ν 3 c 3 kt 1 dν h 6.63 10 34 J s Planck k 1.38 10 23 J K 1 Boltzmann u(ν, T ) T ν e hν c = 3 10 8 m s 1 2. Planck λ = c/ν Rayleigh-Jeans u(ν, T )dν = 8πν2 kt dν c
More informationChap10.dvi
=0. f = 2 +3 { 2 +3 0 2 f = 1 =0 { sin 0 3 f = 1 =0 2 sin 1 0 4 f = 0 =0 { 1 0 5 f = 0 =0 f 3 2 lim = lim 0 0 0 =0 =0. f 0 = 0. 2 =0. 3 4 f 1 lim 0 0 = lim 0 sin 2 cos 1 = lim 0 2 sin = lim =0 0 2 =0.
More information6 2 2 x y x y t P P = P t P = I P P P ( ) ( ) ,, ( ) ( ) cos θ sin θ cos θ sin θ, sin θ cos θ sin θ cos θ y x θ x θ P
6 x x 6.1 t P P = P t P = I P P P 1 0 1 0,, 0 1 0 1 cos θ sin θ cos θ sin θ, sin θ cos θ sin θ cos θ x θ x θ P x P x, P ) = t P x)p ) = t x t P P ) = t x = x, ) 6.1) x = Figure 6.1 Px = x, P=, θ = θ P
More information<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>
電気電子数学入門 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/073471 このサンプルページの内容は, 初版 1 刷発行当時のものです. i 14 (tool) [ ] IT ( ) PC (EXCEL) HP() 1 1 4 15 3 010 9 ii 1... 1 1.1 1 1.
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 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 information1 yousuke.itoh/lecture-notes.html [0, π) f(x) = x π 2. [0, π) f(x) = x 2π 3. [0, π) f(x) = x 2π 1.2. Euler α
1 http://sasuke.hep.osaka-cu.ac.jp/ yousuke.itoh/lecture-notes.html 1.1. 1. [, π) f(x) = x π 2. [, π) f(x) = x 2π 3. [, π) f(x) = x 2π 1.2. Euler dx = 2π, cos mxdx =, sin mxdx =, cos nx cos mxdx = πδ mn,
More information名古屋工業大の数学 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 information