我的代码如下
function x=Problem1()
clear all
%  把函数Problem_1()文件和ProblemData_1数据保存放在桌面或者指定的matlab搜索文件夹(Current Folder设置)下
%  在Command Window里面输入x=Problem_1()即可得到x
clc
load 2003
x=zeros(1,length(P));
%  初始值
r0=0;
for i=1:length(P)
    x(i)=fsolve(@(r)myfun(r,P(i),B(i,:),FROM(i,:)),r0);   
end
x=x';
end

%  求解的r的函数
function y = myfun(r,P,B,FROM)
F=B(1);
for i=1:12
        if i<12
                F=F+(FROM(i+1)-r)*B(i)/((1+r)^i);
        else
            F=F+(FROM(13)-r)*B(12)/(r*((1+r)^11));
        end
end
y=P-F;
%  求解y等于0
end

function x=Problem1()
clear all
%  把函数Problem_1()文件和ProblemData_1数据保存放在桌面或者指定的matlab搜索文件夹(Current Folder设置)下
%  在Command Window里面输入x=Problem_1()即可得到x
clc
load 2003
x=zeros(1,length(P));
%  初始值
r0=0;
for i=1:length(P)
    x(i)=fsolve(@(r)myfun(r,P(i),B(i,:),FROM(i,:)),r0);   
end
x=x';
end

%  求解的r的函数
function y = myfun(r,P,B,FROM)
F=B(1);
for i=1:12
        if i<12
                F=F+(FROM(i+1)-r)*B(i)/((1+r)^i);
        else
            F=F+(FROM(13)-r)*B(12)/(r*((1+r)^11));
        end
end
y=P-F;
%  求解y等于0
end

这是数据和m文件，报错情况如下图

FROM就是公式里的FROE，我没有改过来。