diff --git a/+unittests/+quantity/testDiscrete.m b/+unittests/+quantity/testDiscrete.m
index 68c2ff3beab74dc1c403cf69da42a029e90866c9..cd0a07aa0ad3e3069c0212812684722c6f9046e7 100644
--- a/+unittests/+quantity/testDiscrete.m
+++ b/+unittests/+quantity/testDiscrete.m
@@ -760,17 +760,16 @@ testCase.verifyEqual(azZeta+bzZeta, ABZZETA.on());
 testCase.verifyEqual(a+a, ALrA.on(z), 'RelTol', 1e-3);
 
 %% additional test
-% z = linspace(0, 1, 21);
-% zeta = linspace(0, 1, 41);
-a = quantity.Discrete(ones(numel(z), numel(zeta)), 'grid', {z, zeta}, ...
+a = quantity.Discrete(zGrid, 'grid', {z, zeta}, ...
 	'gridName', {'z', 'zeta'}, 'name', 'a');
-b = quantity.Discrete(ones(numel(zeta), numel(z)), 'grid', {zeta, z}, ...
+b = quantity.Discrete(zetaGrid.', 'grid', {zeta, z}, ...
 	'gridName', {'zeta', 'z'}, 'name', 'b');
 c = a + b;
 cMat = c.on();
+cMatReference = zGrid + zetaGrid;
 %%
-testCase.verifyEqual(numel(cMat), numel(z) * numel(zeta));
-testCase.verifyEqual(cMat(:), 2*ones(numel(z) * numel(zeta), 1));
+testCase.verifyEqual(numel(cMat), numel(cMatReference));
+testCase.verifyEqual(cMat(:), cMatReference(:));
 end
 
 function testInit(testCase)