# Wprawka 4

## Napisz symulację przewodnictwa cieplnego.
Symulacja działa następująco: całość odbywa sie na gridzie 100x100. (Możesz założyć że ma topologię toroidalną, czyli zawija się.)
W symulacji czas jest dyskretny. W chwili T=0, ustawiamy temp[x][y] = 0.
Kilka pól (wybrane dowolnie) jest wybrane jako źródło, i ma **zawsze** temp = 100, a kilka pól jest wbrane jako ujście, i ma **zawsze** temp = -100 (ignorując updejty).

Update wykonujemy równolegle dla **wszystkich pozycji jednocześnie**, i polega na wykonaniu

temp[x][y] = (1-eps) * temp[x][y] + eps*(temp[x-1][y]+temp[x][y-1]+temp[x+1][y]+temp[x][y+1])/4

dla jakiegoś małego eps (np eps = 0.01).

Zadanie polega na tym by zwizualizować ewolucje systemu przez pierwsze kilka kroków: niska temperatura powinna być rysowana np. niebieskim kwadratem, a wysoka np. czerwonym kwadratem, z gradientem pomiędzy. (Dobór palety zostawiam do wyboru.)