#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-1D.dat");
  
  double a=uniform(),b=uniform(),c=rand()%5,d=rand()%5;
  
  for(int i=0;i<101;i++)
  {
    double x = i*0.01;
    output << x << " " << a*x*(x-b) + d*sin(c*Pi*x) << "\n";
  }
}