#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 x_vector("x_vector.dat");
  std::ofstream y_vector("y_vector.dat");
  std::ofstream results_matrix("results_matrix.dat");
  
  double a=uniform(),b=uniform(),c=rand()%5,d=rand()%5;
  
  for(int i=0;i<101;i++)x_vector << i*0.01 << "\n";
  for(int j=0;j<101;j++)y_vector << j*0.01 << "\n";
  
  for(int i=0;i<101;i++)
  {
    double x = i*0.01;
    for(int j=0;j<101;j++)
    {
      double y = j*0.01;
      results_matrix << a*x*x+exp((y-b)*(y-b)) + sin(c*Pi*(x+d*y)) << " ";
    }
    results_matrix << "\n";
  }
}