226
er
1 и
er
2 – абсолютная и относительная точность вычислений (действительные
числа).
Пример 24.2.
Вычислить интеграл из примера 24.1.
-->integrate('(2*x-1)^0.5','x',5,13)
ans =
32.666667
24.3. Интегрирование внешней функции
Наиболее универсальной командой интегрирования в Scilab является:
[I,err]=intg(a,b,name[,er1 [,er2]]),
где
name
– имя функции, задающей подынтегральное выражение; здесь
функция может быть задана в виде набора дискретных точек (как таблица) или
с помощью внешней функции;
a
,
b
– пределы интегрирования;
er
1 и
er
2 – абсолютная и относительная точность вычислений (действительные
числа).
Пример 24.3.
Вычислить интеграл из примера 24.1
-->deff('y=G(x)','y=sqrt(2*x-1)')
-->intg(5,13,G)
ans =
32.666667
Пример 24.4.
Вычислить интеграл
( )
(
)
1
2
0
3
t
dt
sin t
+
∫
-->function y=f(t),y=t^2/sqrt(3+sin(t)),endfunction;
-->[I,er]=intg(0,1,f)
er =
1.933D-15
I =
0.1741192
24. 4. Приближенное дифференцирование, основанное на
интерполяционной формуле Ньютона
Идея численного дифференцирования заключается в том, что функцию
y
(
x
), заданную в равноотстоящих точках
x
i
(
i
=0, 1, ...,
n
) отрезка [
a
,
b
] с
помощью значений
( )
i
i
xf
y
=
приближенно заменяют интерполяционным
