. , , ,

,,,

,

:

.

403

: 頠 . 04-109 /____________/

: .. /____________/

1999.


...

4.24

: = S * cos(x) + q * sin(x),

堠 q - ,

S - .

a=5

EPS.

N ,

X=Xn Dx.


.

:

k=1, 2, ... , n.

,

- .

(k=1, 2, ..., 2n).

. h n, , , h,

( ) :




X(100), Y(100)


: a, q, EXP, Dx, XN, N, ZN, ZK


J = 1

X(J) = XN


XJ = X(J)


S = INTEGR( a, XJ, EPS, ZN, ZK)


Y(J) = S*cos( X(J) )+q*sin( X(J) )

J = J + 1

X(J) = X(J - 1) + Dx


J <= N


: ( X(J), J=1, N ), ( Y(J), J=1, N )


1. X, Y

2. : a, q, EPS, Dx, XN, N, ZN, ZK

3. J, X(J).

4. XJ.

5. S=INTEGR(a, XJ, EPS, ZN, ZK)

6. Y(J), J, X(J).

7. J.

8. ( X(J), J=1, N ), ( Y(J), J=1, N ).



Ϡ S = INTEGR( a, XJ, EPS, ZN, ZK )


I1 = 1

K = 1


I2 = 0

H = ( ZK - ZN ) / K

I = 2


Z2 = ZN + I*H, Z1 = Z2 - H, Z0 = Z1 - H

L2 = ln( XJ + a*Z2 ),

L1 = ln( XJ + a*Z1 ),

L0 = ln( XJ + a*Z0 ),

I2 = I2 + L0 + 4*L1 + L2


I<=K


I = I + 2

| I1 - I2 | < EPS


I1 = I2

K = 2*K INTEGR = I2


INTEGR .

:

a - , .

XJ - = S * cos(x) + q * sin(x), - .

EPS - , .

ZN - , .

ZK - , .

1. I1, K.

2. I2, H, I.

3. Z2, L2, L1, L0, I2 - .

4. I.

5. I.

6. | I1 - I2 | < EPS.

7. I1, K.

8. INTEGR.

: . 403 : 頠 . 04-109 /____________/

 

 

 

! , , , .
. , :