If, as suffices for all practical purposes, we limit the application of the formulae to points in advance of the plane at which the wave is supposed to be broken up, we may use simpler methods of resolution than that above considered. It appears indeed that the purely mathematical question has no definite answer. In illustration of this the analogous problem for sound may be referred to. Imagine a flexible lamina to be introduced so as to coincide with the plane at which resolution is to be effected. The introduction of the lamina (supposed to be devoid of inertia) will make no difference to the propagation of plane parallel sonorous waves through the position which it occupies. At every point the motion of the lamina will be the same as would have occurred in its absence, the pressure of the waves impinging from behind being just what is required to generate the waves in front. Now it is evident that the aerial motion in front of the lamina is determined by what happens at the lamina without regard to the cause of the motion there existing. Whether the necessary forces are due to aerial pressures acting on the rear, or to forces directly impressed from without, is a matter of indifference. The conception of the lamina leads immediately to two schemes, according to which a primary wave may be supposed to be broken up. In the first of these the element dS, the effect of which is to be estimated, is supposed to execute its actual motion, while every other element of the plane lamina is maintained at rest. The resulting aerial motion in front is readily calculated (see Rayleigh, Theory of Sound, § 278); it is symmetrical with respect to the origin, i.e. independent of θ. When the secondary disturbance thus obtained is integrated with respect to dS over the entire plane of the lamina, the result is necessarily the same as would have been obtained had the primary wave been supposed to pass on without resolution, for this is precisely the motion generated when every element of the lamina vibrates with a common motion, equal to that attributed to dS. The only assumption here involved is the evidently legitimate one that, when two systems of variously distributed motion at the lamina are superposed, the corresponding motions in front are superposed also.

The method of resolution just described is the simplest, but it is only one of an indefinite number that might be proposed, and which are all equally legitimate, so long as the question is regarded as a merely mathematical one, without reference to the physical properties of actual screens. If, instead of supposing the motion at dS to be that of the primary wave, and to be zero elsewhere, we suppose the force operative over the element dS of the lamina to be that corresponding to the primary wave, and to vanish elsewhere, we obtain a secondary wave following quite a different law. In this case the motion in different directions varies as cosθ, vanishing at right angles to the direction of propagation of the primary wave. Here again, on integration over the entire lamina, the aggregate effect of the secondary waves is necessarily the same as that of the primary.

In order to apply these ideas to the investigation of the secondary wave of light, we require the solution of a problem, first treated by Stokes, viz. the determination of the motion in an infinitely extended elastic solid due to a locally applied periodic force. If we suppose that the force impressed upon the element of mass D dx dy dz is

DZ dx dy dz,

being everywhere parallel to the axis of Z, the only change required in our equations (1), (2) is the addition of the term Z to the second member of the third equation (2). In the forced vibration, now under consideration, Z, and the quantities ξ, η, ζ, δ expressing the resulting motion, are to be supposed proportional to eint, where i = √(-1), and n = 2π/τ, τ being the periodic time. Under these circumstances the double differentiation with respect to t of any quantity is equivalent to multiplication by the factor -n², and thus our equations take the form

(b²Δ² + n²)ξ + (a² − b²) = 0 }  (7).
dx
(b²Δ² + n²)η + (a² − b²) = 0
dy
(b²Δ² + n²)ζ + (a² − b²) = −Z
dz

It will now be convenient to introduce the quantities.ω1, ω2, ω3 which express the rotations of the elements of the medium round axes parallel to those of co-ordinates, in accordance with the equations

ω3 = , ω1 = , ω2   (8).
dy dx′dz dy′dx dz′

In terms of these we obtain from (7), by differentiation and subtraction,

(b²Δ² + n²) ω3 = 0 }  (9).
(b²Δ² + n²) ω1 = dZ/dy
(b²Δ² + n²) ω2 = −dZ/dx