//2. feladat
	cout << "# 2. feladat: kimenet a sajatertek.dat.\n";
	int N2[7]={10,50,100,250,500,1000,2000};
	int ido[7][2];
	double E=4;
	double D=-2;
	int m;
	for(m=0;m<7;m++)
	{
		int NN=N2[m]*N2[m];
		double * H = NULL;
		H = new (nothrow) double [NN];
//H torlese
		for(k=0;k<N2[m];k++)
		{
			for(l=0;l<N2[m];l++)
			{
				H[k+N2[m]*l]=0;
			}
		}
//H kitoltese
		H[0]=E;
		H[1]=D;
		H[N2[m]-1]=D;
		H[N2[m]*(N2[m]-1)]=D;
		for(k=1;k<N2[m]-1;k++)
		{
			H[(N2[m]+1)*k-1]=D;
			H[(N2[m]+1)*k]=E;
			H[(N2[m]+1)*k+1]=D;
		}
		H[N2[m]*N2[m]-1]=E;
		H[N2[m]*N2[m]-2]=D;

//Ellenorzo kiiratas
		/*for(k=0;k<N2[m];k++)
		{
			for(l=0;l<N2[m];l++)
			{
				cout << H[k+N2[m]*l] << " ";
			}
			cout << endl;
		}*/
//Sajatertekfeladat
		char jobzh = 'N';
		char uploh = 'U';
		int ldh = N2[m];
		int N3= N2[m];
		double WH[N2[m]];
		double * workh = NULL;
		workh = new (nothrow) double [NN];
		int lworkh = -1;
		int infoh; 
		time_t t1,t2;
		//int ido1, ido2;
//Tajekozodas
		dsyev_(&jobzh, &uploh, &N3, H, &ldh, WH, workh, &lworkh, &infoh);
		lworkh=workh[0];
		workh = new (nothrow) double [lworkh];
//Valodi sajatertekfeladat idomeressel
		time(&t1);
		dsyev_(&jobzh, &uploh, &N3, H, &ldh, WH, workh, &lworkh, &infoh);
		//s=t2-t1;
		time(&t2);
		ido[m][0]=N2[m];
		ido[m][1]=difftime(t2,t1);
		//cout << "t= " << difftime(t2,t1) << " s"<< endl;
		/*FILE * sajat;
		sajat=fopen("sajat.dat", "w");
		for(k=0;k<N2[m];k++)
		{
			fprintf(sajat, "%24.16lf\n", WH[k]);
		}
		fclose(sajat);*/
		if(m==6)
		{
//Mivel a sajatertekek a ket szelsoerteket kiveve ketszeresen degeneraltak, igy eleg csak a [0,pi]-be eso k-khoz kiiratni oket.
			double * sajat;
			sajat = new (nothrow) double [3/2*N2[m]];
			int p;
			double pi=4*atan(1);
			//fprintf(stdout, "pi =%24.16lf\n", pi);

			FILE * sajatert;
			sajatert=fopen("sajatertek.dat", "w");
			fprintf(sajatert, "# A 2000×2000-es matrix sajatertekei.\n# Elso oszlop: k, 2.: egzakt E(k), 3.: numerikus E(k). Csak 0<=k<=pi-ben szamolva, k=pi-re szimmetrikus az E(k) fuggveny.\n");

			for(p=0;p<=N2[m]/2;p++)
			{
				sajat[p]=2*pi/N2[m]*p;
				sajat[p+1]=E+2*D*cos(sajat[p]);
				sajat[p+2]=WH[2*p];

				fprintf(sajatert, "%24.16lf %24.16lf %24.16lf\n", sajat[p], sajat[p+1], sajat[p+2]);
			}
			fclose(sajatert);
		}
		delete []H;
		delete []workh;
		//H=NULL;
		//workh=NULL;
	}
	//cout << "m= "<< m << endl;
	cout << "# Beleirtam a programba a C++ nyelv sajat idomerojet, de tul pontatlan. Ezert a masodik feladatot kulon futtattam le es a Linux beepitett time parancsaval mertem az idejet (real time). Ezt abrazoltam N fuggvenyeben a mellekelt futido.ps fajlban.\n";
	FILE * futido;
	futido=fopen("ido.dat", "w");
	fprintf(futido, "# A kulonbozo meretu matrixok sajatertekfeladatanak megoldasi ideje. Elso oszlop: dimenzioszam, masodik: futasi ido masodpercre kerekitve.\n");
	for(m=0;m<7;m++)
	{
		fprintf(futido, "%d %d\n", ido[m][0], ido[m][1]);
	}
	fclose(futido);