Документ взят из кэша поисковой машины. Адрес оригинального документа : http://crydee.sai.msu.ru/ftproot/users/vab/oscil/deter.m Дата изменения: Tue May 3 03:04:12 2005 Дата индексирования: Tue Oct 2 15:27:52 2012 Кодировка:

function DD= deter(om);
global ro go % for boundary condition
global omegao2 lorder
global sol_out sol_in kk_out
global r_in r_cut r_out
omegao2=om;
if lorder>0,
r_tp=r_in^(lorder-1);
Y0 = [omegao2*r_in*r_tp; lorder*r_tp];
else
Y0 = [omegao2; 0];
end
opt_ode = odeset('RelTol',3e-5);
sol_in = ode45(@Diffur, [r_in r_cut], Y0, opt_ode);
%
Yf0 = deval (sol_in, r_cut);
%
% downward
Y0 = [go(1); 1];
sol_out = ode45(@Diffur, [r_out r_cut], Y0, opt_ode);
%
Yf1 = deval(sol_out, r_cut);
kk_out=Yf0(1)/Yf1(1);
%
DD=det([Yf0 Yf1]);
%disp(['try: ' num2str(om) ' ' num2str(DD)])