Plotting the stellar IMF and setting for loop conditions

2 次查看（过去 30 天）

HC98 2020-1-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/502957-plotting-the-stellar-imf-and-setting-for-loop-conditions

Hi all, I'm trying to plot the stellar IMF using MATLAB and want to transfer the following Pyhon code:

m = logspace(-2, 2, 400);
def chabrier03individual(m):
	k = 0.158 * exp(-(-log(0.08))**2/(2 * 0.69**2))
	return numpy.where(m <= 1,\
	        0.158*(1./m) * exp(-(log(m)-log(0.08))**2/(2 * 0.69**2)),\
        	        k*m**-2.3)
                 
            

This is my attempt:

                 %% chabrier03individual
k1=0.158*exp(-(-log(0.08))^2/(2 * 0.69^2))
for j=m
       u=0.158*(1./m).*exp(-(log(m)-log(0.08)).^2/(2 * 0.69^2));
end 
for m=1
        v = k1*m^(-2.3);
 end 