A distribution analyzing device for analyzing a distribution of a field having a property satisfying a Laplace equation, the distribution analyzing device comprising: an obtaining unit configured to obtain measurement data of the field measured, through a sensor sensing area, independently at each of rotation angles of the sensor sensing area and at each of grid coordinate positions of the sensor sensing area on a measurement plane where the field is measured, the sensor sensing area being (i) a finite area which relatively rotates about an axis perpendicular to the measurement plane, (ii) an area which relatively moves such that a predetermined position in the area is brought to each of the grid coordinate positions in a state where the area has relatively rotated, and (iii) an area in which the field is sensed as an aggregate; and a calculation unit configured to calculate the distribution of the field from the measurement data of the field, using an arithmetic expression obtained by deriving a target harmonic function, which indicates the distribution of the field, using a condition that a convolution of the target harmonic function and a shape function, which indicates a shape of a cross section of the finite sensor sensing area along a plane parallel to the measurement plane, is equal to a provisional harmonic function derived by solving the Laplace equation using the measurement data of the field and a size of the sensor sensing area in a direction perpendicular to the measurement plane, characterized in that the arithmetic expression used by the calculation unit to calculate the distribution of the field from the measurement data of the field includes division of an expression obtained by integrating, over the rotation angles, an expression including multiplication of the shape function and the provisional harmonic function, by an expression obtained by integrating, over the rotation angles, an expression including a square of an absolute value of the shape function, the shape function and the provisional harmonic function each having been subjected to a Fourier transform. an obtaining unit configured to obtain measurement data of the field measured, through a sensor sensing area, independently at each of rotation angles of the sensor sensing area and at each of grid coordinate positions of the sensor sensing area on a measurement plane where the field is measured, the sensor sensing area being (i) a finite area which relatively rotates about an axis perpendicular to the measurement plane, (ii) an area which relatively moves such that a predetermined position in the area is brought to each of the grid coordinate positions in a state where the area has relatively rotated, and (iii) an area in which the field is sensed as an aggregate; and a calculation unit configured to calculate the distribution of the field from the measurement data of the field, using an arithmetic expression obtained by deriving a target harmonic function, which indicates the distribution of the field, using a condition that a convolution of the target harmonic function and a shape function, which indicates a shape of a cross section of the finite sensor sensing area along a plane parallel to the measurement plane, is equal to a provisional harmonic function derived by solving the Laplace equation using the measurement data of the field and a size of the sensor sensing area in a direction perpendicular to the measurement plane, characterized in that the arithmetic expression used by the calculation unit to calculate the distribution of the field from the measurement data of the field includes division of an expression obtained by integrating, over the rotation angles, an expression including multiplication of the shape function and the provisional harmonic function, by an expression obtained by integrating, over the rotation angles, an expression including a square of an absolute value of the shape function, the shape function and the provisional harmonic function each having been subjected to a Fourier transform.