Ejemplo n.º 1
0
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)));
}