Quote:
Originally Posted by ceavou22
1. The problem statement, all variables and given/known data
Use a Green's function to solve:
u" + 2u' + u = e-x
with u(0) = 0 and u(1) = 1 on 0<=x<=1
(Relevant information edited out. I really hate the way this is explained in that textbook anyways )
3. The attempt at a solution
Solving for the homogeneous equation first:
u" + 2u' + u = 0
From the characteristic equation,
lambda2 + 2lambda + 1 = 0
lambda= -1 (repeated root)
Characteristic solution:
u1(x) = c1e-x and u2 = c2e-xx
To satisfy boundary conditions,
u(0) = c1e0 + c2e0(0)
c1 = 0
Thus take f(x) = e-xx
and
u(1) = c1e-1 + c2e-1(1)
get c1 = e/2 and c2 =e/2
Thus take g(x) = (1/2)e-x+1(x+1)
evaluating the wroknskian, W = (1/2)e-2x+1
I contruct green's function as per the formula provided above and carry out the integral and get final answer of,
u(x) = xe-x(x2 - 3)
which obviously doesn't satisfy the second boundary condition of u(1) = 1.
I found out the solution from Wolfram alpha:
u = (1/2)e-xx(x+2e-1)
I've also tried many different combinations of f(x) and g(x) but none seen to work
Any help appreciated !
Ok i think there is a "conceptual" error here, and some calculation errors.
Lets fix first the former since its more important.
This is how you do an inhomogenous diff equation in general:
Get a solution for the inhomogeneous equation with some method (for example Green's method). Then get a complete system of solutions for the homogeneous equation. Then you add the right combination of the homogenous solutions to the inhom. solution to satisfy conditions.
So for this exercise:
v := (1/2)*x
2*e
-x is a solution for the inhom. equation (which you should have gotten using the formula i think).
(What you wrote is not a solution for the inhom. equation).
As you found yourself a complete system of solutions for the hom. equation is:
u
1(x) = e
-x
and u
2 = e
-xx
(compl. system here means thatthese span the vspace of hom solutions)
Now you are looking for the solution in the form of:
v + c
1u
1+c
2u
2
Subsitituting this into the boundary conditions provides you with 2 linear equations for c
1 and c
2, which I am sure you can solve.
Now why did you not get this inhom solution v using Green's method?
Well due to calculations errors i suppose. For example your g does not vanish in 1, like it should.
There might be more errors, you will have to check yourself.
Let me know if something is not clear. Good luck with finishing the problem.