int main(void){
printf("2*2 = %d\n",quadrat(2));
return 0;
}
const double Experiment::Intensity(double alpha,const double lambda){
const double k = 2 * M_PI / lambda;
// Normalized: I(alpha) / I_0
// return pow(sin(k/2 * width * sin(alpha)/ (k/2 * width * sin(alpha)) ),2) * pow(cos(k/2*distance*sin(alpha)),2) ; // (k/2 * width * sin(alpha))^2; //* cos(k/2 *alpha *sin(alpha))^2;
return quadrat(sinc(k/2. * slit.Width() *sin(alpha))) * quadrat(cos(k/2.*slit.Distance()*sin(alpha)));
}