Tuesday 12 August 2014

matrix multiplication

#include <stdio.h>
void main()
{
int m[3][3], n[3][3], i, j, h[9];
printf("please enter the elements of matrix A\n");
for (i = 0; i < 3; i++)
for (j = 0; j < 3; j++)
{
printf("enter the element a(%d,%d)", i + 1, j + 1);
scanf("%d", &m[i][j]);
}
printf("entered matrix\n");
for (i=0;i<3;i++)
{
for (j=0;j<3;j++)
{
printf("%d\t",m[i][j]);
}
printf("\n");
}

printf("enter the elements of matrix B \n");
for (i = 0; i < 3; i++)
for (j = 0; j < 3; j++)
{
printf("enter the elements of b(%d,%d)", i + 1, j + 1);
scanf("%d", &n[i][j]);
}
printf("entered matrix\n");
for (i=0;i<3;i++)
{
for (j=0;j<3;j++)
{
printf("%d\t",n[i][j]);
}
printf("\n");
}
printf("\n\n\n");
h[0]=m[0][0]*n[0][0]+m[0][1]*n[1][0]+m[0][2]*n[2][0];
h[1]=m[0][0]*n[0][1]+m[0][1]*n[1][1]+m[0][2]*n[2][1];
h[2]=m[0][0]*n[0][2]+m[0][1]*n[1][2]+m[0][2]*n[2][2];
h[3]=m[1][0]*n[0][0]+m[1][1]*n[1][0]+m[1][2]*n[2][0];
h[4]=m[1][0]*n[0][1]+m[1][1]*n[1][1]+m[1][2]*n[2][1];
h[5]=m[1][0]*n[0][2]+m[1][1]*n[1][2]+m[1][2]*n[2][2];
h[6]=m[2][0]*n[0][0]+m[2][1]*n[1][0]+m[2][2]*n[2][0];
h[7]=m[2][0]*n[0][1]+m[2][1]*n[1][1]+m[2][2]*n[2][1];
h[8]=m[2][0]*n[0][2]+m[2][1]*n[1][2]+m[2][2]*n[2][2];
printf("the multiplication is…………\n");
printf("%d\t%d\t%d\n ",h[0],h[1],h[2]);
printf("%d\t%d\t%d\n ",h[3],h[4],h[5]);
printf("%d\t%d\t%d  ",h[6],h[7],h[8]);
}

Wednesday 6 August 2014

solving 3 variable equations

#include<stdio.h>
void main()
{
    int i,j,a[3][4];
    float det,detx,dety,detz,x,y,z;
    printf("a1x+b1y+cz=d1\na2x+b2y+c2z=d2\na3x+b3y+c3z=d3\n");
    printf("enter value of a1,b1,c1,d1,a2,b2,c2,d2,a3,b3,c3,d3");
    for(i=0;i<3;i++)
    {
        for(j=0;j<4;j++)
        {
            scanf("%d",&a[i][j]);
        }
    }
det=a[0][0]*a[1][1]*a[2][2]-a[0][0]*a[2][1]*a[1][2]-a[0][1]*a[1][0]*a[2][2]+a[0][1]*a[2][0]*a[1][2]+a[0][2]*a[1][0]*a[2][1]-a[0][2]*a[2][0]*a[1][1];
detx=a[0][3]*a[1][1]*a[2][2]-a[0][3]*a[2][1]*a[1][2]-a[0][1]*a[1][3]*a[2][2]+a[0][1]*a[2][3]*a[1][2]+a[0][2]*a[1][3]*a[2][1]-a[0][2]*a[2][3]*a[1][1];
dety=a[0][0]*a[1][3]*a[2][2]-a[0][0]*a[2][3]*a[1][2]-a[0][3]*a[1][0]*a[2][2]+a[0][3]*a[2][0]*a[1][2]+a[0][2]*a[1][0]*a[2][3]-a[0][2]*a[2][0]*a[1][3];
detz=a[0][0]*a[1][1]*a[2][3]-a[0][0]*a[2][1]*a[1][3]-a[0][1]*a[1][0]*a[2][3]+a[0][1]*a[2][0]*a[1][3]+a[0][3]*a[1][0]*a[2][1]-a[0][3]*a[2][0]*a[1][1];
if(det==0&&detx==0&&dety==0&&detz==0)
{
    printf("there are infinite solution");
}
else if(det!=0)
{
    x=detx/det;
    y=dety/det;
    z=detz/det;
    printf("there is a unique solution\n");
    printf("x=%f\n",x);
    printf("y=%f\n",y);
    printf("z=%f\n",z);
}
else
{
    printf("the system is incosistent and has no solution");
}
}