Численное решение уравнений и систем

Для численного решения систем уравнений, когда не удаётся найти аналитическое решение или по другим соображениям, можно использовать функцию вида

fsolve( eqns, vars, options ).

Здесь кроме уравнений eqns и переменных vars задаётся параметр options, управляющий решением. Этот параметр может принимать одно из следующих значений:

complex – ищется один корень (или все корни, для многочленов) на множестве комплексных чисел с плавающей запятой;

fulldigits – это значение параметра предотвращает потерю количества значащих цифр, задающихся параметром среды Digits, при промежуточных вычислениях. Эта опция замедляет работу функции;

maxsols= n – задаёт наименьшее количество разыскиваемых корней. Эта опция значима только для многочленов;

a .. b или x = a .. b или { x = a .. b, y = c .. d,... } – задание интервалов поиска корней. Конечные точки включены.

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

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

Примеры:

> fsolve(tan(sin( x ))= x, x );

1.557085816

> poly:= 23*x^5 + 105*x^4 - 10*x^2 + 17*x:

> fsolve( poly, x, -1..1 );

0, -.6371813185

> fsolve( poly, x, maxsols=3 );

0, -4.536168981, -.6371813185

> q:= 3*x^4 - 16*x^3 - 3*x^2 + 13*x + 16: fsolve( q, x, 1..2 );

1.324717957

> fsolve( q, x, complex);

-.6623589786 -.5622795121 I, -.6623589786 +.5622795121 I,

1.324717957, 5.333333333

> f:= sin (x+y) - exp (x)*y = 0: g:= x^2 - y = 2:

> fsolve( {f,g},{x,y},{x=-1..1,y=-2..0} );

{y = -1.552838698, x = -.6687012050}