Методическое пособие по изучению математического пакета Maxima

Содержание

1. Введение

2. Установка и запуск программы

3. Описание интерфейса wxMaxima

4. Используемые обозначения для ввода команд

5. Решение задач элементарной математики

6. Решение задач линейной алгебры

7. Решение задач математического анализа

8. Решение обыкновенных дифференциальных уравнений

9. Задания для самостоятельной работы

Введение

Важное замечание: данное пособие разработано для Maxima 5.13 (RC2). В предыдущих версиях Maxima могут быть особенности записи выражений, отличные от указанных в данном пособии. Например, в предыдущих версиях Maxima названия встроенных функций были регистронезависимыми, а начиная с версии 5.11 это было изменено. Подробнее см. руководство пользователя для конкретной версии Maxima.

Maxima - программа для выполнения математических вычислений, символьных преобразований и построения графиков. С каждой новой версией в Maxima появляются новые функциональные возможности и виды решаемых задач.

Система аналитических вычислений Maxima идеально подходит как для изучения школьниками старших классов, так и студентами, его могут использовать профессиональные математики для проведения сложных расчетов и исследований.

Основными преимуществами программы являются:

• возможность свободного использования (Maxima относится к классу свободных программ и распространяется на основе лицензии GNU);

• возможность функционирования под управлением различных ОС (в частности Linux и Windows);

• небольшой размер программы (дистрибутив занимает порядка 23 мегабайт, в установленном виде со всеми расширениями потребуется около 80 мегабайт);

• широкий класс решаемых задач;

• возможность работы как в консольной версии программы, так и с использованием одного из графических интерфейсов (xMaxima, wxMaxima или как плагин (plug-in) к редактору TexMacs);

• расширение wxMaxima (входящее в комплект поставки) предоставляет пользователю удобный и понятный интерфейс, избавляет от необходимости изучать особенности ввода команд для решения типовых задач;

• интерфейс программы на русском языке;

• наличие справки и инструкций по работе с программой (русскоязычной версии справки нет, но в сети Интернет присутствует большое количество статей с примерами использования Maxima);

Установка и запуск программы

Скачать последнюю версию программы Вы можете с ее сайта в сети Интернет: http://maxima.sourceforge.net/. Русская локализация сайта: http://maxima.sourceforge.net/ru/

Рассмотрим пример установки программы под управлением операционной системы Windows XP. Не забудьте, что для установки программ, Вам необходимы права Администратора, если Ваш компьютер подключен к сети или по каким-то причинам права вашей учетной записи ограничены.

Полученный файл, например maxima-5.12.99rc2.exe (размер файла около 20 мегабайт), является исполняемым. Для начала установки программы достаточно нажать на него два раза мышкой.

Сразу появится окно выбора локализации (выбираем русский язык).

В следующем окне жмем «далее», внимательно читаем лицензионное соглашение, выбираем «я принимаю условия соглашения» и снова жмем «далее» (два раза).

ИНН

Выбираем путь установки программы (можно оставить его без изменения).

Переходим к выбору устанавливаемых компонент. Из всего перечисленного для нас «лишними» являются Пакеты поддержки языков Maxima.

Установка — Maxima

Выбор компонентов

Какие компоненты должны быть установлены?

Выберите компоненты, которые Вы хотите установить; снимите Флажки с компонентов, устанавливать которые не требуется. Нажмите «Далее*-, когда Вы будете готовы продолжить.

[Выборочная установка	В
0 Ядро Maxima с интерфейсом командной строки	70,4 Мб
0 Графическая оболочка wxMaxinna	5,8 Мб
0 Графическая оболочка ХМ axima	2,3 Мб
1 1 Пакеты поддержки языков Maxima
1 1 Испанский	3,8 Мб
1 1 Португальский	3,5 Мб
1 1 Бразильский Португальский	3,8 Мб

Текущий выбор требует не менее 78,9 Мб на диске.

< Назад

Далее >

Отмена

В следующих окнах нам будет предложено выбрать, куда помещать ярлыки для запуска программы (в меню «Пуск», на рабочий стол и т.д.).

Завершающим этапом будет окно с предложением начать установку. По окончании установки жмем «далее» и «завершить».

Интерфейс wxMaxima

Для удобства работы сразу обратимся к графическому интерфейсу wxMaxima, т.к. он является наиболее дружественным для начинающих пользователей системы.

Достоинствами wxMaxima являются:

• возможность графического вывода формул

Вывод формулы в xMaxima (текстовый режим) Вывод формулы в wxMaxima (графический режим)

• упрощенный ввод наиболее часто используемых функций (через диалоговые окна), а не набор команд, как в классической Maxima.

• разделение окна ввода данных и области вывода результатов (в классической Maxima эти области объединены, и ввод команд происходит в единой рабочей области с полученными результатами).

Рассмотрим рабочее окно программы. Сверху вниз располагаются: текстовое меню программы - доступ к основным функциям и настройкам программы.

В текстовом меню wxMaxima находятся функции для решения большого количества типовых математических задач, разделенные по группам: уравнения, алгебра, анализ, упростить, графики, численные вычисления. Ввод команд через диалоговые окна упрощает работу с программой для новичков. Например, пункт меню Анализ/Интегрировать позволяет вычислить определенный или неопределенный интеграл:

После ввода необходимых параметров, в рабочем окне мы увидим команду и результат вычисления:

(%il) integrate (3*х+5/х, х) ;

(%ol) 5 1од{х)+

2

Ниже располагается графическое меню - наиболее часто используемые функции для работы с файлами: открыть / сохранить / печать данных, а также функции правки - копировать / удалить / вставить текст и другие, затем расположено окно вывода результатов расчетов, поле ввода команд и кнопки

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

При написании подпрограмм и длинных команд нажмите на кнопку (многострочный ввод), которая располагается слева от окна ввода команд.

При использовании интерфейса wxMaxima, Вы можете выделить в окне вывода результатов необходимую формулу и вызвав контекстное меню правой кнопкой мыши скопировать любую формулу в текстовом виде, в формате TEX или в виде графического изображения, для последующей вставки в какой-либо документ.

С%i1) 'integrateк, -5, 5) ;

	5
	X2 +3 х +7 d
_ [		Копировать Копировать текст Копировать ТеХ
		| Копировать изображение |
		Б число с плавающей точкой
		Решить ... Решить численно ...
		Упростить выражение Факторизовать выражение Раскрыть выражение Подставить ,,,
		Интегрировать ... Дифференцировать ...
		График 2D ... График 3D ...

Также в контекстном меню, при выборе результата вычисления, Вам будет предложен ряд операций с выбранным выражением (например, упрощение, раскрытие скобок, интегрирование, дифференцирование и др.).

Ввод простейших команд в wxMaxima

Все команды вводятся в поле ВВОД, разделителем команд является символ ; (точка с запятой). После ввода команды необходимо нажать клавишу Enter для ее обработки и вывода результата. В ранних версиях Maxima и некоторых ее оболочках (например, xMaxima) наличие точки с запятой после каждой команды строго обязательно. Поэтому автор настоятельно рекомендует добавлять ; после каждой команды.

В случае, когда выражение надо отобразить, а не вычислить, перед ним необходимо поставить знак ' (одинарная кавычка). Но этот метод не работает, когда выражение имеет явное значение, например, выражение sin(7r) заменяется на значение равное нулю.

(%i34) integrate (x*x+3_f х) ;

з

(%о34) + 3 х

3

(%i35) ^т integrate (х*х+3,х) ; (%о35)

Обозначение команд и результатов вычислений

(%il) х: 3 ; (%ol) 3

(%i2) у: 10; (%о2) 10

(%i3) х/у; з

(%оЗ) — ю

После ввода, каждой команде присваивается порядковый номер. На приведенном выше рисунке, введенные команды имеют номера 1-3 и обозначаются соответственно (%i1), (%i2), (%i3). Результат вычисления также имеет порядковый номер, например (%o1), (%o2) и т.д. Где "i" - сокращение от англ. input (ввод), а "о" - англ. output (вывод).

х² +3 dх

Этот механизм позволяет избежать в последующих вычисления повторения полной записи уже выполненных команды, например (%i1)+(%i2) будет означать добавление к выражению первой команды - выражения второй и

последующего вычисления результата. Также можно использовать и номера результатов вычислений, например (%o1)*(%o2).

Для последней выполненной команды в Maxima есть специальное обозначение - %.

Пример: Вычислить значение производной функции

2 ¹

у(х) = х +— в точке х=1.

(%i 4) diff(х^Л2+1/х,х) ; 1

(%о4) 2 х

г

к

(%i5) (%о4),к=1; (%о5) 1

Используемые обозначения для ввода команд Ввод числовой информации.

Правила ввода чисел в Maxima точно такие, как и для многих других подобных программ. Целая и дробная часть десятичных дробей разделяются символом точка. Перед отрицательными числами ставится знак минус. Числитель и знаменатель обыкновенных дробей разделяется при помощи символа / (прямой слэш).

Обратите внимание, что если в результате выполнения операции получается некоторое символьное выражение, а необходимо получить конкретное числовое значение в виде десятичной дроби, то решить эту задачу позволит применение оператора numer. В частности он позволяет перейти от обыкновенных дробей к десятичным.

(%i 4 б) 3/7 + 5/3; 44

(%о 4 б) — 21

(%i 47) 3/7 + 5/3,numer; (%о 47) 2.095238095238095

Арифметические операции

Обозначение арифметических операций в Maxima ничем не отличается от классического представления: + , - , * , /. Возведение в степень можно обозначать несколькими способами: ^А , ^АА , **. Извлечение корня степени n записываем, как степень 1/n. Введем еще одну полезную операцию - нахождение факториала числа. Эта операция обозначается восклицательным знаком, например 5!.

Для увеличения приоритета операции, как и в математике, используются круглые скобки: ().

Константы

В Maxima для удобства вычислений есть ряд встроенных констант, самые распространенные из них показаны в следующей таблице:

Название	Обозначение
слева (в отношении пределов)	minus
справа (в отношении пределов)	plus
(плюс бесконечность)	inf
-да (минус бесконечность)	minf
Комплексная бесконечность	infinity
п (число Пи)	%pi
e (экспонента)	%e
Мнимая единица )	%i
Истина	true
Ложь	false
Золотое сечение (1 + sqrt(5))/2	%phi

Переменные

Для хранения результатов промежуточных расчетов применяются переменные. Заметим, что при вводе названий переменных, функций и констант важен регистр букв, так переменные x и X это две разные переменные.

Присваивание значения переменной осуществляется с использованием символа: (двоеточие), например x:5.

Если необходимо удалить значение переменной (очистить ее), то применяется метод kill:

kill(x) - удалить значение переменной x;

kill(all) - удалить значения всех используемых ранее переменных.

(%i18) (%olB)

(%i19) (%ol9)

(%i20) (%o20)

(%i21) (%o21)

Математические функции

В Maxima имеется достаточно большой набор встроенных математических функций. Вот некоторые из них:

тригонометрические функции: sin (синус), cos (косинус), tan(тангенс), cot (котангенс);

• обратные тригонометрические функции: asin (арксинус), acos (арккосинус), atan (арктангенс), acot (арккотангенс);

• sec (секонс, sec x = 1/cos x), csc (косеконс, csc x = 1/sin x);

• sinh (гиперболический синус), cosh (гиперболический косинус), tanh (гиперболический тангенс), coth (гиперболический котангенс), sech (гиперболический секанс), cosh (гиперболический косеканс);

• log (натуральный логарифм);

• sqrt (квадратный корень);

• mod (остаток от деления);

• abs (модуль);

• min(x1,...,xn) и max(x1,...,xn) - нахождение минимального и максимального значения в списке аргументов;

• sign (определяет знак аргумента: pos - положительный, neg - отрицательный, pnz - не определен, zero - значение равно нулю);

Правило записи функций

х: 5 ; 5

х; 5

kill (х) ; done

х;

х

Для записи функции необходимо указать ее название, а затем, в круглых скобках записать через запятую значения аргументов. Если значением аргумента является список, то он заключается в квадратные скобки, а элементы списка также разделяются запятыми.

sin(x);

integrate(sin(x),x,-5,5); plot2d([sin(x)+3,cos(x)],[x,-%pi,%pi],[y,-5,5]); Пользовательские функции

Пользователь может задать собственные функции. Для этого сначала указывается название функции, в скобках перечисляются названия аргументов, после знаков := (двоеточие и равно) следует описание функции. После задания пользовательская функция вызывается точно так, как и встроенные функции Maxima.

Пример:

Нужно помнить, что не следует использовать для функций названия, зарезервированные для встроенных функций Maxima.

Перевод сложных выражений в линейную форму записи

Одним из самых сложных занятий для начинающих пользователей системы Maxima является запись сложных выражений, содержащих степени, дроби и другие конструкции, в линейной форме (в текстовой форме записи, при помощи ASCII символов, в одну строку).

Для облегчения данного процесса можно дать несколько рекомендаций (основаны на многолетнем опыте преподавания Maxima студентам и школьникам):

Математическая запись	Линейная запись выражения
к + 2 у -7	(x+2)/(y-7)
( х + 3 )2 ^	(x+3)**(2*y)

• Функция не существует отдельно от своих аргументов (если таковые имеются). Поэтому, например, при возведении в степень, можно взять всю функцию с аргументами в скобки, а потом уже возводить полученную конструкцию в нужную степень: (sin(x))**2 Очень часто пользователи пытаются возвести в степень только название функции, забывая про аргументы: sin**2(x) - это не правильно!

• Также помните, что аргументы функции записываются в скобках, через запятую. Недопустима запись функции sin(2*x) в виде sin*2*x или sin2x. Сначала напишите название функции, затем поставьте после нее пустые скобки, и только потом, в этих скобках, напишите все ее аргументы, разделяя их запятыми. Никаких конструкций между названием функции и открывающейся скобкой быть не должно!

• В случае записи сложного выражения разбейте его на несколько простых составляющих, введите их по отдельности, а затем объедините, используя рассмотренные ранее обозначения введенных команд.

Пример:

Необходимо ввести следующее выражение:

Данное выражение можно разделить на 3 части, записать их по отдельности и

затем соединив, используя результаты промежуточных расчетов

Решение задач элементарной математики

К таким задачам можно отнести вычисление и преобразование арифметических выражений, построение графиков функций, решение уравнений и систем алгебраических уравнений.

Рассмотрим ряд полезных функций:

assume - ввод ограничений и forget - снятие ограничений

(%i2 0) (%o20)

(%i21) assume (к<0) ; (%o21) [ x< 0 ]

(%i22) sqrt(х^л2) ; (%o22) -x

(%i23) forget(x<0); (%o23) [ x< 0 ]

(%i2 4) sqrt (х^л2) ; (%o2 4) |x|

divide - нахождение частного и остатка от деления одного многочлена на другой

(%i 10) divide (х^Л3-2, х-1) ; (%о10) f х² +х + 1 , - 1 J

Первый элемент полученного массива - частное, второй - остаток от деления.

factor - разложение на множители

(%±19) factor(а*х^Л2+а*х+а); (%о19) а(х² + х + l)

(%i2 0) factor(х^Л2+2*х+1) ; (%о20) (х + I)²

expand - раскрытие скобок

(%i9) expand ( (2 + 3*х)*(3*у+5*х)) ; (%о9) 9 ху + 6у + 15 х² + 10 х

gcd - наибольший общий делитель многочленов

(%il2) а/ (5*х) +b/x-c/x;

с Ъ а (%о12) +— +-

гг к 5 к

(%i 13) ratsimp (а/ (5*х)+b/x-c/x) ; 5 с - 5 ъ - а

(%о13)

5 х

partfrac - преобразовать в простые дроби по заданной переменной Пусть дано выражение

к

(%i43) tan (х) *sec (х) ^л2 + ( (1-sin (х) ^л2) *cos (х) ) / cos (х) ^л2;

1 -sin(к) cos ( к)

. г

(%о 43) sec(x) tan( х.) +-

(%i 4 4) trigsimp(%) ;

sin ( к) + cos ( к)

(%o 4 4)

cos ( к)

trigexpand (тригонометрическое раскрытие скобок) - использует формулы преобразования сумм двух углов для представления введенного выражения в как можно более простом виде - где в качестве аргумента только одна переменная

(%i25) trigexpand(sin(u+v)*cos(u) ^л3) ;

3

(%o25) cos(u) ( cos ( и ) sin( v) + si n ( и ) cos ( v))

trigreduce (приведение тригонометрическое) - преобразует тригонометрическое выражение к сумме элементов, каждый из которых содержит только sin или cos

(%i26) trigreduce (sin (u+v) *cos (u) ^A3) ;

sin{ v + 1 u) + sin( 7 - 2 u) 3 sin{ v + 2 u)+3 sin( r)

(%o2 6) +

8 8

Нахождение корней уравнений и систем алгебраических уравнений

solve - решение алгебраических уравнений и их систем (в качестве параметров в первых квадратных скобках указывается список уравнений через запятую, во-вторых - список переменных, через запятую)

(%i27) solve ( [х+у=0,х-у=3] , [х,у]) ;

(%о27) [ [_X=~_f у = -~] ]

2 2

В случае множества решений у тригонометрических уравнений выдается соответствующее сообщение только и одно из решений

(%i2 8) solve ([sin(x)=0]_f [ х] ) ;

"solve' is using arc-trig functions to get a solution. Some solutions will be lost. (%o2 8) [x = 0 ]

Также Maxima позволяет находить комплексные корни (%i29) solve ([х^Л2 + х+1=0] , [x] ) ;

VT4i + 1 -\ZlT4i - 1

(%o2 9) [x= ,x= ]

2 2

Построение графиков функций

В Maxima имеется несколько альтернативных библиотек для отображения графиков функций, наборов точек, трехмерных тел, градиентов и т.д. Рассмотрим несколько примеров построения графиков. По умолчанию используется библиотека Plot, которая и будет рассмотрена ниже. Также для решения некоторых задач рекомендуется ознакомиться с возможностями библиотеки Draw (см. справочную систему Maxima).

Сразу стоит отметить, что при построении графика стоит обращать внимание на масштаб делений по осям. И, если, вместо ожидаемой окружности, мы увидим эллипс, то, скорее всего это произошло именно из-за несоответствия масштаба делений по осям.

Пример 1. Двумерный график. На графике показаны три зависимости. График построен с использованием функции plot2d. Первый аргумент - список функций, второй и третий - ограничения по осям координат. Третий аргумент является необязательным. Если его не указать - он будет подобран автоматически.

x

Пример 2. Построение графиков функций, заданных параметрически. Приведем пример построения графика функции с зависимостью, заданной параметрически:

(%i31) plot2d ([parametric, cos(t), sin(t), [t,-%pi,%pi], [nticks,80]], [x, -4/3, 4/3]);

Параметр nticks,80 указывает сколько опорных точек используется для построения графика. В данном примере, при значении 4 получим треугольник, при значении 5 - четырехугольник и т.д.

Пример 3. Построение графиков функций в полярных координатах. На первом этапе необходимо задать функцию для построения графиков в полярных координатах.

(%i11) plot_polar(expr,range):=block([theta_var: range[1]],plot2d( ['parametric,cos(theta_var)*expr, sin(theta_var)*expr,range,[nticks,80]]))

Теперь для построения дуги необходимо ввести

cos(t)

Приведем еще один пример построения графика в полярных координатах. На этот раз для функции sin(t)+2*cos(2*t).

-2.5

2.5

1.5

1

0.5

0

-0.5

-1

-1.5

-2

-2

-1

-0.5

0

0.5

1

1.5

-1.5

Пример 4. Построение графика по точкам.

Зададим координаты как два массива значений (в первом координаты х точек, во втором - координаты у):

(%i21) xx:[10, 20, 30, 40, 50];

(%i22) у у : [ . б, .9, 1.1, 1.3, 1.4]; и вызовем функцию построения графика:

Еще раз отметим опорные точки на осях координат (на этот раз укажем Maxima не соединять точки между собой):

(%i27) plot2d([discrete,хх,уу] , [style, points]) ;

1.4 1.3 1.2 1.1 1

0.9 0.8 0.7 0.6 0.5

я

га тз

■ffi ®

10

15

20

25

35

40

45

50

30

Пример 5. Затабулировать функцию и построить ее график.

Сначала создадим список абцисс точек табуляции. Его можно задать как в

виде функции (см. ниже), так и просто записав координаты в квадратных

скобках, через запятую (как было рассмотрено ранее).

Затем задаем функцию и вызываем функцию построения графика.

(%i 9) х:makelist (х/10, х, 0, 10);

113213749

(%о9) [ О — — — , 1 ]

10 5 10 5 2 5 10 5 10 (%i10) f (х) :=х^Л2+5*х+3;

(%о10) JT(x): = x +5х + 3

(%ill) plot2d([dis crete,x,f(x)],[style,points]);

6

8

7

го тз <u <u

5

4

3

0

0.2

0.4

0.6

1

0.8

Пример 6. Построение трехмерных графиков Для начала рассмотрим простой пример:

(%i39) plot3d (2^Л(-u^A2 + v^A2), [u, -3, 3], [v, -2, 2]);

Пример 7. В следующем примере «трехмерность» заменена градиентным переходом цветов.

(%i40) plot3d (cos (-хЛ2 + ул3/4), [x, -4, 4], [y, -4, 4],

Более подробно о построении графиков можно посмотреть в справочной системе Maxima.

Решение задач линейной алгебры

В данном разделе будут рассмотрены основные операции с матрицами. Для задания матрицы используется функция matrix:

(%i 4 8) matrix ([17, 3], [-8, 11]); 17 з

(%о48)

- 8 11

Операции с матрицами

[gnuplot_preamble, "set view map; unset surface"], [gnuplot_pm3d, true], [grid, 150, 150]);

Рассмотрим основные операции с матрицами на примере, для этого зададим две матрицы: x и у.

к: matrix ([17, 3], [-8, 11]); 17 3 - 8 ll

(%i5) y: matrix ([%pi, %e]_f [a, b]); %pi %e

(%o5)

(%о6)

Поэлементное сложение, вычитание, умножение и деление матриц, возведение в степень.

(%i 6) х + у;

%pi + 17 Не + 3

(%о7)

а - 8 to + 11

(%i7) х - у;

(%i4) (%o4)

17 - %pi 3 -

- а - 8 11 - to

(%iВ) х * у;

(%о8)

17 %рi 3 - 8 а 11 to

Пример применения функции к матрице:

(%ill) sqrt (x) ;

л/ IT' Л/з^

<■ 1? о з

(%o14)

(%oll)

^3/Ј%i л[ТГ

(%i13) exp (x) ; (%o13)

о " s . 11

Умножение матриц

(%i14) x.y;

3 a + IT %pi 3 to + IT 11 a - 8 %pi 11 ь - 8

Возведение в степень (в данном примере получается обратная матрица)

Функции для работы с матрицами

determinant - нахождение определителя матрицы:

(%i5 0) determinant(matrix ([17, 3], [-8, 11])); (%о50) 211

eigenvalues - нахождение собственных значений матрицы: (%il7) eigenvalues (х) ;

(%о17) f Г 14 —J~15~4i , -\flS~4i +14 ],[ 1 г 1 ] ]

invert - получение обратной матрицы:

(%i55) invert (matrix ([17, 3], [-8, 11]));

11	3
Zll	Zll
0	17
Zll	Zll

(%o55)

minor - определяет минор матрицы. Первый аргумент - матрица, второй и третий - индексы строки и столбца соответственно:

ll

(%i51) mino г (mat rix ([17, 3], [-8, 11] ) , 1, 1) ;

(%о51)

rank - ранг матрицы:

(%i53) rank(mat rix ([17, 3], [-8, 11])); (%o53) 2

submatrix - возвращает матрицу, полученную из исходной удалением соответствующих строк и (или) столбцов. В качестве параметров следуют номера удаляемых строк, исходная матрица, номера удаляемых столбцов.

(%i22)

(%о22)

(%о23)

(%i5 4) (%o5 4)

М: matrix ([3, 7, аа, bb] , [-1, 8, 5, 2], [9, 2, 11, 4]); 3 1 аа hh

-18 5 2 9 2 11 4

(%i23) submatrix(1,2,М);

9 2 11 4

(%i2 4) submatrix(1,М,2,3) ; - l 2

(%o2 4)

9 4

transpose - транспонирование матрицы:

transpose (matrix ([17, 3], [-8, 11])); 17 - 8

3 11

Решение задач математического анализа

Нахождение производной

Функция diff позволяет найти производные, как первого, так и более высоких порядков. При наличии у функции нескольких переменных можно найти частную производную по одной из них.

Синтаксис: diff(функция, переменная, порядок производной);

Пример: найти первую производную функции y(x)=e^x/x².

Сначала введем функцию: y(x):=exp(x)/x^A2; (обратите внимание, что в отличие от присвоения значения переменной, здесь используется комбинация символов ":=" (двоеточие и равно)), а затем найдем ее производную по переменной х.

Для этого введем команду: diff(y(x),x,1); или diff(y(x),x);. В случае первой производной ее порядок можно не указывать.

у (х) : =ехр (х) /х^Л2;

ехр ( гг)

у(х):=

z

к

(%i56) (%о56)

При вычислении кратных производных по нескольким переменным после указания функции перечисляются переменные дифференцирования с указанием соответствующих кратностей, например, diff(xA8*yA5, x, 4, y, 2);

Интегрирование

Для нахождения неопределенного интеграла в качестве аргументов указывается функция и переменная интегрирования.

Синтаксис: integrate(функция, переменная);

2

integrate(х^Л2+5*х+3, x) ;

3 2

+ + 3 x

3 2

Пример: вычислить интеграл функции x +5x+3 по переменной х.

(%i59) (%о59)

При нахождении значения определенного интеграла помимо рассмотренных параметров указываются пределы интегрирования.

В качестве пределов интегрирования могут фигурировать бесконечность (inf) и минус бесконечность (minf).

Синтаксис: integrate(функция, переменная, левый предел, правый предел);

Пример: вычислить интеграл функции sin(x) по переменной х на отрезке от 0 до п.

(%i60) integrate(sin(x)_fx_fO_f%pi);

В случае, когда интеграл расходится Maxima выдает сообщение "Integral is divergent".

Например integrate(1/x, x, 0, inf); выдаст именно такое сообщение.

В некоторых случаях Maxima может попросить доопределить некоторую переменную, как в случае интегрирования функции xⁿ:

Ts п + 1 zero or nonzero ? nonze ro ; Нахождение пределов

Синтаксис: limit( функция, переменная, значение); или limit( функция, переменная, значение, слева/справа);

Предел слева обозначается minus, а справа - plus. Пример 1:

Разложение в ряд Тейлора

Синтаксис: 1ау1ог(функция, x, a, n). Параметр n определяет, до какой степени (x-a) находить решение.

Пример: Разложить функцию sin(x)+ax+1 в ряд Тейлора по степеням x+1 до третьей степени.

Нахождение суммы ряда

Синтаксис: Биш(функция, переменная, индекс начального члена ряда, индекс конечного члена ряда).

Пример 1:

В данном примере числовой результат получается сразу после записи выражения.

Пример 2:

В данном примере необходимо упростить полученный результат, применив simpsum.

(%i77) sum (l/3^Ai, i, 1, inf) ;

(%i78) sum (l/3^Ai, i, 1, inf) , simpsum;

l

(%o7 8) -

2

Решение обыкновенных дифференциальных уравнений Пример 1. Решение ОДУ первого порядка.

+ 3 ху = sin( х)

/

Записываем ОДУ, обращая внимание на одинарную кавычку (напомним, что она позволяет отображать, но не вычислять некоторое выражение). Затем применяем функцию ode2. Параметры: уравнение, зависимая переменная, независимая переменная. На завершающем этапе подставляем значение начального условия задачи: у(п)=0, для этого воспользуемся функцией ic1 (начальная задача для ОДУ первого порядка).

(%il) х^Л2* ^т diff(у,к) + 3*у*х = sin(x)/x; f _л \ „i

d

(%ol) к'

(%i2) ode2(%, у, х);

- cos ( к ) (%о2) у=

(%i3) icl (%o2_fx=%pi,y=0) ;

X

+ 3 х у =-

cos(к ) + 1

( %o3 ) y=

Пример 2. Решение ОДУ второго порядка.

, г

( d

- о

Ах

-у + у

₄dx у

d

Первые два этапа схожи с предыдущим примером.

(%i2) ^т diff (у_f x_f 2) + y*^Tdiff (у,x)^л3 = 0; , г

( d ^

(%о2)

= 0

d к

-у + у

_vdx у

(%i3) ode2(%,у,х);

у +6 kkl у

(%оЗ)

Для подстановки начального условия используем функцию ic2 (начальное условие для ОДУ второго порядка) и сразу упрощаем полученный ответ при помощи уже рассмотренной функции ratsimp.

(%i5) rat simp (ic2 (%o3, x=0 , y=0 , Miff (y,x)=2) ) ;

2 у ^ — 3 у

(%o5) =x

В случае граничной задачи применяется функция bc2 (в качестве параметров

указываются значения функции в двух точках)

(%i6) bc2 (%оЗ, х=0, у=1, х=1, у=3) ;

у³ - ю у з

(%об) =х

6 2

Пример 3. Решить задачу Коши для дифф. уравнения второго порядка:

на отрезке [1,3], при начальных условиях y'(1)=1, y(1)=1.

На первом этапе записываем дифф. уравнение и решаем его при помощи

функции ode2.

Задаем начальные условия.

(%i14) ic2 (%о12, х=1,у=1, ^т diff (х, у) =1) ;

6 logг (х )

Методическое пособие по изучению математического пакета Maxima 1

Содержание 1

Введение 3

Установка и запуск программы 3

Интерфейс wxMaxima 6

Ввод простейших команд в wxMaxima 13

Обозначение команд и результатов вычислений 13

Используемые обозначения для ввода команд Ввод числовой информации. 15

Арифметические операции 15

Константы 16

Переменные 16

Математические функции 17

Правило записи функций 18

Пример: 20

Перевод сложных выражений в линейную форму записи 20

Пример: 21

Решение задач элементарной математики 22

Нахождение корней уравнений и систем алгебраических уравнений 25

Построение графиков функций 25

Решение задач линейной алгебры 47

Операции с матрицами 47

Функции для работы с матрицами 51

Решение задач математического анализа 54

Нахождение производной 54

Интегрирование 54

Синтаксис: integrate(функция, переменная); 54

Нахождение суммы ряда 56

Пример 1: 56

Пример 2: 56

Решение обыкновенных дифференциальных уравнений Пример 1. Решение ОДУ первого порядка. 56

Ответ построим в виде графика функции на отрезке [1,3]

Пример 4. Нахождение числового значения ОДУ методом Рунге-Кутта 4-го порядка.

Решить численно ОДУ для t от 0 до 8 с шагом 0.1, при х=1 в т. t=0.

clx ₇ = t -х

Подключаем необходимую библиотеку:

(%il8) load("dynamics") ;

Вызываем метод Рунге-Кутта 4-го порядка, результат заносим в переменную sol.

(%i2 4) sol: rk (t-x^A2, x, 1, [t, 0, 8, 0 . 1] ) ;

Отобразим полученный результат графически. Подключаем графическую библиотеку и вызываем метод draw2d, для отображения полученного решения.

(%i23) load("draw");

(%i26) draw2d (points (sol) )

^{dY 2} 2 . ₁ = у - X +1

dt

Подключаем необходимую библиотеку: (%ilB) load("dynamics") ;

Вызываем метод Рунге-Кутта 4-го порядка, результат заносим в переменную

(%i19)

sol: rk( [4-х^Л2-4*у^Л2,у^Л2-х^Л2 + 1] , [х,у], [-1.25,0.75], [t,0,4,0.02]) ;

Подключаем графическую библиотеку и вызываем метод draw3d, для отображения полученного решения.

(%i23) load("draw");

(%i2 4) draw3d (points (sol) )

0.8 0.6 0.4 0.2 0 -0.2 -0.4 -0.6 -0.8

Задачи для самостоятельной работы

Для самопроверки рекомендуется использовать задачи из Сборника типовых расчетов по курсу высшей математики. Учебное пособие. /Под ред. В.Б. Миносцева, -М: МГИУ, 2004. -582с.