gamma function error in calculation

Question

yogeshwari patel 2022-4-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1698005-gamma-function-error-in-calculation

评论： Torsten 2022-4-17

在 MATLAB Online 中打开

% Starting value

The above formula is coded as follows:

syms x a
Y=sym(zeros(1));
Y(1)=0;
a=1/2
for i=1:4
    if i==5
        A=1
    else
        A=0
    end
    if i==4
        B=1
    else
        B=0
    end
    Y(i+1)=simplify((gamma(a*(i-1)+1)/gamma((a*(i-1)+3/2))*(A-Y(i)+((2*B)/gamma(5/2)))));
end
disp(Y)

But it is showing the calculation error Y(5)=1 but the value is shown in MATLAB is as follows:('2535301200456458897054207582575/2535301200456458802993406410752'). Ecxept Y(5) all values are zeros

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

请先登录，再进行评论。

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

Answer 1

VBBV 2022-4-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1698005-gamma-function-error-in-calculation#answer_944090

在 MATLAB Online 中打开

syms x a

Y=sym(zeros(1));

Y(1)=0;

a=1/2

a = 0.5000

for i=1:4

if i==5

A=1;

B = 0;

Y(i+1)=simplify((gamma(a*(i-1)+1)/gamma((a*(i-1)+3/2))*(A-Y(i)+((2*B)/gamma(5/2)))));

else

A=0;

B =1 ;

Y(i+1)=simplify((gamma(a*(i-1)+1)/gamma((a*(i-1)+3/2))*(A-Y(i)+((2*B)/gamma(5/2)))));

end

if i==4

B=1;

A = 0;

Y(i+1)=simplify((gamma(a*(i-1)+1)/gamma((a*(i-1)+3/2))*(A-Y(i)+((2*B)/gamma(5/2)))));

else

B=0;

A = 1;

Y(i+1)=simplify((gamma(a*(i-1)+1)/gamma((a*(i-1)+3/2))*(A-Y(i)+((2*B)/gamma(5/2)))));

end

A = 0

B = 0

A = 1

A = 0

B = 0

A = 1

A = 0

B = 0

A = 1

A = 0

B = 1

disp(vpa(Y,4))

3 个评论
显示 1更早的评论隐藏 1更早的评论

yogeshwari patel 2022-4-17

The correct value are Y(1)=0; Y(2)=0;Y(3)=0, Y(4)=0, Y(5)=1 and rest value s are zero. Now if I work with pen and paper the values are same as I mention above but when I work with MATLAB code the value of Y(5)=2535301200456458897054207582575/2535301200456458802993406410752 .Why that error is appaering