/* --------------------------------- */
/* save as fa.h */
/* --------------------------------- */
double f(
double x,
double y,
double z)
{
return( y*sin(x*y*z)*z );
}
char feq[] = "y*sin(x*y*z)*z";
/* --------------------------------- */
double f_x(
double x,
double y,
double z)
{
return( (y*y*z*z*cos(x*y*z)) );
}
char Dfeq_x[] = "y**2*z**2*cos(x*y*z)";
/* --------------------------------- */
double f_y(
double x,
double y,
double z)
{
return( (x*y*z*z*cos(x*y*z)+z*sin(x*y*z)) );
}
char Dfeq_y[] = "x*y*z*z*cos(x*y*z)+z*sin(x*y*z)";
/* --------------------------------- */
double f_z(
double x,
double y,
double z)
{
return( (x*y*y* z*cos(x*y*z)+y*sin(x*y*z)) );
}
char Dfeq_z[] = "x*y**2*z*cos(x*y*z)+y*sin(x*y*z)";
/* --------------------------------- */
double f_xx(
double x,
double y,
double z)
{
return( (-pow(y,3)*pow(z,3)*sin(x*y*z) ) );
}
char Dfeq_xx[] = "-y**3*z**3*sin(x*y*z)";
/* --------------------------------- */
double f_yy(
double x,
double y,
double z)
{
return((2*x*z*z*cos(x*y*z)-x*x*y*pow(z,3)*sin(x*y*z)));
}
char Dfeq_yy[] = "2*x*z**2*cos(x*y*z)-x**2*y*z**3*sin(x*y*z)";
/* --------------------------------- */
/* --------------------------------- */
double f_zz(
double x,
double y,
double z)
{
return((2*x*y*y*cos(x*y*z)-x*x*z*pow(y,3)*sin(x*y*z)));
}
char Dfeq_zz[] = "2*x*z**2*cos(x*y*z)-x**2*z*y**3*sin(x*y*z)";
/* --------------------------------- */