还是不行，也许是函数调用的问题？

有可能是函数调用的问题。在MATLAB中，可以使用函数句柄来调用函数，而不必使用`feval`命令。请尝试以下修改后的代码：

f=@(t) exp(-(t.^2)./2)./sqrt(2.*pi);
PHI=zeros(1,8);
CumDis=zeros(1,8);
PHI(1) = 0.5;
for i = 2:8
    PHI(i) = PHI(i-1)+0.5;
end
CumDis=0.5+CumDis; % 输出累积分布函数值结果
for i = 1:8
    [R,quad,err,h]=romber(f, 0, PHI(i), i , 1e-8);
    CumDis(i) = CumDis(i) + quad;
end

function [R,quad,err,h]=romber(f,a,b,n,tol)
M = 1;
h = b-a;
err = 1;
J = 0;
R = zeros(4, 4);
R(1, 1) = h*(f(a)+f(b))/2;
while((err>tol)&(J<n))|(J<4)
    J = J + 1;
    h = h/2;
    s = 0;
    for p=1:M
        x=a+h*(2*p-1);
        s=s+f(x);
    end
    R(J+1,1)=R(J,1)/2+h*s;
    M=2*M;
    for K=1:J
        R(J+1,K+1)=R(J+1,K) +(R(J+1,K)-R(J,K))/(4*K-1);
    end
    err=abs(R(J,J)-R(J+1,K+1));
end
quad=R(J+1,J+1);
end

如果这个还是不行的话，您可以提供一下完整的错误提示信息，我再帮您看看。