beta_biaozhunquxian.m

%%%%Betaֵм
%%%%ԵõͬBetaֵı׼ߣֻݵ࣬ٶȱ
clear;clc;
digits(5)                 %%ʾݵλΪǾ
m=sym('m');
A=0.1:0.1:8.5; 
Beta=2;
i=length(Beta);
j=length(A);                  %%betaAֵijֵȷѭ
for k=1:i                    %%betaֵѭ
    for l=1:j                %%Aֵѭ
        fun=exp(-(m+(Beta(k).*Beta(k))/(4.*m)))./m;
        y=int(fun,m,A(l),inf)
        fn=vpa(y)
        fnn(k,l)=double(fn)            %%fnͱתΪֵfnn        
        wc(k,l)=2.*besselk(0,Beta(k))-fnn(k,l)        
    end
    wz=wc(k,:)
loglog(A,wz,'-b*')
hold on 
end