/* abc_la.c  just the functions that need to be computed, not the solver */

#include <math.h>

double a1(double x, double y)
{
  return exp(x/2.0)*exp(y)/2.0;
}

double b1(double x, double y)
{
  return 0.7/(x*x*y*y+0.5);
}

double c1(double x, double y)
{
  return (4.0 - exp(x) - exp(y/2.0))*2.0;
}

double d1(double x, double y)
{
  return x*x+y;
}

double e1(double x, double y)
{
  return x*y*y;
}

double f1(double x, double y)
{
  return 3.0*x + 2.0*y;
}

double u(double x, double y)
{
  return x*x*x + 2.0*y*y*y + 3.0*x*x*y + 4.0*x*y*y + 5.0*x*y +
    6.0*x + 7.0*y + 8.0;
}

double c(double x, double y)
{
  return 0.5*exp(x/2.0)*exp(y)*(6.0*x+6.0*y)             +
    0.7*(6.0*x + 8.0*y + 5.0)/(x*x*y*y+0.5)              +
    (8.0 - 2.0*exp(x) - 2.0*exp(y/2.0))*(12.0*y + 8.0*x) +
    (x*x+y)*(3.0*x*x + 6.0*x*y + 4.0*y*y + 5.0*y + 6.0)  +
    x*y*y*(6.0*y*y + 3.0*x*x + 8.0*x*y + 5.0*x +7.0)     +
    (3.0*x + 2.0*y)*(x*x*x + 2.0*y*y*y + 3.0*x*x*y + 4.0*x*y*y +
     5.0*x*y + 6.0*x + 7.0*y + 8.0);
}