怎麽用 matlab畫出這個公式的圖，求編程過程

clear

cc=3.0e8;

lambda0=[632.8 488.0 441.6 510.5].*1.0e-9;

nu0=cc./lambda0;

Dnu_D=[1.5,2.7,1.1,2.3].*1.0e9;

nu=linspace(0.1,9)*1.0e14;

st=['--';'.-';'o-';'*-'];

figure(1)

for i=1:4

g_D=log(2*sqrt(log(2)/pi)/Dnu_D(i))+(-4*log(2)*(nu/Dnu_D(i)-nu0(i)/Dnu_D(i)).^2);

H(i)=plot(nu,g_D,st(i,:));

hold on

end

xlabel('\nu(nm)')

ylabel('ln(g_D)')

legend(H,'\lambda_0=632.8nm','\lambda_0=488.0nm','\lambda_0=441.6nm','\lambda_0=510.5nm','location','northwest')