diff --git a/atmat/lattice/element_creation/atoctupole.m b/atmat/lattice/element_creation/atoctupole.m
new file mode 100644
index 000000000..e3edaf872
--- /dev/null
+++ b/atmat/lattice/element_creation/atoctupole.m
@@ -0,0 +1,43 @@
+function elem=atoctupole(fname,varargin)
+%ATOCTUPOLE Creates an octupole element
+%
+%  ATOCTUPOLE(FAMNAME,LENGTH,POLYNOMA,POLYNOMB,PASSMETHOD)
+%	
+%  INPUTS
+%  1. FNAME      - Family name 
+%  2. LENGTH     - Length[m]
+%  3. POLYNOMA   - Skew [dipole quad sext oct];	 
+%  4. POLYNOMB   - Normal [dipole quad sext oct]; 
+%  5. PASSMETHOD - Tracking function. Defaults to 'StrMPoleSymplectic4Pass'
+%
+%  OPTIONS (order does not matter)
+%    R1			 -	6 x 6 rotation matrix at the entrance
+%	 R2        	 -	6 x 6 rotation matrix at the entrance
+%	 T1			 - 	6 x 1 translation at entrance 
+%	 T2			 -	6 x 1 translation at exit
+%	 NumIntSteps -   Number of integration steps
+%	 MaxOrder    -    Max Order for multipole (1 up to quadrupole)
+%
+%  OUTPUTS
+%  1. ELEM - Structure with the AT element
+%
+%  EXAMPLES
+%    1. atoctupole(famname,length,polynoma,polynomb,passmethod,'fieldname1',value1,...)
+%   each pair {'fieldname',value} is added to the element
+%
+%  See also atdrift, atquadrupole, atsextupole, atsbend, atrbend, atskewquad,
+%          atthinmultipole, atmarker, atcorrector
+
+% Input parser for option
+[rsrc,L,PolynomA,PolynomB,method] = decodeatargs({0,0,0,'StrMPoleSymplectic4Pass'},varargin);
+[L,rsrc]                          = getoption(rsrc,'Length',L);
+[PolynomA,rsrc]                   = getoption(rsrc,'PolynomA',PolynomA);
+[PolynomB,rsrc]                   = getoption(rsrc,'PolynomB',PolynomB);
+[method,rsrc]                     = getoption(rsrc,'PassMethod',method);
+[cl,rsrc]                         = getoption(rsrc,'Class','Octupole');
+
+
+% Build the element
+elem=atbaselem(fname,method,'Class',cl,'Length',L,...
+    'PolynomA',PolynomA,'PolynomB',PolynomB,'DefaultMaxOrder',3,rsrc{:});
+end