Listing of file='19' on disk='vmedia/glbr_22.wvd.zip'
# Sector 40, program filename = '19'
0001 DIM A(10,14),X(13)
0005 PRINT "NO. OF VARIABLES?"
: INPUT M
0010 PRINT "NO. OF CONSTRAINTS?"
: INPUT N
0015 PRINT "ENTER MATRIX A"
0020 FOR I=2 TO N+1
0025 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)
0030 A(I,N+M+1)=A(I,M+2)
: A(I,M+2)=0
: IF I=2 THEN 40
0035 A(I,M+I-1)=A(I,M+1)
: A(I,M+1)=0
0040 NEXT I
0045 PRINT "ENTER OBJECTIVE FUNCTION"
0050 INPUT A(1,1),A(1,2),A(1,3),A(1,4),A(1,5),A(1,6)
0055 PRINT
: R=1
0060 FOR I=1 TO M
: X(I)=1
: NEXT I
0065 FOR I=2 TO N+1
0070 IF A(I,M+I-1)<>-1 THEN 85
: X(M+I-1)=1
0075 FOR J=1 TO N+M
: A(N+2,J)=A(N+2,J)-A(I,J)
: NEXT J
0080 R=N+2
0085 NEXT I
0090 S,T=1
0095 FOR I=2 TO N+M
: IF A(R,I)<A(R,S) THEN 100
: S=I
0100 IF A(R,I)>=A(R,T) THEN 110
: T=I
0110 NEXT I
0120 IF A(R,T)<0 THEN 140
: IF R=1 THEN 300
0130 IF A(R,S)>1E-4 THEN 280
: R=1
: GOTO 90
0140 S=1
0150 FOR I=2 TO N+1
: IF A(I,T)<=0 THEN 190
0160 Y=A(I,N+M+1)/A(I,T)
: IF S=1 THEN 180
0170 IF Y>=A(S,N+M+1)/A(S,T) THEN 190
0180 S=I
0190 NEXT I
0200 IF S=1 THEN 290
0210 FOR I=1 TO N+M
: IF X(I)=1 THEN 220
: IF A(S,I)=1 THEN 230
0220 NEXT I
0230 X(I)=1
: X(T)=0
: Y=A(S,T)
0240 FOR I=1 TO N+M+1
: A(S,I)=A(S,I)/Y
: NEXT I
0250 FOR I=1 TO N+2
: IF I=S THEN 270
: Y=A(I,T)
0260 FOR J=1 TO N+M+1
: A(I,J)=A(I,J)-Y*A(S,J)
: NEXT J
0270 NEXT I
: GOTO 90
0280 PRINT "INFEASIBLE"
: STOP
0290 PRINT "UNBOUNDED"
: STOP
0300 FOR J=1 TO M
0310 IF X(J)=0 THEN 320
: X(J)=0
: GOTO 340
0320 FOR I=2TO N+1
: IF A(I,J)=1 THEN 330
: NEXT I
0330 X(J)=A(I,N+M+1)
0340 NEXT J
0350 Y=A(1,N+M+1)
0360 PRINT
: PRINT "OBJ. FUNC.=";Y
: PRINT
0370 FOR I=1TO M
: PRINT "X(";I;")=";X(I)
: NEXT I
0380 END