Chapter 22

Jupyter Notebook

The Jupyter Notebook is an open-source web application that allows you to cre-
ate and share documents that contain live code, equations, visualizations and
text.

The Notebook has support for over 40 programming languages, including Python.

Figure 22.1: Jupyter Notebook [16]

Web:
http://jupyter.org

Wikipedia:
https://en.wikipedia.org/wiki/ProjectJ upyter

110
22.1 JupyterHub
JupyterHub is a multi-user version of the notebook designed for companies,
classrooms and research labs [17].

JupyterHub runs in the cloud or on your own hardware.

JupyterHub is open-source and designed to be run on a variety of infrastructure.

This includes commercial cloud providers, virtual machines, or even your own
laptop hardware.

Web:
http://jupyter.org/hub

22.2 Microsoft Azure Notebooks

Microsoft Azure Notebooks is a version of Jupyter Notebook from Microsoft.

The good thing about Microsoft Azure Notebooks is that you have the infras-
tructure and everything up and running ready for you to use. You can use it
for free as well.

Web:
https://notebooks.azure.com

Example 22.2.1. Example Name

Figure 22.2 shows an overview of my Azure Notebook Projects.

Figure 22.2: Azure Notebook Projects

Figure 22.3 shows an overview of my Azure Notebook Project Notebooks.

Figure 22.4 shows an example of a simple Notebook.

[End of Example]

111
Figure 22.3: Azure Notebook Project Notebooks

Figure 22.4: Azure Notebook Example

112
 Part V

Python for Mathematics

Applications

113
 Chapter 23

Mathematics in Python

Python is a powerful tool for mathematical calculations.

If you are looking for similar using MATLAB, please take a look at these re-
sources:
https://www.halvorsen.blog/documents/programming/matlab/

23.1 Basic Math Functions

The Python Standard Library consists of different modules for handling file
I/O, basic mathematics, etc. You don’t need to install these separately, but you
need to important them when you want to use some of these modules or some
of the functions within these modules.

In this chapter we will focus on the math module that is part of the Python
Standard Library.

The math module has all the basic math functions you need, such as: Trigono-
metric functions: sin(x), cos(x), etc. Logarithmic functions: log(), log10(), etc.
Constants like pi, e, inf, nan, etc. etc.

Example 23.1.1. Using the math module

We create some basic examples how to use a Library, a Package or a Module:

If we need only the sin() function we can do like this:

1 from math im por t s i n
2
3 x = 3.14
4 y = sin (x)
5
6 print (y)

If we need a few functions we can do like this

114
1 from math im por t s i n , c o s
2
3 x = 3.14
4 y = sin (x)
5 print (y)
6
7 y = cos (x)
8 print (y)

If we need many functions we can do like this:

1 from math im por t ∗
2
3 x = 3.14
4 y = sin (x)
5 print (y)
6
7 y = cos (x)
8 print (y)

We can also use this alternative:

1 im po rt math
2
3 x = 3.14
4 y = math . s i n ( x )
5
6 print (y)

We can also write it like this:

1 im po rt math a s mt
2
3 x = 3.14
4 y = mt . s i n ( x )
5
6 print (y)

[End of Example]

There are advantages and disadvantages with the different approaches. In your
program you may need to use functions from many different modules or pack-
ages. If you import the whole module instead of just the function(s) you need
you use more of the computer memory.

Very often we also need to import and use multiple libraries where the different
libraries have some functions with the same name but different use.

Other useful modules in the Python Standard Library are statistics (where
you have functions like mean(), stdev(), etc.)

For more information about the functions in the Python Standard Library,
see:
https://docs.python.org/3/library/

115
23.1.1 Exercises
Below you find different self-paced Exercises that you should go through and
solve on your own. The only way to learn Python is to do lots of Exercises!

Exercise 23.1.1. Create Mathematical Expressions in Python

Create a function that calculates the following mathematical expression:
p
z = 3x2 + x2 + y 2 + eln (x) (23.1)
Test with different values for x and y.

[End of Exercise]

Exercise 23.1.2. Create advanced Mathematical Expressions in Python

Create the following expression in Python:

ln (ax2 + bx + c) − sin(ax2 + bx + c)
f (x) = (23.2)
4πx2 + cos(x − 2)(ax2 + bx + c)

Given a = 1, b = 3, c = 5 Find f (9)

(The answer should be f (9) = 0.0044)

Tip! You should split the expressions into different parts, such as:

poly = ax2 + bx + c

num = . . .
den = . . .
f = ...

This makes the expression simpler to read and understand, and you minimize
the risk of making an error while typing the expression in Python.

When you got the correct answer try to change to, e.g., a = 2, b = 8, c = 6

Find f (9)

[End of Exercise]

Exercise 23.1.3. Pythagoras

116
 Figure 23.1: Right-angled triangle

Pythagoras theorem is as follows:

c2 = a2 + b2 (23.3)

Create a function that uses Pythagoras to calculate the hypotenuse of a right-

angled triangle (Figure 23.1), e.g.:

1 def pythagoras (a , b)
2 ...
3 ...
4 return c

[End of Exercise]

Exercise 23.1.4. Albert Einstein

Given the famous equation from Albert Einstein:

E = mc2 (23.4)

The sun radiates 385x1024 J/s of energy.

Calculate how much of the mass on the sun is used to create this energy per day.

How many years will it take to convert all the mass of the sun completely? Do
we need to worry if the sun will be used up in our generation or the next? justify
the answer.

The mass of the sun is 2x1030 kg.

117
 [End of Exercise]

Exercise 23.1.5. Cylinder Surface Area

Create a function that finds the surface area of a cylinder based on the height
(h) and the radius (r) of the cylinder. See Figure ??.

Figure 23.2: cylinder

[End of Exercise]

23.2 Statistics
23.2.1 Introduction to Statistics
Mean or average:
The mean is the sum of the data divided by the number of data points. It is
commonly called “the average”,

Formula for mean:

N
x1 + x2 + x3 + ... + xN 1 X
x̄ = = xi (23.5)
N N i=1

Example 23.2.1. Mean

Given the following dataset: 2.2, 4.5, 6.2, 3.6, 2.6

Mean:
N
1 X 2.2 + 4.5 + 6.2 + 3.6 + 2.6 19.1
x̄ = xi = = = 3.82 (23.6)
N i=1 5 5

118
 [End of Example]

Variance:

Variance is a measure of the variation in a data set.

N
1 X
var(x) = (xi − x̄)2 (23.7)
N i=1
Standard deviation:
The standard deviation is a measure of the spread of the values in a dataset
or the value of a random variable. It is defined as the square root of the variance.

v
u N
√ u1 X
std(x) = σ = var = t (xi − x̄)2 (23.8)
N i=1

We typically use the symbol σ for standard deviation.

We have that σ 2 = var(x)

23.2.2 Statistics functions in Python

Mathematical statistics functions in Python:
https://docs.python.org/3/library/statistics.html

statistics is part of the The Python Standard Library.

For more information about the functions in the Python Standard Library,
see:
https://docs.python.org/3/library/

Example 23.2.2. Statistics using the statistics module in Python Standard

Library
Below you find some examples how to use some of the statistics functions in the
statistics module in Python Standard Library:

1 im po rt s t a t i s t i c s a s s t
2
3 data = [ − 1 . 0 , 2 . 5 , 3 . 2 5 , 5 . 7 5 ]
4
5 #Mean o r Average
6 m = s t . mean ( data )
7 p r i n t (m)
8
9 # Standard D e v i a t i o n
10 s t d e v = s t . s t d e v ( data )

119
11 print ( st dev )
12
13 # Median
14 med = s t . median ( data )
15 p r i n t ( med )
16
17 # Variance
18 v a r = s t . v a r i a n c e ( data )
19 p r i n t ( var )
Listing 23.1: Statistics functions in Python

[End of Example]

IMPORTANT: Do not name your file ”statistics.py” since the import will be
confused and throw the errors of the library not existing and the mean function
not existing.

You can also use the NumPy Library. NumPy is the fundamental package for
scientific computing with Python.

Here you find an overview of the NumPy library:

http://www.numpy.org

Example 23.2.3. Statistics using the NumPy Library

Below you find some examples how to use some of the statistics functions in
NumPy:

1 im po rt numpy a s np
2
3 data = [ − 1 . 0 , 2 . 5 , 3 . 2 5 , 5 . 7 5 ]
4
5 #Mean o r Average
6 m = np . mean ( data )
7 p r i n t (m)
8
9 # Standard D e v i a t i o n
10 s t d e v = np . s t d ( data )
11 print ( st dev )
12
13 # Median
14 med = np . median ( data )
15 p r i n t ( med )
16
17 # Minimum Value
18 minv = np . min ( data )
19 p r i n t ( minv )
20
21 # Maxumum Value
22 maxv = np . max( data )
23 p r i n t ( maxv )
Listing 23.2: Statistics using the NumPy Library

120
 [End of Example]

Exercise 23.2.1. Create your own Statistics Module in Python

Using the built-in functions in the Python Standard Library or the NumPy li-
brary is straightforward.

In order to get a deeper understanding of the mathematics behind these func-

tions and to learn more Python programming, you should create your own
Statistics Module in Python.

Create your own Statistics Module in Python (e.g., ”mystatistics.py) and then
create a Python Script (e.g., ”testmystatistics.py) where you test these func-
tions.

You should at least implement functions for mean, variance, standard deviation,
minimum and maximum.

[End of Exercise]

23.3 Trigonometric Functions

Python offers lots of Trigonometric functions, e.g., sin, cos, tan, etc.

Note! Most of the trigonometric functions require that the angle is expressed in
radians.
Example 23.3.1. Trigonometric Functions in Math module

1 im po rt math a s mt
2
3 x = 2∗mt . p i
4
5 y = mt . s i n ( x )
6 print (y)
7
8 y = mt . c o s ( x )
9 print (y)
10
11 y = mt . tan ( x )
12 print (y)
Listing 23.3: Trigonometric Functions in Math module

Here we have used the Math module in the Python Standard Library.

For more information about the functions in the Python Standard Library,
see:
https://docs.python.org/3/library/index.html

121
 [End of Example]

Example 23.3.2. Plotting Trigonometric Functions

In the example above we used some of the trigonometric functiosn in basic cal-
culations.

Lets see if we are able to plot these functions.

1 im po rt math a s mt
2 im po rt m a t p l o t l i b . p y p l o t a s p l t
3
4 xdata = [ ]
5 ydata = [ ]
6
7 f o r x in range (0 , 10) :
8 xdata . append ( x )
9 y = mt . s i n ( x )
10 ydata . append ( y )
11
12 p l t . p l o t ( xdata , ydata )
13 p l t . show ( )
Listing 23.4: Plotting Trigonometric Functions

In the example we have plotted sin(x), we can easily extend the program to plot
cos(x), etc.

For more information about the functions in the Python Standard Library,
see:
https://docs.python.org/3/library/index.html

[End of Example]

Example 23.3.3. Trigonometric Functions using NumPy

The problem with using the Trigonometric functions in the the Math module
from the Python Standard Library is that they don’t handle an array as input.

We will use the NumPy library instead because they handle arrays, in addition
to all the handy functionality in the NumPy library.

1 im po rt numpy a s np
2 im po rt m a t p l o t l i b . p y p l o t a s p l t
3
4 xstart = 0
5 x s t o p = 2∗ np . p i
6 increment = 0.1
7
8 x = np . a r a n g e ( x s t a r t , xstop , i n c r e m e n t )
9
10 y = np . s i n ( x )

122
11 plt . plot (x , y)
12 plt . t i t l e ( ’ y=s i n ( x ) ’ )
13 plt . xlabel ( ’x ’ )
14 plt . ylabel ( ’y ’ )
15 plt . grid ()
16 plt . a x i s ( [ 0 , 2∗ np . pi , −1, 1 ] )
17 plt . show ( )
18
19 y = np . c o s ( x )
20 plt . plot (x , y)
21 p l t . t i t l e ( ’ y=c o s ( x ) ’ )
22 plt . xlabel ( ’x ’ )
23 plt . ylabel ( ’y ’ )
24 plt . grid ()
25 p l t . a x i s ( [ 0 , 2∗ np . pi , −1, 1 ] )
26 p l t . show ( )
27
28 y = np . tan ( x )
29 plt . plot (x , y)
30 p l t . t i t l e ( ’ y=tan ( x ) ’ )
31 plt . xlabel ( ’x ’ )
32 plt . ylabel ( ’y ’ )
33 plt . grid ()
34 p l t . a x i s ( [ 0 , 2∗ np . pi , −1, 1 ] )
35 p l t . show ( )
Listing 23.5: Trigonometric Functions using NumPy

This Python script gives the plots as shown in Figure 23.3.

[End of Example]

Exercise 23.3.1. Create Python functions for converting between radians an

degrees
Since most of the trigonometric functions require that the angle is expressed in
radians, we will create our own functions in order to convert between radians
and degrees.

It is quite easy to convert from radians to degrees or from degrees to radians.

We have that:

2π[radians] = 360[degrees] (23.9)

This gives:
180
d[degrees] = r[radians] × ( ) (23.10)
π
and
π
r[radians] = d[degrees] × (
) (23.11)
180
Create two functions that convert from radians to degrees (r2d(x)) and from
degrees to radians (d2r(x)) respectively.

123
 These functions should be saved in one Python file .py.

Test the functions to make sure that they work as expected.

[End of Exercise]

Exercise 23.3.2. Trigonometric functions on right triangle

Given right triangle as shown in Figure 23.4.

Create a function that finds the angle A (in degrees) based on input arguments
(a,c), (b,c) and (a,b) respectively.

Use, e.g., a third input “type” to define the different types above.

Use you previous function r2d() to make sure the output of your function is in
degrees and not in radians.

Test the function to make sure it works properly.

Tip! We have that:

a a
sin(A) = → A = arcsin( ) (23.12)
c c
b b
cos(A) = → A = arccos( ) (23.13)
c c
a a
tan(A) = → A = arctan( ) (23.14)
b b

We may also need to use the Pythagoras’ theorem:

c2 = a2 + b2 (23.15)

1 >>> a=5
2 >>> b=8
3 >>> c = s q r t ( a ∗∗2 + b ∗ ∗ 2 )
4
5 >>> A = r i g h t t r i a n g l e ( a , c , ’ s i n ’ )
6 A =
7 32.0054
8
9 >>> A = r i g h t t r i a n g l e ( b , c , ’ c o s ’ )
10 A =
11 32.0054
12 >>> A = r i g h t t r i a n g l e ( a , b , ’ tan ’ )
13 A =
14 32.0054

We also see that the answer in this case is the same, which is expected.

124
 [End of Exercise]

Exercise 23.3.3. Law of Cosines

Given the triangle as shown in Figure 23.5.

Create a function where you find c using the law of cosines.

c2 = a2 + b2 − 2ab cos(C) (23.16)

Test the functions to make sure it works properly.

[End of Exercise]

Exercise 23.3.4. Plotting Trigonometric Functions

Plot sin(θ) and cos(θ) for 0 ≤ θ ≤ 2π in the same plot (both in the same plot
and in 2 different subplots).

Make sure to add labels and a legend and use different line styles and colors for
the plots.

[End of Exercise]

23.4 Polynomials
A polynomial is expressed as:

p(x) = p1 xn + p2 xn−1 + ... + pn x + pn+1 (23.17)

where p1 , p2 , p3 , ... are the coefficients of the polynomial.

We will use the Polynomial Module in the NumPy Package.

Web:
https://docs.scipy.org/doc/numpy-1.13.0/reference/routines.polynomials.polynomial.html

Other Resources:

Python Advanced Course Topics - Polynomials:

https://www.python-course.eu/polynomialc lassi np ython.php

125
 126

Figure 23.3: Trigonometric Functions

Figure 23.4: Right Triangle

Figure 23.5: Law of Cosines

127
 Part VI

Resources

128
Chapter 24

Python Resources

Here you find my Web page with Python resources [1]:

https://www.halvorsen.blog/documents/programming/python/

Python Home Page [6]:

https://www.python.org

Python Standard Library [18]:

https://docs.python.org/3/library/index.html

24.1 Python Distributions

Anaconda:
https://www.anaconda.com

24.2 Python Libraries

NumPy Library:
http://www.numpy.org

SciPy Library:
https://www.scipy.org

Matplotlib Library:
https://matplotlib.org

24.3 Python Editors

Spyder:
https://www.spyder-ide.org

129
Visual studio Code:
https://code.visualstudio.com

Visual Studio:
https://visualstudio.microsoft.com

PyCharm:
https://www.jetbrains.com/pycharm/

Wing:
https://wingware.com

Jupyter Notebook:
http://jupyter.org

24.4 Python Tutorials

Python Tutorial - w3schools.com [13]:
https://www.w3schools.com/python/

The Python Guru [19]:

https://thepythonguru.com

Wikibooks - A Beginner’s Python Tutorial:

https://en.wikibooks.org/wiki/AB eginner

TutorialsPoints - Python Tutorial:

https://www.tutorialspoint.com/python/

The Hitchhiker’s Guide to Python:

https://docs.python-guide.org

Google’s Python Class:

https://developers.google.com/edu/python/

24.5 Python in Visual Studio

Work with Python in Visual Studio
https://docs.microsoft.com/visualstudio/python/

130
Bibliography

[1] H.-P. Halvorsen, “Technology blog - https://www.halvorsen.blog,” 2018.

[2] H.-P. Halvorsen, “Technology blog - https://en.wikipedia.org/wiki/Python( programmingl anguage),00 201
[3] T. . T. P. Languages, “The 2018 top programming languages
- https://spectrum.ieee.org/at-work/innovation/the-2018-top-
programming-languages,” 2018.
[4] S. Overflow, “Stack overflow developer survey 2018 -
https://insights.stackoverflow.com/survey/2018/,” 2018.
[5] stackoverflow.blog, “The incredible growth of python -
https://stackoverflow.blog/2017/09/06/incredible-growth-python/,”
2018.
[6] python.org, “python.org - https://www.python.org,” 2018.
[7] python.org, “The python tutorial - https://docs.python.org/3.7/tutorial/,”
2018.

[8] python.org, “Python 3.7.1 documentation - https://docs.python.org/3.7/,”

2018.
[9] scipy.org, “Scipy - https://www.scipy.org,” 2018.
[10] matplotlib.org, “Matplotlib - https://matplotlib.org,” 2018.

[11] pandas, “pandas - http://pandas.pydata.org,” 2018.

[12] Wingware, “Wingware python ide - https://wingware.com,” 2018.
[13] w3schools.com, “Python tutorial - https://www.w3schools.com/python/,”
2018.

[14] Wikipedia, “Debugging - https://en.wikipedia.org/wiki/Debugging,” 2018.

[15] TechBeamers, “Get the best python ide -
https://www.techbeamers.com/best-python-ide-python-programming/,”
2018.

[16] Jupyter, “Jupyter - https://jupyter.org,” 2018.

[17] JupyterHub, “Jupyterhub - http://jupyter.org/hub,” 2018.

131
[18] python.org, “The python standard library -
https://docs.python.org/3/library/,” 2018.

[19] T. P. Guru, “The python guru - https://thepythonguru.com,” 2018.

132
 Part VII

Solutions to Exercises

133
 Start using Python

Simulation and Plotting of Dynamic System

Given the autonomous system:
ẋ = ax (1)
Where:
1
a=−
T

where T is the time constant.

The solution for the differential equation is:

x(t) = eat x0 (2)

Set T=5 and the initial condition x(0)=1.

Create a Script in Python (.py file) where you plot the solution x(t) in the time
interval:
0 ≤ t ≤ 25

Add Grid, and proper Title and Axis Labels to the plot.

Python Script:
1 im po rt math a s mt
2 im po rt numpy a s np
3 im po rt m a t p l o t l i b . p y p l o t a s p l t
4
5
6 # Model P a r a m e t e r s
7 T = 5
8 a = −1/T
9
10 # Simulation Parameters
11 x0 = 1
12 t = 0
13
14 tstart = 0

134
15 t s t o p = 25
16
17 increment = 1
18
19 x = []
20 x = np . z e r o s ( t s t o p +1)
21
22 t = np . a r a n g e ( t s t a r t , t s t o p +1 , i n c r e m e n t )
23
24
25 # Define the Function
26 f o r k in range ( tstop ) :
27 x [ k ] = mt . exp ( a ∗ t [ k ] ) ∗ x0
28
29
30 # Plot the Simulation Results
31 plt . plot (t , x)
32 p l t . t i t l e ( ’ S i m u l a t i o n o f Dynamic System ’ )
33 plt . xlabel ( ’ t ’ )
34 plt . ylabel ( ’x ’ )
35 plt . grid ()
36 p l t . a x i s ( [ 0 , 25 , 0 , 1 ] )
37 p l t . show ( )

The simulation gives the results as shown in Figure 1.

Figure 1: Simulation of Dynamic System

[End of Exercise]

135
Python Programming
c Hans-Petter Halvorsen

August 12, 2020

ISBN:978-82-691106-4-7

136
Python Programming