f=-pi:0.1:pi;
g1=exp（cos（f））/（2*pi*besseli（01））;
plot（fg1）;
hold?on;
g2=exp（2*cos（f））/（2*pi*besseli（02））;
plot（fg2）;
g3=exp（5*cos（f））/（2*pi*besseli（05））;
plot（fg3）;
g4=exp（10*cos（f））/（2*pi*besseli（010））