小弟求matlab程序

d=100:0.1:2000;%d是一个从100到2000,步进值为0.1的行向量
x1=35+(0.019*(d-100));%x1也是一个行向量
x2=0.15./(4*pi*d);%分子为一个数值,分母是向量,所以用点除法,指每个元素均要除
x3=1000;%x3为常数
p=x1.*x2.*x2*x3*x3;%向量相乘是点乘运算
figure(1);
subplot(1,1,1);
plot(d,p);
grid on;
已经调试好,且有详细说明