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

UptoLike

24
> T:=0.7;
delta:=0.5;
h:=evalf(delta/5);
p:=4;
N:=(h)^(-p);
:= T .7
:= δ .5
:= h .1000000000
:= p 4
:=
N
10000.00000
> unassign('tau');
> simplify((((tau-T)-I*h)^p+((tau-T)+I*h)^p)/
((tau-T)^2+h^2)^p);
8.
+
+
+ 7000.
τ
3
7200. τ
2
3220.
τ
527. 2500.
τ
4
() + 10. τ
2
14. τ 5.
4
> Ker:=tau->8.*(-7000.*tau^3+7200.*tau^2-3220.*tau+527.+
+2500.*tau^4)/(10.*tau^2-14.*tau+5.)^4:
Продолжение
> T:=0.7;
delta:=0.5;
h:=evalf(delta/5);
p:=4;
N:=(h)^(-p);
T := .7
δ := .5
h := .1000000000
p := 4
N := 10000.00000
> unassign('tau');
> simplify((((tau-T)-I*h)^p+((tau-T)+I*h)^p)/
  ((tau-T)^2+h^2)^p);
     −7000. τ 3 + 7200. τ 2 − 3220. τ + 527. + 2500. τ 4
8.                                            4
                   ( 10. τ 2 − 14. τ + 5. )
> Ker:=tau->8.*(-7000.*tau^3+7200.*tau^2-3220.*tau+527.+
+2500.*tau^4)/(10.*tau^2-14.*tau+5.)^4:


                                                           Продолжение


                                              24