# How to solve ODE y'' - 4y' + y - x = 0 using rk4?

asked 2016-09-06 03:29:59 +0100

This post is a wiki. Anyone with karma >750 is welcome to improve it.

How to solve ODE y''+y''-4y'+x = 0 using rk4?

edit retag close merge delete

This looks like homework. If you want some help, you should ask more precise questions related to your research in solving those exercises, especially where you are locked.

( 2016-09-06 09:32:07 +0100 )edit

Sort by ยป oldest newest most voted

answered 2016-09-06 19:08:22 +0100

This post is a wiki. Anyone with karma >750 is welcome to improve it.

I tryed

 # w = y''
# z = y'

var('x, w, z')

desolve_system_rk4( [w , 4*z + x] , [w,z] , ics = [0,1,0], ivar = x, end_points = 10, step = 0.1)


But does no work correty :(

more
more

I got a way to solve the differential equation. If wrong please correct;)

To solve this differential equation we must first make a variable substitution to reduce the differential equation for a first order and thus create a system of ODEs. variable change:

w1 = y w2 = y w3 = y '

Thus, deriving the above variables, we have a system of equations of the form:

w1 '= w2 w2 '= f (x, w1, w2) = 4 * w2 + x

desolve_system_rk4([w2, 4*w2 + x ],[w1 , w2], ics = [0,1,0],ivar = x, step = 0.1, end_points = 2)

( 2016-09-06 23:18:07 +0100 )edit

Please start posting anonymously - your entry will be published after you log in or create a new account.

## Stats

Asked: 2016-09-06 03:29:59 +0100

Seen: 295 times

Last updated: Sep 06 '16