Квадратный корень из матрицы в maple. Основные объекты и команды Maple

В Maple имеется несколько способов представления функции.

Способ 1. Определение функции с помощью оператора присваивания (:= ): какому-то выражению присваивается имя, например:

> f:=sin(x)+cos(x);

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

> x:=Pi/4;

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

Чтобы насовсем не присваивать переменной конкретного значения, удобнее использовать команду подстановки subs({x1=a1, x2=a2,…, },f), где в фигурных скобках указываются переменные хi и их новые значения аi (i =1,2,…), которые следует подставить в функцию f . Например:

> f:=x*exp(-t);

> subs({x=2,t=1},f);

Все вычисления в Maple по умолчанию производятся символьно, то есть результат будет содержать в явном виде иррациональные константы, такие как, и другие. Чтобы получить приближенное значение в виде числа с плавающей запятой, следует использовать команду evalf(expr,t), где expr – выражение, t – точность, выраженная в числах после запятой. Например, в продолжение предыдущего примера, вычислим полученное значение функции приближенно:

> evalf(%);

Здесь использован символ (% ) для вызова предыдущей команды.

Способ 2. Определение функции с помощью функционального оператора, который ставит в соответствие набору переменных (x1,x2,…) одно или несколько выражений (f1,f2,…) . Например, определение функции двух переменных с помощью функционального оператора выглядит следующим образом:

> f:=(x,y)->sin(x+y);

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

Способ 3. С помощью команды unapply(expr,x1,x2,…) , где expr – выражение, x1,x2,… – набор переменных, от которых оно зависит, можно преобразовать выражение expr в функциональный оператор. Например:

> f:=unapply(x^2+y^2,x,y);

В Maple имеется возможность определения неэлементарных функций вида

посредством команды

> piecewise(cond_1,f1, cond_2, f2, …).

Например, функция

записывается следующим образом.

Основные объекты (определение, ввод, действия с ними)

Числа

Maple V работает с числами следующего типа:

целыми десятичными (О, 1, 123, -456 и т.д.),

рациональными в виде отношения целых чисел (7/9, -123/127 и т.д.),

радикалами,

вещественными с мантиссой и порядком (1.23E5, 123.456E-10)

комплексными (2+3*I)

Целые числа задаются в виде последовательности цифр от 0 до 9.

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

Обыкновенные дроби задаются с помощью операции деления лвух целых чисел.

Заметим, что Maple автоматически произвадит сокращение дробей. Над обыкновенными дробями можно выполнять все основные арифметические операции. Если при задании дроби ее знаменатель сокращается, то такая "дробь" трактуется программой Maple как целое число. Для преобразования обыкновенной дроби в десятичную служит команда evalf() . Второй параметр этой команды задает число значащих цифр. Заметим, что десятичное представление всего лишь аппроксимация точной величины, представленной обыкновенной дробью, т.е. дробь и ее десятичное представление не являются идентичными объектами Maple.

Радикалы задаются как результат возведения в дробную степень целых или дробных чисел, или вычисления из них же квадратного корня функцией sqrt() , или корня n -ой степени функцией surd(число, n) .

Числа с плавающей точкой задаются в виде целой и дробной частей, разделенных десятичной точкой. Их можно представтить также, используя так называемую экспоненциальную форму записи (для указания порядка применяется символ e или E ).

Константы

Maple содержит целый ряд предопределенных именованных констант - таких, к значениям котрых можно обращаться по имени. Часть этих констант не может быть изменена. К ним относятся:

Число е задается как exp(1) .

Посмотреть все константы, определенные в Maple, можно, исполнив команду: ?ininame

Кроме перечисленных на странице Cправки констант все прерменные, имена которых начинаются с _Env , по умолчанию являются системными константами Maple.

Строки

Строки - любой набор символов, заключенный в ДВОЙНЫЕ кавычки. Длина строки в Maple практически не ограничена и может достигать на 32-битных компьютерах длины 268 435 439 символов.

Переменные, неизвестные и выражения

Каждая переменная Maple имеет имя, представляющее последовательность латинских символов, начинающихся с буквы, причем строчные и прописные буквы считаются различными. Кроме букв в именах переменных могут использоватьсчя цифры и знак подчеркивания, но ПЕРВЫМ символом имени должна быть БУКВА.

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

неизвестная величина , а выражение, содержащее неизвестные, называется символьным выражением. Именно для работы с такими выражениями прежде всего и разрабатывался Maple.

Важной операцией в Maple, связанной с выражениями, является операция присваивания (:=) . Она имеет следующий синтаксис: переменная:= выражение; Здесь в левой части задается имя переменной, а в правой - любое выражение, которое может быть числовым, символьным или просто другой переменной.

Переменные позволяют хранить и обрабатывать разнообразные типы данных. При этом по умолчанию переменная Maple имеет тип symbol, предсталяющий символьную переменную, и ее значением является ее собственное имя. При присвоении перемнной какого-нибудь значения, ее тип изменяется на тип присвоенного ей значения.

Внутренняя структура объектов Maple

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

Команды, позволяющие выделять части объектов:

rhs(уравн)	Выделение правой части уравнения (или конца диапазона)
lhs(уравн)	Выделение левой части уравнения (или начала диапазона)
numer(дробь)	Выделение числителя числовой или алгебраической дроби
denum(дробь)	Выделение знаменателя числовой или алгебраической дроби
nops(выр)	Определяет количество операндов в выражении
op(выр) op(n,выр)	Выдает операнды выражения в виде списка, Извлекает n-ый операнд выражения
select(б ф, выр)	true
remove(б ф, выр)	Выделяет в выражении операнды, для которых булева функция дает значение false
indets(выр, тип)	Выделяет в выражении подвыражения заданного типа("*", "+" ...)

Познакомимся с этими командами более подробно.

Уравнение представляется в виде двух выражений, соединенных знаком равенства. Его не следует путать с операцией присваивания (:=). Уравнение является объектом Maple и служит для задания действительных уравнений. Его можно использовать в правой части операции присваивания, именуя тем самым уравнение.

В функции has() можно задать несколько подвыражений в виде списка. Ее результатом будет ИСТИНА тогда и только тогда, когда найдено хотя бы одно из подвыражений в списке.

Подстановка и преобразование типов

При выполнении математических преобразований часто необходимо произвести замену переменных в выражении, функции, уравнении и т.д., то есть вместо какой-то переменной подставить ее представление через некоторые другие переменные. А иногда необходимо выполнить преобразование выражения одного типа в другой. (Такое преобразование типов может потребоваться для выполнения некоторых команд, не работающих с исходным тиом выражения). Для этих целей в Maple существуетy несколько команд:

subs(подстановка, ВЫРАЖЕНИЕ)	Синтаксическая подстановка одного выражения вместо другого в ВЫРАЖЕНИЕ
algsubs(подстановка, ВЫРАЖЕНИЕ)	Алгебраическая подстановка одного выражения вместо другого в ВЫРАЖЕНИЕ
subsop(N=новое значение, ВЫРАЖЕНИЕ)	Подстановка нового значения вместо N-го операнда ВЫРАЖЕНИЯ
convert(ВЫРАЖЕНИЕ, тип)	Преобразует ВЫРАЖЕНИЕ в новый тип данных
whattype(ВЫРАЖЕНИЕ)	Определяет тип выражения.

Для подстановки вместо некоторой переменной (выражения) другого выражения служит команда subs() , синтаксис которой имеет следующий вид: subs(старое выражение=новое выражение, ВЫРАЖЕНИЕ) subs(s1, s2, .. sn, ВЫРАЖЕНИЕ) subs(, ВЫРАЖЕНИЕ) где каждое из s1,..sn является уравнением, определяющим подстановку.

Первая форма команды анализирует ВЫРАЖЕНИЕ , определяет в нем все вхождения старое выражение и подставляет вместо них новое выражение .

Вторая форма команды позволяет выполнить серию подстановок в ВЫРАЖЕНИЕ , причем подстановки выполняются последовательно, начиная с s1 . Это означает, что после выполнения первой подстановки, определенной s1 , Maple отыскивает вхождения левой части уравнения s2 во вновь полученном выражении и заменяет каждое такое вхождение на выражение, заданное в правой части уравнения s2 .

То есть вхождения выражений, заданных в левых частях уравнений s1, s2 , определяются в исходном параметре ВЫРАЖЕНИЕ . (см. примеры)

Воспользоваться командой simplify() , указав в качестве параметра требуемую замену (см. след раздел).

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

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

В Maple имеется несколько способов представления функции.

> f:=sin(x)+cos(x);

> x:=Pi/4;

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

> f:=x*exp(-t);

> subs({x=2,t=1},f);

> evalf(%);

Здесь использован символ (% ) для вызова предыдущей команды.

> f:=(x,y)->sin(x+y);

> f:=unapply(x^2+y^2,x,y);

В Maple имеется возможность определения неэлементарных функций вида

посредством команды

> piecewise(cond_1,f1, cond_2, f2, …).

Например, функция

записывается следующим образом.

Кафедра: Информационные Технологии

Лабораторная Работа

На тему: "СИНТАКСИС, ОСНОВНЫЕ ОБЪЕКТЫ И КОМАНДЫ СИСТЕМЫ MAPLE "

Москва, 2008 год

Цели работы :

· знать основные объекты и переменные системы Maple;

· знать и уметь применять команды, используемые при работе с объектами и переменными системы Maple;

· знать синтаксис основных математических функций системы Maple.

Введение

Система аналитических вычислений Maple – интерактивная система. В данном случае это означает, что пользователь вводит команду или оператор языка Maple в области ввода рабочего листа и, нажав клавишу , сразу же передает ее аналитическому анализатору системы, который выполняет ее. При правильном введении команды в области вывода появляется результат выполнения этой команды, если команда содержит синтаксические ошибки или ошибки выполнения, система печатает сообщение об этом. Если ошибку надо исправить, то следует вернуться к оператору, откорректировать его и снова выполнить. Выполнив введенную команду, система ожидает очередной команды от пользователя. Можно вернуться в любой момент к любой команде или оператору на рабочем листе, подкорректировать его и снова выполнить. Однако, если на рабочем листе есть команда, использующая результат вновь вычисленной, то ее следует также снова вычислить, установив на нее курсор, и, нажав клавишу , а если таких команд много, то можно выполнить команду графического интерфейсаEdit ® Execute ® Worksheet для повторного вычисления всех команд рабочего листа.

Каждый оператор или команда обязательно завершаются разделительным знаком. Таких знаков в системе Maple два – точка с запятой (;) и двоеточие (:). Если предложение завершается точкой с запятой, то оно вычисляется, а в области вывода отображается результат. При использовании двоеточия в качестве разделителя команда выполняется, но результаты ее работы не отображаются в области вывода рабочего листа. Это удобно, например, при программировании в Maple, когда нет необходимости в выводе каких-то промежуточных результатов, получаемых из операторов цикла, так как вывод этих результатов может занять много места на рабочем листе, да и может потребоваться значительное количество времени на их отображение.

Здесь и далее для команд Maple используется запись в форме синтаксиса языка Maple. Если при выполнении примеров возникает желание отображать команды в математической нотации, то следует командой Options ® Input Display ® Standard Math Notation установить соответствующий режим отображения.

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

Простейшими объектами , с которыми может работать Maple , являются числа, константы и строки.

Числа

Числа в системе Maple могут быть следующих типов: целые, обыкновенные дроби, радикалы, числа с плавающей точкой и комплексные. Первые три типа чисел позволяют выполнять точные вычисления (без округлений) разнообразных математических выражений, реализуя точную арифметику. Числа с плавающей точкой являются приближенными, в которых число значащих цифр ограничено. Эти числа служат для приближения (или аппроксимации) точных чисел Maple. Комплексные числа могут быть как точными, если действительная и мнимая части представлены точными числами, так и приближенными, если при задании действительной и мнимой частей комплексного числа используются числа с плавающей точкой.

Целые числа задаются в виде последовательности цифр от 0 до 9. Отрицательные числа задаются со знаком минус (–) перед числом, нули перед первой ненулевой цифрой являются не значащими и не влияют на величину целого числа. Система Maple может работать с целыми числами произвольной величины, количество цифр практически ограничено числом 2 28 . Вычисления с целыми числами реализуют четыре арифметических действия (сложение +, вычитание –, умножение *, деление /) и вычисление факториала (!).

Maple представляет большое целое число, которое не помещается в строке области вывода используя символ обратного слэша (\) в качестве символа продолжения вывода на следующей строке. Последняя команда вычисляет количество цифр в результате предыдущего вычисления. В ней в качестве параметра используется операция%, которая является всего лишь удобной формой ссылки на результат выполнения предыдущей операции. В Maple имеются еще две подобные операции, которые идентифицируют результаты предпредыдущей и предпредпредыдущей команд.Их синтаксис выглядит, соответственно, следующим образом:

В Maple имеется достаточно большой набор команд, позволяющих выполнить действия, специфичные при обработке целых чисел: разложение на простые множители (ifactor), вычисление частного (iquo) и остатка (irem) при выполнении операции целого деления, нахождения наибольшего общего делителя двух целых чисел (igcd), выполнение проверки, является ли целое число простым (isprime) и многое другое.

Для проверки вычисления частного и остатка двух целых чисел использованы операции получения результата выполнения предыдущей (вычисление частного) и предпредыдущей (вычисление остатка) команд. Результатом команды isprime () является булева константа true (истина) или false (ложь).

Набрав в области ввода рабочего листа команду? integer, можно получить список всех команд для работы с целыми числами

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

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

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

Дробь и ее десятичное представление не являются идентичными объектами Maple. Десятичное представление всего лишь аппроксимация точной величины, представленной обыкновенной дробью.

Радикалы задаются как результат возведения в дробную степень целых или дробных чисел, или вычисления из них же квадратного корня функцией sqrt(), или вычисления корня n ‑ой степени с помощью функции surd (число, n). Операция возведения в степень задается символом ^ или последовательностью из двух звездочек (**). При возведении в степень дробей их следует заключать в круглые скобки, как, впрочем, и дробный показатель степени. При задании радикалов также производятся возможные упрощения, связанные с вынесением из-под знака радикала максимально возможной величины.

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

Числа с плавающей точкой задаются в виде целой и дробной частей, разделенных десятичной точкой, с предшествующим знаком числа, например, 3.4567, -3.415. Числа с плавающей точкой можно задавать, используя так называемую экспоненциальную форму записи, в которой сразу же после вещественного числа с плавающей точкой или обычного целого, называемого мантиссой, ставится символ е или е, после которого задается целое число со знаком (показатель степени). Такая форма записи означает, что мантиссу следует умножить на десять в степени числа, соответствующего показателю степени, чтобы получить значение числа, записанного в такой экспоненциальной форме. Например, 2.345е4 соответствует числу 23450.0. Таким образом, можно представлять очень большие по абсолютному значению числа (показатель степени положительное число) или очень маленькие (показатель степени отрицательное число).

Из чисел составляются математические выражения с помощью арифметических операций. Символы арифметических операций в Maple перечислены в табл. 1.

Таблица 1. Арифметические операции

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

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

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

Maple умеет работать и с комплексными числами . Для мнимой единицы

в Maple используется константа I . Задание комплексного числа не отличается от его обычного задания в математике.

04. 01 Преобразование уравнений. Команды lhs и rhs

* Entering and Manipulating Equations: The lhs and rhs commands *

Напомним, что уравнению, точно так же как и выражению, можно присвоить имя. В следующей командной строке мы введём уравнение и присвоим ему имя " eq1 " :

> eq1:=x^3-5*x^2+23=2*x^2+4*x-8;

Мы можем вывести на экран левую и правую части уравнения по-отдельности при помощи команд lhs и rhs :

> lhs(eq1);

> rhs(eq1);

Воспользуемся командами lhs и rhs для того, чтобы привести уравнение к стандартному виду, в котором все члены собраны слева, а справа остался только 0:

> eq2:=lhs(eq1)-rhs(eq1)=0;

04. 02 Нахождение точных корней. Команда solve

* Finding Exact Solutions: The solve command *

Рассмотрим вначале рациональные уравнения. Известно, что существуют алгоритмы определения точных корней рациональных корней вплоть до 4-го порядка включительно. В Maple-команду solve и заложены эти алгоритмы.

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

> solve(3*x^3-4*x^2-43*x+84=0,x);

Обратите внимание: в команде мы указываем, относительно какой переменной следует решать уравнение. Хотя в нашем конкретном случае это и не обязательно:

> solve(3*x^3-4*x^2-43*x+84=0);

Maple нашел все 3 действительных корня и вывел их (в неупорядоченном виде ).

Иногда очень важно выбрать конкретный корень, чтобы потом использовать в дальнейших преобразованиях именно его. Для этого заранее следует присвоить имя результату исполнения команды solve . Назовём его X . Тогда конструкция X будет соответствовать первому корню из списка (подчеркнем: это не обязательно меньший корень! ), X - второму корню, и т.д. (Скобки - квадратные! ):

> X:=solve(x^2-5*x+3=0,x);

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

> x=%;

Ещё раз подчеркнём: практика показывает, что уравнению целесообразно присвоить имя. Традиционно в Maple такое имя начинается с букв eq :

> eq1:=7*x^3-11*x^2-27*x-9=0;

(Не путать оператор присваивания " := " со знаком равенства " = " !)

Теперь решим уравнение при помощи команды solve . Множеству корней присвоим имя X :

> X:=solve(eq1,x);

Для убедительности проверим, нет ли среди найденных корней посторонних. Проверку выполним непосредственной подстановкой

> subs(x=X,eq1);

Разумеется, частенько "точные" решения довольно громоздки, если не сказать иначе. Например, это касается уравнения :

> eq1:=x^3-34*x^2+4=0;

> X:=solve(eq1,x);

Теперь Вам понятно, о чем речь? Для себя заметьте, что мнимая единица в Maple обозначается посредством прописной буквы I . Разумеется, в таких случаях не грех найти приближенные значения корней. Имея на руках точное решение, Вы и сами сообразите, как это сделать:

> evalf(X);

В подобных ситуациях хорошей альтернативой команде solve является fsolve , особенности которой будут рассмотрены в следующем параграфе.

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

Решим уравнение :

> solve(5*exp(x/4)=43,x);

Иногда (а в тригонометрии - всегда ) Maple, по умолчанию , не выводит всё множество корней:

> solve(sin(x)=1/2,x);

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

Упражнение 4.1

Решить уравнение Разберитесь, сколько различных корней имеет уравнение. Как Maple поступает при наличии равных корней?

Совет : разложите на множители левую часть уравнения.

> solve(x^3-11*x^2+7*x+147=0,x);

> factor(x^3-11*x^2+7*x+147);

Корень х = 7 является двукратным, а потому у кубического уравнения только два различных корня. Разложение на множители левой части уравнения - тому подтверждение.

04. 03 Поиск приближенных корней. Команда fsolve

* Finding Approximate Solutions: The fsolve command *

Для приближенного решения уравнений используется Maple-команда fsolve . В случае рационального уравнения, fsolve выводит весь список действительных корней (см. Пример 01). Для трансцендентных уравнений эта команда, по умолчанию, выводит только один корень (см. Примеры 02 и 03).

При помощи fsolve найдём приближенные значения сразу всех четырёх действительных корней рационального уравнения :

> eq:=x^4-x^3-17*x^2-6*x+2=0;

> fsolve(eq,x);

Эти четыре корня и составляют исчерпывающее решение исходного рационального уравнения (хотя и приближенное ).

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

> eq:=x^3+1-exp(x)=0;

> fsolve(eq,x);

Maple и вывел только один корень. На этот раз Maple не стал "живописать". Как теперь убедиться в том, что других действительных корней нет? Следующий пример даёт такой инструментарий.

Получить все действительные корни уравнения и убедиться в этом.

Шаг первый ( Основная идея ) : найдём графическое решение уравнения. Для этого построим график функции, стоящей в левой части уравнения. Абсциссы точек пересечения этого графика с осью Ох и будут искомыми корнями.

> plot(x^3+1-exp(x),x=-3..5,y=-5..15);

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

Второй корень очевиден: х = 0. А как поточнее найти остальные?

Шаг второй ( Уточнение ) : применим команду fsolve более "зряче". В Maple предусмотрена возможность указания промежутка, на котором отыскиваются корни. В частности, для определения отрицательного корня нашего уравнения, укажем, что поиски следует вести в "районе" [-1;-0.2]. Об этом красноречиво свидетельствует графическое решение.

> fsolve(eq,x=-1..-.2);

Оставшиеся корни явно принадлежат промежуткам и . Расскажем об этом команде fsolve :

> fsolve(eq,x=1..2);
fsolve(eq,x=4..5);

Ну а что произойдёт, если мы подсунем Maple "пустой участок"? Например, отрезок для нашего уравнения. Там графического решения явно нет:

> fsolve(eq,x=2..4);

Maple выдаёт название команды, само уравнение, имя аргумента и отрезок. Т.е. ничего нового. Мол: "Ищите корни сами, а я не нашел".

Шаг третий ( Дополнительный анализ ) : Как теперь удостовериться в том, что найдены все корни , а не только в видимой области графического решения? Для этого следует расширить интервал поисков:

> plot(x^3+1-exp(x),x=-3..50,y=-10..15);

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

Попробуем в других местах: справа и слева от области найденных корней.

> fsolve(eq,x=5..50);

> fsolve(eq,x=-50..-1);

И здесь ни одного дополнительного корня! Поняв, что с влиянием показательной части уравнения всё ясно, делаем окончательные выводы.

Исчерпывающее решение уравнения состоит из четырёх корней: -.8251554597 , 0 , 1.545007279 , 4.567036837 .

Применим команду fsolve для приближенного решения трансцендентного уравнения .

Как и в предыдущем случае, найдём вначале качественное графическое решение. Для этого ещё нужно угадать, как разбросать по обеим частям уравнения его члены. Но графические возможности Maple настолько великолепны, что почти всегда можно собирать все члены уравнения с одной стороны.

Рассмотрим уравнение, равносильное данному: . Абсциссы точек пересечения графика функции, стоящей в левой части уравнения, с осью Ох и будут искомыми корнями.

> eq:=x^2/20-10*x-15*cos(x+15)=0;

> plot(lhs(eq),x=-10..10);

График указывает область поисков корней: промежуток . Настаёт черёд команды fsolve :

> fsolve(eq,x=1..2);

Корень найден. Но, очевидно, он - не единственный. Расширьте область поисков и ещё раз примените команду fsolve для отыскания второго корня.

Упражнение 4.2

Найти все действительные корни уравнения , начав с графического решения.

Построим график левой части уравнения:

> eq:=x^5-4*x^3+3*x^2+7*x-1=0;

> plot(lhs(eq),x=-5..5,y=-5..5);

В результате находим корни уравнения в первом приближении: -2 ; -1.5 ; 0 . Применим теперь команду fsolve без указания диапазона поиска (оценим возможности Maple ):

> fsolve(eq,x);

С удовлетворением отмечаем, что Maple выводит все три корня (Не будем забывать, что решали рациональное уравнение.)

Упражнение 4.3

Найти все корни уравнения . Воспользуйтесь графическим решением. Проверьте каждый корень непосредственной подстановкой.

Приведём уравнение к стандартному (для этого раздела) виду:

> eq:=x^2-2-ln(x+5)=0;

Теперь построим график левой части уравнения:

> plot(lhs(eq),x=-10..10);

По всей видимости, существует два корня. Один примерно равен -2, а другой, похоже, 2.

Применим команду fsolve , ограничив диапазон поиска:

> x:=fsolve(eq,x=-5..0);

> x:=fsolve(eq,x=1..3);

Непосредственной подстановкой проверим корни:

> evalf(subs(x=x,eq));

Обратите внимание: в обоих случаях истинного равенства нет. С учётом ошибок при округлении, разумное расхождение вполне допустимо.

Убедитесь в отсутствии других корней. Ответ обоснуйте.

Упражнение 4.4

Графики функций и дважды пересекаются на отрезке [-5;5].

а). Постройте в одной системе координат графики обеих функций и при помощи мыши найдите координаты точек пересечения.

b). Составьте уравнение, корнями которого являются абсциссы точек пересечения графиков.

c). Используйте команду fsolve для решения этого уравнения.

d). Используйте результаты из пункта с) для оценки ординат точек пересечения графиков.

e). У Вас не создалось впечатление, что линии могут пересекаться и в третьей точке с координатами (1;9)? Используйте fsolve и графические возможности Maple, чтобы убедиться в противном.

> y1:=10-x^2;

> y2:=4*sin(2*x)+5;

Теперь построим графики функций:

> plot(,x=-5..5);

Приближенные координаты точек пересечения: (-1.8, 6.6) и (2.75, 2) .

b) Составим уравнение:

> eq:= y1=y2;

с) Команда fsolve поможет найти соответствующие корни:

> x1:=fsolve(y1=y2,x=-4..0);

> x2:=fsolve(y1=y2,x=0..4);

d) Используем команду subs для определения соответствующих ординат точек пересечения:

> y:=subs(x=x1,y1);

> y:=subs(x=x2,y1);

Общие точки графиков: (-1.800,6.763) и (2.773,2.311) .

e) Графически исследуем окрестность точки х = 1:

> plot(,x=.5..1.5);

Команда fsolve на этот раз позволит доказать отсутствие корней вблизи точки х = 1:

> fsolve(y1=y2,x=.5..1.5);

04. 04 Решение уравнений в общем виде

* Solving Literal Equations *

Во многих случаях Maple находит решение уравнения в общем (символьном) виде. Речь идёт об уравнении (а не системы!), содержащем несколько переменных. Решение состоит в том, что какая-то из переменных выражается через остальные.

Пусть необходимо решить уравнение относительно переменной g . По привычке, используем команду solve . И она оправдывает наши надежды: