diff --git a/src/generateOmega.m b/src/generateOmega.m
index c0df53cb0863f62d0585ad07f8a4c60f0aa74b39..686e6fb114ebb18deca03da30d3ab96422b95ff5 100644
--- a/src/generateOmega.m
+++ b/src/generateOmega.m
@@ -2,12 +2,19 @@ function [P_Omega, P_OmegaC] = generateOmega(OmegaTypes, nxhat, nd, params)
 
 % Simulates projection operator corresponding to some missing data segment Omega.
 %
-% INPUT:    OmegaTypes  Cell-string-Array; possible strings are 'rectangle',
-%                       'constant_height', 'constant_width', 'random', 'random position constant size'
-%        nxhat -       Number of observation directions/ rows
-%        nd -           Number of illumination directions/ columns
-%        params -       5xlength(OmegaTypes) - array that specifies the
-%                       shapes of the components of Omega
+% INPUT:    OmegaTypes  Cell-string-Array, possible strings are 'rectangle',
+%                       'constant_height', 'constant_width', 'random',
+%                       'random position constant size'.
+%           nxhat       Number of observation directions, i.e. rows.
+%           nd          Number of illumination directions, i.e. columns.
+%           params      Further specification of the shapes of the components of
+%                       Omega, 5*nr_of_OmegaTypes-array.
+% OUTPUT:   P_Omega     Projection operator corresponding to non-observable part Omega,
+%                       nxhat*nd-array, only 0 and 1 entries, zero-matrix for splitting only.
+%           P_OmegaC    Projection operator corresponding to observable part Omega^C,
+%                       nxhat*nd-array, only 0 and 1 entries, one-matrix for splitting only.
+%
+% *************************************************************************************   
     [~,nr_of_OmegaTypes] = size(OmegaTypes);
     P_OmegaC = ones(nxhat,nd);
     for iterk = 1:nr_of_OmegaTypes