Codes: part 33

M:matrix([%e^(%i*a),0,0],[0,%e^(%i*b),0],[0,0,%e^(%i*c)]);
expand(determinant(subst(-b-c,a,M)));

Advertisements