added sumset.py Fixes: #{6563}

DIRECTORY.md

@@ -574,6 +574,7 @@
   * [Sum Of Arithmetic Series](maths/sum_of_arithmetic_series.py)
   * [Sum Of Digits](maths/sum_of_digits.py)
   * [Sum Of Geometric Progression](maths/sum_of_geometric_progression.py)
+  * [Sumset](maths/sumset.py)
   * [Sylvester Sequence](maths/sylvester_sequence.py)
   * [Test Prime Check](maths/test_prime_check.py)
   * [Trapezoidal Rule](maths/trapezoidal_rule.py)
@@ -632,7 +633,7 @@
 
 ## Physics
   * [Horizontal Projectile Motion](physics/horizontal_projectile_motion.py)
-  * [Lorenz Transformation Four Vector](physics/lorenz_transformation_four_vector.py)
+  * [Lorentz Transformation Four Vector](physics/lorentz_transformation_four_vector.py)
   * [N Body Simulation](physics/n_body_simulation.py)
   * [Newtons Second Law Of Motion](physics/newtons_second_law_of_motion.py)
 

maths/sumset.py

@@ -0,0 +1,37 @@
+"""
+
+Calculates the SumSet of two sets of numbers (A and B)
+
+Source:
+    https://en.wikipedia.org/wiki/Sumset
+
+"""
+
+
+def sumset(set_a: set, set_b: set) -> set:
+    """
+    :param first set: a set of numbers
+    :param second set: a set of numbers
+    :return: the nth number in Sylvester's sequence
+
+    >>> sumset({1, 2, 3}, {4, 5, 6})
+    {5, 6, 7, 8, 9}
+
+    >>> sumset({1, 2, 3}, {4, 5, 6, 7})
+    {5, 6, 7, 8, 9, 10}
+
+    >>> sumset({1, 2, 3, 4}, 3)
+    Traceback (most recent call last):
+    ...
+    AssertionError: The input value of [B=3] is not a set
+    """
+    assert isinstance(set_a, set), f"The input value of [set_a={set_a}] is not a set"
+    assert isinstance(set_b, set), f"The input value of [B={set_b}] is not a set"
+
+    return {a + b for a in set_a for b in set_b}
+
+
+if __name__ == "__main__":
+    from doctest import testmod
+
+    testmod()