#include<math.h>
#include<iostream>
using namespace std;

/* Calcul de la densite de probabilite d'un paquet d'ondes */

int main(){
  double a = 1.e5, hbar = 197., t = 1.e5, m = 0.5, x = 1.e5, k0 = 5.e-4;
  double rho, u, v;

  u = 1 + 4 * hbar*hbar * t*t / (m*m * pow(a,4));
  v = x - hbar * k0 * t / m;
  rho = sqrt(2./M_PI) * 1./(a*sqrt(u)) * exp(-2.*v*v/(a*a*u));
  cout << "Densite de probabilite (1/fm) : " << rho << endl;
  return 0;
}