Update quadratic.for

Cleaned up code, made the write statement for the header fortran66 compliant