How to insert legends in the surfdata?

I have calculated z=f(x,y) for x from 0 to 5 and y from 0 to 5 by taking step size 0.5  and stored it in a matrix M. Now I am trying to plot the resultant matrix having order 11x11. Where x and y values are ranging from 0 to 5. But If I use this command
matrixplot(M, heights = histogram, colorscheme = ["Blue", "Green", "Yellow", "Red"])
It is taking 0 to 11 on both x and y axis. It must be 0,0.5,0.1,0.15...5.
How to change the values of x and y axis manually?

Given an almost contact metric manifold M(\phi,\xi,\eta, g), we say
that M is a generalized Sasakian-space-form if there exist three functions f1, f2, f3
on M such that the curvature tensor R is given by

R(X,Y)Z=f_1{g(Y,Z)X-g(X,Z)Y}+f_2{g(X,\phiZ)\phiY-g(Y,\phiZ)\phiX+2g(X,\phiY)\phiZ}+f_3{\eta(X)\eta(Z)Y-\eta(Y)\eta(Z)X+gg(X,Z)\eta(Y)\xi-g(Y,Z)\eta(X)\xi}

In (2n+1) dimensional generalized Sasakian space form M²ⁿ⁺¹(f_1,f_2,f_3), we have the following relations.

S(X,Y)=(2nf_1+3f_2-f_3)g(X,Y)-(3f_2+(2n-1)f_3)\eta(X)\eta(Y)

S(X,\xi)=2n(f_1-f_3)\eta(X)

C\bar(\xi,X)Y=[f_1-f_3-(r/2n(2n-1))][g(X,Y)\xi-\eta(Y)X]

P(X,Y)Z=R(X,Y)Z-(1/(n-1))[S(Y,Z)X-S(X,Z)Y]

R(X,Y)\xi=(f_1-f_3){\eta(Y)X-\eta(X)Y}

R(\xi,X)Y=(f_1-f_3){g(X,Y)\xi-\eta(Y)X}

for any vector fields X, Y on M, where R, S, C\bar, and r denote the Riemannian curvature tensor, Ricci tensor, concircular curvature tensor and scalar curvature of M²ⁿ⁺¹(f1, f2, f3), respectively 

Using above equations I have to evaluate P(C\bar(\xi,X)Y,Z)U.

Manually It is tedious job. Can we find the value by maple? Is there any option to solve these type of problems?

If yes, I can solve many more, which helps a lot in my work.. Thanks in advance.

As I am beginner in maple, how to verify Bianchi Identities?

For Riemannian manifold,

1) R(X, Y)*Z+R(Y, Z)*X+R(Z, X)*Y = 0 

2) ((∇)[X]R)(Y,Z)W+((∇)[Y]R)(Z,X)W+((∇)[Z]R)(X,Y)W=0

I am not getting how to define vector fields. From my previous question I understood defining vector fields particularly. From that post I tried to verify the proof, but I cant. Please tell me about how to define vector fields without taking examples.

Thank in advance.

I am beginner in maple. And my field is Differential geometry. I am trying this calculation through maple.

I have these vector fields e1=z²∗D_x,e2=z²∗D_y,e3=z²∗D_z

Now I have to calculate 2g(∇_e1e3,e1) using maple. where The Riemannian connection ∇ of the metric tensor g is given by Koszul's formula

2g(∇_XY,Z)=X_g(Y,Z)+Y_g(Z,X)−Z_g(X,Y)−g(X,[Y,Z])−g(Y,[X,Z])+g(Z,[X,Y])

By calculated manually I got 2g(∇_e1e3,e1)=2g(−(2/z) e1,e1)

So help me in simplifying tedious calculations using maple. I will try to learn.