У меня есть следующее ДУ и задание к нему: Определить методом коллокации приближенные решения следующих краевых задач. Использовать следующие базовые функции:
Я в этом во всем новичек поэтому нужна помощь! Проблема в реализации...В интернете мало очень информации о реализации этого метода на python или матлабе. Поэтому, уважаемые, прошу у вас помощи!
Еще формулы:
Она должна удовлетворять следующим условиям:
Неоднородные краевый условия:
Будем требовать, чтобы невязка
(8.40) равнялась нулю на
некоторой системе точек
x1,x2, ... ,xn
отрезка[a;b],которые называются точками коллокации, причем количество таких точек должно быть равно числу коэффициентов ci в выражении (8.39)
Тогда для определения c1, c2, .., cn получаем
систему уравнений:
Я начал писать код:
import sympy as sym
x = sym.symbols('x')
y = sym.Function('y')(x)
eq = sym.Eq(sym.diff(y, x, 2) + sym.diff(y, x) - y / x - 8 * x * x + 8 * x - 1.5, 0)
sol = sym.dsolve(eq)
fc = sol.rhs
eq1 = sym.Eq(fc.subs({x: 0}), 0)
dfc = sym.diff(fc, x)
eq2 = sym.Eq(fc.subs({x: 1}), 1)
csol = sym.solve([eq1, eq2], dict=True)
fres = fc.subs(csol[0])
print(fres)
Но я получаю следующую ошибку
raise NotImplementedError("The given ODE " + str(eq) + " cannot be solved by"
NotImplementedError: The given ODE -8.0*x**2 + 8.0*x + Derivative(y(x), x) + Derivative(y(x), (x, 2)) - 1.5 - y(x)/x cannot be solved by the factorable group method