Численное интегрирование

WCRF2007

Кто мог бы помочь с написанием программы численного вычисления двойного интеграла по треугольной области. Я просто написал в Matlabe программу(с применением встроенных функций quad и int но именно этот интеграл он не считает. Подинтегральная функция видимо такая с особенностью. Что вообще можете полезного посоветовать?

seregaohota

график построить

kachokslava

что за функция?

WCRF2007

Код в Matlabe
Функция, которую надо интегрировать Epsilon(в самом низу)
% Высота волновода км
Lz=100;
% Длина волновода км
Lx=500;
% Координаты центра возмущения км
X0=200;
Z0=75;
% Характерные размеры возмущения км
sx=25;
sz=10;
%Частота падающей СДВ-волны Гц
frequency=20*10^3;
% Длина волны в вакууме Lambda=(2pi*c/Omega)

% Частота столкновения электронов(z) Гц
v_electron=@(z5*(10^6)*exp(-0.15*(z-70;
% Концентрация электронов см^(-3)
%N0=@(z30.74*exp(0.24*(z-87;
% EN- усиливающий фактор
%fplot(N0,[0 3*Lz])
EN=1000;
N=@(z,x30.74*exp(0.24*(z-87.*(1+EN*exp(-(z-Z0).^2/sz^2-(x-X0).^2/sx^2;
%N0=@(z30.74*exp(0.24*(z-87;
%N=@(x,z1+EN*exp(-(z-Z0).^2/sz^2-(x-X0).^2/sx^2.*N0(z;
%figure;
%ezsurfc(N,[0 Lz 0 Lx])
Epsilon_0=8.8542*10^(-12); % Диэлектрическая постоянная 8,8542*10^(-12)
image_one=sqrt(-1);
e=1.6*10^(-19); % Заряд электрона
m=9.1*10^(-31); % масса электрона
Omega=2*pi*frequency; %Циклическая частота волны
Omega_plazm=@(x,z30.74*exp(0.24*(z-87.*(1+EN*exp(-(z-Z0).^2/sz^2-(x-X0).^2/sx^2.*e^2/(m*Epsilon_0.^0.5);
%Epsilon=@(x,z1-30.74*exp(0.24*(z-87.*(1+EN*exp(-(z-Z0).^2/sz^2-(x-X0).^2/sx^2.*e^2/(m*Epsilon_0.^0.5)/Omega).^2/(1-image_one*frequency/Omega;
Epsilon=@(x,z1-Omega_plazm(x,z)/Omega).^2)./(1-image_one*v_electron(z)/Omega;

k11122nu

уджос. Тебя не затруднит перевести это в ТеХ или представить в виде картинки?
Оставить комментарий
Имя или ник:
Комментарий: