This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: ACM 95/100b Problem Set 7 Solutions Faisal Amlani March 16, 2009 Problem 1 Consider the ODE d dx bracketleftbigg (1 x 2 ) dy dx bracketrightbigg + y = x 2 over the interval 1 x 1 with boundary conditions y ( 1) finite y (1) finite Solve this problem using a series of Legendre polynomials. That is represent y ( x ) = summationdisplay n =0 A n P n ( x ) and solve for the A n . Does the solution exist for all values of ? SOLUTION : We know that the Legendre polynomials are the complete eigenfunctions of the homogeneous equation with these boundary conditions, each associated with eigenvalues n = n ( n + 1) . It is of Sturm Liouville form with r ( x ) = 1 . So from the beginning of the Week 7 Notes, we have that for y ( x ) = A n n ( x ) = A n P n ( x ) , the solution to the inhomogeneous equation is y ( x ) = summationdisplay n =0 f n  n P n ( x ) where f n are the coefficients of the expansion of the right hand side f ( x ) = x 2 in terms of the eigen functions (Legendre polynomials.) That is, f n = integraltext 1 1 r ( x ) f ( x ) n ( x ) dx integraltext 1 1 r ( x ) 2 n ( x ) dx = integraltext 1 1 x 2 P n ( x ) dx integraltext 1 1 P 2 n ( x ) dx = 2 n + 1 2 integraldisplay 1 1 x 2 P n ( x ) dx from the orthogonality relation of P n . But from the Rodrigues formula, we know that P ( x ) = 1 1 ( x 2 1) = 1 , P 1 ( x ) = 1 2 d dx ( x 2 1) 1 = x, P 2 ( x ) = 1 2 2 2! d 2 dx 2 ( x 2 1) 2 = 1 8 (12 x 2 4) = 1 2 (3 x 2 1) , . . . 1 We easily suspect that we can express x 2 in terms of just P and P 2 . Indeed, from f = 1 2 integraldisplay 1 1 x 2 dx = 1 2 23 = 1 3 , f 1 = 3 2 integraldisplay 1 1 x 2 ( x ) dx = 3 2 (0) = 0 , f 2 = = 5 2 integraldisplay 1 1 x 2 1 2 (3 x 2 1) = 5 4 8 15 = 2 3 , we have that 1 3 P + 2 3 P 2 = 1 3 + 2 3 1 2 (3 x 2 1) = x 2 . So f 3 , f 4 , ... must be zero and we have the solution y ( x ) = summationdisplay n =0 f n  n P n ( x ) = f  P ( x ) + f 2  2 P 2 ( x ) = 1 / 3 + 2 / 3  2(2 + 1) ( 1 2 (3 x 2 1)) = 1 3 + 1 3 1  6 (3 x 2 1) = x 2 2 (  6) . Clearly, the solution doesnt exist for = 0 , 6 . Problem 2 Consider the Chebyshev ODE given by d dx bracketleftbigg radicalbig 1 x 2 dy dx bracketrightbigg + y ( x ) 1 x 2 = 0 1 x 1 (a) Classify the singular points at x = 1 (b) Using a Taylor series about x = 0 generate Taylor series for two linearly independent solutions. Hint: it may be advantageous to rewrite the ODE so that the leading coefficient of y is a poly nomial....
View
Full
Document
This note was uploaded on 12/12/2009 for the course ACM 95b taught by Professor Nilesa.pierce during the Winter '09 term at Caltech.
 Winter '09
 NilesA.Pierce

Click to edit the document details