Интеграл 0 - бесконечность / График

@Jabbar · Регистрация: 28.09.2015

Студворк — интернет-сервис помощи студентам

Помогите с вычислением интеграла и графиком, пожалуйста.
Неопределенный, с пределами 0, +inf.
Подынтегральная функция:
f=(sqrt(E1-Eg-dEg)-sqrt(E1-Eg))/(E1*(E1^2-E^2))
Интегрирование производится по E1.
Eg=13.4; %константа
dEg и E задаются соответственно:
dEg= [-0.0143 -0.0230 -0.0337];
vxN=size(dEg);
lxN=vxN(2);
и
E = 1.2:.001:2;
vE = size(E);
lE = vE(2);

Думал, что нужно использовать либо int(), либо integral(), но то что я пишу, мягко скажем не работает.
К примеру такое:

Matlab M
1
2
3
4
5
6
7
8
9
N=zeros(lxN,lE);
syms E1
for j=1:lxN
    for i=1:lE
        N(j,i) = int((sqrt(E1-Eg-dEg(j))-sqrt(E1-Eg))/(E1*(E1^2-E(i)^2)), 'E1', 0, inf);
    end
end
plot(E1,N)
grid on

выдает:
The following error occurred converting from sym to double:
Error using mupadmex
Error in MuPAD command: Singularity. [arctanh]
Evaluating: symobj::double

Пишу так:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
dabgr = @(E1, E, dEg) (sqrt(E1-Eg-dEg)-sqrt(E1-Eg))./(E1.*(E1.^2-E^2));
 
dnBGR = zeros(lxN,lE);
 
for j=1:lxN
    for i=1:lE
        dnBGR(j,i) = integral(@(E1)dabgr(E1, i, j), 0, Inf);
    end
end
 
figure;
 
grid on;
plot(E, dnBGR, '-');

выдает предупреждения вида:
The integral may not exist, or it may be difficult to approximate numerically to the requested accuracy.

PS. Абсолютно уверен что делаю что-то не так :/
PSS. Стоит ли вообще здесь писать такие for циклы, или можно ограничиться поэлементными операциями?

@VSI · 28.09.2015, 22:52

Jabbar, загляните в самый низ этой страницы в раздел "Похожие темы". Можете найти для себя много полезного...

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11680&d=1772460536 Одним из. . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .
SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .

@VSI Модератор $Эксперт по математике/физике$ 5291 / 4073 / 1392 Регистрация: 30.07.2012 Сообщений: 12,490
	28.09.2015, 22:52
	Jabbar, загляните в самый низ этой страницы в раздел "Похожие темы". Можете найти для себя много полезного... 0