> restart:
> f0:=x1^3 * sin(x2 + 2.5 * x3) - (x1 + x2)^2 / (1.0 + x3^2);
> f1:=exp(-x1 * x2 + x3) / cos(1.0 + x2 * x3 - x1);
> coefficients:=<0.3, 2.9>;
> f:=f0*f1;
> x:=<1.2,2.3,3.4>;
> value:=evalf(subs(x1=x[1],x2=x[2],x3=x[3],f));

                                                     2
                       3                    (x1 + x2)
               f0 := x1  sin(x2 + 2.5 x3) - ----------
                                                    2
                                            1.0 + x3


                              exp(-x1 x2 + x3)
                     f1 := ----------------------
                           cos(-1.0 - x2 x3 + x1)


                                        [0.3]
                        coefficients := [   ]
                                        [2.9]


           /                                2\
           |  3                    (x1 + x2) |
           |x1  sin(x2 + 2.5 x3) - ----------| exp(-x1 x2 + x3)
           |                               2 |
           \                       1.0 + x3  /
      f := ----------------------------------------------------
                          cos(-1.0 - x2 x3 + x1)


                                   [1.2]
                                   [   ]
                              x := [2.3]
                                   [   ]
                                   [3.4]


                        value := -21.84284310

> gradient:=evalf(subs(x1=x[1],x2=x[2],x3=x[3],<diff(f,x1),diff(f,x2),di
> ff(f,x3)>));

                                  [102.6601407 ]
                                  [            ]
                      gradient := [-292.6801235]
                                  [            ]
                                  [-235.1693312]

> hessian:=evalf[5](subs(x1=x[1],x2=x[2],x3=x[3],<<diff(diff(f,x1),x1),d
> iff(diff(f,x1),x2),diff(diff(f,x1),x3)>|<diff(diff(f,x2),x1),diff(diff
> (f,x2),x2),diff(diff(f,x2),x3)>|<diff(diff(f,x3),x1),diff(diff(f,x3),x
> 2),diff(diff(f,x3),x3)>>));

                         [-877.51    2784.0     2029.8 ]
                         [                             ]
              hessian := [2784.0     -8604.0    -6382.5]
                         [                             ]
                         [2029.8     -6382.5    -4598.0]

> 
