if i or j <> 1,3,7,... then    
'perciò tutti gli A1(1,x), A1(x,1),A1(3,x),A1(x,3),A1(7,x),A1(x,7) (con x variabile da 0 alla 'maxdim dell'array) saranno esclusi  
for i=0 to maxdim  
for j=0 to maxdim  
A2(i,j)=A1(i,j)  
next  
next  
end if

     Dim Nums as integer = {1 , 3 ,.................}
     For i=0 to maxdim  
        For z=0 to Nums.Lenght -1
             if i=Nums(z) then goto 20
        Next
        For j=0 to maxdim 
             For K=0 to Nums.Lenght -1
                  if j=Nums(K) then goto 10
             Next
        A2(i,j)=A1(i,j)  
10     Next
20  Next

For i=0 to maxdim  
    For j=0 to maxdim  
        Select Case True
            Case <indici_da_non_copiare)
                ' non copia
            Case Else
                A2(i,j) = A1(i,j)
         End Select
    Next
Next 

' creo gli array e li valorizzo
Dim a1(4, 9) As Integer
Dim a2(4, 9) As Integer
Dim s As String = ""
Dim indiciDaNonCopiare As String = "1,3,7"
For i As Integer = 0 To 4
    For j As Integer = 0 To 9
        a1(i, j) = i + j
    Next
Next

' copia l'array a1 nell'array a2
For i As Integer = 0 To 4
    For j As Integer = 0 To 9
        Select Case True
            Case indiciDaNonCopiare.Contains(CStr(i)), indiciDaNonCopiare.Contains(CStr(j))
                ' non copiare
            Case Else
                a2(i, j) = a1(i, j)
                s = String.Format("{0}{1},{2}{3}", s, CStr(i), CStr(j), Environment.NewLine)
        End Select
    Next
Next

' mostra gli array copiati
MessageBox.Show(s)

private static int[,] Riduci(int[,] originale, int[] indici)
        {
            int[,] ridotto = new int[10, 10]; //uno sufficientemente grande.
            int count = 0;
            for (int i = 0; i < originale.GetLength(0); i++)
                for (int j = 0; j < originale.GetLength(1); j++)
                    if (!indici.Contains(i) && !indici.Contains(j))
                    {
                        int ri = count / (originale.GetLength(1) - indici.GetLength(0));
                        int rj = count % (originale.GetLength(1) - indici.GetLength(0));
                        ridotto[ri, rj] = originale[i, j];
                        count++;
                    }

            return ridotto;
        }

int[,] originale = new int[5, 8];
            //inizializzo sequenzialmente
            for (int i = 0; i < originale.GetLength(0); i++)
            {
                for (int j = 0; j < originale.GetLength(1); j++)
                {
                    originale[i, j] = (i * originale.GetLength(1)) + j;
                    Console.Write(originale[i, j] + " - ");
                }
                Console.WriteLine();
            }

            int[] indici = { 0, 2, 4 }; //elenco indici da eliminare
            int[,] ridotto = Riduci(originale, indici);

            //output
            Console.WriteLine("Output matrice ridotta");
            for (int i = 0; i < originale.GetLength(0) - indici.GetLength(0); i++)
            {
                for (int j = 0; j < originale.GetLength(1) - indici.GetLength(0); j++)
                {
                    Console.Write(ridotto[i, j] + " - ");
                }
                Console.WriteLine();
            }

''' <summary>
''' Rappresenta una matrice a due dimensioni. Fornisce un metodo per restituire un istanza ridotta attravrso gli indici da escludere.
''' </summary>
''' <typeparam name="T">Tipo di elementi contenuti nella matrice</typeparam>
Public Class MatriceRiducibile(Of T)
	'contenitore interno unidimensionale degli elementi
	Private elementi As List(Of T)

	''' <summary>
	''' Crea una nuova istanza della classe MatriceRiducibile
	''' </summary>
	''' <param name="righe">Dimensione 0 (righe)</param>
	''' <param name="colonne">Dimensione 1 (colonne)</param>
	Public Sub New(righe As Integer, colonne As Integer)
		Length0 = righe
		Length1 = colonne

		elementi = New List(Of T)()
	End Sub

	''' <summary>
	''' Restituisce la dimensione 0 (righe)
	''' </summary>
	Public Property Length0() As Integer
		Get
			Return m_Length0
		End Get
		Private Set
			m_Length0 = Value
		End Set
	End Property
	Private m_Length0 As Integer

	''' <summary>
	''' Restituisce la dimensione 1 (colonne)
	''' </summary>
	Public Property Length1() As Integer
		Get
			Return m_Length1
		End Get
		Private Set
			m_Length1 = Value
		End Set
	End Property
	Private m_Length1 As Integer

	''' <summary>
	''' Ottiene o imposta l'elemento agli indici specificati 
	''' </summary>
	''' <param name="r">Indice di riga</param>
	''' <param name="c">Indice di colonna</param>
	''' <exception cref="System.IndexOutOfRangeException"></exception>
	''' <exception cref="System.InvalidOperationException"></exception>
	''' <remarks>L'inserimento di nuovi elementi può avvenire solo in modo sequenziae</remarks>
	Public Default Property Item(r As Integer, c As Integer) As T
		Get
			'check sugli indici
			If r > Length0 OrElse c > Length1 Then
				Throw New IndexOutOfRangeException("indici fuori dimensione")
			End If

			'elemento non ancora inserito
			If r * Length1 + c > elementi.Count Then
				Return Nothing
			End If

			Return elementi.ElementAt(r * Length1 + c)
		End Get
		Set
			'check sugli indici
			If r > Length0 OrElse c > Length1 Then
				Throw New IndexOutOfRangeException("indici fuori dimensione")
			End If

			'elemento non ancora inserito
			If r * Length1 + c > elementi.Count Then
				Throw New InvalidOperationException("Elemento non modificabile")
			End If

			'gli indici corrispondono alla posizione del prossimo elemento quindi aggiungo...
			If r * Length1 + c = elementi.Count Then
				elementi.Add(value)
			Else
				'...altrimenti l'elemento è gia stato inserito quindi modifico
				elementi(r * Length1 + c) = value
			End If
		End Set
	End Property

	''' <summary>
	''' Restituisce una nuova istanza della classe MatriceRiducibile che contiene tutti gli elementi dell'oggetto corrente eccetto quelli i cui indici di riga o di colonna sono presenti nell'insieme di input
	''' </summary>
	''' <param name="indici">Insieme degli indici da escludere</param>
	''' <exception cref="System.ArgumentOutOfRangeException"></exception>
	Public Function Riduci(indici As Integer()) As MatriceRiducibile(Of T)
		If Not IndiciDiRiduzioneValidi(indici) Then
			Throw New ArgumentOutOfRangeException("L'insieme contiene indici non validi")
		End If

		Dim ridotta As New MatriceRiducibile(Of T)(Length0 - indici.GetLength(0), Length1 - indici.GetLength(0))

		Dim c As Integer = 0
		For i As Integer = 0 To Length0 - 1
			For j As Integer = 0 To Length1 - 1
				If Not indici.Contains(i) AndAlso Not indici.Contains(j) Then
					Dim ri As Integer = c \ ridotta.Length1
					Dim rj As Integer = c Mod ridotta.Length1
					ridotta(ri, rj) = Me(i, j)
					c += 1
				End If
			Next
		Next

		Return ridotta
	End Function

	Private Function IndiciDiRiduzioneValidi(indici As Integer()) As Boolean
		Return indici.All(Function(i) i <= Length0 AndAlso i <= Length1)
	End Function
End Class

Public Shared Function MKred(ByRef DTG As DataGridView, ByRef ArrayV() As Integer, ByRef matk As Array) As Array 

      
        Dim dsNl As New DataSet("MK") 'crea il dataset
        Dim dtNl As New DataTable("MKred") ' crea il datatable

        Dim col As DataColumn

        For c As Integer = 0 To dim_tot
            col = New DataColumn(c, GetType(String))
            dtNl.Columns.Add(col)
        Next

        For i As Integer = 0 To dim_tot
            Dim newRow As DataRow = dtNl.NewRow()
            For j As Integer = 0 To dim_tot

                newRow(j) = matk(i, j)
            Next
            dtNl.Rows.Add(newRow)
        Next
        'aggiunge la tabella al dataset
        dsNl.Tables.Add(dtNl)

        DTG.DataSource = dtNl
        DTG.RowHeadersVisible = False
        DTG.AutoSize = False

        Dim v As Integer
        Dim l, di As Integer
        di = Math.Sqrt(matk.Length) - 1
        l = ArrayV.Length - 1
        For v = 0 To l
            If v = 0 Then
                DTG.Columns.RemoveAt(ArrayV(v))
                DTG.Rows.RemoveAt(ArrayV(v))
            Else
                DTG.Columns.RemoveAt(ArrayV(v) - v)
                DTG.Rows.RemoveAt(ArrayV(v) - v)
            End If
        Next

               
        Dim MRKr(maxdim - 1, maxdim - 1) As Double

        For i As Integer = 0 To maxdim - 1
            For j As Integer = 0 To maxdim - 1
                MRKr(i, j) = CDbl(DTG.Item(j, i).Value)
            Next
        Next

        Return MRKr

    End Function

'nb maxdim e dim_tot possono essere passati all'interno della function o messi tra i parametri