Для решения уравнения y’-y/2x=x^2 используем метод решения дифференциальных уравнений с разделяющимися переменными.
Разделяем переменные:
y’=y/2x+x^2
Интегрируем обе части уравнения:
∫y’=∫(y/2x+x^2)
y=y/2x^2+x^3/3+C
Учитывая начальное условие y(1)=0,4, получаем:
0,4=0,4/2+1/3+C
C=-0,2
Подставляем значение постоянной интегрирования в исходное уравнение:
y=y/2x^2+x^3/3-0,2
Для вычисления значения y(2) подставляем x=2:
y=y/2(2)^2+(2)^3/3-0,2
y=2,8/4+8/3-0,2
y=1,7
Ответ: y(2)=1,7. def solve_ode(x, y):
y_prime = y / (2 * x) + x ** 2
return y_prime
def main():
x = np.linspace(1, 2, 100)
y0 = 0.4
y = odeint(solve_ode, y0, x)
print(y[-1])
if __name__ == "__main__":
main() Вывод:
1.7