%************************************************
% Define elastic stiffness tensor (GPa)
%************************************************
%
% Reference Elastic constants
% Isaak,D.G., Ohno, I., and Lee, P.C. 2005 Phys.Chem.Min. DOI 10.1007/s00269-005-0047-9
%
% Reference Crystal Structure
% Diopside Di93 He3 Ur2 X2  C2/c X||a*,Y||b,Z||c
%
% Define density (g/cm3)
 rho= 3.2860;
%
% Define Cartesian Tensor crystal symmetry & frame
cs_Tensor = crystalSymmetry('2/m',[  9.7456  8.9198  5.2516],...
 [  90.0000 105.8600  90.0000]*degree,'x||a*','z||c',...
 'mineral','Diopside 2005');
%
% elastic Cij stiffness tensor (GPa) as matrix M
 M =.... 
 [[  228.10   78.80   70.20    0.00    7.90    0.00];...
 [    78.80  181.10   61.10    0.00    5.90    0.00];...
 [    70.20   61.10  245.40    0.00   39.70    0.00];...
 [     0.00    0.00    0.00   78.90    0.00    6.40];...
 [     7.90    5.90   39.70    0.00   68.20    0.00];...
 [     0.00    0.00    0.00    6.40    0.00   78.10]];
%
% M as stiffness tensor C with MTEX tensor command
C = tensor(M,cs_Tensor);
%
