00001 #include <math.h>
00002 #include <iostream.h>
00003 #include "ReactorFortran.hh"
00004
00005 double U235flux(double enu){
00006
00007
00008
00009 int dim = 15;
00010 float etab[15]={1.0,1.5,2.0,2.5,3.0,3.5,4.0,4.5,5.0,5.5,6.0,6.5,7.0,7.5,8};
00011 float ftab[15]={2.44,1.73,1.28,0.860,0.580,0.410,0.282,0.184,0.118,0.0720,
00012 0.0408,0.0195,0.0104,0.00506,0.00204};
00013
00014 double sum;
00015 for(int n=0;n<dim;n++)sum+=ftab[n];
00016
00017 float ptab[dim];
00018 for(int n=0;n<dim;n++)ptab[n]=ftab[n]/sum;
00019
00020 int kpoly = 2;
00021 float e=enu;
00022
00023 return divdif_(ptab,etab,dim,e,kpoly);
00024
00025 }
00026