We have the following linear ordinary differential equation and we want to solve for
(x-2)y' - y = 2(x-2)^3
Can you help me find where the error is? I am getting incorrect results.
First, express the equation in the standard form of an exact differential equation:
y + 2(x-2)^3 + (-x+2)y' = 0 where
M = y + 2(x-2)^3 and
N = -x+2
We check if the equation is exact. It isn't, so we try to find an integrating factor
We find our
u(x) and multiply our equation with it:
uy + u2(x-2)^3 + u(-x+2)y' = 0
u(x) = 1/(-x+2)^2 and we get:
1/(-x+2)^2 * y + 1/(-x+2)^2 * 2(x-2)^3 + 1/(-x+2)^2 * (-x+2)y' = 0
Clean it up and we get:
1/(-x+2)^2 * y + 2x - 4 + 1/(-x+2) * y' = 0
Let's denote our potential function with F (in some materials they use the phi symbol).
dF/dx = 1/(-x+2)^2 * y + 2x - 4
In order to find
F, let's integrate both sides and add
g(y), which represents the information that we are unable to retrieve from a partial derivative of x. Let's absorb the constant in it.
F = (-x+2)^(-1) * y + x^2 - 4x + g(y)
Let's find out what
dF/dy = (-x+2)^(-1) + g'(y) = uN
uN from earlier:
dF/dy = (-x+2)^(-1) + g'(y) = (-x+2)^(-1)
And we get
g'(y) = 0 so
g(y) = C
F = (-x+2)^(-1) * y + x^2 - 4x + C = 0
y = x^3 - 6x^2 + 8x + Cx - 2C
This is incorrect. Where did I make a mistake?
P.S. I realize there are easier choices for
u available, but let's not go there. I really want to find out what is wrong with this solution. Before posting this I verified that my choice for
u(x) doesn't contain variables other than x and I confirmed that the equation with
u(x) is exact by comparing partials derivatives. I also tried to look for blunders by verifying just about everything with Wolfram Alpha. I can't find the error.