initial commit for testing operator overloading on RHS

Author: AndreyAPopov

src/+otp/RHS.m

@@ -33,6 +33,13 @@
                 obj.(f) = extras.(f);
             end
         end
+
+        function newRHS = plus(obj, other)
+            objF   = obj.F;
+            otherF = other.F;
+            newF = @(t, y) objF(t, y) + otherF(t, y);
+            newRHS = otp.RHS(newF);
+        end
 
         function opts = odeset(obj, varargin)
             opts = odeset( ...