Differential Equations problems

In summary: I'll keep working on the rest.7. {d^2y \over dx^2} +y = {1 \over 1+sin x}First we need to find the general solution of the complimentary equation y''+y=0the characteristic polynomial is m^2+1=0, which gives us m=+/-iso the complimentary solution is y_c=c_1cos x + c_2sin xnow for the particular solution, we'll use the method of undetermined coefficients. Since the right hand side is a trig function, we'll guess a particular solution of y_p=Acos x + Bsin xnow
  • #1
camilus
146
0
  1. Find the eigenvalues and eigenfunctions for the given boundary-value problem.

    [tex]{d \over dx}[x{dy \over dx}] + {\lambda \over x}{y} = 0, [/tex] subject to [tex] y(1)=0, y'(e^\pi)=0[/tex]

  2. Find the eigenvalues and eigenfunctions for the given boundary-value problem.

    [tex]{d \over dx}[{1 \over 3x^2+1}{dy \over dx}] + {\lambda (3x^2+1)}{y} = 0, [/tex] subject to [tex] y(0)=0, y(\pi)=0[/tex]

    hint: let [tex]t=x^3+x[/tex]
  3. Find the general solution of the Cauchy-Euler equation Assume x>0.

    [tex]x^3{d^3y \over dx^3} - 4x^2{d^2y \over dx^2} + 8x{dy \over dx} - 8y = 4ln(x)[/tex]
  4. Use teh variation of parameters to find the general solution of the given differential equation.

    [tex]{d^2y \over dx^2} +y = {1 \over 1+sin x}[/tex]

  5. Given that [tex]y_1(x)=x[/tex] is a solution of the DE [tex](1-x^2)y'' - 2xy' + 2y = 0, -1<x<1,[/tex] use the reduction of order to find a second solution (Do not use the formula.)
  6. Use the method of undetermined coefficients to find the general solution of the give DE.

    [tex]y''+y = x sinx[/tex]

I need a lot of help here! In a little bit, I'll post my work so far, which should lead you in the right direction.
 
Physics news on Phys.org
  • #2
camilus said:
I need a lot of help here! In a little bit, I'll post my work so far, which should lead you in the right direction.
Good. I'll wait for it.
 
  • #3
1.[tex]{d \over dx}[x{dy \over dx}] + {\lambda \over x}{y} = 0, [/tex] subject to [tex] y(1)=0, y'(e^\pi)=0[/tex]

differentiate that first term.

[tex]x{d^2y \over dx^2} +{dy \over dx} +{\lambda \over x}{y} = 0, [/tex]

multiply through by x to make it a Cauchy-Euler equation (making the powers of x and derivatives of y the same).

[tex]x^2{d^2y \over dx^2} +x{dy \over dx} +{\lambda}{y} = 0, [/tex] (1)

now using the substitution [tex]x=e^t [/tex] and [tex]t=ln x[/tex], and taking derivatives of t, you can substitute [tex]x^2y"= {d^2y \over dt^2} - {dy \over dt}[/tex] and [tex]xy'= {dy \over dt} [/tex]

When substituted, the [tex]{dy \over dt} [/tex] cancels out leaving

[tex]{d^2y \over dt^2} + \lambda{y} = 0[/tex] (2)

now getting the auxiliary/characteristic equation of (2) gives us
[tex]m^2+\lambda=0[/tex] or [tex]m=+/- i \sqrt{\lambda}[/tex]

that gives us the complimentary solution [tex]y_c = c_1cos(sqrt{\lambda}t) + c_2 sin (sqrt{\lambda}t)[/tex]

now since we're working with t instead of x, we need to modify our boundaries using t=ln(x). That gives us [tex]y(0)=0, y'(\pi)=0[/tex]

plugging in the first boundary y(0)=0, gives us [tex]c_1 = 0[/tex]

now this is where I'm stuck I think.

we got to get the first derivative of [tex]y_c[/tex], and we could ignore the cosine term because c1=0.

we get [tex]y'_c = c_2\sqrt{\lambda}cos(\sqrt{\lambda}t)[/tex], but when I plug in y(pi)=0, which I am not even sure if its right, I get stuck. All I know is I get to this line [tex]0=c_2\sqrt{\lambda}cos(\sqrt{\lambda}\pi)[/tex] and I am clueless. I know cos(n*pi)=-1 for n=1,3,5... and =1 for n=2,4,6...

now from my notes it says the eigenvalues are [tex]\lambda_n = {n^2\pi^2 \over L^2}[/tex] where L is the second boundary y(L)=y1, and the eigenfunctions are supposed to be [tex]y_n(x) = c_2 sin ({n\pi \over L}{x}) [/tex]for n=1,2,3,4...

****... going to go kill myself...
 
  • #4
number 5... watch this, this is sad...

Given that [tex]y_1(x)=x[/tex] is a solution of the DE [tex](1-x^2)y'' - 2xy' + 2y = 0, -1<x<1,[/tex] use the reduction of order to find a second solution (Do not use the formula.)

attempted solution:

if its given that y1=x, the method of reduction of order says that y2=u(x)y1(x).

so [tex]y_2=ux[/tex] and [tex]y'_2= u + u'x[/tex] and [tex]y''_2=u''x+2u'[/tex]

and putting the DE in standard form.

[tex]y'' - {2x \over 1-x^2}y' + {2 \over 1-x^2}y = 0[/tex]

plugging in [tex]y_2, y'_2, and y''_2[/tex] you get

[tex]u''x+2u' - {2x \over 1-x^2}(u + u'x) + {2 \over 1-x^2}(ux) = 0[/tex]

expanding cancels out the u term, leaving [tex]xu''+(2-{2x^2 \over 1-x^2})u' = 0[/tex], and divide by x, [tex]u''+({2 \over x}-{2x \over 1-x^2})u' = 0[/tex]

Now is where the reduction of order method gets its name. you let w=u', so w'=u", reducing the second-order equation into an eaiser first-order DE.

so [tex]w'+( {2 \over x}-{2x \over 1-x^2})w = 0[/tex] is the result after the substitution. But getting the integrating factor is what got me.

[tex]\mu(x)=e^{\int ({2 \over x}-{2x \over 1-x^2})dx}[/tex], which I would like someone to check this, I get [tex]\mu(x)=e^{ln[x^2(1-x^2)]}=x^2(1-x^2) = x^2-x^4[/tex]

Now if this mu is correct, than multiplying both sides of the DE by the mu will result by definition as

[tex]{d \over dx}{[(x^2-x^4)w] = 0*[x^2(1-x^2)] = 0[/tex]

But this is were I run into problems. Integrating this equation gives me

[tex](x^2-x^4)w= c_1[/tex], and [tex]w= {c_1 \over (x^2-x^4)}[/tex].

Substituting u' for w again, will give us u after integration. But its this integration that got me stuck.

I haven't been able to integrate [tex]u'= {c_1 \over (x^2-x^4)}[/tex]... Assuming I did everything right ealier, and I am supposed to integrate that..
 
  • #5
and I am pretty dumb, I think I am able to integrate that using integration by partial fractions.. going to try right now.
 
  • #6
Well I knocked off #5 I believe, I integrated by partial fractions and got

[tex]u' = c_1\left[ \int {{-1 \over x^2} + {-1/2 \over x+1} + {1/2 \over x-1}} \right]dx[/tex]

so
[tex]u = c_1\left[{1 \over x} + {1 \over 2}{ln(x-1)} - {1 \over 2}{ln(x+1)}\right] +c_2[/tex]

and for convinience, we'll let c1=1 and c2=0 rendering

[tex]u = {1 \over x} + {1 \over 2}{ln(x-1)} - {1 \over 2}{ln(x+1)}[/tex]

and since from the beginning we assumed that [tex]y_2(x) = y_1(x)u(x)[/tex] and y1=x, then [tex]y_2(x) = x \left [ {1 \over x} + {1 \over 2}{ln(x-1)} - {1 \over 2}{ln(x+1)} \right ][/tex]

sweet, 1 down, 5 to go...
 
  • #7
I think the general solution to #6 is

[tex]y(x) = c_1cos x + c_2sin x + xsin x[/tex]

let me know what you get. I am too lazy to type the work right now, I think its right so I rather work on another problem.
 
  • #8
Im almost done with #4 but I am stuck at the integration again. This time its serious.

for u2, I get [tex]u_2=\int {cosx \over 1+sinx}dx[/tex], which simple substitution u=1+sinx can integrate.

I get stuck on u1, where I get [tex]u_1=\int -{sinx \over 1+sinx}dx[/tex]...

Now this is one hell of an integral to compute, I can tell just by putting it into the Wolfram Mathematica Online Integrator. It says that

[tex]\int -{sinx \over 1+sinx}dx = - {{[ cos ({x \over 2}) + sin ({x \over 2}) ] [ xcos({x \over 2}) + (x-2)sin( {x \over 2} ) ]} \over 1+ sin x}[/tex]

how in the world would I do an integral by hand and get this crazy function!
 
  • #9
maybe integration by parts?
 
  • #10
camilus said:
we get [tex]y'_c = c_2\sqrt{\lambda}cos(\sqrt{\lambda}t)[/tex],
I don't think this is right. Remember that y'(e^pi) in the boundary conditions refers to dy/dx. Now you're working in dy/du. The chain rule for dy/dx = dy/du du/dx.

EDIT: This is in reply to qn 1, by the way.
 
  • #11
camilus said:
Well I knocked off #5 I believe, I integrated by partial fractions and got

[tex]u' = c_1\left[ \int {{-1 \over x^2} + {-1/2 \over x+1} + {1/2 \over x-1}} \right]dx[/tex]
That's odd. By partial fractions I got [tex]\frac{1}{x^2} + \frac{1}{2(1-x)} + \frac{1}{2(1+x)}[/tex]. Where did your minus sign come from?

camilus said:
I think the general solution to #6 is
[tex]y(x) = c_1cos x + c_2sin x + xsin x[/tex]
Actually it's a lot more complicated than that. This appears to be a straightforward application of the variation of parameters. The particular solution for y(x) is in the form [tex]y_p(x) = u(x)y_1(x) + v(x)y_2(x)[/tex], where y1 and y2 are individual linearly independent solutions to the homogenous DE.
 
Last edited:
  • #12
camilus said:
I get stuck on u1, where I get [tex]u_1=\int -{sinx \over 1+sinx}dx[/tex]...

Now this is one hell of an integral to compute, I can tell just by putting it into the Wolfram Mathematica Online Integrator. It says that

[tex]\int -{sinx \over 1+sinx}dx = - {{[ cos ({x \over 2}) + sin ({x \over 2}) ] [ xcos({x \over 2}) + (x-2)sin( {x \over 2} ) ]} \over 1+ sin x}[/tex]

how in the world would I do an integral by hand and get this crazy function!
This is for qn 4:

To integrate that you might want to use something known as t-formulae or tangent half-formulae. Or also known by its fancy name Weierstrass substitution.

EDIT: I just tried it out and it seems to work.
 
  • #13
it has nothing to do with the minus sign, what did you get for you integrating factor, mu? EDIT: I am going to make it BOLD on my work so you see what I am talking about.

and I' ll post my work to number 6 tomorrow, I am too lazy to write that all out right now. But the problem says use the method of undetermined coefficients.

you get the homogeneous case out the way, and for the particular solution you let y_p=(Ax+B)cos x (Cx+D)sin x. After the two differentiations and plugging in and cancelling terms, you're left with -2Asinx + 2Bcosx. Since these terms are duplicated in the homogeneous case, you multiply thru by x.

you you got that -2Axsinx+2Bxcos=xsinx

you match up coefficients, making B=0 and A=-1/2...
 
Last edited:
  • #14
camilus said:
it has nothing to do with the minus sign, what did you get for you integrating factor, mu?
My integrating factor is the same as yours.

and I' ll post my work to number 6 tomorrow, I am too lazy to write that all out right now. But the problem says use the method of undetermined coefficients.
Apparently I didn't see that, so yeah, I suppose you should use that instead of variation of parameters.

you get the homogeneous case out the way, and for the particular solution you let y_p=(Ax+B)cos x (Cx+D)sin x. After the two differentiations and plugging in and cancelling terms, you're left with -2Asinx + 2Bcosx. Since these terms are duplicated in the homogeneous case, you multiply thru by x.

you you got that -2Axsinx+2Bxcos=xsinx

you match up coefficients, making B=0 and A=-1/2...
This is with reference to qn 6?

EDIT: Woohoo! 999 posts. Great number. Maybe I should stop posting for a while just to retain that post number.
 
  • #15
camilus said:
number 5... watch this, this is sad...

Given that [tex]y_1(x)=x[/tex] is a solution of the DE [tex](1-x^2)y'' - 2xy' + 2y = 0, -1<x<1,[/tex] use the reduction of order to find a second solution (Do not use the formula.)

attempted solution:

if its given that y1=x, the method of reduction of order says that y2=u(x)y1(x).

so [tex]y_2=ux[/tex] and [tex]y'_2= u + u'x[/tex] and [tex]y''_2=u''x+2u'[/tex]

and putting the DE in standard form.

[tex]y'' - {2x \over 1-x^2}y' + {2 \over 1-x^2}y = 0[/tex]

plugging in [tex]y_2, y'_2, and y''_2[/tex] you get

[tex]u''x+2u' - {2x \over 1-x^2}(u + u'x) + {2 \over 1-x^2}(ux) = 0[/tex]

expanding cancels out the u term, leaving [tex]xu''+(2-{2x^2 \over 1-x^2})u' = 0[/tex], and divide by x, [tex]u''+({2 \over x}-{2x \over 1-x^2})u' = 0[/tex]

Now is where the reduction of order method gets its name. you let w=u', so w'=u", reducing the second-order equation into an eaiser first-order DE.

so [tex]w'+( {2 \over x}-{2x \over 1-x^2})w = 0[/tex] is the result after the substitution. But getting the integrating factor is what got me.

[tex]\mu(x)=e^{\int ({2 \over x}-{2x \over 1-x^2})dx}[/tex], which I would like someone to check this, I get [tex]\mu(x)=e^{ln[x^2(1-x^2)]}=x^2(1-x^2) = x^2-x^4[/tex]

Now if this mu is correct, than multiplying both sides of the DE by the mu will result by definition as

[tex]{d \over dx}{[(x^2-x^4)w] = 0*[x^2(1-x^2)] = 0[/tex]

But this is were I run into problems. Integrating this equation gives me

[tex](x^2-x^4)w= c_1[/tex], and [tex]w= {c_1 \over (x^2-x^4)}[/tex].

Substituting u' for w again, will give us u after integration. But its this integration that got me stuck.

I haven't been able to integrate [tex]u'= {c_1 \over (x^2-x^4)}[/tex]... Assuming I did everything right ealier, and I am supposed to integrate that..

strange, it didnt let me edit.. so I had 2 quote myself..
 
  • #16
Yes, and I did say that I got the same integrating factor as you did.
 
  • #17
My integrating factor is the same as yours.
Okay, then this is my work.
[tex]\mu(x)=e^{\int ({2 \over x}-{2x \over 1-x^2})dx} = e^{ln x^2}*e^{ln(1-x^2)}dx}= x^2(1-x^2) = x^2(x+1)(x-1)[/tex]

doing partials fractions on that is where I got my answer.


Apparently I didn't see that, so yeah, I suppose you should use that instead of variation of parameters.

This is with reference to qn 6?
Yes. question 6 is undetermined coefficients.

EDIT: Woohoo! 999 posts. Great number. Maybe I should stop posting for a while just to retain that post number.

please don't lol. I really need help, and the help is greatly appreciated.
 
  • #18
camilus said:
Okay, then this is my work.
[tex]\mu(x)=e^{\int ({2 \over x}-{2x \over 1-x^2})dx} = e^{ln x^2}*e^{ln(1-x^2)}dx}= x^2(1-x^2) = x^2(x+1)(x-1)[/tex]
Well for one thing [tex]x^2(1-x^2)[/tex] is not the same as [tex]x^2(x^2-1)[/tex].

please don't lol. I really need help, and the help is greatly appreciated.
Oh well I guess I'm past that anyway. Maybe next time when I'm at 1,111 posts.
 
  • #19
Defennder said:
Well for one thing [tex]x^2(1-x^2)[/tex] is not the same as [tex]x^2(x^2-1)[/tex].

Oh well I guess I'm past that anyway. Maybe next time when I'm at 1,111 posts.

yeah, there you go. I used 1-x^2, that's how the problem appears on my paper.

any more help?? especially with the first 2. I could probably get the rest on my own but those first two eigeinvalues/eigenfunctions are killing me!
 
  • #20
Defennder said:
That's odd. By partial fractions I got [tex]\frac{1}{x^2} + \frac{1}{2(1-x)} + \frac{1}{2(1+x)}[/tex]. Where did your minus sign come from?

HAHA I made the same mistake. I used 1-x^2, but when I split it I put (x-1)(x+1) instead of (1-x)(1+x)... so your right in the above quote. Although the work is the same, its just a matter of factoring out the minus sign you're mentioning.

EDIT: nevermind, neither of us is right.

its supposed to be [tex]\frac{1}{2(1+x)} - \frac{1}{x^2} - \frac{1}{2(1-x)}[/tex]
 
  • #21
You used x^2-1 and not 1-x^2. And as for the first two questions, I did reply earlier for qn 1. I haven't looked at 2 yet.
 
  • #22
Defennder said:
I don't think this is right. Remember that y'(e^pi) in the boundary conditions refers to dy/dx. Now you're working in dy/du. The chain rule for dy/dx = dy/du du/dx.

EDIT: This is in reply to qn 1, by the way.

I did this, instead of y(x) boundaries, transformed the equation using t=ln(x), so the boundaries instead of y(1)=0 and y'(e^pi)=0 become y(0)=0 and y'(pi)=0
 
  • #23
Well, with help from a friend on another forum, I think we solved #1.

we got:

CptBork said:
Yeah, you forgot to cancel a [tex]\pi[/tex] there in your answer when you plug in [tex]L=\pi[/tex]. You should just have [tex]\lambda_n=\frac{n^2}{4}[/tex], where n is odd, i.e. [tex]n=2k+1,\ k=0,1,2,\ldots[/tex]. Once you've finished solving and substituted back for [tex]x[/tex], the easiest thing to do if you're not sure is to check if your solutions satisfy the ODE and the boundary conditions.

I had wrote that first, [tex]\lambda_n=\frac{n^2}{4}[/tex], but then I changed it when I noticed we needed a multiple of pi in the eigenfunction. Because the eigenvalues in my notes had a pi..

Although I though I think you're right, and its a pi^2 and not a pi^4 in lambda_n. The pi will come from the initial conditions in this case, and cancel out with pi that's supposed to be there leaving just [tex]\lambda_n=\frac{n^2}{4}[/tex].

So (hopefully) we're done. The eigenvalues are [tex]\lambda_n=\frac{n^2}{4}[/tex], and the eigenfunctions are [tex]y_n = c_2sin({n \over 2}x)[/tex] where n is odd.

HURRAY! lol. Now that we (again, hopefully) knocked this one off, and we have an idea of what these eigenvalues/eigenfunctions problems are asking us to find, let's tackle number two. The differential equation is more complicated, but thankfully, we don't need a derivative, both boundaries are y(x)s, no y'(x)s..
 

1. What are differential equations?

Differential equations are mathematical equations that describe the relationship between a function and its derivatives. They are used to model and solve problems related to change and rates of change.

2. What types of problems can be solved using differential equations?

Differential equations can be used to solve a wide range of problems in various fields such as physics, engineering, economics, and biology. They are particularly useful in modeling systems that involve change over time, such as population growth, motion, and chemical reactions.

3. What is the process for solving a differential equation?

The process for solving a differential equation involves identifying the type of equation, finding its general solution, and then applying initial conditions to find a particular solution. This may involve using various techniques such as separation of variables, integration, and substitution.

4. Can differential equations have multiple solutions?

Yes, differential equations can have multiple solutions. In some cases, there may be an infinite number of solutions, while in others there may be a finite number of solutions. The uniqueness of a solution depends on the initial conditions given.

5. Are there any applications of differential equations in real life?

Yes, differential equations have numerous applications in real life. They are used to model and solve problems in fields such as engineering, physics, economics, and biology. For example, they are used to predict the spread of diseases, design bridges and buildings, and analyze electrical circuits.

Similar threads

  • Calculus and Beyond Homework Help
Replies
6
Views
792
  • Calculus and Beyond Homework Help
Replies
18
Views
1K
  • Calculus and Beyond Homework Help
Replies
5
Views
669
  • Calculus and Beyond Homework Help
Replies
15
Views
979
  • Calculus and Beyond Homework Help
Replies
20
Views
314
  • Calculus and Beyond Homework Help
Replies
2
Views
669
  • Calculus and Beyond Homework Help
Replies
10
Views
282
  • Calculus and Beyond Homework Help
Replies
3
Views
861
  • Calculus and Beyond Homework Help
Replies
8
Views
708
  • Calculus and Beyond Homework Help
Replies
2
Views
401
Back
Top