[原创]自己写的二分法求方程的根FORTRAN源程序(练手儿用的)

houston

[watermark]
PROGRAM ERFENFA
IMPLICIT NONE
REAL*8 F
REAL*8 A,B,C
REAL*8 EMISO,E
EMISO=0.5E-5
PRINT*,';[ERFENFA V.1.0 -(C) BAICANG LIU 10 2004]';
WRITE(*,*)"LEASE INPUT A B"
READ(*,*)A,B
IF(F(A)*F(B).GT.0.0)THEN
WRITE(*,*)"LEASE INPUT A B AGAIN"
WRITE(*,*)"NOW THE ROOT IS NOT IN THE DOMAIN"
ELSE
E=1.0
DO WHILE(E.GT.EMISO)
C=(A+B)/2.0
IF(F(A)*F(C).GE.0.0)THEN
IF(F(A)*F(C).EQ.0.0)THEN
WRITE(*,*)C
CONTINUE
ELSE
A=C
ENDIF
ELSE
B=C
ENDIF
E=ABS(B-A)
ENDDO
WRITE(*,*)';THE ROOT OF THE EQUATION IS:';,C
ENDIF
END
REAL*8 FUNCTION F(X)
REAL*8 X
F=SIN(X)-X**2.0/2.0
END FUNCTION
!REAL*8 FUNCTION F(X)
!REAL*8 X
!F=X**3.0-X-1.0
!END FUNCTION[/watermark]

dut

二分法收敛速度比较慢，你的程序中对于多个根的情况没有考虑！另外，误差以及区间长度的设定是一个很重要的问题。

houston

写该段程序的目的,主要是让大家了解一下FOR90写的程序,其实根FOR77差不太多,有一些新的功能,如MODULE等很不错.

taibaisha

挺好,等我会了,些个好的给你
我脸皮挺厚吧,哈哈!