怎样用MATLAB实现改变两束光振幅后的光强分布图

I_0=1:10;
d=0.02;D1=50;lamda=5e-7;
x=-67:67;
k=pid/(lamdaD1);

I=4I_0cosd(k*x).cosd(kx)

plot([0,0],[0,45],'k')
hold on
plot(x,I,'b')
hold off
text(40,42,'光强分布曲线I')
text(2,43,'y轴')现象和发生背景

报错内容：错误使用 *
内部矩阵维度必须一致。

想要改一下他人提供的代码，想要在I—0变化下得到的I。
我将I第一行I-0=10；改成I-0=1：10；就发生了报错。
请问有人知道怎么修改才正确？