Решение бесконечного интеграла

[Ayamy]

все правильно.

[type_Oleg]

Сейчас начал считать, но порядки - уж слишком..
Вот значения функции на нижнем пределе интегрирования

(c1)exp(-m*x)/(L^5*(exp(c2/L)-1) =
= (3.75e-16)exp(-250*1.18)/((0.6e-6)^5*(exp(2.56e-6/0.9e-6)-1) =
= 2.865e-144 / 9.56e-30 = 3.00e-115
10 в минус 115. Что-то не то. Интеграл соответственно - примерно 1.5e-121

Наверное, иза exp(-m*x) . e в степени минус 295

[type_Oleg]

Ну, вообщем разбирайесь с формулой и параметрами. А я сделал интеграл " до бесконечности ".

program Integral;

var y,x,L,c1,c2,m, tin,dL : Double;
const 
  NINT=2000;  // число трапеций разбиения
  EPS=0.0001;  // относительная погрешность

// функция, которую нужно проинтегрировать 
function fun(ll:Double):Double;
begin
 fun:=c1*Exp(-m*x)/(ll*ll*ll*ll*ll*(Exp(c2/ll)-1));  
end;

// функция для расчета интеграла от L1 до L2 методом трапеций
function IntegrTrap(L1,L2:Double;n:Integer):Double;
var i:Integer;
    dl,s:Double;
begin
  s:=0;
  dl:=(L2-L1)/n;
  for i:=1 to n do
     s:=s+dl*(fun(L1+i*dl)+fun(L1+(i-1)*dl))/2;
  IntegrTrap:=s;
end;

// сама программа

begin
 x:=1.18;
 c1:=3.75e-16;
 c2:=2.56e-6;
 m:=250;
 L:=0.9e-6;
 y:=0;
 dL:=L;  // шаги приближения к  " бесконечности "
 repeat 
  tin:=IntegrTrap(L,L+dL,NINT);	
  y:=y+tin;
  L:=L+dL;
 until Abs(tin)<Abs((EPS*y));
 Writeln(y);
 Readln;	
end .

Принцип такой. Задаемся относительной погрешностью EPS (в данном случае 0,0001) .Выбыраем шаг приближения к бесконечности ( в данном случае L ).
Вычисляем интеграл от L до 2L, суммируем, затем от 2L до 3L, суммируем , и т.д. До тех пор, пока очередной " маленький " интеграл не станет меньше 0,0001 от уже накопленного.
Проверил на других функциях, нормально интегрирует " до бесконечности ".

[Ayamy]

там единственная проблема с х-толщиной поглощающего слоя..мне просто дано что это толщина атмосферы (из-за нее экспонента становится нулем).

Спасибо большое Вам за помощь))) Очень помогла Ваша программа.