// La biblioteca del sistema que utilice es
#include <stdio.h>
// Función main
 // Recibe: void
 // Devuelve: int
int main(void)
{
	float k,a,det=0;
	float m1[5][5],m2[5][5],m3[5][5]={0,0,0,0,0};
	int r1,c1,r2,c2,r3,c3,i,j,l=0;
	int x,n,m=0;
	
	printf("\tMenu de operaciones\n");
	printf("\n 1.- Suma de matrices");
	printf("\n 2.- Resta de matrices");
	printf("\n 3.- Multiplicacion de un escalar por una matriz");
	printf("\n 4.- Inversa de una matriz con determinante");
	
	printf("\n\nEscoga la operacion\n");
	scanf("%d",&x);
	
	switch("%d",x)
	{
		case 1:
			{
			printf("Escoja el tamaño de la matriz\n");
			scanf("%d",&n);
			printf("Escoja los valores de la matriz 1\n");
				for(r1=0;r1<n;r1++)
				{
					for(c1=0;c1<n;c1++)
					{
						scanf("%f",&m1[r1][c1]);
					}
				}
			
			printf("\nLa matriz 1 es:\n\n");
    		for(i=0;i<n;i++)
    			{
        		for(j=0;j<n;j++) printf("%10f",m1[i][j]);
        		printf("\n");
    			}
			
			printf("\nEscoja los valores de la matriz 2\n");
				for(r2=0;r2<n;r2++)
				{
					for(c2=0;c2<n;c2++)
					{
						scanf("%f",&m2[r2][c2]);
					}
				}
			
			printf("\n\nLa matriz 2 es:\n\n");
    		for(i=0;i<n;i++)
    			{
        		for(j=0;j<n;j++) printf("%10f",m2[i][j]);
        		printf("\n");
    			}
			
				
				for(r3=0;r3<n;r3++)
				{
					for(c3=0;c3<n;c3++)
					{
						m3[r3][c3]=m1[r3][c3]+m2[r3][c3];
					}
				}
			
    			printf("\n\nLa suma de la matriz 1 y la matriz 2 es:\n\n");
    			for(i=0;i<n;i++)
    			{
        		for(j=0;j<n;j++) printf("%10f",m3[i][j]);
        		printf("\n");
    			}
			
			
			getchar();
			break;
			}
			
		case 2:
			{
				printf("Escoja el tamaño de la matriz\n");
				scanf("%d",&n);
				printf("Escoja los valores de la matriz 1\n");
				for(r1=0;r1<n;r1++)
				{
					for(c1=0;c1<n;c1++)
					{
						scanf("%f",&m1[r1][c1]);
					}
				}
			
				printf("\nLa matriz 1 es:\n\n");
    			for(i=0;i<n;i++)
    				{
        				for(j=0;j<n;j++) printf("%10f",m1[i][j]);
        				printf("\n");
    				}
			
				printf("\nEscoja los valores de la matriz 2\n");
				for(r2=0;r2<n;r2++)
				{
					for(c2=0;c2<n;c2++)
					{
						scanf("%f",&m2[r2][c2]);
					}
				}
			
				printf("\n\nLa matriz 2 es:\n\n");
    			for(i=0;i<n;i++)
    			{
        			for(j=0;j<n;j++) printf("%10f",m2[i][j]);
        			printf("\n");
    			}
			
				
				for(r3=0;r3<n;r3++)
				{
					for(c3=0;c3<n;c3++)
					{
						m3[r3][c3]=m1[r3][c3]-m2[r3][c3];
					}
				}
			
    			printf("\n\nLa resta de la matriz 1 y la matriz 2 es:\n\n");
    			for(i=0;i<n;i++)
    			{
        			for(j=0;j<n;j++) printf("%10f",m3[i][j]);
        			printf("\n");
    			}
			
			
				getchar();
				break;
			}	
		
		case 3:
			{
				printf("\nEscriba las dimensiones de la matriz\n");
				scanf("%d",&n);
				printf("\nAsigne valores a la matriz\n");
					for(r1=0;r1<n;r1++)
					{
						for(c1=0;c1<n;c1++)
						{
							scanf("%f",&m1[r1][c1]);
						}
					}
				printf("\nLa matriz  es:\n\n");
    				for(i=0;i<n;i++)
    				{
        				for(j=0;j<n;j++) printf("%10f",m1[i][j]);
        				printf("\n");
    				}
				printf("\nEscoja el valor del escalar\n");
				scanf("%f",&k);
				
				for(r2=0;r2<n;r2++)
				{
					for(c2=0;c2<n;c2++)
					{
						m2[r2][c2]=m1[r2][c2]*k;
					}
				}
				
				printf("\nLa matriz multiplicada por el escalar es:\n\n");
    				for(i=0;i<n;i++)
    				{
        				for(j=0;j<n;j++) printf("%10f",m2[i][j]);
        				printf("\n");
    				}
				
				getchar();
				break;
			
			}
			
		case 4:
			{
				printf("Introduce la dimension de la matriz\n");
				scanf("%d",&n);
				m=n-1;
				printf("\nAsigne valores a la matriz\n");
					for(r1=0;r1<n;r1++)
					{
						for(c1=0;c1<n;c1++)
						{
							scanf("%f",&m1[r1][c1]);
						}
					}
				printf("\nLa matriz  es:\n\n");
    				for(i=0;i<n;i++)
    				{
        				for(j=0;j<n;j++) printf("%10f",m1[i][j]);
        				printf("\n");
    				}
    			
    			/*Calculando el determinante*/
    			
    			det=m1[0][0];
    				for(i=0;i<m;i++)
    					{
    						a=i+1;
    							for(j=a;j<n;j++)
    								{
    									for(l=a;l<n;l++)
    										{
    											m1[j][l] = ( m1[i][i]*m1[j][l]-m1[i][l]*m1[j][i] )/m1[i][i]; 
    										}
    								}
							det=det*m1[i+1][i+1];
    					}
			
				printf("\nEl determinante de la matriz es: %f\n",det);
				
					if(det==0)
						{
							printf("\nNo tiene inversa");
						}
					else
						{	
							printf("\nTiene inversa");
						}
					
				getchar();
				break;
					
		
				
			}
				
			default:
				{
					printf("No es una opcion valida");
					break;
				}
	}
	getchar();
	return 0;

 // se devuelve un 0.
   //que en este caso quiere decir que la salida se ha efectuado con éxito.
}