close all; clear; clc;

model = Dipole_AxialSym_Eigenvalue;
mphplot(model,'pg1');
figure; mphmesh(model);

MK = mphmatrix(model, 'sol1', 'Out', {'Kc','Ec','Dc'});
K1 = MK.Kc;
D1 = MK.Dc;
save('test2.mat','K1','D1');
[ eigvec , omega]= eigs (MK.Kc, MK.Dc, 10,-8e7*2*pi*1i) ;
eigfreq = sort(real(-diag(omega)/2/pi/1i))