原创 高斯積分
!// 注: 參考來源http://fcode.cn/algorithm-73-1.html Module Gauss_Legendre !//高斯—勒讓德積分高斯點及權重的求解模塊 Implicit none Intege
2020-02-25 23:14:06
原创 拉格朗日插值fortran程序
Program LagrangianInterpolation Implicit none Integer :: i, fileid Integer, parameter :: n1 = 51 !// 插值基節點數
2020-02-25 23:14:06
1
原创 切比雪夫插值fortran程序
!// Chebyshev interpolation Module interpolation Contains Subroutine newtdd( x,y,c ) !// 計算牛頓插值係數 Implicit none
2020-02-25 23:14:06
原创 使用置換矩陣求解線性方程組2
!// 說明：至此，第二章有關方程求解的部分就此結束。下節我們正式進入第二章的迭代部分，下期見！！！ !// -----------原方程----------- !// 2x + 1y + 5z = 5 !// 4x
2020-02-25 23:14:06
原创 迭代法求解非線性問題
!// 方程組：-u^3 + v = 0 !// u^2 + v^2 - 1 = 0 !// 即：f1(u,v) = -u^3 + v f2(u,v) = u^2 + v^2 - 1 !// 多元牛頓法的求解適用於小型
2020-02-25 23:14:06
原创 python函數文檔小記
python可以對函數進行函數文檔說明，這樣可以大大增加代碼的可閱讀性。 def func1(x, y): "my first function!" z = x + y return z print(func
2020-02-25 23:14:06
2
原创 完全QR分解
Program QR_complete Implicit none Integer :: i, j Integer, parameter :: m = 4, n = 2 Real(kind=8) :: A(m,n), Q
2020-02-25 23:14:06
原创 預條件共軛梯度法求解線性方程組2
!// 在病態矩陣問題上，共軛梯度法比部分主元的高斯消去法還要差 !// 這種情況下，可通過預條件得到緩解，主要是將問題轉化爲良態矩陣系統，再實施CG方法 !// 此代碼的預條件因子M= ( D + wL ) * Inv(D) * (
2020-02-25 23:14:06
原创 高斯消去法求解線性方程組
!//---------------原方程------------- !// x + y + z + w = 10 !// 2x + 3y + z + w = 15 !// 3x - y + 2z - w = 3 !
2020-02-25 23:14:06
原创 使用切線法求解方程零點
Program SecantMethod Implicit none Real(kind=8) :: x0 = 0.d0, x1 = 1.d0 Real(kind=8) :: xtmp Real(kind=8), pa
2020-02-25 23:14:06
原创 使用置換矩陣求解線性方程組1
!// 說明：至此，第二章有關方程求解的部分就此結束。下節我們正式進入第二章的迭代部分，下期見！！！ !// -----------原方程----------- !// 2x + 1y + 5z = 5 !// 4x
2020-02-25 23:13:56
原创 梯形積分與辛普森積分
Program main !// 辛普森積分 Implicit none Real(kind=8), parameter :: a = 0.d0, b = 1.d0 !//積分區間爲[a,b], 被積函數爲f(x) = x
2020-02-25 23:13:56
原创 matplotlib繪圖：調整座標刻度朝內/外
import numpy as np import matplotlib.pyplot as plt X = np.linspace(-np.pi, np.pi, 256, endpoint = True) C, S = np.cos
2020-02-25 23:13:56
6
原创 複合體型辛普森積分
Program main !// 複合梯形積分和辛普森積分 Implicit none Real(kind=8), parameter :: a = 0.d0, b = 1.d0 !//積分區間爲[a,b], 被積函數爲f
2020-02-25 23:13:56
原创 牛頓插值fortran程序
!// Newton interpolation Module interpolation Contains Subroutine newtdd( x,y,c ) !// 計算牛頓插值係數 Implicit none Rea
2020-02-25 23:13:56