In [1]:
import math
e = math.exp(1.0)
In [2]:
e
Out[2]:
All of the functions we have written so far are void; they print something but their return value is None
.
Now, we are (finally) going to write fruitful functions (Function with return value).
In [3]:
# Function that returns the area of cirle
def area(radius):
temp = math.pi * radius**2
return temp
We have seen the return
statement before, but in a fruitful function the return
statement includes an expression. This statement means: “Return immediately from this function and use the following expression as a return value.” The expression can be arbitrarily complicated, so we could have written this function more concisely:
In [4]:
def area(radius):
return math.pi * radius**2
On the other hand, temporary variables like temp
often make debugging easier.
Sometimes it is useful to have multiple return statements, one in each branch of a conditional:
In [6]:
def absolute_value(x):
if x < 0:
return -x
else:
return x
Since these return
statements are in an alternative conditional, only one will be executed.As soon as a return
statement executes, the function terminates without executing any subsequent statements. Code that appears after a return statement, or any other place the flow of execution can never reach, is called dead code.
In [ ]:
As you write larger functions, you might find yourself spending more time debugging.
To deal with increasingly complex programs, you might want to try a process called incremental development. The goal of incremental development is to avoid long debugging sessions by adding and testing only a small amount of code at a time.
As an example, suppose you want to find the distance between two points, given by the coordinates $(x_1, y_1)$ and $(x_2, y_2)$.
$$ distance = \sqrt{(x_2 - x_1)^2 + (y_2 - y_1)^2} $$The first step is to consider what a distance
function should look like in Python. In other words, what are the inputs (parameters) and what is the output (return value)?
In this case, the inputs are two points, which you can represent using four numbers. The return value is the distance, which is a floating-point value.
Already you can write an outline of the function:
def distance(x1, y1, x2, y2):
return 0.0
Obviously, this version doesn’t compute distances; it always returns zero. But it is syntactically correct, and it runs, which means that you can test it before you make it more complicated.
We can start adding code to the body. A reasonable next step is to find the differences $x_2-x_1$ and $y_2-y_1$. The next version stores those values in temporary variables and prints them.
In [10]:
def distance(x1, y1, x2, y2):
dx = x2 - x1
dy = y2 - y1
print('dx is', dx)
print('dy is', dy)
return 0.0
In [11]:
distance(1,2,4,6)
Out[11]:
If the function is working, it should display 'dx is 3' and 'dy is 4' for values (1,2,4,6) which will result (3,4,5). If so, we know that the function is getting the right arguments and performing the first computation correctly. If not, there are only a few lines to check.
Next we compute the sum of squares of dx and dy:
In [13]:
def distance(x1, y1, x2, y2):
dx = x2 - x1
dy = y2 - y1
dsquared = dx**2 + dy**2
print('dsquared is: ', dsquared)
return 0.0
In [14]:
distance(1,2,4,6)
Out[14]:
Again, you would run the program at this stage and check the output (which should be 25). Finally, you can use math.sqrt
to compute and return the result:
In [15]:
def distance(x1, y1, x2, y2):
dx = x2 - x1
dy = y2 - y1
dsquared = dx**2 + dy**2
result = math.sqrt(dsquared)
return result
In [16]:
distance(1,2,4,6)
Out[16]:
If that works correctly, you are done. Otherwise, you might want to print the value of result
before the return statement.
The final version of the function doesn’t display anything when it runs; it only returns a value. The print
statements we wrote are useful for debugging, but once you get the function working, you should remove them. Code like that is called scaffolding because it is helpful for building the program but is not part of the final product.
When you start out, you should add only a line or two of code at a time. As you gain more experience, you might find yourself writing and debugging bigger chunks. Either way, incremental development can save you a lot of debugging time.
The key aspects of the process are:
Start with a working program and make small incremental changes. At any point, if there is an error, you should have a good idea where it is.
Use temporary variables to hold intermediate values so you can display and check them.
Once the program is working, you might want to remove some of the scaffolding or consolidate multiple statements into compound expressions, but only if it does not make the program difficult to read.
In [ ]:
As you should expect by now, you can call one function from within another. This ability is called composition.
As an example, we’ll write a function that takes two points, the center of the circle and a point on the perimeter, and computes the area of the circle.
Assume that the center point is stored in the variables xc
and yc
, and the perimeter point is in xp
and yp
. The first step is to find the radius of the circle, which is the distance between the two points. We just wrote a function, distance
, that does that:
radius = distance(xc, yc, xp, yp)
The next step is to find the area of a circle with that radius; we just wrote that, too:
result = area(radius)
Encapsulating these steps in a function, we get:
def circle_area(xc, yc, xp, yp):
radius = distance(xc, yc, xp, yp)
result = area(radius)
return result
The temporary variables radius
and result
are useful for development and debugging, but once the program is working, we can make it more concise by composing the function calls:
In [17]:
def circle_area(xc, yc, xp, yp):
return area(distance(xc, yc, xp, yp))
In [18]:
def is_divisible(x, y):
if x % y == 0:
return True
else:
return False
It is common to give boolean functions names that sound like yes/no questions; is_divisible
returns either True
or False
to indicate whether x is divisible by y.
In [19]:
is_divisible(6,4)
Out[19]:
In [20]:
is_divisible(6,3)
Out[20]:
The result of the ==
operator is a boolean, so we can write the function more concisely by returning it directly:
In [21]:
def is_divisible(x, y):
return x % y == 0
Boolean functions are often used in conditional statements:
In [24]:
x = 9
y = 14
if is_divisible(x, y):
print('x is divisible by y')
else:
print('x is not divisible by y')
In [ ]:
To give you an idea of what you can do with the tools you have learned so far, we’ll evaluate a few recursively defined mathematical functions. A recursive definition is similar to a circular definition, in the sense that the definition contains a reference to the thing being defined. A truly circular definition is not very useful:
vorpal: An adjective used to describe something that is vorpal
If you saw that definition in the dictionary, you might be annoyed. On the other hand, if you looked up the definition of the factorial function, denoted with the symbol !, you might get something like this:
$$ 0! = 1 $$$$ n! = n(n-1)! $$This definition says that the factorial of 0 is 1, and the factorial of any other value, n, is n multiplied by the factorial of n-1.
So 3! is 3 times 2!, which is 2 times 1!, which is 1 times 0!. Putting it all together, 3! equals 3 times 2 times 1 times 1, which is 6.
If you can write a recursive definition of something, you can usually write a Python program to evaluate it. The first step is to decide what the parameters should be. In this case it should be clear that factorial
takes an integer:
def factorial(n):
If the argument happens to be 0, all we have to do is return 1:
def factorial(n):
if n == 0:
return 1
Otherwise, and this is the interesting part, we have to make a recursive call to find the factorial of n-1 and then multiply it by n:
def factorial(n):
if n == 0:
return 1
else:
recurse = factorial(n-1)
result = n * recurse
return result
In [25]:
def factorial(n):
if n == 0:
return 1
else:
recurse = factorial(n-1)
result = n * recurse
return result
factorial(6)
Out[25]:
Let's define the fibonacci which is recursively defined mathematical function and one of the most common recursive introduction example:
In [27]:
def fibonacci (n):
if n == 0:
return 0
elif n == 1:
return 1
else:
return fibonacci(n-1) + fibonacci(n-2)
fibonacci(12)
Out[27]:
If you try to follow the flow of execution here, even for fairly small values of n, your head explodes. But according to the leap of faith, if you assume that the two recursive calls work correctly, then it is clear that you get the right result by adding them together.
What happens if we call factorial
and give it 1.5 as an argument?
In [28]:
fibonacci(1.5)
It looks like an infinite recursion. But how can that be? There is a base case—when n == 0
. But if n is not an integer, we can miss the base case and recurse forever.
In the first recursive call, the value of n is 0.5. In the next, it is -0.5. From there, it gets smaller (more negative), but it will never be 0.
We have two choices. We can try to generalize the factorial function to work with floating-point numbers, or we can make factorial
check the type of its argument. The first option is called the gamma function and it’s a little beyond the scope of this book. So we’ll go for the second.
We can use the built-in function isinstance
to verify the type of the argument. While we’re at it, we can also make sure the argument is positive:
In [29]:
def factorial (n):
if not isinstance(n, int):
print('Factorial is only defined for integers.')
return None
elif n < 0:
print('Factorial is not defined for negative integers.')
return None
elif n == 0:
return 1
else:
return n * factorial(n-1)
The first base case handles nonintegers; the second catches negative integers. In both cases, the program prints an error message and returns None
to indicate that something went wrong:
In [30]:
factorial("fred")
In [31]:
factorial(-2)
This program demonstrates a pattern sometimes called a guardian. The first two conditionals act as guardians, protecting the code that follows from values that might cause an error. The guardians make it possible to prove the correctness of the code.
Exercise 1: The Ackermann function, $A(m, n)$, is defined:
$$ A(m,n) =\begin{cases}n+1 & if m = 0 \\ A(m 1, 1) & if m > 0 and n = 0\\A(m 1, A(m, n 1))& if m > 0 and n > 0. \end{cases} $$Write a function named ack
that evaluates Ackermann’s function. Use your function to evaluate ack(3, 4)
, which should be 125. What happens for larger values of m and n?
In [ ]:
Exercise 2: A palindrome is a word that is spelled the same backward and forward, like “noon” and “redivider”. Recursively, a word is a palindrome if the first and last letters are the same and the middle is a palindrome.
The following are functions that take a string argument and return the first, last, and middle letters:
def first(word):
return word[0]
def last(word):
return word[-1]
def middle(word):
return word[1:-1]
palindrome.py
and test them out. What happens if you call middle
with a string with two letters? One letter? What about the empty string, which is written '' and contains no letters?is_palindrome
that takes a string argument and returns True
if it is a palindrome and False
otherwise. Remember that you can use the built-in function len
to check the length of a string.
In [ ]:
Exercise 3: A number, a, is a power of b if it is divisible by b and a/b is a power of b. Write a function called is_power
that takes parameters a and b and returns True
if a is a power of b. Note:you will have to think about the base case.
In [ ]:
Exercise 4: The greatest common divisor (GCD) of a and b is the largest number that divides both of them with no remainder.
One way to find the GCD of two numbers is based on the observation that if r is the remainder when a is divided by b, then gcd(a, b) = gcd(b, r)
. As a base case, we can use gcd(a, 0) = a
.
Write a function called gcd
that takes parameters a and b and returns their greatest common divisor.
In [ ]: