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