That's exactly what I am doing, in fact I have about 6 or 7 sucessful programs made in Fortran 95, here a copy of one of them:
! New Quadratic formula calculator
program new_quad_solve
implicit none
real :: a, b, c
real :: disc, root1, root2
print *, ' Enter values for a, b...