% This script is for investigation of deter function
%Загрузка модели S
global c2o N2o go ro lorder
%global Xfin_i Yfin_i Xfin_o Yfin_o  % results of current solutions
mods=load('model-s.txt'); 
%Константы (Книга такая-то)
Grav = 6.670e-8;
Msun = 1.98e33;
Rsun = 6.9599e10;
%Передача переменным значений из модели 
r = mods(:,1);
p = mods(:,4);
rho = mods(:,5);
g1 = mods(:,10);
A = mods(:,15);
lq = mods(:,2);
clear mods
%Степень моды
lorder = 1;
% Other variables from model
c2 = g1.*p./rho;
N2 = (A.*exp(lq)*Grav*Msun)./(r.^3);
g = (1./r.^2).*(Grav.*(exp(lq)*Msun));
%Обезразмеривание величин
ro = r/Rsun;
c2o = c2*(Rsun/(Grav*Msun));
N2o = N2*((Rsun^3)/(Grav*Msun));
rhoo = rho*(4*pi*(Rsun^3)/Msun);
go = g*(Rsun*Rsun)/(Grav*Msun);
%omega in mHz
%list - probe values of frequencies
l_list=length(list);
dd=[];
for ii=1:l_list,
omega = list(ii);
omega = 2*pi*omega/1e3;
omegao2 = (omega^2)*(Rsun^3)/(Grav*Msun);
dd=[dd; deter(omegao2)];
disp([list(ii) dd(end)]);
end
plot(list,dd);
