550 REM OBLICZENIE ZANURZENIA 560 I = 1: SM 1 = SM* 1.025/GAMA 570 1 = 1 + 1
580 IF D(I)<SM1 THEN GOTO 570 590 A = 2*D(I -1) - 4*D(I)+2*D(I +1)
600 B=4*D(I)-3*D(I-1)-D(1 + 1)
610 DELTA = B*B-4*A*(D(I-1)-SM l)
620 DELTA = SQRT(DELTA)
630 DT = (—B + DELTA)/2/A 640 T=2+0.5*I + DT
650 U = F(I— 1): V = F(I): W=F(I +1): GOSUB 1000
660 XF=Y: U = S(1 — 1): V = S(I): W = S(I +1): GOSUB 1000
670 XS = Y: U = K(I— 1): V= K(I): W = K(I +1):GOSUB 1000
680 ZM = Y: U = M(I— 1): V = M(I): W = M(I +1): GOSUB 1000
690 MJ = Y
700 GM=ZM —ZG1
710 TRYM = SM*(XG — XF)/MJ
720 TTR = -TRYM*XS/LPP: TTD = TRYM +TTR
730 TD = T+TTD:TR=T + TTR
740 I = INT((SM 1 -10000)/2500 +1): DEL=(SM 1 - 7500 - 2500*I)/2500 750 ZGK = R(I) + DEL*(R(I +1)- R(I))
770 PRINT " - WYNIKI OBLICZEŃ 780 PRINT
790 PRINT "WYPORNOŚĆ = GM = ”;GM
800 PRINT "POŁOŻENIE ZG = ”;ZG1 "KRYTYCZNE POŁOŻENIE ZG = ”;ZGK 810 PRINT "ZANURZENIE ŚREDNIE ODPOWIADAJĄCE RÓWNEJ STEPCE=” 820 PRINT "ZANURZENIE NA DZIOBIE = ”;TD 830 PRINT "ZANURZENIE NA RUFIE = ”;TR 840 STOP
1000 REM PODPROGRAM INTERPOLACJI KRZYWOLINIOWEJ
1010 A = 2*U-4*V + 2*W
1020 B = 4*V —3*U —W
1030 Y = A*DT*DT+B*DT+U
1040 RETURN
1050 END
madę by voyteck (dec'2004)
250