NUMERICAL IMPLEMENTATION IN 1 D AND 2 D OF A SHAPE OPTIMIZATION PROBLEM WITH ROBIN BOUNDARY CONDITIONS