index must be a positive integer or logical,这个问题出在哪

运行的时候出现：Attempted to access f(0.01,8,0.62,0.000492407); index must be a positive integer or logical.
不知道问题出在哪里，用clear和clc清屏之后，单独运行f(0.01,8,0.62,0.000492407)可以得到答案-0.5787，请各位大神帮忙看看这问题是出在哪里了
主程序
b1=0.006;
r=0.62
R=0.5;
a=8
n=0.25;
m=0.03;
f=0.997;
faiA=fai(R,a,r);
rouA=rou(R,a,r);
beitaA=beita(R,rouA,faiA);
a1=0.001;
A=a1*faiA;
a2=a1*rouA;
b2=a2;
rouB=b2/b1；
RB=0.01
f1=f(RB,a,r,rouB)
下面是我的m文件中的函数
1.function beita1=beita(R,rou,fai)
beita1=(1+R*rou)/fai;
2.function fai1=fai(R,a,r)
fai1=((R^a+1+(1-R)^a)/(r+1))^(1/a);
3.function rou1=rou(R,a,r)
rou1=(R^(a-1)-r*(1-R)^(a-1))/(1+r*(1-R)^(a-1));
4.function f1=f(R,a,r,rouB)
f1=R^(a-1)-r*(1-R)^(a-1)-rouB*(1+r*(1-R)^(a-1));