A numerical analysis was developed to evaluate the complete load-deflection and moment-curvature relationships for square and L-shaped slender reinforced concrete columns subjected to biaxial bending and axial load. This computer model can be used for any cross section geometry and material properties of normal concrete. The analysis was based on a deformation control and both the ascending and descending branches of curves can be studied. The finite difference method was introduced to calculate the deflections which satisfy the compatibility equations. Six square slender columns and eight L-shaped slender columns were tested to compare their experimental load-deformation results with the analytical results derived from the theoretical studies. A satisfactory agreement was achieved for the present study. The results of present study can be used for future design reference.