{ "id": "5c2e193ff592f822bf36d18d52a5936b4a6de042", "text": "IV. On the Equation of Laplace's Functions, &c. By W. F. Donkin, M.A., F.R.S., F.R.A.S., Savilian Professor of Astronomy in the University of Oxford.\n\nReceived December 3,—Read December 11, 1856.\n\nSection I.\n\nThe equation\n\n$$\\frac{d^2u}{dx^2} + \\frac{d^2u}{dy^2} + \\frac{d^2u}{dz^2} = 0,$$\n\ntransformed to polar coordinates by putting\n\n$$x = r \\sin \\theta \\cos \\phi, \\quad y = r \\sin \\theta \\sin \\phi, \\quad z = r \\cos \\theta,$$\n\nbecomes, as is well known,\n\n$$\\frac{d^2u}{d\\theta^2} + \\cot \\theta \\frac{du}{d\\theta} + \\frac{1}{(\\sin \\theta)^2} \\frac{d^2u}{d\\phi^2} + r^2 \\frac{d^2u}{dr^2} + 2r \\frac{du}{dr} = 0,$$\n\nwhich may be written in the form\n\n$$\\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + \\left( \\frac{d}{d\\phi} \\right)^2 + (\\sin \\theta)^2 r \\frac{d}{dr} \\left( r \\frac{d}{dr} + 1 \\right) u = 0; \\ldots \\ldots \\ldots (1.)$$\n\nand if it be assumed that\n\n$$u = u_0 + u_1 r + u_2 r^2 + \\ldots \\ldots,$$\n\nthen by substituting this value of $u$ in either of the two equations last written, we find that $u_n$ is a function of $\\theta$ and $\\phi$ satisfying the equation\n\n$$\\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + \\left( \\frac{d}{d\\phi} \\right)^2 + n(n+1)(\\sin \\theta)^2 u_n = 0, \\ldots \\ldots \\ldots \\ldots (2.)$$\n\nwhich, under a slightly different form, is commonly called the Equation of Laplace's Functions.\n\nThis equation was first solved in finite terms by Mr. Hargreave*, but in a form very inconvenient for practical applications. A solution free from this objection was afterwards obtained by Professor Boole†, by a method explained in his memoir \"On a General Method in Analysis,\" Philosophical Transactions, 1844. Lastly, in the second volume of the Journal referred to, the same mathematician gave two more solutions, one of which however is reducible, as he states, to Mr. Hargreave's form; the other, though much more convenient than this, is still for most purposes probably less useful than that given in the first volume, from which it differs essentially in form, as well as in the method by which it is deduced.\n\n* Philosophical Transactions, 1841.\n† Cambridge and Dublin Mathematical Journal, vol. i. p. 10.\nIn the following pages I have treated the equation (2.) by a very simple method*. The result bears a general resemblance to Professor Boole's first solution, and I conceive that the two forms must be capable of being identified by the assumption of a proper relation between the arbitrary functions; but I am not able at present to show this identity.\n\nSome further investigations are added, which it is unnecessary to notice beforehand.\n\n1. Putting $k$ for $\\frac{d}{d\\phi}$ in equation (2.), we have\n\n$$\\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + k^2 + n(n+1)(\\sin \\theta)^2 u_n = 0. \\quad \\ldots \\ldots \\ldots \\ldots \\quad (3.)$$\n\nIn the case of $n=0$, this becomes\n\n$$\\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + k^2 u_0 = 0, \\text{ or } \\frac{d^2 u_0}{dt^2} + k^2 u_0 = 0$$\n\n(where $t=\\log \\tan \\frac{\\theta}{2}$). Hence\n\n$$u_0 = C_1 \\left( \\tan \\frac{\\theta}{2} \\right)^{k\\sqrt{-1}} + C_2 \\left( \\tan \\frac{\\theta}{2} \\right)^{-k\\sqrt{-1}}. \\quad \\ldots \\ldots \\ldots \\ldots \\quad (4.)$$\n\nIf, however, $k=0$, we should have\n\n$$u_0 = C_1 + C_2 \\log \\tan \\frac{\\theta}{2}, \\quad \\ldots \\ldots \\ldots \\ldots \\quad (5.)$$\n\n* It may be worth while to give a preliminary illustration of this method (which contains no novelty except in detail), by applying a similar process to the well-known equation\n\n$$\\frac{d^2 u}{dx^2} + \\left( k^2 - \\frac{n(n+1)}{x^2} \\right) u = 0.$$\n\nThis may be written in the form\n\n$$\\left( D - \\frac{n}{x} \\right) \\left( D + \\frac{n}{x} \\right) + k^2 \\right) u = 0, \\quad \\left( D = \\frac{d}{dx} \\right);$$\n\nnow\n\n$$\\left( D + \\frac{n}{x} \\right) \\left( D - \\frac{n}{x} \\right) = \\left( D - \\frac{n-1}{x} \\right) \\left( D + \\frac{n-1}{x} \\right);$$\n\nhence if we put $u = \\left( D - \\frac{n}{x} \\right) v$, and then operate on each side with $\\left( D - \\frac{n}{x} \\right)^{-1}$, we have\n\n$$\\left( D - \\frac{n-1}{x} \\right) \\left( D + \\frac{n-1}{x} \\right) + k^2 \\right) v = \\left( D - \\frac{n}{x} \\right)^{-1} 0.$$\n\nAssuming 0 as the value of the right-hand member, and then putting $v = \\left( D - \\frac{n-1}{x} \\right) w$, and so on successively, it is obvious that we shall ultimately have\n\n$$u = \\left( D - \\frac{n}{x} \\right) \\left( D - \\frac{n-1}{x} \\right) \\cdots \\left( D - \\frac{1}{x} \\right) z,$$\n\n$$\\left( D^2 + k^2 \\right) z = 0;$$\n\nhence, observing that $D - \\frac{i}{x} = x^i D x^{-i}$, we obtain finally,\n\n$$u = x^n \\left( D - \\frac{1}{x} \\right)^n \\left( C_1 \\sin kx + C_2 \\cos kx \\right),$$\n\nwhich agrees with the result obtained by Professor Boole's general method (Philosophical Transactions, 1844).\nwhich is the complete solution of (2.) in the particular case in which \\( u_n \\) does not contain \\( \\varphi \\), and \\( n = 0 \\).\n\nBut putting \\( \\frac{d}{d\\varphi} \\) for \\( k \\) in (4.), and arbitrary functions of \\( \\varphi \\) instead of \\( C_1, C_2 \\), we have\n\n\\[\nu_0 = f\\left(\\varphi + \\sqrt{-1} \\log \\tan \\frac{\\theta}{2}\\right) + F\\left(\\varphi - \\sqrt{-1} \\log \\tan \\frac{\\theta}{2}\\right); \\quad \\ldots \\quad (6.)\n\\]\n\ninstead of which we may evidently write\n\n\\[\nu_0 = f\\left(e^{\\varphi \\sqrt{-1}} \\tan \\frac{\\theta}{2}\\right) + F\\left(e^{-\\varphi \\sqrt{-1}} \\tan \\frac{\\theta}{2}\\right), \\quad \\ldots \\quad (7.)\n\\]\n\nas well as other forms. Either of these is the complete value of \\( u_n \\) in the case of \\( n = 0 \\).\n\n2. Returning to the general case, if we put for shortness\n\n\\[\n(\\sin \\theta \\frac{d}{d\\theta})^2 + n(n+1)(\\sin \\theta)^2 = \\xi_n,\n\\]\n\nit is evident that\n\n\\[\n\\xi_{-n} = \\xi_{n-1}; \\quad \\ldots \\quad (8.)\n\\]\n\nand also, since\n\n\\[\n(\\sin \\theta \\frac{d}{d\\theta} + n \\cos \\theta)(\\sin \\theta \\frac{d}{d\\theta} - n \\cos \\theta) = (\\sin \\theta \\frac{d}{d\\theta})^2 + n(n+1)(\\sin \\theta)^2 - n^2,\n\\]\n\nif we put\n\n\\[\nw_n = \\sin \\theta \\frac{d}{d\\theta} + n \\cos \\theta,\n\\]\n\nwe have\n\n\\[\nw_n w_{-n} + n^2 = \\xi_n, \\quad \\ldots \\quad (9.)\n\\]\n\nfrom which, changing \\( n \\) into \\(-n\\), and observing (8.), we get\n\n\\[\nw_{-n} w_n + n^2 = \\xi_{n-1};\n\\]\n\non the other hand, changing \\( n \\) into \\( n-1 \\) in (9.), we have\n\n\\[\nw_{n-1} w_{-(n-1)} + (n-1)^2 = \\xi_{n-1}.\n\\]\n\nComparing the last two results, we see that the operation \\( w_n \\) possesses the following property, namely,\n\n\\[\nw_{-n} w_n + n^2 = w_{n-1} w_{-(n-1)} + (n-1)^2. \\quad \\ldots \\quad (10.)\n\\]\n\nNow the equation to be solved (3.), art. 1, is\n\n\\[\n(w_n w_{-n} + n^2 + k^2)u_n = 0. \\quad \\ldots \\quad (11.)\n\\]\n\nLet \\( u_n = w_n v \\), and let the operation \\( w_n^{-1} \\) be performed on each side; assuming for the present that \\( w_n^{-1} \\) may be put \\( = 0 \\) without ultimate loss of generality, we obtain\n\n\\[\n(w_{-n} w_n + n^2 + k^2)v = 0;\n\\]\n\nand this, by virtue of (10.), becomes\n\n\\[\n(w_{n-1} w_{-(n-1)} + (n-1)^2 + k^2)v = 0;\n\\]\n\nwhich, compared with (11.), gives the relation \\( v = u_{n-1} \\), and consequently, since \\( u_n = w_n v \\),\n\n\\[\nu_n = w_n u_{n-1},\n\\]\n\nor\n\n\\[\nu_n = (\\sin \\theta \\frac{d}{d\\theta} + n \\cos \\theta)u_{n-1}.\n\\]\nThis relation is true whatever be the value of \\( n \\). But as we are now supposing \\( n \\) a positive integer, it is evident that we have\n\n\\[\nu_n = \\omega_n \\omega_{n-1} \\ldots \\omega_2 \\omega_1 u_0, \\quad \\ldots \\ldots \\ldots \\ldots . (12.)\n\\]\n\nwhere \\( u_0 \\) is given by (6.) or (7.) of the last article in the general case, and by (5.) in the case in which \\( u_n \\) is a function of \\( \\theta \\) without \\( \\phi \\).\n\nThis solution contains two arbitrary functions in the former case, and two arbitrary constants in the latter. It is therefore, in general, the complete solution. (See, however, arts. 8, 9, 17, &c.)\n\n3. The operative symbol \\( \\omega_n \\omega_{n-1} \\ldots \\omega_2 \\omega_1 \\), written at length, is\n\n\\[\n(\\sin \\theta \\frac{d}{d\\theta} + n \\cos \\theta)(\\sin \\theta \\frac{d}{d\\theta} + (n-1) \\cos \\theta) \\ldots (\\sin \\theta \\frac{d}{d\\theta} + \\cos \\theta). \\quad \\ldots \\ldots . (13.)\n\\]\n\nNow it is evident that\n\n\\[\n\\sin \\theta \\frac{d}{d\\theta} + n \\cos \\theta = \\frac{1}{(\\sin \\theta)^{n-1}} \\frac{d}{d\\theta} (\\sin \\theta)^n\n\\]\n\n(the subject of operation being of course omitted on both sides); and if this be put in the form\n\n\\[\n\\frac{1}{(\\sin \\theta)^n} \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta (\\sin \\theta)^{n-1},\n\\]\n\nit will be easily seen that the expression (13.) is equivalent to\n\n\\[\n\\frac{1}{(\\sin \\theta)^n} (\\sin \\theta \\frac{d}{d\\theta} \\sin \\theta)^n, \\quad \\ldots \\ldots \\ldots \\ldots \\ldots . (14.)\n\\]\n\nwhich may also be put in the less symmetrical form\n\n\\[\n\\frac{1}{(\\sin \\theta)^{n-1}} \\left( (\\sin \\theta)^2 \\frac{d}{d\\theta} \\right)^n \\sin \\theta.\n\\]\n\nEither of these forms has the peculiarity, as compared with (13.), of being intelligible without supposing \\( n \\) a positive integer (setting aside the difficulties of general differentiation).\n\n4. The results of the preceding articles may be summed up as follows:—\n\nThe complete solution of the equation\n\n\\[\n\\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + n(n+1)(\\sin \\theta)^2 u_n = 0 \\text{ is}\n\\]\n\n\\[\nu_n = (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\left( C_1 + C_2 \\log \\tan \\frac{\\theta}{2} \\right); \\quad \\ldots \\ldots . (15.)\n\\]\n\nand the complete solution of the equation\n\n\\[\n\\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + \\left( \\frac{d}{d\\phi} \\right)^2 + n(n+1)(\\sin \\theta)^2 u_n = 0 \\text{ is}\n\\]\n\n\\[\nu_n = (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\left[ f(e^{\\sqrt{-1} \\tan \\frac{\\theta}{2}}) + F(e^{-\\sqrt{-1} \\tan \\frac{\\theta}{2}}) \\right]. \\quad \\ldots . (16.)\n\\]\n\nIt may be observed that (15.) is deducible from (16.) by taking \\( f(x) = C_1 + \\frac{1}{2} C_2 \\log x, \\)\n\n\\( F(x) = \\frac{1}{2} C_2 \\log x. \\)\n5. The expression (16.) may be compared with the following result of Professor Boole's first solution*:\n\n\\[ u_n = F\\left(\\mu, \\frac{\\sqrt{(1-\\mu^2)}}{\\mu} e^{\\theta \\sqrt{-1}}\\right), \\]\n\nwhere \\( F(\\mu, e^{\\theta \\sqrt{-1}}) \\)\n\n\\[ = \\mu^n \\left( \\frac{d}{d\\mu} \\frac{1}{\\mu} \\right)^n \\left\\{ (\\mu + \\mu^2)^n \\psi \\left( \\frac{\\mu e^{\\theta \\sqrt{-1}}}{1+\\mu} \\right) + (\\mu - \\mu^2)^n \\chi \\left( \\frac{\\mu e^{\\theta \\sqrt{-1}}}{1-\\mu} \\right) \\right\\} \\]\n\nand \\( \\mu = \\cos \\theta \\). I have stated above my impression that this form must be capable of transformation into identity with (16.); but it is probable that this can only be effected by means of some not obvious theorem of differentiation. Examples of such theorems will be seen below (arts. 14 and 16).\n\nSection II.\n\n6. I proceed to discuss some particular applications of the formulæ obtained above.\n\nIf the expression \\((1-2r \\cos \\theta + r^2)^{-\\frac{1}{2}}\\) be developed in the form\n\n\\[ p_0 + p_1 r + p_2 r^2 + \\ldots, \\]\n\nthen \\( p_n \\) satisfies the equation of Laplace's functions, and does not contain \\( \\varphi \\). Hence \\( p_n \\) must be given by the expression (15.), art. 4, if the constants be properly assumed; and since it is evident that \\( C_2 \\) must be \\(=0\\) in this case, we have\n\n\\[ p_n = c_n (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n .1, \\]\n\nwhere \\( c_n \\) is a function of \\( n \\).\n\nThis form of \\( p_n \\) may be independently verified, and \\( c_n \\) determined at the same time, as follows:\n\nPut \\( r(\\sin \\theta)^{-1} = \\varepsilon \\), then\n\n\\[ (1-2r \\cos \\theta + r^2)^{-\\frac{1}{2}} = \\frac{1}{\\sin \\theta} \\left( 1 + (\\cot \\theta - \\varepsilon)^2 \\right)^{-\\frac{1}{2}}, \\]\n\nand if we put \\( \\cot \\theta = t \\), this may be written\n\n\\[ \\sqrt{1+t^2} e^{-t \\frac{d}{dt}} \\frac{1}{\\sqrt{1+t^2}}, \\]\n\nwhere \\( \\varepsilon \\) is to be treated as constant till after all operations. Now in general, if \\( P, Q \\) be any symbols whatever, quantitative or operative, we know that \\( Pf(Q) \\frac{1}{P} = f(PQ \\frac{1}{P}) \\).\n\nHence it is evident that since \\( \\sqrt{1+t^2} = \\frac{1}{\\sin \\theta} \\) and \\( -\\frac{d}{dt} = (\\sin \\theta)^2 \\frac{d}{d\\theta} \\), the above expression is equivalent to\n\n\\[ e^{\\varphi \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta} .1, \\]\n\nand therefore this is a symbolical form of \\((1-2r \\cos \\theta + r^2)^{-\\frac{1}{2}}\\). Restoring the value of \\( \\varepsilon \\) after development, we find that the coefficient of \\( r^n \\) is\n\n* Cambridge and Dublin Journal, vol. i. p. 18.\n\\[\n\\frac{1}{1 \\cdot 2 \\cdot 3 \\ldots n} (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\cdot 1.\n\\]\n\nIn other words, if we put \\( v_n = (\\sin \\theta)^{-n} (\\sin \\theta \\frac{d}{d\\theta} \\sin \\theta)^n \\cdot 1 \\), then\n\n\\[\n(1 - 2r \\cos \\theta + r^2)^{-\\frac{i}{2}} = v_0 + v_1 \\frac{r}{1} + v_2 \\frac{r^2}{1 \\cdot 2} + \\ldots + v_n \\frac{r^n}{1 \\cdot 2 \\ldots n} + \\ldots \\ldots .\n\\]\n\n(From this it is easy to deduce the relation \\( v_n = (2n-1) \\cos \\theta \\cdot v_{n-1} - (n-1)^2 v_{n-2} \\) by means of the known relation between the coefficients in the above series and those in the development of \\( (1 - r^2)(1 - 2r \\cos \\theta + r^2)^{-\\frac{3}{2}} \\).)\n\n7. Adopting the notation of the last article, we have\n\n\\[\ne^{\\rho \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta} (\\sin \\theta)^i = \\frac{1}{\\sin \\theta} e^{-\\rho \\frac{d}{d\\theta} (\\sin \\theta)^{i+1}}\n\\]\n\n\\[\n= \\frac{1}{\\sin \\theta} \\left( 1 + (t - \\varphi)^2 \\right)^{-\\frac{i+1}{2}}\n\\]\n\n\\[\n= (\\sin \\theta)(1 - 2r \\cos \\theta + r^2)^{-\\frac{i+1}{2}}.\n\\]\n\nHence this theorem, which will be useful afterwards. The coefficient of \\( \\frac{r^n}{1 \\cdot 2 \\ldots n} \\) in the development of \\( (1 - 2r \\cos \\theta + r^2)^{-\\frac{i+1}{2}} \\) is\n\n\\[\n(\\sin \\theta)^{-(n+i)} (\\sin \\theta \\frac{d}{d\\theta} \\sin \\theta)^n (\\sin \\theta)^i. \\ldots \\ldots \\ldots . \\quad (17.)\n\\]\n\n8. Let us next consider the development of the expression\n\n\\[\n\\{1 - 2r (\\cos \\theta \\cos \\varphi + \\sin \\theta \\sin \\varphi \\cos \\varphi) + r^2\\}^{-\\frac{i}{2}}. \\ldots \\ldots \\ldots . \\quad (18.)\n\\]\n\nSince the coefficient of \\( r^n \\) is necessarily expressible in the form\n\n\\[\nq_0 + q_1 \\cos \\varphi + q_2 \\cos 2\\varphi + \\ldots + q_n \\cos n\\varphi,\n\\]\n\nand must also satisfy the equation (2.), art. 1, we find, on substitution in that equation, that \\( q_i \\) must be a solution of the equation\n\n\\[\n\\left\\{ \\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + n(n+1)(\\sin \\theta)^2 - i^2 \\right\\} u_n = 0. \\ldots \\ldots \\ldots . \\quad (19.)\n\\]\n\nI proceed to consider this equation in a general manner, without reference, in the first instance, to the problem immediately in hand.\n\nSince the equation (19.) only differs from (3.), art. 1, in having \\(-i^2\\) instead of \\(k^2\\), we may apply the process used in the solution of that equation. We thus get\n\n\\[\nu_n = (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\left( C_1 \\left( \\tan \\frac{\\theta}{2} \\right)^i + C_2 \\left( \\cot \\frac{\\theta}{2} \\right)^i \\right). \\ldots \\ldots \\quad (20.)\n\\]\n\nNow, although this expression appears to contain two arbitrary constants, it will be shown presently, that when \\( i \\) is an integer and not greater than \\( n \\), it really contains only one. This is seen immediately to be the case when \\( i = 0 \\), and may be easily verified in other simple cases, such as \\( i = 1, n = 2, \\&c. \\)\n9. In fact, the process of art. 2 gives\n\n\\[ u_n = \\varpi_n u_{n-1} = \\varpi_n \\varpi_{n-1} u_{n-2} = \\ldots = \\varpi_n \\varpi_{n-1} \\cdots \\varpi_{n-m+1} u_{n-m}, \\]\n\nwhere \\( u_{n-m} \\) satisfies the equation\n\n\\[ (\\varpi_{n-m} \\varpi_{-(n-m)} + (n-m)^2 - r^2) u_{n-m} = 0; \\]\n\nand in general we can only arrive at an integrable form by taking \\( m = n \\). But when \\( i \\) is an integer not greater than \\( n \\), we shall get an integrable form by taking \\( m \\) so as to satisfy the condition \\((n-m)^2 = i^2\\). Suppose \\( i \\) positive and less than \\( n \\), and take \\( m = n - i \\); we then have for the solution of (19.), \n\n\\[ u_n = \\varpi_n \\varpi_{n-1} \\cdots \\varpi_{i+1} u_i, \\]\n\nwhere \\( u_i \\) satisfies the equation \\( \\varpi_i \\varpi_{-i} u_i = 0 \\), or\n\n\\[ (\\sin \\theta \\frac{d}{d\\theta} + i \\cos \\theta)(\\sin \\theta \\frac{d}{d\\theta} - i \\cos \\theta) u_i = 0 : \\]\n\nthe integration of this is easy, and the result is\n\n\\[ u_i = C_1 (\\sin \\theta)^i + C_2 (\\sin \\theta)^i \\int \\frac{a\\theta}{(\\sin \\theta)^{2i+1}}; \\]\n\nso that, in the case considered, the integral of (19.) may be expressed in the form\n\n\\[ u_n = (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^{n-i} (\\sin \\theta)^{2i} \\left( C_1 + C_2 \\int d\\theta (\\sin \\theta)^{-2i-1} \\right) . . . (21.) \\]\n\n(for \\( \\varpi_n \\varpi_{n-1} \\cdots \\varpi_{i+1} \\) is equivalent to \\( (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^{n-i} (\\sin \\theta)^{i} \\) (see art. 3)).\n\nIt is evident that the expression (21.) contains two independent arbitrary constants, whatever be the value of \\( i \\). And \\( i \\) being supposed a positive integer, the second term will always involve a logarithmic function; hence the form (20.), which involves no such function, cannot be the complete solution of the equation (19.), and therefore we shall not really limit its generality by putting \\( C_2 = 0 \\). Moreover, since the coefficient of \\( r^n \\cos i \\varphi \\) in the development of (18.) cannot contain any logarithmic function, it follows that, for the purpose of that development, we must put \\( C_2 = 0 \\) in the expression (21.).\n\nComparing the two forms thus obtained for the coefficient of \\( r^n \\cos i \\varphi \\), we obtain incidentally the following theorem, namely,\n\n\\[ \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^i \\left( \\tan \\frac{\\theta}{2} \\right)^i = c (\\sin \\theta)^{2i}, \\ldots \\ldots \\ldots \\ldots \\ldots (22.) \\]\n\n\\( c \\) being a constant, of which the value will appear afterwards.\n\nThe preceding investigation shows that the assumption \\( \\varpi_n^{-1} 0 = 0 \\) (instead of the general value \\( C (\\sin \\theta)^{-n} \\)) is liable in certain cases to limit the generality of the result. I shall return to this point presently, but proceed now to complete the development of the expression (18.).\n\n10. Since the coefficient of \\( r^n \\cos i \\varphi \\) must contain \\( \\theta \\) and \\( \\theta' \\) symmetrically, and satisfy both the equation (19.) and a similar equation in \\( \\theta' \\), it is evident that it must be of the form \\( f(\\theta)f(\\theta') \\). Hence, by means of the conclusions of art. 9, we arrive at the following\nresult: put, for convenience, \\( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta = \\Theta, \\sin \\theta' \\frac{d}{d\\theta'} \\sin \\theta' = \\Theta' \\); then the coefficient of \\( r^n \\) in the development of \\( (1 - 2r(\\cos \\theta \\cos \\theta' + \\sin \\theta \\sin \\theta' \\cos \\phi) + r^2)^{-\\frac{1}{2}} \\) may be expressed in either of the two forms,\n\n\\[\n\\frac{1}{(\\sin \\theta \\sin \\theta')^n} \\Theta^n \\Theta'^n \\left\\{ a_0 + a_1 \\tan \\frac{\\theta}{2} \\tan \\frac{\\theta'}{2} \\cos \\phi + a_2 \\left( \\tan \\frac{\\theta}{2} \\tan \\frac{\\theta'}{2} \\right)^2 \\cos 2\\phi + \\ldots + a_n \\left( \\tan \\frac{\\theta}{2} \\tan \\frac{\\theta'}{2} \\right)^n \\cos n\\phi \\right\\}\n\\]\n\n\\[\n+ \\frac{1}{(\\sin \\theta \\sin \\theta')^n} \\left\\{ c_0 \\Theta^n \\Theta'^n + c_1 \\Theta^{n-1} \\Theta'^{n-1} (\\sin \\theta \\sin \\theta')^3 \\cos \\phi + \\ldots + c_n (\\sin \\theta \\sin \\theta')^{2n} \\cos n\\phi \\right\\}\n\\]\n\nand it only remains to determine the constants.\n\n11. For this purpose we may adopt a process of which the first part is taken from Laplace. Let \\( \\cos \\theta = \\mu, \\cos \\theta' = \\mu' \\), and suppose the coefficient of \\( r^n \\cos i\\phi \\) in the development of \\( (1 - 2r(\\mu \\mu' + \\sqrt{1 - \\mu^2} \\sqrt{1 - \\mu'^2} \\cos \\phi) + r^2)^{-\\frac{1}{2}} \\) to be itself developed in a series of powers and products of \\( \\mu \\) and \\( \\mu' \\). It is only necessary to ascertain the term independent of \\( \\mu, \\mu' \\), and the term containing the product \\( \\mu \\mu' \\); and for this end we may write the above expression thus, \\( (1 - 2r(\\mu \\mu' + \\cos \\phi) + r^2)^{-\\frac{1}{2}} \\), since the terms \\( \\mu^2, \\mu'^2 \\) cannot affect the required result; and this again may (for the present purpose) be considered equal to \\( (1 - 2r \\cos \\phi + r^2)^{-\\frac{1}{2}} + \\mu \\mu' r (1 - 2r \\cos \\phi + r^2)^{-\\frac{3}{2}} \\); and if we put \\( 2 \\cos \\phi = x + \\frac{1}{x} \\), so that \\( 1 - 2r \\cos \\phi + r^2 = (1 - rx)(1 + \\frac{r}{x}) \\), and develope in the usual manner, we easily find the following results:\n\n1st, if \\( n-i \\) be even, the term independent of \\( \\mu, \\mu' \\) in the coefficient of \\( r^n \\cos i\\phi \\) is\n\n\\[\n2 \\cdot \\frac{1^2 \\cdot 3^2 \\cdot 5^2 \\cdots (n+i-1)^2 \\cdot 1^2 \\cdot 3^2 \\cdot 5^2 \\cdots (n-i-1)^2}{1 \\cdot 2 \\cdot 3 \\cdots (n+i) \\cdot 1 \\cdot 2 \\cdot 3 \\cdots (n-i)}\n\\]\n\nand the term containing \\( \\mu \\mu' \\) vanishes.\n\n2nd, if \\( n-i \\) be odd, the term independent of \\( \\mu, \\mu' \\) vanishes, and the term containing \\( \\mu \\mu' \\) is\n\n\\[\n2 \\cdot \\frac{1^2 \\cdot 3^2 \\cdots (n+i)^2 \\cdot 1^2 \\cdot 3^2 \\cdots (n-i)^2}{1 \\cdot 2 \\cdot 3 \\cdots (n+i) \\cdot 1 \\cdot 2 \\cdot 3 \\cdots (n-i)} \\mu \\mu',\n\\]\n\nthe factor 2 to be omitted in each case when \\( i = 0 \\).\n\n12. On the other hand, the coefficient of \\( r^n \\cos i\\phi \\) is (art. 10)\n\n\\[\nc_i (\\sin \\theta \\sin \\theta')^{-i} \\Theta^{n-i} \\Theta'^{n-i} (\\sin \\theta)^{2i} (\\sin \\theta')^{2i};\n\\]\n\nand we may find the term independent of \\( \\mu, \\mu' \\), and the term involving \\( \\mu \\mu' \\) in the development of this, as follows:\n\nIt is easily found, from the theorem of art. 7, that\n\n\\[\n(\\sin \\theta)^{-i} \\Theta^{n-i} (\\sin \\theta)^{2i} = 1 \\cdot 2 \\cdots (n-i) (\\sin \\theta)^n \\times (\\text{coefficient of } r^{n-i} \\text{ in the development of } (1 - 2r \\cos \\theta + r^2)^{-(i+\\frac{1}{2}))};\n\\]\n\nand for our present purpose we may put \\( \\sin \\theta = 1 \\), and\n\n\\[\n(1 - 2r \\cos \\theta + r^2)^{-(i+\\frac{1}{2})} = (1 + r^2)^{-(i+\\frac{1}{2})} + (2i+1)\\mu \\mu r (1 + r^2)^{-(i+\\frac{3}{2})}.\n\\]\nHence it is easily found that when \\( n-i \\) is even, the term independent of \\( \\mu, \\mu' \\) in (25.) is\n\n\\[\nc_i \\cdot 1^2 \\cdot 3^2 \\cdot 5^2 \\cdots (n-i-1)^2 \\cdot (2i+1)^2 \\cdot (2i+3)^2 \\cdots (n+i-1)^2;\n\\]\n\nand when \\( n-i \\) is odd, the term involving \\( \\mu \\mu' \\) is\n\n\\[\nc_i \\cdot \\mu \\mu' \\cdot 1^2 \\cdot 3^2 \\cdot 5^2 \\cdots (n-i)^2 \\cdot (2i+1)^2 \\cdot (2i+3)^2 \\cdots (n+i)^2.\n\\]\n\nIf these expressions be equated to those obtained in art. 11, we find that whether \\( n-i \\) be even or odd the value of \\( c_i \\) is\n\n\\[\n2 \\cdot \\frac{1^2 \\cdot 3^2 \\cdot 5^2 \\cdots (2i-1)^2}{1 \\cdot 2 \\cdot 3 \\cdots (n-i) \\cdot 1 \\cdot 2 \\cdot 3 \\cdots (n+i)},\n\\]\n\nthe factor 2 being omitted when \\( i=0 \\).\n\n13. We have then the following result:\n\nThe coefficient of \\( r^n \\cos i \\phi \\) in the development of\n\n\\[\n(1-2r(\\cos \\theta \\cos \\theta' + \\sin \\theta \\sin \\theta' \\cos \\phi)+r^2)^{-\\frac{1}{2}} \\text{ is }\n\\]\n\n\\[\n2 \\cdot \\frac{1^2 \\cdot 3^2 \\cdot 5^2 \\cdots (2i-1)^2}{1 \\cdot 2 \\cdot 3 \\cdots (n-i) \\cdot 1 \\cdot 2 \\cdot 3 \\cdots (n+i)} (\\sin \\theta \\sin \\theta')^{-n} \\Theta^{n-i} \\Theta'^{n-i} (\\sin \\theta)^{2i} (\\sin \\theta')^{2i},\n\\]\n\nwhere \\( \\Theta = \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta, \\Theta' = \\sin \\theta' \\frac{d}{d\\theta'} \\sin \\theta' \\); and with respect to the numerical coefficient, it is to be observed that when \\( i=0 \\) the factor 2 is to be omitted, and the numerator considered to become unity. The extreme values of \\( i \\) are of course 0 and \\( n \\), and the term \\( 1 \\cdot 2 \\cdot 3 \\cdots (n-i) \\) is to be taken =1 when \\( i=n \\).\n\n14. It is now easy to ascertain the values of the coefficients \\( a_0, a_1, \\&c. \\) in the form (23.), art. 10. But first it is necessary to prove that\n\n\\[\n(\\sin \\theta \\frac{d}{d\\theta} \\sin \\theta)^n (\\tan \\frac{\\theta}{2})^n = 1 \\cdot 3 \\cdot 5 \\cdots (2n-1)(\\sin \\theta)^n.\n\\]\n\nWe have proved this already (art. 9, equation (22.)), except in so far as the numerical coefficient was left undetermined, so that we have only to establish its value. Now the expression on the left of the above equation may be written\n\n\\[\n(\\sin \\theta)^{-1} \\left( (\\sin \\theta)^2 \\frac{d}{d\\theta} \\right)^n \\sin \\theta \\left( \\tan \\frac{\\theta}{2} \\right)^n;\n\\]\n\nand if we put \\( \\cot \\theta=t \\), so that \\( \\tan \\frac{\\theta}{2} = \\sqrt{1+t^2}-t \\), and \\( (\\sin \\theta)^2 \\frac{d}{d\\theta} = -\\frac{d}{dt} \\), the equation (27.) becomes\n\n\\[\n(-)^n \\left( \\frac{d}{dt} \\right)^n \\left( \\frac{\\sqrt{1+t^2}-t}{\\sqrt{1+t^2}} \\right)^n = \\frac{1 \\cdot 3 \\cdot 5 \\cdots (2n-1)}{(1+t^2)^{\\frac{2n+1}{2}}};\n\\]\n\nand since it is only the value of the coefficient that is in question, it will be enough to prove this when \\( t=0 \\), i.e. to show that the coefficient of \\( t^n \\) in the development of\n\n\\[\n\\frac{(t-\\sqrt{1+t^2})^n}{\\sqrt{1+t^2}}\n\\]\n\nis \\( \\frac{1 \\cdot 3 \\cdot 5 \\cdots (2n-1)}{1 \\cdot 2 \\cdot 3 \\cdots n} \\). Now, expanding the numerator by the binomial theorem, the above\nexpression becomes\n\n\\[\n\\frac{t^n}{\\sqrt{1+t^2}} - nt^{n-1} + \\frac{n(n-1)}{1.2} t^{n-2} \\sqrt{1+t^2} - \\frac{n(n-1)(n-2)}{1.2.3} t^{n-3}(1+t^2) + \\ldots + (-)^n(1+t^2)^{\\frac{n-1}{2}};\n\\]\n\nand if we develope each term and collect the coefficient of \\( t^n \\), we find that it is\n\n\\[\n1 + \\frac{n(n-1)}{2^2} + \\frac{n(n-1)(n-2)(n-3)}{2^2.4^2} + \\ldots + \\frac{n(n-1)\\ldots(n-2i+1)}{2^2.4^2.6^2\\ldots(2i)^2} + \\ldots,\n\\]\n\na series of which the sum * is \\( \\frac{1.3.5\\ldots(2n-1)}{1.2.3\\ldots n} \\). The truth of (27.) is therefore established.\n\n15. It follows that\n\n\\[\n\\Theta^n \\left( \\tan \\frac{\\theta}{2} \\right)^i = \\Theta^{n-i} \\Theta^i \\left( \\tan \\frac{\\theta}{2} \\right)^i = 1.3\\ldots(2i-1)\\Theta^{n-i}(\\sin \\theta)^{2i},\n\\]\n\nand consequently the coefficient of \\( r^n \\cos i\\varphi \\) may also be expressed in the form\n\n\\[\n\\frac{2}{1.2.3\\ldots(n-i).1.2.3\\ldots(n+i)} (\\sin \\theta \\sin \\theta')^{-n}\\Theta^n \\Theta'^n \\left( \\tan \\frac{\\theta}{2} \\right)^i \\left( \\tan \\frac{\\theta'}{2} \\right)^i,\n\\]\n\nthe factor 2 being omitted as before when \\( i=0 \\).\n\nThe law of the numerical coefficients in this case is remarkable. For if the expression \\( (1+\\cos \\varphi)^n \\) be developed in the form \\( A_0 + A_1 \\cos \\varphi + A_2 \\cos 2\\varphi + \\ldots + A_n \\cos n\\varphi \\), it will be found that\n\n\\[\nA_i = 1.2.3\\ldots n.1.3.5\\ldots(2n-1) \\cdot \\frac{2}{1.2.3\\ldots(n-i).1.2.3\\ldots(n+i)} \\quad \\text{(omitting the 2 when } i=0),\n\\]\n\nso that we may express the law of the development of\n\n\\[\n(1-2r(\\cos \\theta \\cos \\theta' + \\sin \\theta \\sin \\theta' \\cos \\varphi) + r^2)^{-}\n\\]\n\nas follows:—the coefficient of \\( r^n \\) is\n\n\\[\n\\frac{1}{1.2.3\\ldots n.1.3.5\\ldots(2n-1)} (\\sin \\theta)^{-n}(\\sin \\theta')^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\left( \\sin \\theta' \\frac{d}{d\\theta'} \\sin \\theta' \\right)^n U_n,\n\\]\n\nwhere\n\n\\[\nU_n = A_0 + A_1 \\cos \\varphi \\tan \\frac{\\theta}{2} \\tan \\frac{\\theta'}{2} + A_2 \\cos 2\\varphi \\left( \\tan \\frac{\\theta}{2} \\right)^2 \\left( \\tan \\frac{\\theta'}{2} \\right)^2 + \\ldots + A_n \\cos n\\varphi \\left( \\tan \\frac{\\theta}{2} \\right)^n \\left( \\tan \\frac{\\theta'}{2} \\right)^n,\n\\]\n\nand \\( A_i \\) is the coefficient of \\( \\cos i\\varphi \\) in the development of \\( (1+\\cos \\varphi)^n \\).\n\nAlthough this form of the development follows a more simple and remarkable law than that of art. 13, it is evident that for actual calculation it would be much less convenient, since the number of operations is much greater, and the differentiations more complex. But as the complete explicit form of the development is known independently, this consideration is not of much practical importance.\n\n* \\( \\frac{1.3\\ldots(2n-1)}{1.2\\ldots n} \\) is the coefficient of \\( x^n \\) in the development of \\( (1-2x)^{-\\frac{1}{2}} \\); now\n\n\\[\n(1-2x)^{-\\frac{1}{2}} = ((1-x)^2-x^2)^{-\\frac{1}{2}} = (1-x)^{-1} + \\frac{1}{2}(1-x)^{-3}x^2 + \\frac{1.3}{2.4}(1-x)^{-5}x^4 + \\ldots;\n\\]\n\nand if each term be developed and the coefficient of \\( x^n \\) collected, the result is the series in the text.\n\nThe equation (27.) probably admits of some simpler demonstration, which I have failed to perceive.\n16. Professor Boole, in the paper above referred to, obtains the following expression for the coefficient of \\( r^n \\cos i\\varphi \\); namely,\n\n\\[\n\\frac{2f(\\mu)f(\\mu')}{1.2....(n-i).1.2...(n+i)},\n\\]\n\nwhere\n\n\\[\nf(\\mu) = (1-\\mu^2)^i \\mu^{n-i} \\left( \\frac{d}{d\\mu} \\frac{1}{\\mu} \\right)^n \\mu^{n+i}(1+\\mu)^{n-i}.\n\\] \n\n(29.)\n\nHere the numerical coefficients are the same as in the form just given. It follows, therefore, that the expression on the right of (29.) must be equivalent to\n\n\\[\n(\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\left( \\tan \\frac{\\theta}{2} \\right)^i.\n\\]\n\nThis equivalence, and that expressed by equation (27.), art. 14, are instances of theorems by no means obvious or easy to verify directly. (See art. 5.)\n\nFor actual calculation, Professor Boole's form would be preferable to either of mine; for though the number of operations is much greater than in the form (26.), their result may be assigned with much greater facility. But considered merely as an analytical expression, the form (26.) is the simplest of the three.\n\nSection III.\n\n17. It was shown above (art. 9), that the expression\n\n\\[\nu_n = (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\left( C_1 \\left( \\tan \\frac{\\theta}{2} \\right)^i + C_2 \\left( \\tan \\frac{\\theta}{2} \\right)^{-i} \\right)\n\\]\n\nfails in certain cases to represent the complete solution of the equation\n\n\\[\n\\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + n(n+1)(\\sin \\theta)^2 - i^2 \\right) u_n = 0,\n\\]\n\nnamely, when \\( i \\) is an integer, and not greater than \\( n \\). (Similarly, the expression \\( u = C_1 e^{ix} + C_2 e^{-ix} \\) fails to give the complete solution of \\( \\left( \\frac{d}{dx} \\right)^2 - i^2 \\) \\( u = 0 \\), when \\( i = 0 \\).)\n\nIt is desirable, therefore, to investigate the solution of the equation in a more general manner. Putting, as before, \\( \\sin \\theta \\frac{d}{d\\theta} + n \\cos \\theta = \\omega_n \\), the equation to be integrated is\n\n\\[\n(\\omega_n \\omega_{-n} + n^2 - i^2) u_n = 0.\n\\] \n\n(30.)\n\nLet \\( u_n = \\omega_n v \\), then\n\n\\[\n(\\omega_n \\omega_{-n} \\omega_n + (n^2 - i^2) \\omega_n) v = 0,\n\\]\n\nand any value of \\( v \\) which satisfies this, will give a value of \\( u_n \\) satisfying (30.). Now \\( \\omega_n \\omega_{-n} + n^2 \\) is identically the same as \\( \\omega_{n-1} \\omega_{-(n-1)} + (n-1)^2 \\) (10.), art. 2; hence the equation becomes\n\n\\[\n\\omega_n (\\omega_{n-1} \\omega_{-(n-1)} + (n-1)^2 - i^2) v = 0;\n\\]\n\nsimilarly, putting \\( v = \\omega_{n-1} \\omega \\), we get\n\n\\[\n\\omega_n \\omega_{n-1} (\\omega_{n-2} \\omega_{-(n-2)} + (n-2)^2 - i^2) \\omega = 0,\n\\]\nand so on, till we have, finally,\n\n\\[ u_n = \\omega_n \\omega_{n-1} \\cdots \\omega_1 z, \\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots (31.) \\]\n\nand\n\n\\[ \\omega_n \\omega_{n-1} \\cdots \\omega_1 (\\omega_0^2 - i^2) z = 0. \\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots (32.) \\]\n\nThe point to be observed is, that every value of \\( z \\) which satisfies (32.) will give a value of \\( u_n \\) satisfying (30.); hence, although the complete value of \\( z \\) will contain \\( n+2 \\) constants, it is certain that \\( n \\) of them will be destroyed by the \\( n \\) direct operations of the expression (31.). In general, the two constants left are those introduced by the inverse operation \\((\\omega_0^2 - i^2)^{-1}\\); but in the exceptional cases noticed above, one of these disappears, and one of those introduced by the other inverse operations remains instead.\n\n18. The complete value of \\( u_n \\), in its most general form, is therefore\n\n\\[ u_n = \\omega_n \\omega_{n-1} \\cdots \\omega_1 (\\omega_0^2 - i^2)^{-1} \\omega_1^{-1} \\omega_2^{-1} \\cdots \\omega_n^{-1}.0. \\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots (33.) \\]\n\nIt was shown before that \\( \\omega_n \\omega_{n-1} \\cdots \\omega_1 = (\\sin \\theta)^{-n} \\left( \\frac{\\sin \\theta}{d \\theta} \\right)^n \\); in like manner we have \\( \\omega_1^{-1} \\omega_2^{-1} \\cdots \\omega_n^{-1} = \\left( \\frac{1}{\\sin \\theta} \\left( \\frac{d}{d \\theta} \\right)^{-1} \\frac{1}{\\sin \\theta} \\right)^n (\\sin \\theta)^n \\), the latter operation being the inverse of the former.\n\nThe expression \\((\\omega_0^2 - i^2)^{-1}\\) is equivalent to \\(\\left( \\frac{d}{dt} \\right)^2 - i^2 \\right)^{-1}\\), if \\( t = \\log \\tan \\frac{\\theta}{2} \\); and this may be resolved into either of the two following forms; namely,\n\n\\[ \\frac{1}{2} \\left[ \\left( \\frac{d}{dt} - i \\right)^{-1} - \\left( \\frac{d}{dt} + i \\right)^{-1} \\right], \\]\n\nor\n\n\\[ \\left( \\frac{d}{dt} \\pm i \\right)^{-1} \\left( \\frac{d}{dt} \\mp i \\right)^{-1}. \\]\n\nThe former, omitting the useless factor \\( \\frac{1}{2} \\), gives\n\n\\[ \\left( \\tan \\frac{\\theta}{2} \\right)^i \\left( \\frac{d}{d \\theta} \\right)^{-1} \\frac{1}{\\sin \\theta} \\left( \\cot \\frac{\\theta}{2} \\right)^i - \\left( \\cot \\frac{\\theta}{2} \\right)^i \\left( \\frac{d}{d \\theta} \\right)^{-1} \\frac{1}{\\sin \\theta} \\left( \\tan \\frac{\\theta}{2} \\right)^i; \\]\n\nthe latter, taking the lower signs, gives\n\n\\[ \\left( \\tan \\frac{\\theta}{2} \\right)^i \\left( \\frac{d}{d \\theta} \\right)^{-1} \\frac{1}{\\sin \\theta} \\left( \\cot \\frac{\\theta}{2} \\right)^{2i} \\left( \\frac{d}{d \\theta} \\right)^{-1} \\frac{1}{\\sin \\theta} \\left( \\tan \\frac{\\theta}{2} \\right)^i. \\]\n\nEither of these expressions is to be substituted in (33.), and thus the complete value of \\( u_n \\) is obtained by means of \\( n+2 \\) integrations, and \\( n \\) differentiations. In general the \\( n \\) integrations implied in the expression \\(\\left( \\frac{1}{\\sin \\theta} \\left( \\frac{d}{d \\theta} \\right)^{-1} \\frac{1}{\\sin \\theta} \\right)^n 0\\) are useless, so that we may substitute 0 for this expression, as was done in the process first given, without ultimate loss of completeness.\n\n19. If we suppose \\( i \\) an integer, and less than \\( n \\), we may stop the process of art. 17 at an earlier stage, thus\n\n\\[ u_n = \\omega_n \\omega_{n-1} \\cdots \\omega_{i+1} z, \\]\n\nand\n\n\\[ \\omega_n \\omega_{n-1} \\cdots \\omega_i \\omega_{i+1} z = 0, \\]\n\nso that we now have\n\n\\[ u_n = \\omega_n \\omega_{n-1} \\cdots \\omega_{i+1} \\omega_{i+1}^{-1} \\omega_{i+1}^{-1} \\cdots \\omega_n^{-1}.0; \\]\nor, reducing as before, and observing that \\( \\varphi^{-i} \\) is equivalent to\n\n\\[\n(\\sin \\theta)^i \\left( \\frac{d}{d\\theta} \\right)^{-1} (\\sin \\theta)^{-i-1},\n\\]\n\nwe have, finally,\n\n\\[\nu_n = (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^{-i} (\\sin \\theta)^{2i} \\left( \\frac{d}{d\\theta} \\right)^{-1} (\\sin \\theta)^{-2i} \\left( \\frac{1}{\\sin \\theta} \\left( \\frac{d}{d\\theta} \\right)^{-1} \\frac{1}{\\sin \\theta} \\right)^{n-i+1}. . . (34.)\n\\]\n\nIf in this expression we put 0 for the result of all but one of the operations indicated in the last term, we reproduce the form (21.) obtained in art. 9, which appears in fact to be always complete. The form (34.), however, as before remarked, will never contain superfluous constants.\n\n20. It is an interesting question whether the forms thus obtained on the supposition that \\( i \\) is an integer, are not really general, inasmuch as they are expressible without reference to that supposition. I believe that they are; but I doubt whether, in the existing state of analysis, it can be proved either that they are or are not (see § IV.).\n\n21. It is easy to obtain other forms analogous to those given in the preceding articles, containing inverse in place of direct operations. But such forms are objectionable, because they necessarily give rise to superfluous constants, the discrimination of which may be a matter of difficulty.\n\nSection IV.\n\n22. It was observed at the beginning of this paper that the equation\n\n\\[\n\\frac{d^2u}{dx^2} + \\frac{d^2u}{dy^2} + \\frac{d^2u}{dz^2} = 0,\n\\]\n\nwhen transformed to polar coordinates, may be written in the form\n\n\\[\n\\left\\{ \\left( \\sin \\theta \\frac{d}{d\\theta} \\right)^2 + \\left( \\frac{d}{d\\phi} \\right)^2 + (\\sin \\theta)^2 r \\frac{d}{dr} \\left( r \\frac{d}{dr} + 1 \\right) \\right\\} u = 0. . . . . (35.)\n\\]\n\nI do not know whether this form has been noticed; but it is remarkable in this respect, that if the symbol \\( r \\frac{d}{dr} \\) be replaced by \\( n \\), it becomes identical with the equation discussed in the first section, and of which the solution was shown to be\n\n\\[\n(\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\left\\{ f \\left( e^{\\varphi \\sqrt{-1}} \\tan \\frac{\\theta}{2} \\right) + F \\left( e^{-\\varphi \\sqrt{-1}} \\tan \\frac{\\theta}{2} \\right) \\right\\}. . . . . (36.)\n\\]\n\nIt is not easy to verify, by direct substitution, that this expression satisfies the differential equation of which it is the solution; at least I have not yet succeeded in doing so; and even if this were accomplished, it is most likely that difficulties would remain to be overcome before we could establish the legitimacy of extending this result to the case in which \\( n \\) is a symbol of unlimited signification. Still it seems worth while to try, at least as an experiment, the consequence of assuming that the form (36.) will give the solution of (35.) on putting \\( r \\frac{d}{dr} \\) instead of \\( n \\).\nOmitting the second arbitrary function to save space, we should have (observing that \\( r \\) must be considered to enter the arbitrary functions in an arbitrary manner)\n\n\\[\nu = (\\sin \\theta)^{-r} \\frac{d}{d\\theta} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^r f \\left( r, e^{\\sqrt{-1}} \\tan \\frac{\\theta}{2} \\right).\n\\] \n\n(37.)\n\n23. There is no difficulty in the interpretation of this expression. It is only necessary to recollect that \\( e^x \\frac{d}{dx} \\varphi(x) = \\varphi(cx) \\), and we have the following result:—let \\( \\frac{r}{\\sin \\theta} = \\varepsilon \\), then\n\n\\[\nu = f \\left( \\varepsilon \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta, e^{\\sqrt{-1}} \\tan \\frac{\\theta}{2} \\right),\n\\]\n\nwhere \\( \\varepsilon \\) is to be put \\( = \\frac{r}{\\sin \\theta} \\) after all other operations, and the function must be supposed to be developed as if the two symbols \\( \\varepsilon \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta, e^{\\sqrt{-1}} \\tan \\frac{\\theta}{2} \\) were commutative, the powers of the former being always prefixed to those of the latter.\n\nIt is easy to derive from the form (38.), particular expressions which we already know to be solutions of (35.). For example, supplying the other arbitrary function, it is evident that if we take \\( f(x, y) = \\frac{1}{2} x^n y^i \\), and \\( F(x, y) \\) the same, we get\n\n\\[\nu = r^n (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\left( \\tan \\frac{\\theta}{2} \\right)^i \\cos i \\varphi,\n\\]\n\nwhich, we know, satisfies (35.).\n\nAgain, we may take \\( f(x, y) = e^x \\), or\n\n\\[\nu = e^{\\varepsilon \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta} \\cdot 1,\n\\]\n\nwhich has been already shown (art. 6) to represent\n\n\\[\n(1 - 2r \\cos \\theta + r^2)^{-\\frac{1}{2}}.\n\\]\n\nThese verifications afford, I think, ground to believe that the form (37.) is a true solution of the equation (35.). And here I leave the subject, hoping that it may attract the attention of some mathematician more able than myself to bring it to a satisfactory conclusion.\n\nNote on the last paragraph.\n\nReceived December 6, 1856.\n\nIf the expression \\( f \\left( \\varepsilon \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta, e^{\\sqrt{-1}} \\tan \\frac{\\theta}{2} \\right) \\) be supposed to be developed and arranged according to powers of the term \\( \\varepsilon \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\), each term in the result will be of the form\n\n\\[\nr^n (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\psi_n \\left( e^{\\sqrt{-1}} \\tan \\frac{\\theta}{2} \\right),\n\\]\n$\\psi_n$ representing a function of which the form may be considered to depend in an arbitrary manner upon the value of $n$. Now we know that every such term as the above satisfies the differential equation, and it may be asked, what more is required to prove the correctness of the solution? I answer, that the point in question is whether\n\n$$\\sum r^n (\\sin \\theta)^{-n} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^n \\psi_n \\left( e^{\\theta \\sqrt{-1}} \\tan \\frac{\\theta}{2} \\right) . . . . . . (A.)$$\n\nis equivalent to the most general interpretation of\n\n$$(\\sin \\theta)^{-r} \\frac{d}{dr} \\left( \\sin \\theta \\frac{d}{d\\theta} \\sin \\theta \\right)^r f(r, e^{\\theta \\sqrt{-1}} \\tan \\frac{\\theta}{2});$$\n\nand, if not, whether the most general interpretation would satisfy the Differential Equation (35.).\n\nThe expression (A.) gives, in fact, no new result, being an evident consequence of the conclusions of the former sections.\n\nPostscript (Added May 29, 1857).\n\nMr. Cayley has been kind enough to communicate to me direct verifications of the equation (27.), art. 14, and of the identity referred to in art. 16. Assuming a formula established in Mr. Cayley's paper \"On certain Formulae for Differentiation,\" &c.*, the former of the two theorems just mentioned is easily obtained, the latter not without a good deal of trouble. The investigations of Section II. may be compared with some of those in the Introduction to Murphy's 'Electricity;' but I have not at present examined the latter with a view to such comparison, and therefore merely mention the subject.\n\n* Cambridge and Dublin Journal, vol. ii. p. 124, equation (2.).", "source": "olmocr", "added": "2026-01-12", "created": "2026-01-12", "metadata": { "Source-File": "/home/jic823/projects/def-jic823/royalsociety/pdfs/108610.pdf", "olmocr-version": "0.3.4", "pdf-total-pages": 16, "total-input-tokens": 26652, "total-output-tokens": 15125, "total-fallback-pages": 0 }, "attributes": { "pdf_page_numbers": [ [ 0, 0, 1 ], [ 0, 2142, 2 ], [ 2142, 4698, 3 ], [ 4698, 7023, 4 ], [ 7023, 9725, 5 ], [ 9725, 12258, 6 ], [ 12258, 15120, 7 ], [ 15120, 18295, 8 ], [ 18295, 21766, 9 ], [ 21766, 24710, 10 ], [ 24710, 27833, 11 ], [ 27833, 30436, 12 ], [ 30436, 33883, 13 ], [ 33883, 36821, 14 ], [ 36821, 39347, 15 ], [ 39347, 40986, 16 ] ], "primary_language": [ "en", "en", "en", "en", "en", "en", "en", "en", "en", "en", "en", "en", "en", "en", "en", "en" ], "is_rotation_valid": [ true, true, true, true, true, true, true, true, true, true, true, true, true, true, true, true ], "rotation_correction": [ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ], "is_table": [ false, false, false, false, false, false, false, false, false, false, false, false, false, false, false, false ], "is_diagram": [ false, false, false, false, false, false, false, false, false, false, false, false, false, false, false, false ] }, "jstor_metadata": { "identifier": "jstor-108610", "title": "On the Equation of Laplace's Functions, &c", "authors": "W. F. Donkin", "year": 1857, "volume": "147", "journal": "Philosophical Transactions of the Royal Society of London", "page_count": 16, "jstor_url": "https://www.jstor.org/stable/108610" } }