Jacobi polynomials in MATLAB

Question

R shah 2016-7-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/293704-jacobi-polynomials-in-matlab

评论： R shah 2016-7-7

i want to execute jacobi polynomials for N = j-3 where j= 0, 1,...,4, a = b = 3 and x = -0.655, the code according to my knowledge will be,

syms j

a= 3;

b=3;

for j = 0:4

    display (jacobiP(j-3,3,3,-0.655))
 end

but this code isn't working kindly guide me.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Karan Gill 2016-7-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/293704-jacobi-polynomials-in-matlab#answer_227887

在 MATLAB Online 中打开

Your first argument "j-3" is negative for j=0,1,2. That's what the error message says: A nonnegative integer or a symbol is expected.

If you replace "j-3" with "j" it works.

syms j
a= 3;
b=3;
for j = 0:4
jacobiP(j,3,3,-0.655)
end
ans =
     1
ans =
   -2.6200
ans =
    3.5765
ans =
   -2.8153
ans =
    0.4398

2 个评论
显示无隐藏无

R shah 2016-7-6

thx sir, but i know this information i also know that jacobi polynomials are defined for n >= 0 but when we take the third derivative of jacobi polynomial we got the jacobi P(j-3, a+3, b+3, -0.655) i've taken the case when a = 0 and b = 0 so i must have to evaluate jacobiP(j-3, 3,3,-0.655). should i consider the values:

0 for j = 0

0 for j = 1

0 for j = 2

-2.,815 for j = 3

0.4398 for j = 4

R shah 2016-7-7

在 MATLAB Online 中打开

whats wrong if i take this code for above problem

    clc; clear all 
  syms j
  a = 3 ;
  b = 3 ;
for j = 0:4
  if  j-3 == -3:-1
     display(jacobiP(j-3,3,3,0.655) == 0)
  else  
   display(jacobiP(j-3,3,3,0.655))
  end
end