数値解析

【numpy】行列操作【その1】

numpyの行列操作で面白い行列の作り方を知ったので書く. a = np.arange(1,11) b = np.arange(0,10) #b: [0,1,2,3,4,5,6,7,8,9] c = b[:,np.newaxis] #c: #[[0] # [1] # [2] # [3] # [4] # [5] # [6] # [7] # [8] # [9]] d = a + c #[[ 1 2 3 4 5 6 7 8 9 10]…

Newton法で平方根を求める

(define (sqrt-iter guess x) (if (good-enough? guess x) guess (sqrt-iter (imporve guess x) x))) (define (improve guess x) (average guess (/ x guess))) (define (average x y) (/ (+ x y) 2)) (define (good-enough? guess x) (< (abs (- (square gu…

ガウスジョルダン法(2)

// 行列Aを上三角行列にして解く方法 上三角行列に出来れば、ベクトルX=[x,y,z]についてz→y→xの順に一つずつ解が求まり簡単だ gist.github.com

ガウスジョルダン法(1)

「数値解析」の講義で示されたプログラムをC++で書きました. // gist.github.com

桁落ちについて

C言語の講義で「桁落ち」という現象を扱いました. 浮動小数点において入力された順に足していくと不確かな値が出力されてしまう場合があります. ↓ ↓ ↓ 情報落ち、桁落ち、丸め誤差、打切り誤差の違い ↑ ↑ ↑ こちらに詳しく載っていたので解説はこちらで~ 桁…