sabit nokta çözümü(fortran90)

[Ahmet EFE]

PROGRAM sabitnokta
IMPLICIT NONE
REAL::X,BETA,X0
INTEGER::N,I
WRITE(*,*)"KÖK İÇİN İLK TAHMİN DEĞERİNİ GİRİN"
READ(*,*)X0
WRITE(*,*)"HATA ÜST SINIRINI GİR"
READ(*,*)BETA
WRITE(*,*)"MAKSİMUM İTERASYON SAYISINI GİR"
READ(*,*)N
DO I=1,N
X=G(X0)
X0=X
WRITE(*,*)X
IF(ABS(X-X0).LE.BETA)THEN
WRITE(*,*)"KÖK=",X
STOP
END IF
X0=X
END DO
WRITE(*,*)"KÖK YAKINSAMIYOR" 
STOP
CONTAINS
FUNCTION G(X)
REAL::G
REAL,INTENT(IN)::X
G=(X**2+2)/3
RETURN
END FUNCTION G
END PROGRAM sabitnokta

ikinci dereceden denklem çözümü için...

denklem: (x**2+2)/3 için yapılmıstır program G=(X**2+2)/3 Burayı değiştirerek programı farklı denklemler için kullanabilirsiniz...

[yusufcelebi]

Eline sağlık kanka yarında sınav vardı )) ))