I was fiddling with numerical range of two by two matrices so I modified my root testing python program to do this. The numerical range of is the set of all values
for all nonzero vectors . This set is a compact set (it can be seen as the image via the mapping of vectors on the unit sphere which is a compact set). It’s convex which is harder to show. For two by two it is an elliptic disc (could be degenerate).
See the result here, it plugs in random vectors and shows the result. Here’s an example plot for the matrix .
The code is really inefficient and eats up all your cpu. There’s no effort to optimize this.