Matrix in a program

10262020, 10:21 PM
Post: #1




Matrix in a program
Hello
I have no practice in programming and when I tried to develop some algebra with the HP prime to find implicit equations the program fails: //Implicit equations by Gaus EXPORT ImpGa(Ma1) BEGIN //ddimension , ccolumn dimension //nactual column , rrow search LOCAL d,c,n,r,Ma2; c:=colDim(Ma1); n:=1; FOR n FROM n=1 TO c DO //Search for 1 in col n r:=col(Ma1,n); // RETURN r; r:=CAS.POS(r,1); IF r==0 THEN r:=n; ELSE IF r≠n THEN SWAPROW(Ma1,r,n); END; END; Ma1:=pivot(Ma1,n,n); END; RETURN Ma1; END; When I tried to run it it gives me "Error: Invalid input". I put in the command line: ImpGa([[1 2 x1] [2 1 x2] [0 1 x3] [1 2 x4] [3 1 x5]] If I delete the columna [x1 x2 x3 x4 x5] it works. Seems that does not accept variables in the program. Doing the same operations manually in CAS mode it works. Can you give some clue about how to introduce this kind of matrix in this program? Thanks very much Toni Garçon 

10272020, 01:03 AM
(This post was last modified: 10272020 01:04 AM by Major.)
Post: #2




RE: Matrix in a program
If you are wanting to use variables then I think it needs to be a CAS program instead of a Home program.
Just add #cas at the very top and #end at the very bottom of the program. Change "EXPORT ImpGa(Ma1)" to "ImpGa(Ma1):=" And then change CAS.POS to just POS. Control System Programs 

10272020, 02:28 AM
Post: #3




RE: Matrix in a program
(10262020 10:21 PM)Tonig00 Wrote: When I tried to run it it gives me "Error: Invalid input". It does not work even with variables deleted. It just seems it did. Anyway, you should not do pivot on the variables. It make no sense. Some Issues: 1. SWAPROW(Ma1,r,n) does nothing. You have to save the returns. 2. pivot on 1 is preferred. But more importantly, pivot must be nonzero. 3. pivot row/col can only be used once. see Gauss Jordan Elimination Through Pivoting This is XCas implementation of revised code. note: XCas is setup 0based, but POS returns 1based (0 signals not found) Code: ImpGa(M, n) := { // pivot upto n columns XCas> M := [[1,2,x1],[2,1,x2],[0,1,x3],[1,2,x4],[3,1,x5]] XCas> ImpGa(M, 2) \(\left(\begin{array}{ccc} 1 & 0 & \mathrm{x1}2\cdot \mathrm{x3} \\ 0 & 1 & \mathrm{x3} \\ 0 & 0 & 2\cdot \mathrm{x1}+\mathrm{x2}+3\cdot \mathrm{x3} \\ 0 & 0 & \mathrm{x1}+\mathrm{x4} \\ 0 & 0 & 3\cdot \mathrm{x1}+7\cdot \mathrm{x3}+\mathrm{x5} \end{array}\right) \) (10272020 01:03 AM)Major Wrote: If you are wanting to use variables then I think it needs to be a CAS program instead of a Home program. Another way is replacing variables with an identity matrix. XCas> M := [[1,2],[2,1],[0,1],[1,2],[3,1]] XCas> ImpGa(augment(M, identity(5)), 2) \(\left(\begin{array}{ccccccc} 1 & 0 & 1 & 0 & 2 & 0 & 0 \\ 0 & 1 & 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 2 & 1 & 3 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 & 1 & 0 \\ 0 & 0 & 3 & 0 & 7 & 0 & 1 \end{array}\right) \) 

10272020, 07:37 PM
Post: #4




RE: Matrix in a program
Hello
Thanks very much for your help. First I tried with the system #cas / #end, and ImpGa(matrix) gives the same as in XCas: [[1,0,x12*x3],[0,1,x3],[0,0,2*x1+x2+3*x3],[0,0,x1+x4],[0,0,3*x3+x5]] Following also the suggested corrections (not including variable column which really has not sense at all): #cas //Implicit ecuations by Gaus ImpGa(Ma1):= BEGIN //ddimension , ccolumn dimension //nactual column , rrow search LOCAL d,c,n,r,Ma2; c:=colDim(Ma1); n:=1; FOR n FROM 1 TO c1 DO //Search for 1 in col n r:=col(Ma1,n); // RETURN r; r:=POS(r,1); IF r==0 THEN r:=n; ELSE IF r≠n THEN SWAPROW(Ma1,r,n); END; END; Ma1:=pivot(Ma1,n,n); END; RETURN Ma1; END; #end I do not know if this is my interpretation, but the ImpGa() function has disappear from the "User" menu and has gone to the "Catlg" menu. Possibly all CAS functions go to general catalog. This raises the question if it is possible to group CAS functions. I have to practice much more. Thanks again to all Toni 

10272020, 07:59 PM
Post: #5




RE: Matrix in a program
Sorry it does not work with any matrix, so I have still errors.
I will correct and send it again. Toni 

« Next Oldest  Next Newest »

User(s) browsing this thread: 1 Guest(s)