
                    2
cs1(1,1,3) :=  - a*u

                                   2  2
cs1(1,2,2) := cos(x(1))*sin(x(1))*a *u

                 2
cs1(1,3,1) := a*u

                                   2  2
cs1(2,1,2) := cos(x(1))*sin(x(1))*a *u

                                      2  2
cs1(2,2,1) :=  - cos(x(1))*sin(x(1))*a *u

                          2    2
cs1(2,2,3) :=  - sin(x(1)) *a*u

                       2    2
cs1(2,3,2) := sin(x(1)) *a*u

                 2
cs1(3,1,1) := a*u

                       2    2
cs1(3,2,2) := sin(x(1)) *a*u

