/* Test function for 2d plots */

Float foit_f (Float x, Float y)
  {
    ROUND_NEAR;
    return (x*x + y*y + x*y - Quarter);
  }

Interval foit_fv (Interval x, Interval y)
  {
    Interval x2 = iv_sqr(x);
    Interval y2 = iv_sqr(y);
    Interval xy = iv_mul(x, y);
    Interval sum = iv_add(iv_add(x2, y2), xy);
    Interval res = iv_affine (sum, One, -Quarter, Zero);
    return (res);
  }

FOIP foit_ff (FOIP x, FOIP y)
  {
    MemP frame = foi_top();
    FOIP x2  = foi_sqr(x);
    FOIP y2  = foi_sqr(y);
    FOIP xy  = foi_mul(x, y);
    FOIP sum = foi_add(foi_add(x2, y2), xy);
    FOIP res = foi_affine(sum, One, -Quarter, Zero);
    return (foi_return(frame, res));
  }

char *foit_fname = "f = x^2 + y^2 + x*y - 1/4";

Interval foit_fxd = {-One, One};
Interval foit_fyd = {-One, One};

int foit_fn = 32;