DO 100 I=1,N
DO 100 J=1,N
YY(1)=X(J)-X(I)
YY(2)=Y(J)-Y(I)
R=DSQRT(YY(1)**2+YY(2)**2)
DO 50 L1=1,2
I1=2*(I-1)+L1
DO 50 L2=1,2
I2=2*(J-1)+L2
I3=NN*(I2-1)+I1
DELT=0.0
IF(L1==L2) DELT=1.0
XK(I3)=-RO*((C1*R-C2*CAPR)*DELT*R*R-C3*R*YY(L1)*YY(L2))/GE
50 CONTINUE
100 CONTINUE
RETURN
END