能不能给个matlab程序，大概编个一百多条左右。

a=[0.150,0.025,0.150,0.025,0.300,0.025,0.300,0.025,0.200,0.025,0.200,0.025,0.130;
　　0,0.025,0.025,0.025,0.025,0.025,0.025,0.025,0.025,0.025,0.025,0.025,0.025];
　　m=[153.15264,0,0,0,0];
　　n=[7.85398*10^9/0.2,0,0,0,0];
　　for i=1:13
　　b(1,i)=polyval(m,a(1,i));
　　b(2,i)=polyval(n,a(2,i));
　　end

　　for w=1:1:600;
　　c=[1;-0.077494*w^2];
　　e(1,1)=c(1,1);
　　e(2,1)=c(2,1);
　　for i=2:13
　　j=b(1,i);
　　k=b(2,i);
　　c=[1,1/k;-j*w^2,1-j*w^2/k]*c;
　　end
　　e(1,w)=w;
　　e(2,w)=c(2,1);
　　end

　　j=1;
　　for w=1:1:599
　　u=w+1;
　　r=e(2,w)*e(2,u);
　　if r<0
　　h(j)=w;
　　j=j+1;
　　end
　　end
　　h

　　for j=1:1:3
　　d=-1;f=1;w=h(j);
　　while f>=0
　　c=[1;-0.077494*w^2];
　　for i=2:13
　　j=b(1,i);
　　k=b(2,i);
　　c=[1,1/k;-j*w^2,1-j*w^2/k]*c;
　　end
　　e=c(2,1);
　　f=e*d;
　　d=c(2,1);
　　w=w+0.001;