Cześć, dostałam gotowy program napisany w Fortranie i moim zadaniem jest napisanie go w Matlabie. Niestety nie znam się na tym pierwszym języku i chociaż mniej więcej wiem jak działa zadany program mam kłopoty z zamianą formatu programu.
Program dotyczy zastosowania metody strzałów przy rozwiązywaniu równania różniczkowego drugiego stopnia (znalezienia wartości własnych)
Bardzo proszę o wszelkie wskazówki i pomoc.
Program w Fortranie:
REAL K
K=1.
DK=1.
TOLK=1.E-05
CALL INTGRT(K,PHIP)
PHIOLD=PHIP
10 CONTINUE
K=K+DK
CALL INTGRT(K,PHIP)
IF (PHIP+PHIOLD .LT. 0) THEN
K=K-DK
DK=DK/2
ENDIF
IF (ABS(DK) .GT. TOLK) GOTO 10
EXACT=4.*ATAN (1.)
PRINT *,' EIGENVALUE, ERROR = ', K, EXACT-K
STOP
END
SUBROUTINE INTGRT(K, PHIP)
REAL K
DATA NSTEP/100/
H=1./NSTEP
PHIM=0.
PHIZ=.01
CONST=(K*H)**2/12
DO 10 IX=1, NSTEP-1
PHIP=2*(1.-5.*CONST)*PHIZ-(1.+CONST)*PHIM
PHIP=PHIP/(1+CONST)
PHIM=PHIZ
PHIZ=PHIP END
CONTINUE
PRINT *,K,PHIP
RETURN
END
Pozdrawiam