The Wronskian 
is defined by
![\[ W(x) = u_1(x)u_2'(x) - u_2(x)u_1'(x) \]](https://www.stumblingrobot.com/wp-content/ql-cache/quicklatex.com-5573e53737d709f43b735c483be5fdaa_l3.png "Rendered by QuickLaTeX.com")
for given functions 
and 
.
Let 
and 
be two solutions of the second-order linear differential equation
![\[ y'' + ay' + by = 0, \]](https://www.stumblingrobot.com/wp-content/ql-cache/quicklatex.com-bc7dc287cfc50d544e76ec281bdcc813_l3.png "Rendered by QuickLaTeX.com")
such that 
is not constant.
- Let
be any solution of the given differential equation. Use the properties of the Wronskian proved in the previous two exercises (here and here) to prove that there exist constants 
and 
such that
![\[ c_1 v_1(0) + c_2 v_2(0) = f(0), \qquad c_1 v_1'(0)+ c_2 v_2'(0) = f'(0). \]](https://www.stumblingrobot.com/wp-content/ql-cache/quicklatex.com-cd82886d97de74a86d8a06bdeb2ba282_l3.png "Rendered by QuickLaTeX.com")
- Prove that every solution of the differential equation has the form
.
- Proof. By the previous exercise, since is not constant we know
. So we may define,
![\[ c_1 = \frac{f(0) v_2'(0) - v_2(0) f'(0)}{v_1(0)v_2'(0) - v_1'(0) v_2(0)}, \qquad c_2 = \frac{-f(0)v_2'(0) + v_1(0) f'(0)}{v_1(0)v_2'(0) - v_2(0) v_1'(0)}. \]](https://www.stumblingrobot.com/wp-content/ql-cache/quicklatex.com-676d7e52fd74f1c13b145cb9f7d5c9cf_l3.png "Rendered by QuickLaTeX.com")
This implies
and
- Proof. Let
be the Wronskian of and . If is any solution then the Wronskian of any pair of 
is a solution of 
. Hence,
![\[ W(x) = W(0) e^{-ax} \]](https://www.stumblingrobot.com/wp-content/ql-cache/quicklatex.com-32a662c8c54d07a270e7b22ef2960c99_l3.png "Rendered by QuickLaTeX.com")
This implies,
![\[ fv_1' - f'v_1 = (f(0)v_1'(0) - v_1(0) f'(0))e^{-ax} \]](https://www.stumblingrobot.com/wp-content/ql-cache/quicklatex.com-08ca39261846706ca39a47322b593d7c_l3.png "Rendered by QuickLaTeX.com")
and
![\[ fv_2' - f'v_2 = (f(0) v_2'(0) - f'(0) v_2(0)) e^{-ax}. \]](https://www.stumblingrobot.com/wp-content/ql-cache/quicklatex.com-5d51e10ee34a2aabf824bbb2e96df9b8_l3.png "Rendered by QuickLaTeX.com")
Since
is constant we know 
for any 
. Hence, 
. Therefore, by part (a),![\[ f = c_1 v_1 + c_2 v_2. \qquad \blacksquare \]](https://www.stumblingrobot.com/wp-content/ql-cache/quicklatex.com-982392706b2d13995dbc69f3ca715ba8_l3.png "Rendered by QuickLaTeX.com")
You skipped lots of nontrivial details in (b). Here is my proof using minimal symbols:
Define two new functions u1 and u2 just like in theorem 8.4, u1 for f(x) and u2 for the linear combination. If we manage to prove that u1=u2 then we are done. To do that we use part (a) to show that u1 and u2 have the same value and derivative at 0 and then apply theorem 8.5
In a, uniqueness of solutions looks to be assumed given that their value and derivative at some point (0) is equal. As far as I know Apostol only proved uniqueness for $\ddot{y} + by = 0$. It looks like its easy to extend that uniqueness prrof for the case $\ddot{y}+ a\dot{y} + by = 0$, but it wasn’t mentioned.
For (b), I solved as follows: Let v1 = u1*e^(-ax/2). Let v2 = u2*e^(-ax/2). Let f=u3*e^(-ax/2). Set y” + ay’ + by = (u” + (4b-a^2)u/4)*e^(-ax/2). From part (a) derive c1u1(0) + c2u2(0) = u3(0) nad c1u1′(0) + c2u2′(0) = u3′(0). Then by using theorem 8.5 you can conclude that c1u1 + c2u2 = u3. Therefore, c1v1 + c2v2 = f
Okay, one final attempt (sorry for all the posts). We started out needing to prove the existence of constants so that two solutions to a diff equation, of which the ratio is not constant, are in fact the basis for all solutions. These constants give us an equality at f(0). (I can’t explain why Apostol had us also find f'(0)). Then, using the properties of the Wronskian we find the existence of constants for arbitrary x, plug these into the equation in part a) for f and find a characterization for all solutions. I hope I haven’t missed something because then I would need to leave another post..
In the first part of the problem, what Apostol wanted was to solve the system of equations algebraically. It’s quite a bit of symbolic manipulation, but you don’t necessarily need any techniques from linear algebra here. Hence the reason for the two systems (one for f and one for f’), you solve for c1 in the first equation in terms of c2, plug it into the second, solve for c2, and then take c2 and plug it back into the equation for c1. It’s just an extension of the initial value problems that we would do earlier in the chapter. The constants c1 and c2 don’t change with x, the values f(0) and f'(0) were probably chosen for simplicity. It just as easily could have been f(a) and f'(a), where a is some arbitrary point on the interval we’re working with.
You had it right when you said that the Wronskian of v1 and v2 is not 0 since v2/v1 is not constant (the solution of a previous exercise), which is why constants c1 and c2 exist, since f, f’, v1, v1′, v2, and v2′ are all assumed to be continuous and differentiable on the interval of choice (i.e., c1 and c2 exist since there exist values for f(0), f'(0), v1(0), v1′(0), v2(0), and v2′(0), and that the denominator is never zero)
As for the second part of the problem, I did it a different (possibly wrong) way, as follows
And since y, v_1 and v_2 all satisfy the equation
Could be wrong, but there it is.
We can show (b) by using (a) and the uniqueness theorem for the solution of y’+by=0. Let g= c1v1 + c2v2 with the property of (a), namely g(0)=f(0) and g'(0)=f'(0), by linearity it fulfills the property g”+ag’+bg=0, we can let g(x)=u1(x)*e^(-ax/2) and f(x)=u2(x)*e^(-ax/2) for some u1 and u2. Then proceeding as in the theorem of uniqueness, we see that g”(x)+ag'(x)+bg(x)=e^(-ax/2) * (u1″(x)- D*u1(x)) and similarly with f where u2 is replaced with u1. Then since the expression above is equal to 0, we know that it is not the case that e^(-ax/2)=0, so the (uj”(x)- D*uj(x)) must equal to 0, for j=1,2, moreover, in order to be able to appeal to the theorem, we need to make sure that u1(0)=u2(0) and u1′(0)=u2′(0), which is direct, we just need to plug in 0 to g and f and use the fact that g(0)=f(0) and g'(0)=f'(0). Hence we have u1=u2, but this just means that g=f. I don’t know if I made an error or not in this proof, if there’s any fault, please point it out. I hope this helps.
Some more thoughts: The problem asks to prove f(0) and f'(0) are equal to the corresponding values at zero of a linearly independent combination of solutions. Since f is a solution of y” + by = 0, it seems the problem then becomes only to prove the basis is also a solution of y” + by = 0. So the fact that we are asked to find f'(0) seems to be relevant yet your solution doesn’t seem to require that information. Also, since W(x) is guaranteed to not be zero, could we have not skipped part a) altogether and just found the constants with W(x) instead of W(0)? I have only briefly looked at linear algebra so I can wait for those studies to get a better proof of these matters, or perhaps Apostol will give me some clues in the next chapter. But I would very much appreciate it if you could elucidate the strategy you took; ie how the basis is proven to be a solution.
Correction: y = e^(-ax/2)*u, where u is a solution of y” + by = 0. And so the purpose of calculating y'(0) in part a) is even more confusing now.
Might need a hint for the general part where you stated ‘Therefore, by part (a),’. It appears you equated the constants c1 & c2, which use values at zero, with similar versions of themselves but with arbitrary x’s, ie W(x) in the denominator rather than W(0). Is this the crux of the proof?