This is the code for calculating solid angle C, surface pressure ps, and field pressure pf coming out from a pulsating cylinder with radius of r and normal velocity vn in an unbounded two dimensional acoustic domain using the solution of Helmholtz Integral Equation.

For verifying the calculation result, the code will also compare the calculation result with the analytical solution.

In case of non-uniqueness problem, for example at r = 1 m and f = 131 Hz, up to 10 CHIEF points can be implemented.

The code is written in a very simple and traditional way so that any beginner in MATLAB programming would be able to understand it easily.