#include <iostream>
#include <cmath>

using namespace std;

int main()
{
  const double E=10. ; //10 eV
  const double PI=3.141592;
  const double k0=sqrt(2.*5.11e5*E)/1.97e2;
  const double delk0=0.5/0.1;
  const double lambda=2.*PI*1.97e2/sqrt(2.*5.11e5*E);
  //cout <<"lambda: "<< lambda << "k0: " << k0 << "delk0: " << delk0 << endl;
  // p85 Iwanami QM text book. delx=0.1nm.
  double RealPhi;
  double k=k0;
  for(double x=-2.; x< 2.; x+=0.02){
    RealPhi=0;
    for(int i=0; i < 21; i++){
      k=k0+(double)(i-10)*delk0*0.2;
      double Ak=exp(-(k-k0)*(k-k0)/(4.*delk0*delk0));
      RealPhi+=Ak*cos(k*x);
      //cout << doremi[i] << " " << freq << endl;
    }
    cout << x << " "<< RealPhi << endl;
  }
  
  return 0;
}