#include <iostream>
#include <fstream>
#include <stdlib.h>
#include <cmath>
#include <ctime>

double uniform()
{
  return 2.*double(rand())/RAND_MAX - 1.;
}

const double Pi=4.*atan(1.);

int main()
{

  srand(time(NULL));
  
  std::ofstream output("results.dat");
  
  double a=uniform(),b=uniform(),c=rand()%5,d=rand()%5;
  
  for(int i=0;i<101;i++)
  {
    double x = i*0.01;
    for(int j=0;j<101;j++)
    {
      double y = j*0.01;
      output << x << " " << y << " " << a*x+exp((y+b)*(y+b)) + cos(c*Pi*(x+d*y));
      output << " " << a*x*x+exp((y-b)*(y-b)) + sin(c*Pi*(x+d*y)) << "\n";
    }
    output << "\n";
  }
}