code: def show_plane(orig, n, scale=1.0, **kwargs):
      
      b1 = orthogonal_vector(n)
      b1 /= la.norm(b1)
      b2 = np.cross(b1, n)
      b2 /= la.norm(b2)
      verts = [orig + scale*(-b1 - b2),
               orig + scale*(b1 - b2),
               orig + scale*(b1 + b2),
               orig + scale*(-b1 + b2)]
      faces = [(0, 1, 2), (0, 2, 3)]
      trimesh3d(np.array(verts), faces, **kwargs)
docs: Show the plane with the given origin and normal. scale give its size