# Can you please check my work?

## Homework Statement

Use the Runge-Kutta method to find approximate values of the solution of the initial-value problem y'+(x^2)y=sin xy, y(1)=pi; h=0.2 at the points xi=x0+ih, where x0 is the point where the initial condition is imposed and I=1, 2.

## Homework Equations

yn+1=yn+hf(xn+1/2h, yn+(1/2)hf(xn, yn))
f(x, y)=sin xy-(x^2)y
h=0.2, x0=1, y0=pi.

## The Attempt at a Solution

y1=pi+0.2f(1+0.1, pi+0.1f(1, pi))
=pi+0.2f(1.1, 2.83291)
=2.4669

y2=2.4669+0.2f(1.3, 2.4669+0.1f(1.2, 2.4669))
=2.4669+0.2f(1.3, 2.11683)
=1.76101

y1=2.475605264, y2=1.825992433.
I got y1=2.4669 and y2=1.76101.

LCKurtz
Homework Helper
Gold Member

## Homework Statement

Use the Runge-Kutta method to find approximate values of the solution of the initial-value problem y'+(x^2)y=sin xy, y(1)=pi; h=0.2 at the points xi=x0+ih, where x0 is the point where the initial condition is imposed and I=1, 2.

## Homework Equations

yn+1=yn+hf(xn+1/2h, yn+(1/2)hf(xn, yn))

Does that formula mean$$y_{n+1} = y_n + hf(x_n + \frac h 2,y_n + \frac h 2f(x_n,y_n))\text{?}$$

haruspex
Homework Helper
Gold Member
Using RK2 and OpenOffice Calc, I get y1=2.4636, y2 =1.8181, which are closer to your answers than to the book's.
However, if I cut the step size to 0.1 and, correspondingly, look at y2, y4 I get nearly the book answers: 2.4734, 1.8247.
This suggests the book is using a more accurate method than RK2 with step size 0.2 Maybe it uses RK4 with step size 0.2 - I haven't tried that.

@LCKurtz , yes, that's the right formula that I'm using.

LCKurtz
Homework Helper
Gold Member
Does that formula mean$$y_{n+1} = y_n + hf(x_n + \frac h 2,y_n + \frac h 2f(x_n,y_n))\text{?}$$

@LCKurtz , yes, that's the right formula that I'm using.

OK. I put that formula in Maple and I get: ##(x_1,y_1)= (1.2, 2.463635926),~(x_2,y_2) = (1.4, 1.818074242)## which agree with Haruspex's results.

• Math10

By the way, do I need to set the mode to radians or degrees on my calc?

What's RK2 and RK4?

LCKurtz
Homework Helper
Gold Member
All I can tell you for sure, is that using the formula you gave us, both Haruspex and I get identical answers independently. I doubt we are both wrong. And, speaking for myself, I am assuming the formula you gave is the appropriate one. If it isn't, all bets are off.

• Math10
LCKurtz
Homework Helper
Gold Member
By the way, do I need to set the mode to radians or degrees on my calc?

• Math10
LCKurtz
Homework Helper
Gold Member
• Math10
haruspex
• 