/* ---------------------------------- */
/* save as fe.h                      */
/* --------------------------------- */
double f(
double x,
double y,
double z)
{
 return( (1) );
}
char  feq[] = "1";
/* --------------------------------- */
double uz(
double x,
double y)
{
 return(0);
}
char  ueq[] = "0";
/* --------------------------------- */
double vz(
double x,
double y)
{
 return((4-y));
}
char  veq[] = "4-y";
/* --------------------------------- */
double sx(
double y)
{
 return((-sqrt(y)));
}
char  seq[] = "-sqrt(y)";
/* --------------------------------- */
double tx(
double y)
{
 return((sqrt(y)));
}
char  teq[] = "sqrt(y)";
/* --------------------------------- */