vector<double> A,B;
double Okres=4.0;
double Fourier( double t)
{
  double y=0;
  if(A.size()!=B.size())
    throw "Blad: Rozne rozmiary tablic A i B.";
  if(A.size()>0) 
    y=0.5*A[0];
  for(int i=1; i<A.size(); i++)
    y+=A[i]*cos(i*2.0*M_PI*t/Okres)+B[i]*sin(i*2.0*M_PI*t/Okres);
  return y;
}