0

```
#include<stdio.h>
#include<conio.h>
void main()
{
float a[5][5],x[3];float t,s;
int i,j,k;
printf("enter a matrix of order 3*4");
for(i=0;i<3;i++)
{
for(j=0;j<4;j++)
{
scanf("%f",&a[i][j]);
}
}
for(i=0;i<3;i++)
{
for(j=0;j<4;j++)
printf("%f ",a[i][j]);
printf("\n");
}
for(i=0;i<=1;i++)
{
for(j=i+1;j<=2;j++)
{
t=a[j][i]/a[i][i];
for(k=i;k<=3;k++)
a[j][k]=a[j][k]-t*a[i][k];
}
}
printf("The upper triangular matrix is as\n");
for(i=0;i<3;i++)
{
for(j=0;j<4;j++)
printf(" %f",a[i][j]);
printf("\n");
}
for(i=2;i>=0;i--)
{
s=0.0;
for(j=i+1;j<3;j++)
{
s=s+a[i][j]*x[j];
x[i]=(a[i][3]-s)/a[i][i];
}
}
printf("\nRoots of the equation are");
for(i=0;i<=2;i++)
printf(" %f",x[i]);
getch();
}
```