Program główny i nagłówek podprogramu kwadrat powinny mieć następującą strukturę:
program rownanie
c Program rozwiazuje rownanie kwadratowe
implicit real*8 (a-h,o-z)
c Czytanie danych
read (*,*) a,b,c
c Obliczanie pierwiastkow
call kwadrat(a,b,c,x1,x2,*10,*20)
print *,"Pierwiastki rownania:",x1,x2
stop
10 print *,"Rownanie nie jest rownaniem kwadratowym"
stop
20 print *,"Rownanie nie ma pierwiastkow rzeczywistych"
stop
end
subroutine kwadrat(a,b,c,x1,x2,*,*)
c Podprogram rozwiazuje rownanie kwadratowe a*x**2+b*x+c dla danych
c wspolczynnikow a, b, c. Rozwiazania sa zapamietane w x1 i x2
implicit real*8 (a-h,o-z)
c
c Sprawdzanie wspolczynnikow
c Jezeli wspolczynniki nieprawidlowe to nastepuje return 1
c Wpisz tutaj odpowiednia instrukcje
c
c Obliczanie delty
c Jezeli delta ujemna to nastepuje return 2
c Wpisz tutaj odpowiednia instrukcje
c
c Obliczanie pierwiastkow x1 i x2
c
return
end