Неопределенный интеграл

Форум пользователей пакета Maple

Модератор: Admin

Кли
Сообщения: 156
Зарегистрирован: Пт ноя 04, 2016 4:54 pm

Неопределенный интеграл

Сообщение Кли » Вс мар 17, 2019 9:51 pm

У меня есть несобственный интеграл второго рода, я вычислил с(особая точка), далее использую метод рунге-ромберга для формулы Симпсона, но че то у меня неправильно вычисляет

Код: Выделить всё

restart;
n:=4;
a:=-0.1*n;
b:=0.1*n;
f:=unapply(sin(x-0.1*n)/sqrt((x+0.1*n)*(x+5)),x);
eps:=evalf(10^(-3)):
ss:=evalf(int(1/sqrt((x+0.1*n)*(x+5)),x));
ñ:=solve(evalf(subs(x=0.4,ss)-subs(x=g,ss)));
assign(c):
N:=5: m:=4:
h:=(b-a)/N:
x[0]:=a:
for k from 1 to N do
x[k]:=x[k-1]+h:
end do:
Y[1]:=0:
for k from 1 to N do 
Y[1]:=evalf(Y[1]+h*(f(x[k-1])+4*f((x[k]+x[k-1])/(2))+f(x[k]))/6):
end do:
print('Y[1]' = Y[1]);
y[0]:=a:
for k from 1 to 2*N do
y[k]:=y[k-1]+(h/2):
end do:
Y[2]:=0:
for k from 1 to 2*N do
Y[2]:=evalf(Y[2]+(h/2)*((f(y[k-1])+4*f((y[k]+y[k-1])*(1/2))+f(y[k]))*(1/6))):
end do:
print('Y[2]' = Y[2]);
while(abs(((Y[2]-Y[1])/((2^m)-1)))) > eps do
N:=N+1:
h:=(b-a)/N:
x[0]:=a:
for k from 1 to N do
x[k]:=x[k-1]+h:
end do:
Y[1]:=0:
for k from 1 to N do 
Y[1]:=evalf(Y[1]+h*(f(x[k-1])+4*f((x[k]+x[k-1])/(2))+f(x[k]))/6):
end do:
y[0]:=a:
for k from 1 to 2*N do
y[k]:=y[k-1]+(h/2):
end do:
Y[2]:=0:
for k from 1 to 2*N do
Y[2]:=evalf(Y[2]+(h/2)*((f(y[k-1])+4*f((y[k]+y[k-1])*(1/2))+f(y[k]))*(1/6))):
end do:
end do:
print('N' = N);
Y:=Y[2]+((Y[2]-Y[1])/((2^m)-1));
int(f(x),x=a..b);
G:=evalf(int(f(x),x=a..b));

Кли
Сообщения: 156
Зарегистрирован: Пт ноя 04, 2016 4:54 pm

Re: Неопределенный интеграл

Сообщение Кли » Вт мар 19, 2019 6:56 pm

там в точке особенность происходит, а точка с это новая точка, один интеграл я оцениваю, другую численно вычисляю