ONLY DO WHAT ONLY YOU CAN DO

こけたら立ちなはれ 立ったら歩きなはれ

連立一次方程式

Clojure で連立一次方程式を解く(ヤコビの反復法)

例として を考える. この方程式を上から順に対角線上の変数について解くと となる. に適当な値を入れて右辺を計算し、 得られた値を新たなとして、計算を繰り返す. 漸化式で書くと (def N 4) (def a [[9.0 2.0 1.0 1.0] [2.0 8.0 -2.0 1.0] [-1.0 -2.0 7.0 -…

Haskell で連立一次方程式を解く(ガウスの消去法)

連立一次方程式 が与えられた場合, まず(1)式を9で割って(1)'とする. (2)式から(1)'式に2を掛けたものを引くとの項が消える. (3)式に(1)'式を足すとの項が消える. (4)式から(1)'式を引くとの項が消える.同様にして, 2番目の式を利用して, 3番目以降の式のyを…

Haskell で連立一次方程式を解く(ヤコビの反復法)

例として を考える. この方程式を上から順に対角線上の変数について解くと となる. に適当な値を入れて右辺を計算し、 得られた値を新たなとして、計算を繰り返す. 漸化式で書くと import Text.Printf import Debug.Trace import Control.Monad n = 4::Int d…

VB で連立一次方程式を解く(ヤコビの反復法)

例として を考える. この方程式を上から順に対角線上の変数について解くと となる. に適当な値を入れて右辺を計算し、 得られた値を新たなとして、計算を繰り返す. 漸化式で書くと Option Explicit Module VB1001 Private Const N As Integer = 3 Public Sub…

Java で連立一次方程式を解く(ヤコビの反復法)

例として を考える. この方程式を上から順に対角線上の変数について解くと となる. に適当な値を入れて右辺を計算し、 得られた値を新たなとして、計算を繰り返す. 漸化式で書くと import java.lang.*; public class Java1001 { private static final int N …

Go言語 で連立一次方程式を解く(ヤコビの反復法)

例として を考える. この方程式を上から順に対角線上の変数について解くと となる. に適当な値を入れて右辺を計算し、 得られた値を新たなとして、計算を繰り返す. 漸化式で書くと package main import "fmt" import "math" const N = 4 func main() { var a…

D言語 で連立一次方程式を解く(ヤコビの反復法)

例として を考える. この方程式を上から順に対角線上の変数について解くと となる. に適当な値を入れて右辺を計算し、 得られた値を新たなとして、計算を繰り返す. 漸化式で書くと import std.stdio; import std.math; const int N = 4; void main(string[] …

C++ で連立一次方程式を解く(修正コレスキー法)

行と列とを入れ替えても(転置行列)一致する行列を対称行列と言う. 修正コレスキー法では, 対称行列 を下三角行列 , 対角行列 , の転置行列 の積に分解し( ), から を求め, から を求める. #include <iostream> #include <iomanip> #include <math.h> using namespace std; void decomp(d</math.h></iomanip></iostream>…

C++ で連立一次方程式を解く(コレスキー法)

行と列とを入れ替えても(転置行列)一致する行列を対称行列と言う. コレスキー法では, 対称行列 を下三角行列 と, の転置行列 との積に分解し( ), から を求め, から を求める. #include <iostream> #include <iomanip> #include <math.h> using namespace std; void decomp(double a[4][</math.h></iomanip></iostream>…

C++ で連立一次方程式を解く(LU分解法)

連立一次方程式 を行列で表すと, こんな感じ このとき とすると, 最初の連立方程式は と表すことができる. このとき となる 上三角行列 , 下三角行列 を考えると, であり, とおくと となる.LU分解法では, 係数行列 を上三角行列 , 下三角行列 に分解し, から…

C++ で連立一次方程式を解く(ガウス・ジョルダン法)

考え方は, ガウスの消去法とほぼ同じ.連立一次方程式 が与えられた場合, まず(1)式を9で割って(1)'とする. (2)式から(1)'式に2を掛けたものを引くとの項が消える. (3)式に(1)'式を足すとの項が消える. (4)式から(1)'式を引くとの項が消える.ここまでは, ガウ…

C++ で連立一次方程式を解く(ガウスの消去法)

連立一次方程式 が与えられた場合, まず(1)式を9で割って(1)'とする. (2)式から(1)'式に2を掛けたものを引くとの項が消える. (3)式に(1)'式を足すとの項が消える. (4)式から(1)'式を引くとの項が消える.同様にして, 2番目の式を利用して, 3番目以降の式のyを…

C++ で連立一次方程式を解く(ガウス・ザイデル法)

考え方は、ヤコビ法とほぼ同じ. ヤコビ法の場合、漸化式で書くと だったが、2番目以降の式のの代わりに, 1番目の式で計算して得られたを使う. 同様に、3番目以降の式のの代わりに, 2番目の式で計算して得られたを使う. まとめると、 #include <iostream> #include <iomanip> #in</iomanip></iostream>…

C++ で連立一次方程式を解く(ヤコビの反復法)

例として を考える. この方程式を上から順に対角線上の変数について解くと となる. に適当な値を入れて右辺を計算し、 得られた値を新たなとして、計算を繰り返す. 漸化式で書くと #include <iostream> #include <iomanip> #include <math.h> using namespace std; void jacobi(double a</math.h></iomanip></iostream>…