clear?al?;
clc;
x=0:0.01:5;
v=[-0.500.511.5];
p=2./（0.5*gamma（v+1））.*（（0.2/（2*0.5））.^（v+1））.*besselk（v0.2/0.5）;
figure（1）
p1=2./（0.5*gamma（v（1）+1））.*（（x/（2*0.5））.^（v（1）+1））.*besselk（v（1）x/0.5）;
plot（xp1）
text（0.2p（1）‘\leftarrowσv=-0.5‘）;
hold?on;
p2=2./（0.5*gamma（v（2）+1））.*（（x/（2*0.5））.^（v（2）+1））.*besselk（v（2）x/0.5）;
plot（xp2）
text（0.2p（2）‘\leftarrowσv=0‘）;
p3=2./（0.5*gamma