Korzystanie z istniejących podprogramów

W pliku /big/staff/adam/FORTRAN/Cw6/jacobi.f na komputerze etoh.chem.univ.gda.pl znajduje się procedura diagonalizacji symetrycznej macierzy kwadratowej metodą Jacobiego. Przeczytaj komentarze na początku podprogramu objaśniające znaczenie parametrów formalnych i następnie napisz program obliczający poziomy energetyczne oraz orbitale molekularne węglowodoru o wiązaniach podwójnych sprzężonych metodą Hückla. Tablica sąsiedztwa ma być wczytywaną z pliku lub klawiatury. Działanie programu sprawdź na przykładzie etylenu, 1,3-butadienu oraz benzenu.

Tablica sąsiedztwa węglowodoru o wiązaniach podwójnych ma wymiar (liczba atomów)$\times$(liczba atomów); jeżeli między atomami i i j jest wiązanie podwójne to wartość elementu i,j tablicy sąsiedztwa wynosi 1, w przeciwnym wypadku 0. Elementy diagonalne tablicy sąsiedztwa są równe 0. Energie poszczególnych poziomów to wartości własne a współczynniki w orbitalach molekularnych to wektory własne tablicy sąsiedztwa.