diff --git a/+quantity/Symbolic.m b/+quantity/Symbolic.m index 7541bf257b760c77d4f91c37056c5993bdad73e5..9b65518b7270779da3ca25741fe51611f6c156b1 100644 --- a/+quantity/Symbolic.m +++ b/+quantity/Symbolic.m @@ -744,36 +744,7 @@ classdef Symbolic < quantity.Function 'name', "{" + obj(1).name + "}^{H}", ... 'domain', obj(1).domain); end % expm() - - function C = int(obj, z, a, b) - % from help sym/int: - % "int(S,v,a,b) is the definite integral of S with respect to v - % from a to b." - - if nargin == 2 - - grdIdx = obj(1).domain.gridIndex(z); - a = obj(1).grid{grdIdx}(1); - b = obj(1).grid{grdIdx}(end); - end - - if nargin == 3 - b = a; - a = z; - z = num2cell( obj(1).variable ); - end - if nargin == 1 - assert(obj.nargin == 1, 'Not Implemented') - z = obj(1).grid{1}; - a = z(1); - b = obj.variable; - z = num2cell( obj(1).variable ); - end - - C = obj.cumInt(z, a, b); - end % int() - function C = cumInt(obj, domain, lowerBound, upperBound) % CUMINT cumulative integration % result = cumInt(obj, domain, lowerBound, upperBound) diff --git a/+signals/SignalModel.m b/+signals/SignalModel.m index 4a6ab5c7a86fb307a6c67ec9b94809e738a264e6..5a71e133d6cfffe095727b27872857be8cb0682e 100644 --- a/+signals/SignalModel.m +++ b/+signals/SignalModel.m @@ -112,7 +112,7 @@ classdef SignalModel end methods (Static) - function sm = constant(varargin) + function sm = constant(varargin) S = 0; p = 1; diff --git a/+unittests/+quantity/testSymbolic.m b/+unittests/+quantity/testSymbolic.m index 38ea50f92254c4aee045afe7aa49161db6520359..68c6373acb4b4bae287a8b390d5fdf5e12da1d4b 100644 --- a/+unittests/+quantity/testSymbolic.m +++ b/+unittests/+quantity/testSymbolic.m @@ -132,6 +132,18 @@ fFlipped2 = f.flipGrid({'z', 'zeta'}); tc.verifyEqual(fReference2.on, fFlipped2.on, 'AbsTol', 10*eps); end % testFlipGrid(); +function testInt(testCase) +t = sym('t'); +T = quantity.Domain.defaultDomain(10, 't'); + +f = quantity.Symbolic(sin(t), 'domain', T); +F = int(f); + +testCase.verifyEqual(F, double( int( sin(t), t, 0, 1 ) ), 'AbsTol', 6e-17); + + +end + function testCumInt(testCase) tGrid = linspace(pi, 1.1*pi, 51)'; sGrid = tGrid;