0

В общем вопрос такой, я написал метод эйлера и модифицированный эйлера как проверить теперь свою функциию и вывести на общий график с моими методами Вот код

import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate as spint
import scipy as sc
def f(x, y):
    return (x**2*np.sin(0.5*y))+2

x0 = 0; y0 = 1; X = 1;
n = 3



def Eyler_mod(x0, y0, f, X, n):#модифицированный метод
    h = X/n                    #Эйлера - с средней точкой
    x = np.zeros(n+1)
    y = np.zeros(x.shape)
    x[0] = x0
    y[0] = y0
    for i in range (n):
        x[i+1] = x[i]+(i+1)*h/2
        y[i+1] = y[i] + (h/2)*f(x[i], y[i])
    return x, y

def Eyler(x0, y0, f, X, n):
    h = X/n
    x = np.zeros(n+1)
    y = np.zeros(x.shape)
    x[0] = x0
    y[0] = y0
    for i in range (n):
        x[i+1] = x[i]+(i+1)*h
        y[i+1] = y[i] + h*f(x[i], y[i])
    return x, y
print(Eyler(x0, y0, f, X, n))
x, y = Eyler(x0, y0, f, X, n)
plt.plot(x, y)
plt.grid()
x, y = Eyler_mod(x0, y0, f, X, n)
plt.plot(x, y,marker='o',color="red")
plt.grid()
print(Eyler_mod(x0, y0, f, X, n))
2
  • 1
    Вывести так же, как и любую другую функцию. А что значит "проверить"? 30 апр 2019 в 11:49
  • @Эникейщик я это сделал таким образом t = np.linspace(0, 1, 11) y0 = [1.0,1.0] sol = spint.odeint(f, y0, t) plt.plot(t, sol[:, 0],marker='o',color="green") Но я не понимаю правильно ли я задал t и y0
    – Dobrotiu
    30 апр 2019 в 11:59

0

Ваш ответ

Нажимая на кнопку «Отправить ответ», вы соглашаетесь с нашими условиями использования и подтверждаете, что прочитали и поняли наши политику конфиденциальности и нормы поведения.

Посмотрите другие вопросы с метками или задайте свой вопрос.