ВУЗ:
Составители:
Рубрика:
Ответы: 1.
1
a
3
ln
x
x+a
+
2x−a
2a
2
x
2
;
1
4b
5
arctg
x
b
+
x(3x
4
+8b
2
x
2
+9b
4
)
12b
4
(x
2
+b
2
)
3
;
5
16
x −
15
64
sin 2x +
+
3
64
sin 4x −
1
192
sin 6x; 2.
π
12a
;
ln 2
2
≈ 0.3465; 1.1827.
§ 26. Аналитическое решение дифференциальных
уравнений и систем
26.1. Решение дифференциального уравнения первого поряд-
ка. По умолчанию все пере менные в Maxima являются независимыми. Поэтому,
перед тем как приступить к з аданию и решению диффере н циального уравнения
y
′
= f(x, y), необходимо указать, что пе ременная y зависит от x:
-->
depends(y,x);
(%)
[y (x)]
Решим дифференциальное уравнение y
′
= 2 − y. Запишем его под именем
eqn:
-->
eqn: d iff(y,x)=2-y;
(%)
d
d x
y = 2 − y
Для решения дифференциального уравнения используется кома н д а ode2().
Решение запишем под именем sol:
-->
sol: o de2(eqn, y, x );
(%)
y = e
−x
(2 e
x
+ %c)
Постоянную c можно найти, если д аны начальные условия. Для этого есть
команда ic1(). Решение с начальными условиями (x = 0, y = 0 ) запишем под
тем же именем sol :
-->
sol: i c1(sol, x=0, y =0);
(%)
y = e
−x
(2 e
x
− 2)
Построим график полученной функции на отрезке x ∈ [0, 5]. Команда
rhs(sol) выдает только правую час ть вы ражения sol (т.е. отбрасыва ет “y =” ):
-->
wxplot2d(rhs(sol), [x,0,5])$
71
1 x 2x−a 1 x(3x4 +8b2 x2 +9b4 ) 5
Ответы: 1. a3 ln x+a + 2a2 x2 ; 4b5 arctg xb + 12b4(x2 +b2 )3 ; 16 x − 15
64 sin 2x +
3 1 π
+ 64
sin 4x − 192
sin 6x; 2. 12a
; ln22 ≈ 0.3465; 1.1827.
§ 26. Аналитическое решение дифференциальных
уравнений и систем
26.1. Решение дифференциального уравнения первого поряд-
ка. По умолчанию все переменные в Maxima являются независимыми. Поэтому,
перед тем как приступить к заданию и решению дифференциального уравнения
y ′ = f (x, y), необходимо указать, что переменная y зависит от x:
--> depends(y,x);
(%) [y (x)]
Решим дифференциальное уравнение y ′ = 2 − y. Запишем его под именем
eqn:
--> eqn: diff(y,x)=2-y;
d
(%) y =2−y
dx
Для решения дифференциального уравнения используется команда ode2().
Решение запишем под именем sol:
--> sol: ode2(eqn, y, x);
(%) y = e−x (2 ex + %c)
Постоянную c можно найти, если даны начальные условия. Для этого есть
команда ic1(). Решение с начальными условиями (x = 0, y = 0) запишем под
тем же именем sol:
--> sol: ic1(sol, x=0, y=0);
(%) y = e−x (2 ex − 2)
Построим график полученной функции на отрезке x ∈ [0, 5]. Команда
rhs(sol) выдает только правую часть выражения sol (т.е. отбрасывает “y =”):
--> wxplot2d(rhs(sol), [x,0,5])$
71
Страницы
- « первая
- ‹ предыдущая
- …
- 69
- 70
- 71
- 72
- 73
- …
- следующая ›
- последняя »
