MATLAB: There is something wrong in R.

formula

clc;

clear all;

close all;

c=3*10.^8;

h=6.625*10.^-34;

k= 1.38*10.^-23;

T=500;

f=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]

R= ((2.*h).*(f.^3))./((c.^2).*(exp.^(h.*f/(k.*T))-1));

plot(f,R)

Best Answer

The problem is in the way you call the exp function.

This works:

R = ((2.*h).*(f.^3))./((c.^2).*(exp(h.*f/(k.*T))-1));

Related Solutions

MATLAB: How to store output of looping integral in an array

I'm guessing you have some kind of units issue, such that something you are calculating is many orders of magnitude out of whack, or that you have something wrong in your formula.

You didn't tell us what all of those constants are (e.g. k), so I just set them all to 1, just to see if I could get your code to work at all.

I got what you did ... all the values the same. But I noticed that you are integrating over a tiny interval of something around 10^-18. So, it's not too surprising that your functions are all effectively equal over that tiny range, and give the same result.

When I integrate over [Ec, Ec+5000], I see variation in the output.

So, I would just carefully check all your units and the formula itself.

MATLAB: Evaluate function handle (determine the value of function handle using the code

Defining lots and lots of function handles is not an efficient approach. One function is much simpler:

myfun(30)
ans = 5×5
         0    0.1513         0         0    0.0008
         0         0    0.0066    0.0013    0.0008
         0    0.1619         0         0    0.0007
         0    0.4139         0         0    0.0008
         0         0         0         0         0
function hf = myfun(age)
hf = zeros(5,5);
hf(1,2) = exp(-0.0625.*age-0.0134); % exp(age effect+time effect) 
hf(1,5) = exp(-9.65573+0.01844+0.08218*age+0.02246); % exp(intercept+ age effect+time effect) 

hf(2,3) = exp(-1.6660-0.1116.*age-0.0025); % exp(intercept+ age effect+time effect) 
hf(2,4) = exp(-8.96236+0.07691.*age + 0.00978); % assuming the death rate of male of same age(Hubener et al.)
hf(2,5) = exp(-9.65573+0.08218.*age+0.02246); % self-mortality 
hf(3,2) = exp(-0.0625.*age-0.0134+0.0676); %exp(intercept+ age effect+time effect+marriage once before)
hf(3,5) = exp(-9.65573+0.08218.*age+0.02246-0.11853);
hf(4,2) = exp(-0.4176-0.0625-0.0134.*age);
hf(4,5) = exp(-9.65573+0.08218.*age+0.02246-0.00415);
end

Best Answer

Related Solutions

MATLAB: How to store output of looping integral in an array

MATLAB: Evaluate function handle (determine the value of function handle using the code

Related Question