Alpha as property

[aquanumerical]

In Depomod.Stats.SpanningEllipse an Ellipse is fitted which has parameters accessed as:
SE.Centre;
SE.MajorSemiAxis;
SE.MinorSemiAxis;

I believe the missing rotation angle parameter can be obtained from:
[eigVec eigVal] = eig(SE.A);
sigmaSE = atan2(eigVec(2,1), eigVec(1,1));

I would suggest exposing it directly as SE.Alpha

    properties
        EnclosingPoints = [];
        Centre = [];
        MajorSemiAxis@double;
        MinorSemiAxis@double;
        **Alpha@double;**
....
        function SE = SpanningEllipse(points)
            SE.EnclosingPoints = points;
            
            [SE.A, SE.Centre] = Depomod.Stats.Utils.MinVolEllipse(points', .01);
            [SE.U SE.D SE.V] = svd(SE.A);

            SE.MajorSemiAxis = 1/sqrt(SE.D(1,1));
            SE.MinorSemiAxis = 1/sqrt(SE.D(2,2));
            **[eigVec eigVal] = eig(SE.A);
            SE.Alpha = atan2(eigVec(2,1), eigVec(1,1));**
        end

test case:
X = [1,1;-0.5,0.5;-1,-1;0.5,-0.5];
SE = Depomod.Stats.SpanningEllipse(X);
[eigVec eigVal] = eig(SE.A);
alphaSE = atan2(eigVec(2,1), eigVec(1,1));
% although they are equivalent the rotation angle is reported as -3*pi/4 we might expect pi/4?