https://gcc.gnu.org/bugzilla/show_bug.cgi?id=98701
--- Comment #6 from Bahram.shahrooz at uc dot edu --- I get wrong answers doing this: gfortran FiberAnalysis.f -o FiberAnalysis I get correct answers doing this: gfortran FiberAnalysis.f -fcheck=all -ffpe-trap=zero,invalid -o FiberAnalysis Moment Axial_Load 3244.08910 -215.38970 3432.22199 -190.65966 3594.93500 -168.86785 3734.50302 -149.77519 3842.77168 -134.14558 3986.24344 -114.80934 4115.25069 -97.21038 4230.37015 -81.27311 4323.74191 -67.84354 4442.09009 -51.57463 4551.67144 -36.41379 4653.40529 -22.24901 4748.46682 -8.91886 4844.11684 4.77310 4933.55651 17.64818 5017.35901 29.77892 5097.44257 41.41472 5193.49624 55.17913 5275.53074 67.17330 5351.16344 78.34124