I think I found my problem. I was choosing a value of a, then getting the values of c, the compression area (Acomp), the distance to the compression area centroid (y'), the stress in the steel (fs = (0.003/c)*(d-c)), then the phi factor (relating to the c / d ratio), then phi*Pn =...