Программирование в математическом пакете Mathcad. Воскобойников Ю.Е - 25 стр.

UptoLike

25
Simpson f a
,
b
,
N
,
()h
ba
N
Sfa() fb()
SS4faih
.
h
2
.
i0N 1
..
for
SS2faih
.
()
.
i1N 1
..
for
h
6
S
.
Adapt f a
,
b
,
( ) eps 10
8
I1 Simpson f a
,
b
,
5
,
()
I2 Simpson f a
,
b
,
10
,
()
I2 I1 I2 eps
<
if
Adapt f a
,
ab
2
,
Adapt f
ab
2
,
b
,
otherwise
Используя эти программы - функции вычислим определен-
ный интеграл от функции f(x) = x
2
на отрезке [0,1]. Точное зна-
чение интеграла равно 1/3=0.33333333333333... Обращение к про-
грамме-функции Adapt дает результат
Adapt f 0
,
1
,
( ) 0.333333333333333
=
Перед обращением к программе-функции Adapt необходимо
описать функцию пользователя f(x) в виде
f(x) : = x
2
,
так как имя функции f(x) используется в качестве фактического
параметра.
3.2. Модульное программирование
в нескольких документах MathCAD
                                        b       a
Simpson ( f , a , b , N )      h
                                            N
                               S        f( a )      f( b )
                                for i ∈ 0 .. N               1
                                                                            h
                                   S        S       4 .f a       i. h
                                                                            2
                                for i ∈ 1 .. N 1
                                  S S 2 . f ( a i. h )
                               h.
                                    S
                               6

                                   8
Adapt( f , a , b )    eps    10
                      I1    Simpson( f , a , b , 5 )
                      I2    Simpson( f , a , b , 10 )
                      I2 if I1          I2 < eps
                                        a       b                       a       b
                      Adapt f , a ,                   Adapt f ,                     ,b   otherwise
                                            2                               2

    Используя эти программы - функции вычислим определен-
ный интеграл от функции f(x) = x2 на отрезке [0,1]. Точное зна-
чение интеграла равно 1/3=0.33333333333333... Обращение к про-
грамме-функции Adapt дает результат
Adapt ( f , 0 , 1 ) = 0.333333333333333
Перед обращением к программе-функции Adapt необходимо
описать функцию пользователя f(x) в виде
                            f(x) : = x2,
так как имя функции f(x) используется в качестве фактического
параметра.

                      3.2. Модульное программирование
                     в нескольких документах MathCAD



                                                     25