Καλώς ήλθες στο φόρουμ.
Για το ερώτημά σου, δεν καταλαβαίνω πολλά πράγματα από αυτά που γράφεις. Για παράδειγμα
1) Γράφεις το σύστημα στην μορφή
![[u_1,..., u_4]= [u_1,..., u_4]A [u_1,..., u_4]= [u_1,..., u_4]A](/forum/ext/geomar/texintegr/latexrender/pictures/11f62066437a1577ee5a5f02f08c4cfd.png)
ενώ ο ορθόδοξος τρόπος είναι
![[u_1,..., u_4]^T = A[u_1,..., u_4]^T [u_1,..., u_4]^T = A[u_1,..., u_4]^T](/forum/ext/geomar/texintegr/latexrender/pictures/3cae91cf3d41dab914fad4c976451b1f.png)
. Το κρατάμε.
2) Από κάτω γράφεις τις εξισώσεις αναλυτικά. Όμως αυτό που γράφεις είναι συμβατό με αυτό που αναφέρω στο α) ως ορθόδοξη γραφή, και όχι με αυτό που γράφεις ο ίδιος.
3) Αναφέρεις ότι η μοναδική λύση είναι η

και λοιπά. Να όμως που αν αντικατστήσουμε τις τιμές αυτές στο σύστημα είτε με την δική σου γραφή, είτε με την δική μου, οι εξισώσεις
δεν ικανοποιούνται (ικανοποιούν μόνο την

, που όμως δεν δίνεται από το σύστημα). Οπότε για τι λύση μιλάς;
Για να συνοψίσω: Αν δεν αποφασίσεις ποιος είναι ο σωστός συμβολισμός, θα προσκρούσεις σε προβλήματα. Η μέθοδος Cramer και Gauss που πας να εφαρμόσεις απαιτούν, εννοείται, τα σωστά σύμβολα.
Τέλος, επειδή η εξίσωση

που γράφεις, δεν είναι μία από τις τέσσερις του συστήματος, αλλά μία προσθήκη, τελικά σημαίνει ότι έχουμε ένα σύστημα 5 εξισώσεων με 4 άγνωστους. Αλλά τότε δεν είναι βέβαιο ότι οι εξισώσεις είναι συμβατές. Γι' αυτό βγάζεις αντιφατικά αποτελέσματα.
Παρακαλώ κάνε άλλη μία προσπάθεια, αλλά με σωστούς συμβολισμούς. Αν υπάρξει ανάγκη, μην διστάσεις να επανέλθεις.