1 DIM A(10,14),X(13)
5 PRINT "NO. OF VARIABLES?": INPUT M
10 PRINT "NO. OF CONSTRAINTS?": INPUT N 
15 PRINT "ENTER MATRIX A" 
20 FOR I=2 TO N+1 
25 INPUT A(I,1),A(I,2),A(I,3),A(I,4),A(I,5),A(I,6),A(I,7),A(I,8),            A(I,9) 
30 A(I,N+M+1)=A(I,M+2): A(I,M+2)=0: IF I=2 THEN 40 
35 A(I,M+I-1)=A(I,M+1): A(I,M+1)=0 
40 NEXT I 
45 PRINT "ENTER OBJECTIVE FUNCTION" 
50 INPUT A(1,1),A(1,2),A(1,3),A(1,4),A(1,5),A(1,6)
55 PRINT : R=1
60 FOR I=1 TO M: X(I)=1: NEXT I
65 FOR I=2 TO N+1 
70 IF A(I,M+I-1)<>-1 THEN 85: X(M+I-1)=1 
75 FOR J=1 TO N+M: A(N+2,J)=A(N+2,J)-A(I,J): NEXT J
80 R=N+2 
85 NEXT I 
90 S,T=1 
95 FOR I=2 TO N+M: IF A(R,I)<A(R,S) THEN 100: S=I 
100 IF A(R,I)>=A(R,T) THEN 110: T=I 
110 NEXT I  
120 IF A(R,T)<0 THEN 140: IF R=1 THEN 300
130 IF A(R,S)>1E-4 THEN 280: R=1: GOTO 90
140 S=1 
150 FOR I=2 TO N+1: IF A(I,T)<=0 THEN 190 
160 Y=A(I,N+M+1)/A(I,T): IF S=1 THEN 180 
170 IF Y>=A(S,N+M+1)/A(S,T) THEN 190  
180 S=I
190 NEXT I  
200 IF S=1 THEN 290
210 FOR I=1 TO N+M: IF X(I)=1 THEN 220: IF A(S,I)=1 THEN 230  
220 NEXT I  
230 X(I)=1: X(T)=0: Y=A(S,T) 
240 FOR I=1 TO N+M+1: A(S,I)=A(S,I)/Y: NEXT I
250 FOR I=1 TO N+2: IF I=S THEN 270: Y=A(I,T)
260 FOR J=1 TO N+M+1: A(I,J)=A(I,J)-Y*A(S,J): NEXT J
270 NEXT I: GOTO 90
280 PRINT "INFEASIBLE": STOP 
290 PRINT "UNBOUNDED": STOP 
300 FOR J=1 TO M  
310 IF X(J)=0 THEN 320: X(J)=0: GOTO 340 
320 FOR I=2TO N+1: IF A(I,J)=1 THEN 330:NEXT I
330 X(J)=A(I,N+M+1)
340 NEXT J
350 Y=A(1,N+M+1)
360 PRINT :PRINT "OBJ. FUNC.=";Y:PRINT 
370 FOR I=1TO M: PRINT "X(";I;")=";X(I): NEXT I
380 END 

