Maple 9.5/10 в математике, физике и образовании

Дьяконов Владимир Павлович

NewtonTolerance=positive — задает погрешность при выполнении ньютоновских итераций в ходе построений.

Функция plot_real_curve вычисляет и строит алгебраическую кривую по точкам. Применение функции plot_real_curve показывает рис. 4.38.

Рис. 4.38. Примеры применения функции plot_real_curve

4.11. Векторные вычисления и функции теории поля

4.11.1. Пакет векторных вычислений VectorCalculus

В Maple 8

были существенно расширены возможности вычислений над векторами (пространственными объектами) и поверхностями. Для этого введен пакет VectorCalculus, который, при вызове, открывает доступ ко многим командам и функция векторного анализа, теории поля и приложений дифференциального исчисления [67, 68] (файл vc):

> restart; with(VectorCalculus); interface(showassumed=0);

Warning, the assigned names <,> and <|> now have a global binding

Warning, these protected names have been redefined and unprotected:

 *, +, Vector, diff, int, limit, series

[&x, *, +, ., <, >, <|>, AddCoordinates, ArcLength, BasisFormat, Binormal, CrossProduct, Curl, Curvature, Del, DirectionalDiff, Divergence, DotProduct, Flux, GetCoordinateParameters, GetCoordinates, Gradient, Hessian, Jacobian, Laplacian, LineInt, MapToBasis, Nabla, PathInt, PrincipalNormal, RadiusOfCurvature, ScalarPotential, SetCoordinateParameters, SetCoordinates, SurfaceInt, TNBFrame, Tangent, TangentLine, TangentPlane, TangentVector, Torsion, Vector, VectorField, VectorPotential, Wronskian, diff, evalVF, int, limit, series]

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

Пакет VectorCalculus ориентирован в первую очередь на решение задач математической физики, использующих методы теории поля и приложения дифференциального исчисления. Он оперирует такими привычными для физиков (разумеется, и для математиков) понятиями, как поток векторного поля, градиент, тор-сион, векторный потенциал и др. Приведенный ниже материал поясняет применение большинства функций этого пакета. Полезно просмотреть и файл VectorCalculus.mws, содержащий примеры его применения. В Интернете можно найти целую серию уроков по векторному анализу и теории поля в виде пакета Calculus IV или V (разработчик проф. J. Wagner).

4.11.2. Объекты векторных вычислений

Вектор в геометрическом представлении в данном пакете по умолчанию задается в прямоугольной системе координат:

> v := Vector( [x,y,z]);

v := хе_х + ye_y + ze_z

Здесь е_х, е_у и е_z — проекции единичного вектора е на оси координат х, у и z. Тип координатной системы (по умолчанию — прямоугольная) можно определить следующим образом:

> attributes(v);

coords = cartesian

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

VectorField(v, с)

где v — вектор и с — опционально заданный параметр в форме name[name, name,...], задающий тип координатной системы.

Можно изменить систему координат, например, задав (с помощью функции установки координат SetCoordinates) полярную систему координат:

> SetCoordinates(polar);

polar

> w := <r,theta>;

w: = r e_r + θ e₀

> attributes(w);

coords = polar

Аналогично можно задать вектор в сферической системе координат:

> SetCoordinates(spherical[r,phi,theta]);

spherical_r,φ,θ

> F := VectorField(<r,0,0>);

F.= rē_r

> attributes(F);

vectorfield, coords = spherical_r,φ,θ

Можно также сменить формат представления вектора и выполнить с ним некоторые простейшие векторные операции:

> BasisForrnat(false);

true

> v := <a,b,c>;

> BasisFormat(true);

false

> v;

ae_r + bе_φ +ce_θ

> SetCoordinates(polar);

polar

> MapToBasis(<r,theta>, 'cartesian');

r cos(θ)e_x + r sin(θ)e_y

> SetCoordinates(spherical);

spherical

> MapToBasis(<r,phi,theta>, 'cartesian');

r sin(φ)cos(θ)e_x + r sin(φ)sin(θ)e_y + r cos(φ)е_z

> SetCoordinates(spherical[r,phi,theta]);

spherical_r,φ,θ

> MapToBasis(VectorField(<r,0,0>), 'cartesian'[x,y,z]);

хē_х + yē_y + zē_z

Пакет VectorCalculus предусматривает возможность задания новой системы координат с помощью команды:

AddCoordinates(newsys, eqns, owrite)

где newsys — спецификация новой системы координат в виде symbol[name, name, …]; eqns — соотношения между координатами новой системы и прямоугольной системы координат, представленные в виде list(algebraic); owrite — заданное опционально равенство.