Reorient Ellipsoid Gained From Ellipsoid Function

I have an ellipsoid I would like to plot. I know that I can do so using the ellipsoid function (link) as follows:

[X,Y,Z] = ellipsoid(xc,yc,zc,xr,yr,zr)
surf(X,Y,Z);

This gives an ellipsoid with mean equal to [xc,yc,zc] and semiaxes with length equal to xr, yr, and zr respectively. However, the semiaxes are necessarily the canonical basis for R3: [1,0,0], [0,1,0], and [0,0,1].

Given vectors three 3-dimensional orthonormal vectors u1, u2, u3, how can I orient the ellipsoid above so that the ith semiaxis is oriented in the direction of ui?

Thank you.