为什么我这个程序无法画出图形？

[sazor]

function jrl = Jr_l(eps)
if length(eps)>1
error('epsilon must be s scale,not a vector!');
end
if eps>1
Fai_L=pi;
else
Fai_L = acos(1-2*eps); %积分限。
end
syms f
int_f = ((1-(1-cos(f))/(2*eps))^(1.1))*cos(f)/(2*pi); %被积函数
jrl =double(int(int_f,f,-Fai_L,Fai_L)); %积分
jrl=real(jrl);
上面这个程序是可运行的，输入一个变量就会出一个结果。
我想画它的曲线图
eps = 0.1:0.1:5;
J = zeros(size(eps));
for i=1:length(eps)
J(i) = Jr_l(eps(i));
end
clear i
plot(eps,J);
xlabel 'eps'
ylabel 'Jr'
却总是显示错误，请各位帮帮忙:'(

[laosam280]

函数文件的变量有作用域，函数运行以后就没有作用了。所以解决的办法有两个，
1. 设置全局变量。但是这种办法在一般的程序里不可取，因为只要你不清除变量，它一直存在；
2. 将函数文件改写为脚本文件，这样就可以了
下面的代码是我稍微修改的，有错误，但是运行的顺序没有问题。你对着你问题再把循环改写好就可以了。

for eps = 0.1:0.1:5;
J = zeros(size(eps));
if length(eps)>1
error('epsilon must be s scale,not a vector!');

if eps>1
Fai_L=pi;
else
Fai_L = acos(1-2*eps); %积分限。
end
end
syms f
int_f = ((1-(1-cos(f))/(2*eps))^(1.1))*cos(f)/(2*pi); %被积函数
jrl =double(int(int_f,f,-Fai_L,Fai_L)); %积分
jrl=real(jrl);



plot(eps,J); hold on;
end
hold off;
xlabel 'eps'
ylabel 'Jr'

[扬帆竞航]

那个function文件里面的变量是不能在其他的函数段里面用的，它就仅仅是一个代号，没有实际意义的；你可以运行你那个function文件，结果出来肯定不是jrl=……，而是ans=……
你可以将画图的函数段也集成到function函数中去，直接对数组作用，画出图像。
fuction Jr_l(eps)
J = zeros(size(eps));
if length(eps)>1
error('epsilon must be s scale,not a vector!');

if eps>1
Fai_L=pi;
else
Fai_L = acos(1-2*eps); %积分限。
end
syms f
int_f = ((1-(1-cos(f))/(2*eps))^(1.1))*cos(f)/(2*pi); %被积函数
jrl =double(int(int_f,f,-Fai_L,Fai_L)); %积分
jrl=real(jrl);

plot(eps,J); hold on;
end
hold off;
xlabel 'eps'
ylabel 'Jr'


end