@@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这里的括号内的代码怎么理解