msci Code

#include "ndisc.h"

#define DTOR PI/180

int main(int argc, char *argv[]) {
  Ndisc *system;
  FILE *fint;		//file for intermediate results
  float b;		//initial conditions
  char noise_type;	//noise type, 'n'=none, 'g'=gauss., 'w'=white
  char std;		//noise standard deviation
  float t;		//time delay
  long n;		//number of deflections

  if (argc < 2) {
    printf("scatter b\n");
    printf("\n");
    printf("where:\n");
    printf("	b     =	impact parameter\n");
    return 1;
  }

  system=new Ndisc(0);

  sscanf(argv[1], "%lf", &b);

  system->scatter_b_corr(b, &t, &n);

  printf("%f %d\n", t, n);

  delete system;

  return 0;

}