[laosam280]

这样就可以了，其实程序没有什么大问题，就是少了几个end和多了几处分号

%%IGvmin is global minimum
IGvmin=1000000000000000;
%%exmperimantal values
M=1000000;
for c=0:1:100000
for omiga=0:1:100
for k=100000:100:100000000
%%igvlmax is local maximum;
IGvlmax=0;
%%IGv’s calculational equation
IGv=sqrt(16*c.^4*omiga.^4+16*k.^4+4*M.^2*k.^2*omiga.^4+32*c.^2*omiga.^2*k.^2-16*M*k.^3*omiga.^2+4*M.^2*c.^2*omiga.^6)/((-M*omiga.^2+2*k).^2+4*c.^2*omiga.^2);
if IGvmin>=IGvlmax
IGvlmax=IGvmin;
end
end
if IGvlmax<=IGvmin
IGvmin=IGvlmax;
Kmin=k; cmin=c;
end
end
end

if IGvlmin>=IGvmax
k=kmin
c=Cmin
end
for omiga=0:1:100
IGv(i)=sqrt(16*c.^4*omiga.^4+16*k.^4+4*M.^2*k.^2*omiga.^4+32*c.^2*omiga.^2*k.^2-16*M*k.^3*omiga.^2+4*M.^2*c.^2*omiga.^6)/((-M*omiga.^2+2*k).^2+4*c.^2*omiga.^2);
end
for i=0:1:100
plot(i,IGv(i))
hold on;
end