home *** CD-ROM | disk | FTP | other *** search
GW-BASIC | 1994-02-04 | 1.6 KB | 91 lines |
- 0 CLS:PRINT
- 10 PRINT "PLOT OF POLAR EQUATION"
- 20 PRINT
- 30 DIM X(91),Y(90)
- 40 N = 90
- 50 PRINT "ABSOLUTE VALUE OF ENDPOINTS";
- 60 INPUT Z
- 70 PRINT
- 80 PRINT "INCREMENT OF X-AXIS = ";Z / 30
- 90 PRINT "INCREMENT OF Y-AXIS = ";Z / 18
- 100 PRINT
- 110 FOR I = 1 TO N
- 120 D = 0.0698132 * I
- 130 REM - ENTER FUNCTION HERE (F AS A FUNCTION OF D)
- 140 X(I) = INT (((F * COS (D) / Z + 1) * 30) + 0.5)
- 150 Y(I) = INT ((( - F * SIN (D) / Z + 1) * 18) + 0.5)
- 160 NEXT I
- 170 FOR J = 1 TO N
- 180 FOR I = 1 TO N - J
- 190 A = X(I)
- 200 B = Y(I)
- 210 IF B < = Y(I + 1) THEN 260
- 220 X(I) = X(I + 1)
- 230 Y(I) = Y(I + 1)
- 240 X(I + 1) = A
- 250 Y(I + 1) = B
- 260 NEXT I
- 270 NEXT J
- 280 T = 1
- 290 FOR P = 0 TO N - 1
- 300 IF Y(P + 1) > = 0 THEN 320
- 310 NEXT P
- 320 FOR I = 0 TO 36
- 330 T = T + P
- 340 P = 0
- 350 IF T > N THEN 370
- 360 IF Y(T) = I THEN 420
- 370 IF I = 18 THEN 400
- 380 PRINT TAB( 30);"*";
- 390 GOTO 860
- 400 S = N + 1
- 410 GOTO 740
- 420 FOR L = T TO N
- 430 IF Y(L) > Y(T) THEN 460
- 440 P = P + 1
- 450 NEXT L
- 460 IF P = 1 THEN 560
- 470 FOR J = 1 TO P
- 480 FOR L = 1 TO P - J
- 490 C = X(T + L - 1)
- 500 A = X(T + L)
- 510 IF C < = A + 0 THEN 540
- 520 X(T + L - 1) = A
- 530 X(T + L) = C
- 540 NEXT L
- 550 NEXT J
- 560 IF I = 18 THEN 730
- 570 L = - 1
- 580 S = 0
- 590 FOR K = 0 TO P - 1
- 600 IF X(T + K) = L THEN 690
- 610 L = X(T + K)
- 620 IF L = 30 THEN 660
- 630 IF L < 30 THEN 670
- 640 IF S = 1 THEN 670
- 650 PRINT TAB( 30);"*";
- 660 S = 1
- 670 IF L > 60 THEN 860
- 680 PRINT TAB( L);"+";
- 690 NEXT K
- 700 IF S = 1 THEN 860
- 710 PRINT TAB( 30);"*";
- 720 GOTO 860
- 730 S = T
- 740 FOR J = 0 TO 60
- 750 IF X(S) < > J THEN 830
- 760 PRINT "+";
- 770 FOR K = S TO T + P - 1
- 780 IF X(K) = X(S) THEN 810
- 790 S = K
- 800 GOTO 840
- 810 NEXT K
- 820 GOTO 840
- 830 PRINT "*";
- 840 NEXT J
- 850 PRINT "X";
- 860 PRINT
- 870 NEXT I
- 880 PRINT TAB( 30);"Y"
- 890 RUN "MATH"
-