Have you seen the proof that the set of all solutions to an nth order homogeneous linear equation form an n dimensional vector space? That is the fundamental concept behind all such problems!
We can think of a 1 dimensional vector space as a line through the origin of a coordinate system, a 2 dimensional vector space as a plane containing the origin, etc. A "linear manifold" would be a line, plane, etc. that does NOT contain the origin. Such a thing is NOT a vector space because the sum of two "objects" in it is not again in it- addition is not "closed". For example, the graph of y= 3x is a line that contains the origin. Two points on that line are (a, 3a) and (b, 3b). The sum of those two is (a+ b, 3a+ 3b)= (a+ b, 3(a+ b)), again on the line so this set is "closed under addition". Multiplying (a, 3a) by the number p gives (pa, p(3a))= (pa, 3(pa)), again on the line. The points on this line form a "vector space" with that addition and scalar multiplication.
The graph of y= 3x+ 2 is a line that does not contain the origin. The points (a, 3a+ 2) and (b, 3b+ 2) are on that line. Their sum is (a+ b, 3a+ 3b+ 2)= (a+ b. 3(a+ b)+ 2) which is NOT on that line. Similarly, p time (3a, 3a+ 2) gives (p(3a), p(3a)+ 3p) NOT on that line. But if, instead of adding like this, we define sum of (a, 3a+ 1) and (b, 3b+ 1) to be (a+ b, 3(a+ b)+ 1) (geometrically this would be "go from (a, 3a+ 1) and (b, 3b+ 1) to (a, 3a), (b, 3b) on the parallel line through the origin, moving along the mutual perpendicular to those parallel lines, add those, then go back to the original line") we can define an addition of such things (and a similar "scalar multiplication").
We do the same kind of thing with linear differential equations. The set of all solutions to a linear, nth order, homogeneous differential equation form a vector space with the usual sum of functions and multiplication by numbers. If y1 and y2 both satisfy y''+ p(x)y'+ y= 0 so does ay1+ by2 for any numbers a and b: (ay1+ by2)''+ p(x)(ay1+ by2)'+ q(x)(ay1+ by2)= a(y1''+ p(x)y1'+ q(x)y1)+ b(y2''+ p(x)y2'+ q(x)y2)= a(0)+ b(0)= 0.
If, rather, y1 and y2 satisfy y''+ p(x)y'+ q(x)y= f(x), with f(x) NOT 0, ay1+ by2 does not satisfy the equation: (ay1+ by2)''+ p(x)(ay1+ by2)'+ q(x)(ay1+ by2)= a(y1''+ p(x)y1'+ q(x)y1)+ b(y2''+ p(x)y2'+ q(x)y2)= a(f(x)+ b(f(x))= (a+ b)f(x) NOT f(x). Instead suppose y1 and y2 satisfy y''+ p(x)y'+ y= 0 and y3 satisfies y''+ p(x)y'+ y= f(x). Then ay1+ by2+ y3 satisfies y''+ p(x)y'+ q(x)y= f(x) for any a, b: (ay1+ by2+ y3)"+ p(x)(ay1+ by2+ y3)'+ q(x)(y1+ y2+ y3)= (ay1''+ ap(x)y1'+ aq(x)y1)+ (by2''+ bp(x)y2'+ bq(x)y2)+ (y3''+ p(x)y3'+ q(x)y3)= a(0)+ b(0)+ f(x)= f(x).
Further if y1 and y2(x) are solutions to the entire equation, y''+ p(x)y'+ q(x)y= f(x), then y3= y1- y2 satisfies y''+ p(x)y'+ q(x)y= 0 (I will leave that "as an exercise for the reader"). So y1= y3+ y2- a function that satisfies the "associated homogenous equation" plus a function that satisfies the entire equation.