Ćwiczenia

1. Korzystając z programu na diagonalizację macierzy opracowanego w ćwiczeniu nr 6 oraz z załączonych plików źródłowych napisz program huckel obliczający energie $\pi$-elektronowe oraz orbitale molekularne cząsteczki o wiązaniach sprzężonych w przybliżeniu $\pi$-elektronowym metodą Hückla. Moduł głowny ma być w pliku huckel.f a program ma działać następująco:
   1. Wczytywać hamiltonian $\bf H$, danej cząsteczki, który dla węglowodorów sprzężonych jest równy macierzy sąsiedztwa. Element macierzy sąsiedztwa jest równy 1, jeżeli atomy są połączone wiązaniem i 0 jeżeli nie. Macierz sąsiedztwa należy wczytać podając współrzędne i wartości niezerowych elementów. Przykładowo, dla cyklobutadienu będzie to wyglądało następująco:

      4          ! liczba atomow
      1 2 1.0    ! 1 atom polaczony z 2-gim; wartosc elementu 1
      2 3 1.0    
      3 4 1.0
      1 4 1.0
      

      Program ma sam uzupełniać elementy położone symetrycznie, czyli dla powyższego przykładu ma automatycznie nadać wartości 1 elementom (2,1), (3,2), (4,3), i (4,1). Uwaga! Jeżeli zamiast 1.0 będzie co innego to mamy inną siłę wiązania $\pi$-elektronowego i może to np. odpowiadać związkowi heterocyklicznemu. W tym przypadku należy też podać explicite elementy diagonalne równe odpowiednim całkom rdzeniowym.

   2. Drukować macierz $\bf H$ w postaci dolnego trójkąta; do tego służy procedura matout.f

   3. Przeprowadzać diagonalizację macierzy $\bf H$; do tego służy procedura jacobi.f.

   4. Drukować energie orbitalne (wartości własne) i orbitale molekularne (wektory własne macierzy $\bf H$); do tego służy procedura eigout.f.

2. Napisz Makefile najpierw w formie najprostszej a potem wykorzystując reguły przyrostkowe (suffix rules).