Question

1. Write a code to calculate nas a function of ga, when NA l300, 600, 3000, 6000, NB (200, 400, 2000, 4000), and g 100, plot them and try to find some tendency when N increases (hint: 4 plots). 2. Write a code to calculate the probability of S2(ga), when Na 1200, 300, NB (2000, 3000), for gr (100, 1000), plot them and try to explain the differences. (hint: 2 plots)

Can someone help me figure out how to write code on Matlab for either of these two questions? My biggest problem is working with the large factorials which Matlab can't seem to handle

Answer 1

Us the following function to compute the factorials of very large numbers.

function [m,e] = mylargefact(n)

f = gammaln(N+1)/log(10);
e = floor(f);
m = 10.^(f-e);

The following method also work for large factorials

syms n
subs(n*factorial(n-1),1000)