Tengo problemas para escribir un programa que debe calcular la temperatura de un punto central de una región rectangular tras ciertos pasos de tiempo.
Para el cual cree una matriz de (N+1)x(M+1) donde N = 2M y el punto central de la región es (N/2)+1 y (M/2)+1. Las celdas superiores tienen una temperatura inicial 20 ºC, el resto de celdas tienen una temperatura inicial de 0 ºC, ademas en las celdas del contorno la temperatura es invariante.
La temperatura en cada celda y paso de tiempo (Tk) se calcula con la temperatura en el paso de tiempo anterior (Tk-1) de las celdas adyacentes, siguiendo la siguiente formula:
(Siendo (i,j) el valor de una celda)
T = (4*[(i+1,j)+(i-1,j)+(i,j+1)+(i,j-1)]+[(i-1,j-1)+(i-1,j+1)+(i+1,j+1)+(i+1,j-1)])/10
Hasta aquí lo tengo todo claro y no tengo mayor problema. El problema empieza cuando para calcular la temperatura tras cada paso de tiempo suma la temperatura ya calculada en el propio salto de tiempo en vez de utilizar la temperatura del anterior salto de tiempo, es decir:
Para calcular la temperatura de la celda (1,6) utiliza las celdas adyacentes, pero anteriormente ya calculó la celda (1,5), y en vez de utilizar la temperatura de (1,5) antes del calculo, la suma después de dicho calculo dando un resultado erróneo.
Hasta ahora tengo escrito lo siguiente:
int main( ) {
int N = 20;
int M = 10;
int A = 2 * N + 1;
int B = 2 * M + 1;
int matriz[A][B];
for( int i = 1; i < A; i++ ) {
for( int j = 0; j < B; j++ ) {
matriz[0][j] = 20;
matriz[i][j] = 0;
}
}
for( int i = 1; i < ( A - 1 ); i++ ) {
for( int j = 1; j < ( B - 1 ); j++ ) {
int n = i - 1;
int e = j + 1;
int s = i + 1;
int o = j - 1;
int T = ( 4 * ( matriz[n][j] + matriz[s][j] + matriz[i][o] + matriz[i][e] ) + ( matriz[n][e] + matriz[n][o] + matriz[s][e] + matriz[s][o] ) ) / 10;
matriz[i][j] += matriz[i][j] + T;
}
}
std::cout << matriz[N + 1][M + 1];
return 0;
}
1
atamaño - 1
.