Документ взят из кэша поисковой машины. Адрес оригинального документа : http://foroff.phys.msu.ru/illposed/programs/Fig4_5.txt Дата изменения: Mon Jul 7 00:04:43 2008 Дата индексирования: Mon Oct 1 22:04:18 2012 Кодировка:

Fig 4.5
=======


00001 IMPLICIT REAL*8(A-H,O-Z)
00002 IMPLICIT INTEGER*4(I-N)
00003 REAL*8 KERN,L1,L2
00004 DIMENSION U(400),U0(64),Z0(64),Z(64),
*KERN(64)
00005 EXTERNAL AK
00006 C=0.0
00007 D=2.0
00008 L1=0.0
00009 L2=1.0
00010 ALPHA=1.0
00011 C1=1.000001
00012 N=64
00013 IMAX=500
00014 ST=(D-C)/N
00015 DO 5 I=1,N
00016 S=(I-0.5)*ST-0.5
00017 Z0(I)=0.
00018 IF(S.LT.0.0.OR.S.GT.1.0)GOTO 5
00019 Z0(I)=((DEXP(-(S-0.3)**2/0.03)+
+DEXP(-(S-0.7)**2/0.03))/0.955040800-
-0.052130913)*1.4*S
00020 5 CONTINUE
00021 DO 554 K=1,N
00022 S=0.5*(L1+L2)+ST*(K-N/2-1)
00023 KERN(K)=0.
00024 IF(S.LT.L1.OR.S.GT.L2)GOTO 554
00025 KERN(K)=AK(S)
00026 554 CONTINUE
00027 HH=1.E-9
00028 DO 771 I=1,N
00029 S=0.0
00030 DO 770 J=1,N
00031 IND=I-J+N/2+1
00032 IF(IND.LT.1.OR.IND.GT.N)GOTO 770
00033 S=S+KERN(IND)*Z0(J)
00034 770 CONTINUE
00035 U0(I)=S*ST
00036 771 CONTINUE
00037 DL=1.E-8
00038 PRINT 999,(Z0(I),I=1,N)
00039 999 FORMAT(' Exact solution:'/(5F11.8))
00040 PRINT 553,DL,HH
00041 553 FORMAT(' Errors:'/
*' in right side:',D16.9/
*' in operator :',D16.9)
00042 CALL PTIKR(AK,U0,A,B,C,D,L1,L2,N,Z,AN,
*DL,HH,C1,IMAX,ALPHA,U,400,IERR)
00043 PRINT 501,IERR,AN,ALPHA,
*A,B,(Z(II),II=1,N)
00044 501 FORMAT(' .'/
*' Return code (IERR) :',
* I5/
*' Discrepancy :',
* D16.9/
*' Regularization parameter :',
* D16.9/
*' Solution found on interval : (',
* F4.1,',',F3.1,')'/
*' .'/' Approximate solution:'/
* (5F11.8))
00045 STOP
00046 END


00001 FUNCTION AK(Y)
00002 IMPLICIT REAL*8(A-H,O-Z)
00003 AK=DEXP(-80.*(Y-.5)**2)
00004 RETURN
00005 END





Exact solution:
.00000000 .00000000 .00000000 .00000000 .00000000
.00000000 .00000000 .00000000 .00000000 .00000000
.00000000 .00000000 .00000000 .00000000 .00000000
.00000000 .00040570 .00469802 .01649271 .03976953
.07814444 .13325059 .20303119 .28077126 .35567888
.41531472 .44932413 .45318938 .43054079 .39307592
.35804191 .34403504 .36623085 .43211954 .53865959
.67164363 .80785932 .92004465 .98364014 .98334749
.91718612 .79650698 .64202556 .47754936 .32383759
.19461398 .09552648 .02555887 .00000000 .00000000
.00000000 .00000000 .00000000 .00000000 .00000000
.00000000 .00000000 .00000000 .00000000 .00000000
.00000000 .00000000 .00000000 .00000000
Errors:
in right side: .999999994D-08
in operator : .999999972D-09
.
Return code (IERR) : 0
Discrepancy : .344029501D-07
Regularization parameter : .573252726D-06
Solution found on interval : ( -.5,1.5)
.
Approximate solution:
-.00017321 -.00009008 .00000358 .00006029 .00008013
.00007742 .00004968 -.00001771 -.00011536 -.00018162
-.00012615 .00008648 .00034550 .00039636 .00003549
-.00048053 .00007864 .00429205 .01602516 .03971215
.07899454 .13503747 .20511177 .28210682 .35544513
.41342703 .44649610 .45048758 .42878790 .39255440
.35866264 .34570400 .36896886 .43575007 .54236307
.67392962 .80724469 .91610707 .97765767 .97787558
.91464292 .79762609 .64528134 .48012463 .32370988
.19255030 .09573080 .03446567 .00344908 -.00649811
-.00528917 -.00088940 .00212787 .00252328 .00120195
-.00032238 -.00105252 -.00086547 -.00023138 .00029757
.00044660 .00027443 .00000811 -.00015972