%specify data
index = Data(:,1); %All Property

%specify model
mod22 = arima( 'ARLags', [1,2], 'MALags', [1,2], 'Distribution', 'gaussian',            'Variance', garch(1,1));

%fit model
[fit22,EstParamCov22,logL(9)] = estimate(mod22,index,'options', options);