Scribd
Upload
26 views

Codigos Python

The document shows code for numerically solving differential equations using the Runge-Kutta method of order 4. It imports NumPy and matplotlib, defines functions for the differential equation and Runge-Kutta steps, runs the simulation, and plots the results.

Uploaded by

Carlos Daniel Dias Soler
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
26 views

Codigos Python

The document shows code for numerically solving differential equations using the Runge-Kutta method of order 4. It imports NumPy and matplotlib, defines functions for the differential equation and Runge-Kutta steps, runs the simulation, and plots the results.

Uploaded by

Carlos Daniel Dias Soler
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

import numpy as np

import matplotlib.pyplot as plt

n = 120
h = 0.25
def f(x,y):
dy = (0.2-2*10**(-6)*y)*y
return dy

y=[]
x=[]
y.append(1000)
x.append(0)

for k in range(0,n):
F1 = f(x[k],y[k])
F2 = f(x[k]+h/2,y[k]+h*F1/2)
F3 = f(x[k]+h/2,y[k]+h*F2/2)
F4 = f(x[k]+h,y[k]+h*F3)

x.append(x[k]+h)
y.append(y[k]+h*(F1+2*F2+2*F3+F4)/6)

for k in range(0,n+1):
print("%12.2f & %12.12f \n"%(x[k],y[k]) )

fig = plt.figure(figsize=(10, 8))


ax = fig.add_subplot(111)

ax.plot(x,y)
fig.savefig('graph_rk4.png')
import numpy as np

import matplotlib.pyplot as plt

from matplotlib import style

style.use('ggplot')

% matplotlib inline

% config Inlinebackend.figure_format='svg'

def dydx(x, y):

return -20*y+20*x+21

n= 10

x0= 0

xf= 1

y0= 0.6

h= (xf-x0)/n

x=[x0]

y=[y0]

for i in range(n):

y0=y0 + h*dydx(x0, y0)

y.append(y0)

x0=x0 + h # x0 +=h

x.append(x0)

print(x)

print(y)

You might also like