Support our Blind Hand Charity by donating at 9818291723. Your contribution has the potential to bring about significant positive change for humanity.

Search This Blog

Wednesday, September 27, 2023

CRAMMER'S RULE

print("Enter Coefficients in ax+by+cz=d format\n")

l=[[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0]]

for i in range(0,3):

    for j in range(0,4):

        x=int(input("Enter the value = "))

        l[j][i]=x


d=[[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0]]

for i in range(0,3):

    for j in range(0,3):

        d[j][i]=l[j][i]


dx=[[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0]]

for i in range(0,3):

    for j in range(0,3):

        if i==0:

            i1=3

        else:

            i1=i

        dx[j][i]=l[i1][j]


dy=[[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0]]

for i in range(0,3):

    for j in range(0,3):

        if i==1:

            i1=3

        else:

            i1=i

        dy[j][i]=l[i1][j]


dz=[[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0]]

for i in range(0,3):

    for j in range(0,3):

        if i==2:

            i1=3

        else:

            i1=i

        dz[j][i]=l[i1][j]


c=0

m=[d,dx,dy,dz]

for k in range(0,4):

    p=[]

    for i in range(0,3):

        for j in range(0,3):

            if j!=i and i!=2:

                x=m[k][1][j]*m[k][2][j+i+1]

                y=m[k][1][j+1+i]*m[k][2][j]

                if i%2==0:

                    p.append(m[k][0][i]*(x-y))

                else:

                    p.append(m[k][0][i]*(y-x))

            if j!=i and i==2 and j!=1:

                x=m[k][1][j]*m[k][2][j+i-1]

                y=m[k][1][j-1+i]*m[k][2][j]

                p.append(m[k][0][i]*(x-y))

    c1=0

    if k==0:

        for i in range(0,len(p)):

            c+=p[i]

    else:

        for i in range(0,len(p)):

            c1+=p[i]

    if k==1:

        print('X =',c1/c)

    elif k==2:

        print('Y =',c1/c)

    elif k==3:

        print('Z =',c1/c)

OUTPUT