Решение дифференциальных уравнений в приложении MathCAD

x₀:=0 y₀:=1

x_i+1=x_i+h

y_i+1=y_i+h×(0.2×y_i+x_i)

Для получения численных значений записываются выражения:

x= y=

Имеются и встроенные функции для решения дифференциальных уравнений, например, встроенная функция rkfixed.

Чтобы решить систему дифференциальных уравнений, приведенную выше, можно записать:

R:=rkfixed(y, x1, x2, Np, D)

Здесь x1, x2 – левая и правая границы интервала, на котором ищется решение; Np – число точек, в которых определяется решение; y_0, y₁ – начальные условия; D(x,y) – вектор правых частей системы. Для определения матрицы с решениями надо набрать R=

Можно также построить графики решения для различных значений i, характеризующие зависимость R(1,i) от R(0,i) и зависимость R(2,i) от R(0,i), то есть зависимость y от x и зависимость z от x.