. , , ,

,,,

,


:

.


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 /____________/ : .. /____________/ 1999. ...

 

 

 

! , , , .
. , :