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

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

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

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

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

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

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

}