如何用以下matlab程序作出相似图像

clear all;clc;figure;
k=1;c=3*10^8;
for lamada=0.15:0.001:11;
    n1=1+0.643356*lamada^2/(lamada^2-0.057789^2)+0.506762*lamada^2/(lamada^2-0.10968^2)+3.8261*lamada^2/(lamada^2-46.3864^2);% 
     n(k)=sqrt(n1);
     k=k+1;
end
lamada=0.15:0.001:11;
plot(lamada,n)
 dlmwrite('Y.dat',n','delimiter','','newline','pc','precision',8);
 dlmwrite('X.dat',lamada','delimiter','','newline','pc','precision',8);
 xlabel('v(100THz)');  
ylabel('nF2');  
title('\lambda(μm) ');