For , the Bernoulli differential equation becomes linear under . Multiplying the equation by gives
On an interval with , the integrating factor is , hence
Therefore the complete set of solution branches is
The original equation requires ; each nonzero branch additionally requires . The zero solution is restored explicitly because the reciprocal substitution excludes it.
On a branch with , the given derivative is nonzero wherever the original equation is defined, so use inverse-function linearization of a first-order equation and regard as a function of :
The integrating factor gives
Thus the nonzero solutions are given implicitly by
Each branch must also retain ; at those excluded points the original right-hand side is undefined and the inverse graph cannot represent a regular . Differentiating the implicit equation gives , equal to the required quotient wherever these restrictions hold.
The additional solution , on any interval avoiding , was lost when dividing by and must be included. It is not defined through the point in the original equation.
For the first-order linear ordinary differential equation, an integrating factor is
Multiplication converts the left side to a derivative:
Integrating gives . The initial value fixes , so
Rewrite the complex position as . Differentiating gives the complex velocity and acceleration
For real , . Their requested magnitudes are therefore
There is no pole for real time. The magnitude is the Euclidean length of the two Cartesian components, not the modulus of either component separately.
Real coefficients imply that the conjugate is also a polynomial root. Write the remaining root as . Then
The coefficient of in this product is . That coefficient is zero in the specified polynomial, so . Hence the two other roots are
As a check, the full polynomial is , corresponding to , and .
The original PDF asks for cube roots of . Their complex modulus is , and their arguments are , . Converting to Cartesian form gives
Cubing any one gives the positive real number .
The square roots simplify to the absolute values and . Put . A nontrivial solution of the homogeneous linear system exists exactly when its coefficient matrix is singular:
The second factor vanishes precisely when , or , giving both coordinate axes. The first factor vanishes when is a nonnegative integer. Hence the complete locus is
These are the two axes together with concentric circles of positive integer radii. The zero-radius case is the origin, already included on the axes; at that point both equations vanish and any nonzero pair works.
Figure 1.
The parameter locus consists of both axes and circles with integer radii
. Details.
Collect the coefficients of the three unknowns in rows. The linear system has matrix form
The next two parts solve this same system by the two requested methods.
Subtract the lower curve from the upper curve on the intersection interval. The definite integral is
The result is in square coordinate units.
The parabola opens upward, with its minimum at and zeros at . The quartic graph is symmetric about the vertical axis, has its flat maximum at , and also has zeros at .
Their intersections solve , so the only real intersections are and . Between them, is the upper curve and the lower curve.
Figure 1.
The parabola, the inverted quartic, and the finite region between their intersections
. Details.
The defining equation gives . Equality is attainable at , with the positive branch giving
Insert the linear initial profile into the transformed integral from (b):
The constant part integrates to . The other part is an absolutely integrable odd function, so its integral vanishes. Hence
As a direct check, and . Thus a linear temperature profile is stationary on the whole line. More generally, every affine function is preserved by this heat-kernel solution.
The indicated integration by substitution gives
Pointwise, the argument of tends to . Since the Gaussian integral is , the expected limit is . To justify it for unbounded data, use the approximate identity property rather than assuming a uniform bound on .
In the original variable, normalization gives
For , continuity makes the difference smaller than any prescribed if is small enough, and the kernel mass is at most one. On the complement the Gaussian mass tends to zero; the data-weighted tail tends to zero by heat-kernel convergence for Gaussian-integrable initial data. This proves
For heat-kernel convergence from a conditionally convergent Gaussian integral, the same tail conclusion follows using the bounded primitive of . The factor has its maximum at , approaching . For sufficiently small , it is monotone on each of the two tails outside , and its values at the tail endpoints tend to zero exponentially. Integration by parts bounds each tail by a constant times its endpoint value, proving the same limit without an absolute-convergence assumption. The approach is from positive time.
Write the kernel as , with . Its partial derivatives are
Thus the heat kernel itself satisfies . The proposed convolution consequently gives
This verifies the heat equation. The parameter is positive time; the real integral formula is not defined at negative time.
The passage of derivatives through the integral needs justification. A sufficient condition is the Gaussian-weighted absolute integrability in heat-kernel convergence for Gaussian-integrable initial data; it holds, in particular, for data of polynomial growth. On compact positive-time intervals strictly below an integrable later time, the differentiated kernels have an integrable Gaussian bound, so differentiation under the integral sign applies.
Even if the PDF's assumed convergence is read as ordinary conditional convergence, one can justify the same calculation. Choose a later time and let be a primitive of . Convergence at spatial point zero and time makes bounded on both tails. Put . This function, its spatial and time derivatives, and their -derivatives decay as Gaussians when . Integration by parts rewrites the convolution as , with vanishing endpoint terms. Boundedness of now permits differentiation under this absolutely convergent integral; since , the same equation follows. Thus one need not silently equate conditional and absolute convergence.
For ,
In particular its absolute value is bounded by , proving convergence when ; the apparent endpoint singularity at one is removable, with value . On any compact parameter interval with lower bound , the differentiated integrand is , dominated by . The dominated convergence theorem therefore justifies differentiation under the integral sign:
Hence . Since the integrand is identically zero at , fixes . Thus
This proves the logarithmic integral representation of the natural logarithm, including its parameter range.
The Leibniz integral rule accounts for both the moving upper endpoint and the integrand's dependence on :
For the specified integrand this gives . When , integration by parts gives
so
For the requested independent check, direct integration gives at . Its derivative is , agreeing with the result above. At zero the apparent singularity is removable: , so .

There are unlisted articles, also show them or only show them.