как сделать чтобы функция не изменяла переданный в нее двумерный массив в главной программе? - C++

#include <iostream>
#include <conio.h>
#include <stdlib.h>
 
using namespace std;
 
void input_matr(float a[][10], int n)
{
    cout<<"Введите коэффициенты и свободные члены системы: "<<endl;
    for ( int i = 0; i < n; i++)
        for ( int j = 0; j < n+1; j++ )
            if (j==n)
            {
                cout<<" = ";
                cin >> a[ i ][ j ];
            }
            else
            {
                cout<<"A["<<i<<","<<j<<"]= ";
                cin >> a[ i ][ j ];
            }
}
 
void out_m_ex( float a[][10], int n)
{
    cout<<"Система уравнений:"<<endl;
    for ( int i = 0; i < n; i++)
        for ( int j = 0; j < n+1; j++ )
            if (j==n)
                cout<<" = "<<a[i][j]<<endl;
            else
                cout<<"  "<<a[i][j];
    cout<<endl;
 }
 
 
void j_gauss(float a[10][10], int n)
{
    int i,j,k;
     for (k=0; k<n; k++)
    {
      for (i=0; i<n; i++)
        for (j=0; j<n+1; j++)
            if ((i!=k)&&(j!=k))
                a[i][j]=a[i][j]-a[i][k]*a[k][j]/a[k][k];
     
      for (i=0; i<n+1; i++)
        if (i!=k)
            a[k][i]=a[k][i]/a[k][k];
     
      for (i=0; i<n; i++)
        if (i==k) 
            a[i][k]=1;
     else 
         a[i][k]=0;
     };
 
     cout<<"Решение системы:"<<endl;
     for (i=0; i<n; i++)
        {
            cout<<a[i][n];
            cout<<endl;
        }
}
 
void out_check(float a[][10], int n)
{
    float h[10];
    int i,j;
 
    for(i=0;i<n;i++)
        for(j=0;i<n;i++)
            h[i]+=a[i][n]*a[i][j];
    cout<<"Проверка:"<<endl;
    for(i=0;i<n;i++)
        cout<<h[i]<<endl;
}
 
void main()
{
    locale loc("russian_Russia.1251"); 
    locale::global(loc); 
 
    float a[10][10];
    int n;
    cout<<"введите количество уравнений: ";
    cin>>n;
    input_matr(a,n);
    
    system("cls");
      cout<<"      Лабораторная работа ЧМ №8"<<endl;
      cout<<"[решение СЛАУ методом jordana-gaussa]"<<endl;
      cout<<endl;
    out_m_ex(a,n);
    j_gauss(a,n);
    out_m_ex(a,n);
    _getch();
}