logo search
mathcad

Гладкие системы

Когда известно, что решение является гладкой функцией, лучше использовать функцию Bulstoer вместо функции rkfixed. Функция Bulstoer использует метод Bulirsch-Stoer в отличие от метода Рунге-Кутты, используемого функцией rkfixed. В этом случае найденное приближенное решение будет немного более точным, чем то, которое возвратила бы функция rkfixed.

Список аргументов и матрица, возвращаемая функцией Bulstoer, те же, что и для функции rkfixed.

Bulstoer ( y, x1, x2, npoints, D)

y =

Вектор начальных условий размерности n, где n есть порядок дифференциального уравнения или число уравнений в системе (если решается система уравнений).

x1, x2 =

Граничные точки интервала, на котором ищется решение дифференциального уравнения. Начальные условия, заданные в векторе y, — это значение решения в точке x1.

npoints =

Число точек (не считая начальной точки), в которых ищется приближенное решение. При помощи этого аргумента определяется число строк (1 + npoints) в матрице, возвращаемой функцией Bulstoer.

D (x, y) =

Функция, возвращающая значение в виде вектора из n элементов, содержащих первые производные неизвестных функций.