# Can you please check my work?

1. Jan 4, 2015

### Math10

1. The problem statement, all variables and given/known data
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.

2. Relevant 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.

3. 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.

2. Jan 4, 2015

### LCKurtz

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{?}$$

3. Jan 4, 2015

### haruspex

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.

4. Jan 5, 2015

### Math10

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

5. Jan 5, 2015

### LCKurtz

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.

6. Jan 5, 2015

7. Jan 5, 2015

### Math10

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

8. Jan 5, 2015

### Math10

What's RK2 and RK4?

9. Jan 5, 2015

### LCKurtz

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.

10. Jan 5, 2015

### LCKurtz

11. Jan 5, 2015

### LCKurtz

12. Jan 5, 2015

### haruspex

It is RK2.

13. Jan 5, 2015

### Math10

Thank you guys for the help!