The infinite square well potential in one dimension has a smooth supersymmetric partner potential which is shape invariant. In this paper, we study the generalization of this to two dimensions by constructing the supersymmetric partner of the disk billiard. We find that the property of shape invariance is lost in this case. Nevertheless, the WKB results are significantly improved when SWKB calculations are performed with the square of the superpotential. We also study the effect of inserting a singular flux line through the center of the disk.
