S1 = 6*cosh(r)*sinh(r)*cos(θ）+3+6 *sinh(r)^2 ,
r和θ 取值范围【-4，4】

绘出这个方程的图

求编程过程

最近在学习MATLAB，老师布置了这个问题，实在是搞不懂啊，希望得到您的解答

x=-4:0.1:4;
y=-4:0.1:4;
[X,Y]=meshgrid(x,y);
Z=6*cos(X)*sin(X)*cos(Y)+3+6 *sin(X)^2;
mesh(X,Y,Z)