Skip to main content

Factorial of a number using Stirling's(Gosper) Approximation using Matlab

Problem

Find the factorial of a number using Modified Stirling's formula(Gospers). Also print the difference upto 5 decimals between the result from the factorial function in Matlab and the result you get from approximations.The relation to find the factorial of a number using the Stirling's formula is given by:
Modified Stirlings approximation(Gospers formula) using Matlab
Gospers formula(Modified Stirling formula)
To know more about Stirling's formula or Gospers formula then go to: Stirling's Approximation - Math.Wolfram

Solution

We will solve this problem using Matlab functions. This is similar to our previous post Velocity of a moving fluid using Matlab. But the little difference between the previous post and this post is that we are required to print the difference upto five decimals of the difference between the factorial function in matlab and the above formula.
It should be remembered that we are asked indirectly to print the modulus of the difference. To know about modulus of a number then please do check out mathcenter website
Program is written in a simple manner so that it is easy to understand, also explanation section is written to make a detailed explanation of the program. Please do check out explanation section if you want a detailed explanation of the program. To solve this problem some of the concepts you need to know are as follows:
1) Matlab functions
2) Matlab comments
3) Matlab Variables
4) factorial of a number
5) abs function
6) fprintf in Matlab

Program

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
function gosper = stirlings(n)
%This function will calculate the factorial of a
%number using the modified stirlings formula or 
%gospers formula
%Use this function as follows:
%>> stirlings(n)
%Where n is the number you want the value of the factorial

gosper = sqrt((2*n+(1/3))*pi)*((n^n)*exp(-n));

real = factorial(n);

difference = abs(gosper - real);

fprintf('The difference between real and gosper approximation is: \n');
fprintf('%.5f\n',difference);
end
I have saved this file as stirlings.m on my local computer.

Explanation

Line 1: We have initiated the function with the return value gosper and the name as stirlings whose parameter(s) are only a number.
Line 2-7: They are matlab comments and they can also be called as the documentation of the function. You can get the documentation of the function by typing the following in the command window:

>> help stirlings
And the result will be as follows:
factorial of a number using stirlings approximation using Matlab
Stirlings approximation using Matlab
Line 9: We are creating a variable named gosper and storing the value of the matlab equivalent of the given equation in the question. To know why I have made many parenthesis around each number then you need to understand about Matlab operator precedence.
Line11: I am calculating the value of the factorial using the factorial function of the Matlab.
Line 13: The difference between the value of the factorial function and the value of the value calculated from the gospers formula is calculated. I have also added the modulus of to the number so that we only get a positive value of the number. Read about the abs function in Matlab to know more about the modulus of a number using matlab.
Line15-16 the difference is printed using fprintf function. 

Input/Output

Difference between the factorial function and stirlings formula output
Modified Stirlings approximation using Matlab

Try it yourself

1)Write a program to ask the user to give two options. Option 1 stating that the value of the factorial is calculated using unmodified stirlings formula and Option 2 using modified stirlings formula. According to the user input calculate the same. Stirling formula. Stirlings formula is as follows:
Stirlings approximation formula
Stirlings formula for factorial of a number
Comment your program so that it will be useful for others.

Final Note

I have tried to explain every part in this post so that it is easy for everyone to understand. If you have any doubt or couldn't understand anything then please do comment below or do contact me. Contact me from here: contact me.
The above program was high lighted using hilite.me.
Keywords: stirlings, gosper, approximation, factorial, abs function
Labels:

Popular posts from this blog

Project Euler Problem 67 Solution with Python

Maximum path sum II By starting at the top of the triangle below and moving to adjacent numbers on the row below, the maximum total from top to bottom is 23. 3 7 4 2 4 6 8 5 9 3 That is, 3 + 7 + 4 + 9 = 23. Find the maximum total from top to bottom in triangle.txt (right click and 'Save Link/Target As...'), a 15K text file containing a triangle with one-hundred rows.
Read more

Problem 43 Project Euler Solution with python

Sub-string divisibility The number, 1406357289, is a 0 to 9 pandigital number because it is made up of each of the digits 0 to 9 in some order, but it also has a rather interesting sub-string divisibility property. Let d 1 be the 1 st digit, d 2 be the 2 nd digit, and so on. In this way, we note the following: d 2 d 3 d 4 =406 is divisible by 2 d 3 d 4 d 5 =063 is divisible by 3 d 4 d 5 d 6 =635 is divisible by 5 d 5 d 6 d 7 =357 is divisible by 7 d 6 d 7 d 8 =572 is divisible by 11 d 7 d 8 d 9 =728 is divisible by 13 d 8 d 9 d 10 =289 is divisible by 17 Find the sum of all 0 to 9 pandigital numbers with this property. One might write a simple solution using direct if else statements for this problem. Using if else statements, the execution time may be a few seconds. But this is not a very good approach. I too had written a program with if else statement which will take each and every permutation of the 0-9 Pandigital and check for the conditions given in the qu...
Read more

Add/Embed SVG to Blogger website

In this post I will tell you my method(trick) of adding SVG images in a blogger website or blog. Before starting , the first thin g I am assu m ing is that you are aware of SVG if you are here. If not please see S calable V ec tor G raphics Recently when I tried to embed a SVG image for a post on pygal, I tried uploading the SVG file and blogger Image uploader came up with an error, because of which I had to find some other way.  SVG File upload Error in Blogger  I started sea rc hing Google " Embed SVG in Blogger " . I found blogorrhea , w h ich gave some i nformatio n on add ing SVG directly as a markup , which worked , but I faced another problem using this . Also th is guy has used lot of Javascript which was confusin g for me, being new to using SVG.   So I first t houg ht of learning on h ow to embed SVG in HTML and t his on e worked out. Actually we can embed SVG in HTML i n following ways: Using Object tag Using Iframe tag Using embed...
Read more