Your code starts with "PROGRAM MAIN", but this main routine ends with a RETURN statement, that is not correct.
You assign double precision values to ZERO,ONE,TWO,and THREE, but double precision has not been declared for these or any other variables.
Many variables like C10, C01, C11, C20, C30...