想知道用matlab写黑体辐射3个公式哪里出错了

clc

c=2.997925e8;

lambda1=[2:0.15:12,13:4:30];

lambda=lambda1*1e-6;

T=1000;

k=1.38054e-23;

h=6.626e-34;

c1=2*pi*h*c^2;

c2=h*c/k;

m1=c1*lambda.^(-5).*exp(-c2./(lambda*T));

figure

plot(lambda1,m1,'LineWidth',2)

xlabel("波长长度(um)")

ylabel("单色辐射度")

%维恩公示

lambda1=7:0.15:30;

lambda=lambda1*1e-6;

m2=2*pi*c*k*T./(lambda.^4);

hold on

plot(lambda1,m2,'LineWidth',2)

%瑞利-金斯公式

lambda1=[0.5:0.15:5,5.5:1:30];

lambda=lambda1*1e-6;

m3=2*pi*h*c^2./lambda.^5./exp(h*c./(lambda*k*T))-1;

hold on

plot(lambda1,m3,'LineWidth',2)

%普朗克公式

legend('维恩公式','瑞利-金斯公式','普朗克公式')

%用

%用频率表示

lambda=0.5:0.1:50;

mu=c./lambda;

%维恩公式

m4=2*pi*h*mu.^3/c^3.*exp(-h*mu/k/T);

%瑞利-金斯公式

m5=2*pi*k*T*mu.^2/c^3;

%普朗克公式

m6=2*pi*h*mu.^3/c^3./(exp(h*mu/k/T)-1);

figure

plot(mu,m4,'LineWidth',2)

hold on

plot(mu,m5,'LineWidth',2)

plot(mu,m6,'LineWidth',2)