1000 REM ***** SIMEQN *** MATHEMATICS PROGRAM*** 1010 REM ***** VERSION 1**UPDATED 10/1/70 1020 REM REVISED MAY 1972 ANDY ROTH C.C.G. 1030 REM SOLVES SIMULTANEOUS LINEAR EQUATIONS 1040 DIM A(25,25),B(25,25),X(25,1),Y(1,25),Z(25,1) 1050 PRINT"LINEAR SYSTEM SOLVER" 1060 INPUT"DO YOU WANT INSTRUCTIONS";A$ 1070 IF A$<>"YES" THEN 5000 1080 PRINT 1090 PRINT" THIS PROGRAM SOLVES SYSTEMS CONSISTING OF N LINEAR" 1100 PRINT"EQUATIONS IN N UNKNOWNS. TO USE, COEFFICIENTS IN THE" 1110 PRINT"EQUATIONS ARE ENTERED WHEN CALLED FOR. STARTING WITH" 1120 PRINT"THE FIRST COEFFICIENT OF THE FIRST EQUATION AND ENDING" 1130 PRINT"WITH THE N-TH COEFFICIENT OF THE N-TH EQUATION. ALL ZERO" 1140 PRINT"COEFFICIENTS MUST BE ENTERED IN THEIR PROPER PLACE" 1150 PRINT"THE RIGHT SIDE CONSTANT TERMS OF THE EQUATIONS ARE THEN" 1160 PRINT"ENTERED WHEN CALLED FOR. IF ADDITIONAL CASES WITH THE" 1170 PRINT"SAME COEFFICIENT MATRIX BUT DIFFERENT RIGHT SIDES" 1180 PRINT"ARE TO BE RUN, THEY MAY BE ENTERED WHEN CALLED FOR" 1190 PRINT" FOR EXAMPLE:" 1200 PRINT" ------------" 1210 PRINT" 3X+5Y-2Z=9 3X+5Y-2Z=19" 1220 PRINT" 7X+ Y =-3 7X+ Y =-3" 1230 PRINT" X-7Y+9Z=14 X-7Y+9Z= 8" 1240 PRINT"WOULD BE ENTERED IN THE FOLLOWING FASHION" 1250 PRINT"# SYSTEMS,# VARIABLES: ? 2,3" 1260 PRINT"COEFFICIENTS: ? 3,5,-2,7,1,0,1,-7,9" 1270 PRINT"RIGHT SIDE FOR CASE 1 ? 9,-3,14" 1280 PRINT"RIGHT SIDE FOR CASE 2 ? 19,-3,8" 1290PRINT 5000 PRINT 5010 INPUT"# SYSTEMS, # VARIABLES";L,M 5020 PRINT"COEFFICIENTS:"; 5030 MAT INPUT A(M,M) 5040 PRINT 5050 PRINT"SOLUTION FOR LINEAR SYSTEM OF ORDER";M 5060 PRINT"COEFFICIENT MATRIX" 5070 MAT PRINT A, 5080 PRINT 5090 MAT B=ZER(M,M) 5100 MAT B=INV(A) 5110 MAT Y=ZER(1,M) 5120 MAT Z=ZER(M,1) 5130 Q=0 5140 Q=Q+1 5150 PRINT"RIGHT SIDES FOR CASE";Q; 5160 MAT INPUT X(M,1) 5170 PRINT 5180 PRINT"DEPENDENT VARIABLE VECTOR FOR CASE";Q 5190 PRINT 5200 MAT Y=TRN(X) 5210 MAT PRINT Y, 5220 MAT Z=B*X 5230 MAT Y=TRN(Z) 5240 PRINT 5250 PRINT"SOLUTION VECTOR FOR CASE";Q 5260 PRINT TAB(17);"VARIABLES" 5270 FOR K=1 TO M 5280 PRINT" X";K, 5290 NEXT K 5300 PRINT 5310 MAT PRINT Y, 5320 PRINT"PROOF OF SOLUTION FOR CASE";Q 5330 PRINT 5340 MAT X= A*Z 5350 MAT Y=TRN(X) 5360 MAT PRINT Y, 5370 IF Q