#include <iostream>
#include <ginac/ginac.h>
using namespace GiNaC;
using namespace std;

//================================================//

ex Herm (int n, symbol& x){
  ex H = exp(-power(x,2));
  return normal (power(-1,n) * diff(H, x, n) / H);
}

//================================================//

double Z(int i){
  double z[50];
  z[0] = 0.1107958724;
  z[1] = 0.3324146923;
  z[2] = 0.5541148236;
  z[3] = 0.7759507615;
  z[4] = 0.9979774361;
  z[5] = 1.220250391;
  z[6] = 1.442825970;
  z[7] = 1.665761509;
  z[8] = 1.889115537;
  z[9] = 2.112947996;
  z[10] = 2.337320464;
  z[11] = 2.562296402;
  z[12] = 2.787941424;
  z[13] = 3.014323580;
  z[14] = 3.241513680;
  z[15] = 3.469585636;
  z[16] = 3.698616859;
  z[17] = 3.928688683;
  z[18] = 4.159886855;
  z[19] = 4.392302079;
  z[20] = 4.626030636;
  z[21] = 4.861175092;
  z[22] = 5.097845105;
  z[23] = 5.336158360;
  z[24] = 5.576241649;
  z[25] = 5.818232135;
  z[26] = 6.062278833;
  z[27] = 6.308544361;
  z[28] = 6.557207032;
  z[29] = 6.808463353;
  z[30] = 7.062531060;
  z[31] = 7.319652822;
  z[32] = 7.580100808;
  z[33] = 7.844182384;
  z[34] = 8.112247311;
  z[35] = 8.384696940;
  z[36] = 8.661996168;
  z[37] = 8.944689217;
  z[38] = 9.233420890;
  z[39] = 9.528965823;
  z[40] = 9.832269808;
  z[41] = 10.14450994;
  z[42] = 10.46718542;
  z[43] = 10.80226075;
  z[44] = 11.15240439;
  z[45] = 11.52141540;
  z[46] = 11.91506194;
  z[47] = 12.34296422;
  z[48] = 12.82379975;
  z[49] = 13.40648734;
  return z[i];
}

//================================================//

int main (){
  symbol x("x");
  int i;
  ex H = Herm(99,x);
  for (i=0; i<50; i++)
    cout << H.subs(x==Z(i)) << endl;
  return 0;
}