boundary_types = p,p
boundary_values = 2,3
initial_gradients =         0, -1.3,  1.3, 
reflection_path = (2,0),(3,1),(2,2),(3,0),(4,0)
grid_size = 5
max_radius = 50%