jinbeihuang
小新爱
柒指
来讨论下高斯消元法和牛顿迭代法
展开Biu
@@18!!来讨论下高斯消元法和牛顿迭代法
下面这个是用牛顿迭代法计算三元一次方程 a,b,c,d都是自己输入
PROGRAM MAIN_ROOT
IMPLICIT NONE
REAL::ROOT,A,B,C,D,X0,X
write( *,* ) "PLEASE INPUT A,B,C,D,X0"
READ (*,*) A,B,C,D,X0
X=ROOT(A,B,C,D,X0)
write(*,*) X
END PROGRAM MAIN_ROOT
FUNCTION ROOT(A,B,C,D,X0) RESULT(ROOT_RESULT)
IMPLICIT NONE
REAL,INTENT(IN)::A,B,C,D,X0
REAL::ROOT_RESULT
reAL::X1,F,F1
X1 = X0
F=A*X1**3+B*X1**2+C*X1+D
write(*,*) F
DO WHILE(ABS(F).GT.1E-6)
F=A*X1**3+B*X1**2+C*X1+D
F1=3*A*X1*X1+2*B*X1+C
X1 = X1-F/F1
write(*,*) X1
ROOT_RESULT = X1
END DO
END FUNCTION ROOT
DO WHILE(ABS(F).GT.1E-6)想问下LZ这里的括号内的代码怎么理解
[查看全文]
manmiray
jinbeihuang
Utend
终于有个人回了
展开Biu
jinbeihuang 发表于 2014-4-3 23:15终于有个人回了,可能Fortran用在科学计算上比较多,所以接触少些吧
嘿嘿,技术区真的挺冷清的~
对啊,而且现在学校里也都是些热门语言
[查看全文]
jinbeihuang
Utend
jinbeihuang