Can you solve this second-order differential equation with given conditions?

  • Thread starter Mervin
  • Start date
In summary, the given equation can be solved through a transformation and then transformed into a second order differential equation with constant coefficients. The final solution is a combination of two functions, with arbitrary constants, which can be obtained using a specific value for K. This solution is also consistent with another solution provided by using a different argument. However, more information is needed on the origin and use of this equation.
  • #1
Hi everyone,

Does anyone know this equation and how to solve it:

(x-a)^2 * (x-b)^2 * y'' = c*y

Physics news on
  • #2
I think the solution is

y = L*(x - A)^((a+K)/(a-b)) . (x - B)^(-(b+K)/(a-b))

Where L, K are arbitrary constants.

Please check (which may be not exactly a doddle).
Last edited:
  • #3
Whoa there, epenguin. Please don't give out answers. This site exists in part to help people learn. Telling someone an answer straight out is harmful, not helpful. There is no learning that way.
  • #4
Thank you very much for the solution. I will check it my self. Don't worry.

Thanks once again!
  • #5
D H said:
Whoa there, epenguin. Please don't give out answers. This site exists in part to help people learn. Telling someone an answer straight out is harmful, not helpful. There is no learning that way.

1. I didn't tell him how I got it. (I don't know that this comes into a standard variety of problems. Mervin asked 'does anyone know...' and in a month no one seemed to.) Perhaps it is and it would be good hear from someone, perhaps Mervin, about what the significant generalisations are, and whether there are significant applications.

2. For actually checking it he will have to do an efficiently organised calculation. It took me some time to get right, well it seemed right last night - I would not bet my house on it, so presumably that is a not quite trivial excercise that gives some insights. Mervin?

3. In a sense the only way to solve a differential equation or integration is to know the answer, or at least enough of its shape to be able to hammer it till it fits. I admit to subversive questionings about whether therefore it is a valid thing to teach or give so much emphasis on it and integrations, given that it will be used only by real professionals, and for understanding physics you only need those answers, which you can check. :uhh:
However I will try to limit my subversion at this site, according to what you say. :smile:

4. OK Mervyn, if yo are presented with a d.e. like that, should we say in general f(x).y'' = y , how would you go about solving it, with this one as example?
  • #6
Hello Mervin,

The solution provided by epenguin seems not to be the right one. The factor c is not present, which is an indication. However it can be solved by a transformation. If you use the following:

[tex]y=(x-b)\cdot \eta[/tex]

You can thransform the given equation into one with constant coefficients, which is easily solved. Let's do this step by step. First use the transformation for y, you get as intermediate results:


The equation becomes:


From the other part of the transformation we have:

  • #7
part two of the solution, I get a database error...

[tex]\frac{d\eta}{dx}=\frac{a-b}{(x-a)(x-b)} \frac{d\eta}{d\xi}[/tex]
[tex]\frac{d^2\eta}{dx^2}=-\frac{(a-b)(2x-(a+b))}{(x-a)^2(x-b)^2}\frac{d\eta}{d\xi} +\frac{(a-b)^2}{(x-a)^2(x-b)^2} \frac{d^2\eta}{d\xi^2}[/tex]

The equation becomes now:


This equation is a second order differential equation with constant coefficients which can be solved in a standard way. The final solution is:

[tex]y(x)=A|x-a|^{\frac{1}{2}(1+\lambda)}|x-b|^{\frac{1}{2}(1-\lambda)} + B|x-a|^{\frac{1}{2}(1-\lambda)}|x-b|^{\frac{1}{2}(1+\lambda)}[/tex]



After all these weeks I finally found it in one of my recently obtained books. The book "Handbook of Exact Solutions for Ordinary Differential Equations" written by Polyanin has this one described. The method is given, but not how the transformation was obtained. However I suspect that it is closely related to Lie point symmetries, which is a theory (new for me) I am starting to learn. I can't use it at this equation yet, but maybe in the future I will be able to show that it is the correct way. Can you give some more information on where this equation is used or coming from?
  • #8
Well I got the hard part right!


I had another look at this equation

(x-a)^2 * (x-b)^2 * y'' = c*y. (1)

For the solution, obtained without a book :tongue2: by a different argument than coomast’s, I gave

y = L*(x - a)^((a+K)/(a-b)) . (x - b)^(-(b+K)/(a-b)) (2)

It is true of this, as can be checked by twice differentiating, that

(x-a)^2 * (x-b)^2 * y'' = C*y (3)

where you can get any C you want by choosing the K. So, late at night I thought therefore 'I've done it'. Not quite right or complete. (It was subliminally disquieting that my formulation after that was so unlike the linear constant coefficient case.)

The way to think of it rather is, if you twice differentiate w.r.t. x, the function that appears in eq. 2 above without the constant, i.e. define

g(x, K) = (x - a)^((a+K)/(a-b)) . (x - b)^(-(b+K)/(a-b)) (4)

you find

(x-a)^2 * (x-b)^2 *g’’(x, K) = (a + K)*(b + K)*g(x, K) (5)

(More explicitly

(x-a)^2 * (x-b)^2 * [(x - a)^((a+K)/(a-b)) . (x - b)^(-(b+K)/(a-b))]’’

= (a + K)*(b + K)*[(x - a)^((a+K)/(a-b)) . (x - b)^(-(b+K)/(a-b))] ) (6)

So this function g is a solution of

(x-a)^2 * (x-b)^2 * y'' = C*y

if (a + K)*(b + K) = C (7)

There are then two values of K, K1, K2 say, that satisfy (7). So

y = g(x, K1) and y = g(x, K2) are solutions and so is

y = A*g(x, K1) + B*g(x, K2) (8)

where A, B are arbitrary constants in the usual way remembered from lde’s.

When we incorporate the K’s which are solutions of (7), the formula is exactly the same as that given by coomast.

As this has turned out not all that easy maybe it is OK to give my way of finding the solution and rewrite in LaTex presently?
Last edited:
  • #9
Hello. I thought after all this time and as a solution, I think identical to mine has been posted, I would indicate how I got it. The idea is basically that the form eq. (2) below is a plausible one natural to try, even if with a bit more insight (or hindsight) another form leads you there too.

I have had trouble with latex posting (database error and I have to trick the system by building up a bit at a time) so I now have to do this in two posts. Hopefully, cough cough, here goes:


[tex] (x - a)^2(x-b)^2y^\prime^\prime = cy\ (1)[/tex]

assume a solution of form

[tex]y = e^p[/tex]. (2)

:rolleyes:(to be continued...)
  • #10

[tex]y^\prime^\prime = (e^p)^\prime^\prime[/tex]

[tex] \ = (p^\prime^\prime + p\prime^2)e^p[/tex]

Combining with (1) we obtain

[tex] p^\prime^\prime + p^\prime^2 = \frac{c}{(x - a)^2(x-b)^2}[/tex] (3)

In order to get such a demominator it seems necessary or probable we need a denominator (x - a)(x-b) in [tex]p^\prime[/tex]. So let

[tex]p^\prime = \frac{f(x)}{(x - a)(x-b)}[/tex] (4)


[tex]p^\prime^\prime + p^\prime^2 = \frac{(x-a)(x-b)f^\prime - [(x-a) + (x-b)]f + f^2}{(x - a)^2(x-b)^2}[/tex] (5)

Assume f is a polynomial of degree n. Then the first two numerator terms are of degree (n + 1) whilst the last is of degree 2n. For a cancellation to leave only a constant as in eq.(3) to be possible, it is necessary that n=1. So let

[tex]f = lx + K[/tex] (6)

The numerator is then found to be

[tex]x^2(l^2-l) + 2xK(l-1) + (abl + aK + bK + K^2)[/tex]

so to give form (3)

[tex]l = 1 [/tex] and

[tex]ab + (a+b)K +K^2 = c[/tex] (7)

and (4) becomes

[tex]p^\prime = \frac{x+K}{(x - a)(x-b)} [/tex].


[tex]p = \int\frac{x+K}{(x - a)(x-b)}\ dx[/tex]

[tex]= \frac{(a+K)\ln(x-a) - (b+K)\ln (x-b)}{(a-b)} + K^\prime[/tex]

So [tex]y = e^p[/tex] works out as

[tex]y = (x-a)^\frac{a+K}{a-b} (x-b)^-^\frac{b+K}{a-b}[/tex]

This y multiplied by any constant is also a solution, and is so for either of the values of K which are solutions of eq. (7), so similarly to the way of solving lde's, and inserting the solutions for K, we find the general solution

[tex]y = A(x-a)^r(x-b)^s + B(x-a)^s(x-b)^r[/tex]

where [tex]r,s = \frac{1}{2}[1 \pm\surd(1+\frac{4c}{(a-b)})][/tex]

and A, B arbitrary constants.

1. What is the meaning of the equation (x-a)^2 * (x-b)^2 * y'' = c*y?

This equation represents a type of second-order linear differential equation, where the second derivative of the function y is equal to a constant multiple of the function itself. The constants a and b represent the roots of the quadratic polynomial (x-a)^2 * (x-b)^2, which determine the behavior of the solution to the equation.

2. How do you solve the equation (x-a)^2 * (x-b)^2 * y'' = c*y?

To solve this equation, you can use the method of undetermined coefficients, where you assume a solution of the form y = e^(rx) and solve for the constant r. This will give you two solutions, and the general solution can be expressed as a linear combination of these solutions. Alternatively, you can also use the method of variation of parameters to find a particular solution.

3. What are the applications of the equation (x-a)^2 * (x-b)^2 * y'' = c*y?

This type of differential equation can be used to model a variety of physical phenomena, such as the motion of a spring-mass system or the oscillations of a pendulum. It can also be used in engineering applications, such as in the design of circuits or control systems.

4. How does the value of the constant c affect the behavior of the solution to (x-a)^2 * (x-b)^2 * y'' = c*y?

The value of the constant c determines the type of solution to the equation. If c is positive, the solution will exhibit exponential growth or decay. If c is negative, the solution will oscillate between positive and negative values. If c is zero, the solution will be a constant function.

5. Can the equation (x-a)^2 * (x-b)^2 * y'' = c*y have complex solutions?

Yes, the equation can have complex solutions, especially when the value of c is negative. In this case, the general solution will involve complex numbers, and the real and imaginary parts of the solution will represent oscillatory behavior. This can be seen in the solutions to the harmonic oscillator equation, which is a special case of this equation.

Similar threads

  • Differential Equations
  • Differential Equations
  • Differential Equations
  • Differential Equations
  • Differential Equations