论坛 › 函数论与泛函分析 › 特殊函数

yaoyuan4102 发表于 2009-4-28 12:45

特殊函数

function BF = Beta(x,y)
format long;
BF = exp(gammaln(x)+gammaln(y)-gammaln(x+y));

function =conjgrad(A,b,x0)
if(nargin == 3)
    eps = 1.0e-6;
end

r1  = b-A*x0;
p1  = r1;
d   = dot(r1,r1)/dot(p1,A*p1);
x   = x0+d*p1;
r2  = r1-d*A*p1;
f   = dot(r2,r2)/dot(r1,r1);
p2  = r2+f*p1;
n   = 1;

for(i=1:(rank(A)-1))
    x0 = x;
    p1 = p2;
    r1 = r2;
    d  = dot(r1,r1)/dot(p1,A*p1);
    x  = x0+d*p1;
    r2 = r1-d*A*p1;
    f  = dot(r2,r2)/dot(r1,r1);
    p2 = r2+f*p1;
    n  = n + 1;
end

d  = dot(r2,r2)/dot(p2,A*p2);
x  = x+d*p2;
n  = n + 1;

羽仙 发表于 2009-5-9 23:29

楼主辛苦了
谢谢分享

xxgzftj 发表于 2011-12-28 15:46

{:soso_e103:}

孪生素数 发表于 2012-2-26 10:43

这个是····？？？

竹下夜月 发表于 2012-2-27 08:48

运行一下，看看是啥东西:lol

645464 发表于 2012-9-10 23:22

哎，，，，，，，

查看完整版本: 特殊函数