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

main (int argc, char *argv[]) {

  int i;
  double bin[4096], blah[4096], val1[4096], val2[4096], resid[4096];
  FILE *pre, *post, *out;

  pre = fopen("55352_411.giantpulse.profile","r");
  post = fopen("55532_437.giantpulse.profile","r");
  out = fopen("gp_diff.txt","w");

  for(i = 0; i < 4096; i++) {
    fscanf(pre,"%lf %lf\n",&bin[i],&val1[i]);
  }

  for(i = 0; i < 4096; i++) {
    fscanf(post,"%lf %lf\n",&blah[i],&val2[i]);
  }

  for(i = 0; i < 4096; i++) {
    resid[i] = val1[i]-val2[i];
    fprintf(out,"%lf %lf\n",bin[i],resid[i]);
  }

}