Danke fuer Deine Bemuehungen.
Aber leider ist es nicht das, was ich brauche.
Zum einen geht es um ein n*n-Feld und zum anderen um einen
Loesungsansatz per linearer Optimierung ( Simplex und so ) 
Richard
Also: Nehme Variablen 1>=x(ij)>=0, x(ij))=1: Dame auf Feld ij, sonst nicht.
Zielfunktion maximiere sum_i,j x(ij)
Nebenbedingungen:
Zeilen: sum_i=1^n x(ij)=1 fuer alle j
Spalten sum_j=1^n x(ij)=1 fuer alle i
Diagonalen
sum_{i+j=c} x(ij)=1 fuer c=2…16
sum_(i-j)=c x(ij)=1 fuer c=-7,…,7
Ohne die Diagonalengleichungen ist das die Bedingung fuer stochastische Matritzen-> Polyeder hat ganzzahlige Ecken.
Ich wuerde tippen, dass das mit den DGs nicht mehr stimmt, dann
musst du das ganzzahlige Problem Loesen, was mit Limplex nicht geht. Dann nimm z.B. branch&bound.
Um evtl zu zeigen, dass das Teil doch ganzzahlig ist, kanst Du versuchen die Diagonalgleichungen in die Zielfunktion aufzunehmen.
MFG
Martin