Эффективный метод вычисления интеграла Адамара на конечном интервале. Добрынина Н.Ф - 23 стр.

UptoLike

ВУЗ: 

ПГУ | Пенза

Составители: 

Рубрика: 

23
.00002241296235( ) − t
k
1.
9
− ) d
⌠
⌡
⎮
⎮
⎮
⎮
t
k
t
+ k 1
+
1
() − τ + .5 .01 I
2
1
() − τ − .5 .01 I
2
τ
⎞
⎠
⎟
⎟
⎟
⎟
⎟
⎟
⎟
.2660992120 d
⌠
⌡
⎮
⎮
⎮
⎮
.499800000
.500400000
+
1
() − τ + .5 .01 I
2
1
() − τ − .5 .01 I
2
τ
+
> A:=evalf(
1/2*sum( Fi(t[k])*(Integ(t[k+1])-Integ(t[k])),k=0..j-2)
+1/2*sum(Fi(t[k])*(Integ(t[k+1])-Integ(t[k])),k=j+2..N-1)
+1/2* Fi( t_1[j] )*(Integ(t[j+2])-Integ(t[j-1])) ) ;
:=
A
-.173478876
Вычисление второго интеграла
> restart;
> z:=tau->exp(-tau^2)*sin(tau);
:= z → τ e
()−τ
2
()sin τ
> plot(z(tau),tau=-2..2);
Продолжение 
                                        t
                                         k+1                                                ⎞
                                                                                           ⎟
                            9 ⌠                         1                    1             ⎟
− .00002241296235( tk − 1. ) ) ⎮
                               ⎮                                   +                    dτ ⎟
                               ⎮
                               ⎮               ( τ − .5 + .01 I ) 2 ( τ − .5 − .01 I ) 2 ⎟⎟
                               ⌡t                                                          ⎟⎟
                                        k                                                   ⎠
                           .500400000
              ⌠                    1                      1
+ .2660992120 ⎮
              ⎮                                +                      dτ
              ⎮
              ⎮           ( τ − .5 + .01 I ) 2
                                                 ( τ − .5 − .01 I ) 2
              ⌡.499800000
> A:=evalf(
1/2*sum( Fi(t[k])*(Integ(t[k+1])-Integ(t[k])),k=0..j-2)
+1/2*sum(Fi(t[k])*(Integ(t[k+1])-Integ(t[k])),k=j+2..N-1)
+1/2* Fi( t_1[j] )*(Integ(t[j+2])-Integ(t[j-1])) ) ;
A := -.173478876




                              Вычисление второго интеграла

> restart;
> z:=tau->exp(-tau^2)*sin(tau);
             ( −τ2 )
z := τ → e             sin( τ )
> plot(z(tau),tau=-2..2);



                                                                            Продолжение




                                               23

Страницы