r1=0:0.001:1; p=(3.*r1.^2)./(3.*r1.^2+2*0.3*(1-r1)); plot(r1,p)这样就对了你的程序分子少了个点，所以算出来的p只是一个数，而不是一个一维矩阵。
