покритикуйте код

[romaf91]

напишите что вы бы переделали
я начинающий
поэтому прошу у вас совета

[Murashov]

хз.. не мне судить конечно, но я бы оформил бы константы заглавными буквами, так более удобно оформлять код. и не стал бы заливать исходник в файлы док))

[Stilet]

for(int i=0;i<=n;i++)
{
dpdR[i]=0;
R[i]=0;
HR[i]=0;
HRx[i]=0;
HRxx[i]=0;
HRxxx[i]=0;
tetta[i]=0;
tettaR[i]=0;
tettaRx[i]=0;
tettaRxx[i]=0;
Z[i]=0;
tettaZx[i]=0;
tettaZRxx[i]=0;
PRx[i]=0;
A[i]=0;
B[i]=0;
dHdt[i]=0;
Ht[i]=0;
otnoshenie[i]=0;
dJvolnadR[i]=0;
Jvolna[i]=0;
}

Вместо этого вроде можно сразу инициализировать нулями:

double dpdR[n+1]={0}; 
double R[n+1]={0},HR[n+1]={0},HRx[n+1]={0},HRxx[n+1]={0},HRxxx[n+1]={0};

Тогда этот цикель не нужен

for(int i=1;i<=n;i++)
{
//cout<<"i="<<i<<"j="<<j<<endl;
//cout<<"HRx= "<<HRx[i]<<"\t\tHRi= "<<HR[i]<<endl;
}

Он конечно скорее всего проигнорируется, но код становится корявочитабельным.
Два даже три цикла можно слить в один

		for (int i=0;i<=n-1;i++)
		{
			
			HRx[i] = (HR[i+1]-HR[i])/dR1;
			if(i>0)			HRxx[i] = (HR[i+1]-2*HR[i]+HR[i-1])/(dR2);
			if(i>1)			HRxxx[i] = (HR[i+1]-3.*HR[i]+3.*HR[i-1]-HR[i-2])/(dR3);
		}

И остальные циклы тоже лишние. Все можно в одном делать.

Короче - низачот, но пойдет.