%这里开始读入数据
q=load('e:\0.1.txt');
for ii=1:length(q(:,1))

v=q(ii,1)
u=q(ii,3)

%下面开始求值
ffv=0;
for i=1:24
    eval(['ffv=ffv+eval(fv',num2str(i),');'])
end

ffu=0;
for i=1:24
    eval(['ffu=ffu+eval(fu',num2str(i),');'])
end

ff1=-ffv/ffu;
ff2=1+(u/v)^2*(1-2*v/u*ff1);

ff=eval(ff2)

hold on

plot(v,ff)
end

%这里开始读入数据
q=load('e:\0.1.txt');

ff=zeros(length(q(:,1)));
for ii=1:length(q(:,1))

v=q(ii,1)
u=q(ii,3)

%下面开始求值
ffv=0;
for i=1:24
    eval(['ffv=ffv+eval(fv',num2str(i),');'])
end

ffu=0;
for i=1:24
    eval(['ffu=ffu+eval(fu',num2str(i),');'])
end

ff1=-ffv/ffu;
ff2=1+(u/v)^2*(1-2*v/u*ff1);

ff(ii)=eval(ff2);

end

v=q(:,1);

plot(v,ff);