Где ошибка?Помогите, пожалуйста

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

Модератор: Admin

Evgeniy455
Сообщения: 1
Зарегистрирован: Ср янв 24, 2018 12:51 am

Где ошибка?Помогите, пожалуйста

Сообщение Evgeniy455 » Ср янв 24, 2018 1:02 am

Здравствуйте, есть код и там где то ошибки. Помогите!
BesJS:=proc(n,x) -BesselJ(n+1,x)+(n*BesselJ(n,x))/(x); end proc; n:=1; y0:=1.5; step:=0.1; ye:=20; for y from y0 to ye by step do `y;`

`if BesJS`(n,y-step)*BesJS(n,y) <0 then do aa:=fsolve(BesJS(n,x),x=y-step..y); print(aa) end do; `end if;` `end do;`

нужно найти все нули производной функции Бесселякоторіе лежат в интервале

Kitonum
Сообщения: 2032
Зарегистрирован: Ср дек 31, 2008 1:55 pm
Откуда: г. Пенза

Re: Где ошибка?Помогите, пожалуйста

Сообщение Kitonum » Ср янв 24, 2018 7:33 am

restart;
BesJS:=proc(n,x)
-BesselJ(n+1,x)+(n*BesselJ(n,x))/x;
end proc:

n:=1; y0:=1.5; step:=0.1; ye:=20.;
for y from y0 to ye by step do
if BesJS(n,y-step)*BesJS(n,y)<0 then
r:=y-step..y; aa:=fsolve(BesJS(n,x), x=r);
print(r, aa) end if;
end do:

plot(BesJS(1,x), x=1.4..20); # График

Markiyan Hirnyk
Сообщения: 1210
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Re: Где ошибка?Помогите, пожалуйста

Сообщение Markiyan Hirnyk » Вс янв 28, 2018 10:27 am

Для нахождения всех нулей аналитической функции в заданном интервале целесообразно применить команду RootFinding:-NextZero в условном цикле.