求教

魅力薄膜

本人根据矩阵传输理论 匆忙学习了下MATLAB 想能够快速设计CLENS 但是在三维矩阵和循环和画图 那里被卡住了 求高手指点。源程序如下clc;clear;
l=input('请输入CLENS长度L：');
n=input('请输入CLENS折射率N：');
r=input('请输入CLENS曲率R：');
t=[1,l/n;(1-n)/r,[(1-n)/r*l+n]/n];
f=-1/t(2,1)
sw=-t(2,2)/t(2,1)
sx=-1/t(2,1)
w0=input('请输入入射光斑半径W0(uM)：');
w0=w0/1000;
lambda=input('请输入波长uM：');
lambda=lambda/1000;
q0=(pi*w0^2/lambda)*i
endcap=input('请输入ENDCAP的长度(mm)：');
b=input('请输入ENDCAP与clens之间的距离(mm):');

d=input('请输入离透镜中心点距离(MM)：');
if endcap==0
    m=[1,d;0,1]*t*[1,b;0,1]
else
     m=[1,d;0,1]*t*[1,b;0,1]*[1,0;0,1.45]*[1,endcap;0,1]
end

q=(m(1,1)*q0+m(1,2))/(m(2,1)*q0+m(2,2))
p=imag(1/q)
w=1/sqrt(-pi*p/lambda)
希望高手能指点下 怎样让D从1:1:100 让所有的D都能对应一个W   B从0:0.05:0.5都能对应一个d