What is an Acceptable Three-Phase (Two Winding) Saturation Curve?

When providing a saturation curve, there seems to be two limitations: 1. The slope between two data points, and 2. The number of datapoints. The real saturation curve would result in an error during energization of the transformer (i.e. zero to nominal voltage). Altering it to the originalSaturation data by elininating the last data point and increasing the rise of the slope of the second to last data point fixed it. In order to obtain a working saturation curve for simulations wiht initial flux of 1pu, multiple data points had to be removed from the linear region.

In short, why do these limmitations occur?

realSaturation = [0 0; 0.007278 0.4386 0.020380 0.7859 0.036390 0.9745 0.082970 1.1760 0.147000 1.2650 0.326100 1.3370 0.751100 1.3740 1.278000 1.3860 2.964000 1.3960 6.204000 1.4000 9.279000 1.4010 20 1.4011 100 1.4012];

originalSaturation = [0 0; 0.007278 0.4386 0.020380 0.7859 0.036390 0.9745 0.082970 1.1760 0.147000 1.2650 0.326100 1.3370 0.751100 1.3740 1.278000 1.3860 2.964000 1.3960 6.204000 1.4000 9.279000 1.401 20 1.4025];

remnantSaturation = [0 0; 0.00 0.4386 0.0 1.2650 0.3261 1.3370 0.751100 1.3740 1.278000 1.3860 2.964000 1.3960 6.204000 1.4000 9.279000 1.401 20 1.4025];