这个函数怎么实现呢？

MiniTime

function y = C1( gama )
%C1 此处显示有关此函数的摘要
%   此处显示详细说明
syms k;
f1 = @(k)dfac(2*k-3)/(factorial(k)*factorial(k+1)) * power(gama/2,k);
y = 3*sqrt(pi)/4 * (1 - gama/4 - symsum(f1,k,1,inf));
end

function y = dfac( x )
%DFAC 双阶乘
%   此处显示详细说明
y = 1;
if x > 0
    for i = x:-2:1
        y = y * i;
    end
end
if x == 0
    y = 1;
end
if x < 0
    x = abs(x);
    for i = x:-2:1
        y = y * i;
    end
    y = 1 / y;
end

//---------------------------------命令行窗口-------------------------------------------//
y = C1(2)无法从 sym 转换为 logical。
出错 dfac (line 5)if x > 0
出错 C1 (line 5)f1 = dfac(2*k-3)/(factorial(k)*factorial(k+1)) * power(gama/2,k)

二爷名声在外

把表达式从符号表达式换成普通表达式试试，这里是函数传递的时候类型出错了……

MiniTime

二爷名声在外 发表于 2015-1-24 20:51
把表达式从符号表达式换成普通表达式试试，这里是函数传递的时候类型出错了……

去掉@（k）吗？

二爷名声在外

如果这样，定义符号那一句也应该去掉

MiniTime

二爷名声在外 发表于 2015-1-24 21:08
如果这样，定义符号那一句也应该去掉

应该是dfac的问题，提示：无法从 sym 转换为 logica