How to find factorial(c) of an array of numbers?

Respected all

I have successfully got the values for j and v. I have got stuck in solving t due to the problem in finding factorial(c). My question is how to code for factorial term in the expression? Is it possible or I am on the wrong path? As without solving factorial, t cannot be determined. What should I do?

Code attached:

A = linspace(1,10)

C = linspace(1, 11)

[k, c] = ndgrid(A, C)

j=gammainc(c,(0.85));

v=gamma(c.^(-1));

t=0.85.*(1-(((0.85).^5.*(1-(((0.85))./(1-(((0.85)).^6)))))).*1.*(1-((0.85).^5.*(factorial(c)).^-1.*(c.^(5-c).^-1)./(exp(0.85).*k.*v+(0.85).^c.*(1-(0.85./c).^(5-c+1)))/factorial(c).*(1-(0.85./c)))));

The following error is being encountered

Error using factorial

N must be an array of real non-negative integers.

How can this error be removed?

Please help

Regards

Surabhi

Answer 1

[Adam Danz]

N must be an array of real non-negative integers.

To troubleshoot that error message, let's go through the list of possibilities.

% Is c an array or real numbers?
all(isreal(c(:)))  % = 1, true

% Is c all non-negative?
all(c(:)>=0)  % = 1, true

% Is c all integers?
all(isinteger(c(:)))  % = 0!  FALSE

There's the problem. C contains non-integers. It makes sense that this would cause an error, right? Factorials are functions of intergers.

factorial(3) = 24

4 * 3 * 2 * 1 = 24

How could that be done with non-integers?

[surabhi sachdeva]

Yes sir you are right But request you to please suggest me can something be done to evaluate fact(c) because factorial c is required to solve the equation

Please help

Regards

[Adam Danz]

There is no way to evaluate factorial(c) if c contains non-integers. It's not mathematically possible. It doesn't even make sense.

If you're certain that the equation you wrote is correct, then your inputs are wrong.

c must contain only integers > 0.

C = linspace(1, 11)
[k, c] = ndgrid(A, C)
%   ^ This must contain ONLY integers >0

Note that little-c is based on the values of big-C. So, make sure big-C is only integers.

[surabhi sachdeva]

Sir,

Actually my code was for 2D representation as below.

c=2:1:11; [It denotes c value to be from 2 to 11 but t is not using all the different values of c]

j=gammainc(c,(0.85));

v=gamma(c.^(-1));

t=0.85.*(1-(((0.85).^5.*(1-(((0.85))./(1-(((0.85)).^6)))))).*1.*(1-((0.85).^5.*(factorial(c)).^-1.*(c.^(5-c).^-1)./(exp(0.85).*k.*v+(0.85).^c.*(1-(0.85./c).^(5-c+1)))/factorial(c).*(1-(0.85./c)))));

Also, I was trying to change the above in 3D format, for which all the below lines worked out well except for fact(c)

A = linspace(1,10)

C = linspace(1, 11)

[k, c] = ndgrid(A, C)

j=gammainc(c,(0.85));

v=gamma(c.^(-1));

t=0.85.*(1-(((0.85).^5.*(1-(((0.85))./(1-(((0.85)).^6)))))).*1.*(1-((0.85).^5.*(factorial(c)).^-1.*(c.^(5-c).^-1)./(exp(0.85).*k.*v+(0.85).^c.*(1-(0.85./c).^(5-c+1)))/factorial(c).*(1-(0.85./c)))));

Is there a way out for this? If you could help me in solving t.

Regrds

[Adam Danz]

I didn't understand, "but t is not using all the different values of c"

What's wrong with this?

A = 1:10;
C = 2:1:11;
[k, c] = ndgrid(A, C)

[surabhi sachdeva]

It is showing NaN result.

[Adam Danz]

There are many moving parts in that very long equation. The NaNs are not being caused by the factorial() output.

There are some key points to understand about this new problem.

Understand these points before solving the new problem.

1) There's probably several ways to change the equation so that non-NaN values are returned but that doesn't mean the changes that were made are correct. For example, (1+inf)/(1-inf) returns a NaN. I could replace either of the 2 inf values with a non-inf number or I could replace both of them and the result will not be nan but I don't know what the equation is supposed to be.

2) The problem could be that the t equation is fine but the inputs are the wrong size or shape. For example, if the inputs are scalar instead of arrays, the currect equation works!

k = 1;
c = 2; 
j = 0.896; 
v = 1.7725; 
t=0.85.*(1-(((0.85).^5.*(1-(((0.85))./(1-(((0.85)).^6)))))).*1.*(1-((0.85).^5.*(factorial(c)).^-1.*(c.^(5-c).^-1)./(exp(0.85).*k.*v+(0.85).^c.*(1-(0.85./c).^(5-c+1)))/factorial(c).*(1-(0.85./c)))))

% RESULT:
% t =  0.98732

3) The problem could be that the t equation was coded incorrectly. I have a feeling this is the problem. Specifically, I think some of the array operations should be element-wise operations instead of matrix/vector operations. To understand this difference, read this page of the documentation. <-- don't skip this part, it's the most important part of the entire answer.

For example, adding this period results in non-NaN output when inputs are arrays.

t=0.85.*(1-(((0.85).^5.*(1-(((0.85))./(1-(((0.85)).^6)))))).*1.*(1-((0.85).^5.*(factorial(c)).^-1.*(c.^(5-c).^-1)./(exp(0.85).*k.*v+(0.85).^c.*(1-(0.85./c).^(5-c+1)))./factorial(c).*(1-(0.85./c)))))
% LOOK HERE ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> ---> --> ^ 

But that doesn't mean the this is the correct solution. That's up to you to think about.

[surabhi sachdeva]

you explained so well

thank you so much sir for your guidance

Regards

Surabhi

[Adam Danz]

Glad I could help. If you get stuck again with this topic, leave a comment to loop me back in.

[surabhi sachdeva]

Thank you sir

Regards