Question: How should Maple 2016 replace evalm in this procedure?

plot3d of procedure Sievert correctly displays the constant curvature Sievert surface, but the procedure uses the deprecated command evalm.

What Maple 2016 statement(s) would create the same value of X in Sievert?

Sievert := proc (B)

local a, b, denom, m, X;

a := sinh(B)*u; b := cosh(B)*v;

denom := sinh(B)*((cosh(2*a)-cos(2*b))*cosh(2*B)+2+cosh(2*a)+cos(2*b));

m := cosh(B)*[sinh(a), sin(b)*cos(v), sin(b)*sin(v)]+[0, -cos(b)*sin(v), cos(b)*cos(v)];

X := evalm([u, 0, 0]-8*cosh(B)*cosh(a)*m/denom);

end proc:

plot3d(Sievert(.75), u = -2.5 .. 2.5, v = -10.5 .. 10.5, scaling = constrained, grid = [30, 100], style = patch, shading = xy, lightmodel = light3, orientation = [-3, 140], title = "Sievert's surface", titlefont = [Courier, bold, 14]);