#include <stdio.h>
#include <math.h>
#include <complex.h>

#define N 100

int main(){
complex psi[N];
double dx=10./N;
double ax;
FILE *fp;
int i,k=4,px=N*.4;
  for(i=0;i<N;i++){
    ax=exp(-(i-px)*(i-px)*dx*dx*5.);
    psi[i]=ax*cexp(I*k*(i-px)*dx);
    printf("%lf %lf %lf\n",i*dx,creal(psi[i]),cimag(psi[i]));
    }

fp=fopen("sch.dat","w");
fwrite(psi,sizeof(complex),N,fp);
fclose(fp);
}