diff --git a/maths/relu.py b/neural_network/activation_functions/relu.py
similarity index 100%
rename from maths/relu.py
rename to neural_network/activation_functions/relu.py
diff --git a/neural_network/activation_functions/scaled_exponential_linear_unit.py b/neural_network/activation_functions/scaled_exponential_linear_unit.py
new file mode 100644
index 000000000000..11e28e1dd470
--- /dev/null
+++ b/neural_network/activation_functions/scaled_exponential_linear_unit.py
@@ -0,0 +1,44 @@
+"""
+Implements the Scaled Exponential Linear Unit or SELU function.
+
+The function takes a vector of K real numbers and two real numbers
+alpha(default = 1.6732) & lambda (default = 1.0507) as input and
+then applies the SELU function to each element of the vector.
+SELU is a self-normalizing activation function. It is a variant
+of the ELU. The main advantage of SELU is that we can be sure
+that the output will always be standardized due to its
+self-normalizing behavior. That means there is no need to
+include Batch-Normalization layers.
+
+References :
+https://iq.opengenus.org/scaled-exponential-linear-unit/
+"""
+
+import numpy as np
+
+
+def scaled_exponential_linear_unit(
+    vector: np.ndarray, alpha: float = 1.6732, _lambda: float = 1.0507
+) -> np.ndarray:
+    """
+    Applies the Scaled Exponential Linear Unit function to each element of the vector.
+    Parameters : vector : np.ndarray
+                 alpha : float (default = 1.6732)
+                 _lambda : float (default = 1.0507)
+    Returns : np.ndarray
+    Formula : f(x) = _lambda * x if x > 0
+                     _lambda * alpha * (e**x - 1) if x <= 0
+    Examples :
+    >>> scaled_exponential_linear_unit(np.array([1.3, 3.7, 2.4]))
+    Output : np.array([1.36591, 3.88759, 2.52168])
+
+    >>> scaled_exponential_linear_unit(np.array([2.342, -3.455, -7.2116, 0.0, -4.532]))
+    Output : np.array([2.4607394, -1.70249977, -1.75673386,  0., -1.73911634])
+    """
+    return _lambda * np.where(vector > 0, vector, alpha * (np.exp(vector) - 1))
+
+
+if __name__ == "__main__":
+    import doctest
+
+    doctest.testmod()
diff --git a/maths/sigmoid.py b/neural_network/activation_functions/sigmoid.py
similarity index 100%
rename from maths/sigmoid.py
rename to neural_network/activation_functions/sigmoid.py
diff --git a/maths/sigmoid_linear_unit.py b/neural_network/activation_functions/sigmoid_linear_unit.py
similarity index 100%
rename from maths/sigmoid_linear_unit.py
rename to neural_network/activation_functions/sigmoid_linear_unit.py
diff --git a/maths/tanh.py b/neural_network/activation_functions/tanh.py
similarity index 100%
rename from maths/tanh.py
rename to neural_network/activation_functions/tanh.py