diff --git a/examples/integration.jl b/examples/integration.jl
index 6908480..1ad3098 100644
--- a/examples/integration.jl
+++ b/examples/integration.jl
@@ -16,7 +16,6 @@ function setindex(x::TaylorArray{T, N, A, P}, index, d) where {T, N, A, P}
     ntuple(i -> i == index + 1 ? d : v[i], Val(P + 1)) |> TaylorArray
 end
 
-
 """
 Computes the taylor integration of order P
 
diff --git a/src/primitive.jl b/src/primitive.jl
index 7966255..bc25b08 100644
--- a/src/primitive.jl
+++ b/src/primitive.jl
@@ -3,7 +3,7 @@ import Base: exp, exp2, exp10, expm1, log, log2, log10, log1p, inv, sqrt, cbrt
 import Base: sin, cos, tan, cot, sec, csc, sinh, cosh, tanh, coth, sech, csch, sinpi, cospi
 import Base: asin, acos, atan, acot, asec, acsc, asinh, acosh, atanh, acoth, asech, acsch
 import Base: +, -, *, /, \, ^, >, <, >=, <=, ==
-import Base: sinc, cosc, hypot, max, min, literal_pow
+import Base: sinc, cosc, hypot, max, min, literal_pow, sincos
 
 Taylor = Union{TaylorScalar, TaylorArray}
 
@@ -72,6 +72,8 @@ for func in (:sin, :cos)
     end
 end
 
+sincos(t::TaylorScalar) = (sin(t), cos(t))
+
 # Binary
 
 ## Easy case