CLS
DIM A(10, 10)
DIM B(10, 10)
DIM P(10)
DIM C1(10)
PRINT "CHARACTERISTIC EQUATION OF A MATRIX"
PRINT
PRINT "F(LAMBDA) = (-1)^N * (L^N-P(1)*L^(N-1)-P(2)*L^(N-2) ... -P(N))"
PRINT "where N is the order"
PRINT
INPUT "KEY IN ORDER:    "; N
PRINT "KEY IN ELEMENTS AS INDICATED:"
PRINT
    FOR I = 1 TO N
        FOR J = 1 TO N
            PRINT "A ( " ; I ; " , " ; J ; ") = ";
            INPUT A(I,J)
            B(I,J) = A(I,J)
        NEXT J
        PRINT
    NEXT I
M = N -1
FOR K = 1 TO M
    T1 = 0
    FOR I= 1 TO N
        T1=T1+B(I,I)
    NEXT I
    A1 = K
    P(K) = T1 / A1
    FOR I = 1 TO N
        B(I, I) = B(I, I) -P(K)
    NEXT I
    FOR J = 1 TO N
        FOR I= 1 TO N
            C1(I)= B(I,J)
        NEXT I
        FOR I = 1 TO N
            B(I,J) = 0
            FOR L = 1 TO N
                B(I,J) = B(I,J) + A(I,L) * C1(L)
            NEXT L
        NEXT I
    NEXT J
NEXT K
P(N) = B(N,N)
PRINT "COEFFICIENTS ARE: "
PRINT
FOR K = 1 TO N
    PRINT "P(";K;") = ";P(K)
NEXT K
END