I need help with the code of the following questions please using JupyterHub notebook

Suppose a plane wave of wavelength , such as light or a sound wave, is blocked by an object with a straight edge, represented by the solid line at the bottom of this figure: The wave will be diffracted at the edge and the resulting intensity at the position (x,z) marked by the dot is given by near-field diffraction theory to be I=8I0([2C(u)+1]2+[2S(u)+1]2) where I0 is the intensity of the wave before diffraction and u=xz2,C(u)=0ucos21t2dt,S(u)=0usin21t2dt. Write a program to calculate I0I and make a plot of it as a function of x in the range 5m to 5m for the case of a sound wave with wavelength =1m, measured z=3m past the straight edge. Calculate the integrals using Gaussian quadrature with N=50 points. You should find significant variation in the intensity of the diffracted sound--enough that you could easily hear the effect if sound were diffracted, say, at the edge of a tall buliding. [20 points] psuedocode