m=mphload('amb_new_L1')
pdom = m.probe.create('pdom', 'DomainPoint');
pdom.model('comp1');
pdom.setIndex('coords3','2699e-4',0,0);
pdom.setIndex('coords3','869e-4 ',0,1);
pdom.setIndex('coords3','2000e-4',0,2);
M = mphstate(m,'sol1','out',{'MA','MB','C','D','Mc'},'input','current','output','comp1.ppb1')      %,'sparse','on')
%str = mphmatrix(m,'sol1','out',{'K','L','M','N'},'initmethod','init')

%save statespace2.mat;