Book Image

NumPy Beginner's Guide - Second Edition

By : Ivan Idris
Book Image

NumPy Beginner's Guide - Second Edition

By: Ivan Idris

Overview of this book

NumPy is an extension to, and the fundamental package for scientific computing with Python. In today's world of science and technology, it is all about speed and flexibility. When it comes to scientific computing, NumPy is on the top of the list. NumPy Beginner's Guide will teach you about NumPy, a leading scientific computing library. NumPy replaces a lot of the functionality of Matlab and Mathematica, but in contrast to those products, is free and open source. Write readable, efficient, and fast code, which is as close to the language of mathematics as is currently possible with the cutting edge open source NumPy software library. Learn all the ins and outs of NumPy that requires you to know basic Python only. Save thousands of dollars on expensive software, while keeping all the flexibility and power of your favourite programming language.You will learn about installing and using NumPy and related concepts. At the end of the book we will explore some related scientific computing projects. This book will give you a solid foundation in NumPy arrays and universal functions. Through examples, you will also learn about plotting with Matplotlib and the related SciPy project. NumPy Beginner's Guide will help you be productive with NumPy and have you writing clean and fast code in no time at all.
Table of Contents (19 chapters)
Numpy Beginner's Guide Second Edition
Credits
About the Author
About the Reviewers
www.PacktPub.com
Preface
Index

Time for action – inverting matrices


The inverse of a matrix A in linear algebra is the matrix A -1, which when multiplied with the original matrix, is equal to the identity matrix I. This can be written, as A* A -1 = I.

The inv function in the numpy.linalg package can do this for us. Let's invert an example matrix. To invert matrices, perform the following steps:

  1. We will create the example matrix with the mat function that we used in the previous chapters.

    A = np.mat("0 1 2;1 0 3;4 -3 8")
    print "A\n", A

    The A matrix is printed as follows:

    A
    [[ 0  1  2]
     [ 1  0  3]
     [ 4 -3  8]]
    
  2. Now, we can see the inv function in action, using which we will invert the matrix.

    inverse = np.linalg.inv(A)
    print "inverse of A\n", inverse

    The inverse matrix is shown as follows:

    inverse of A
    [[-4.5  7.  -1.5]
     [-2.   4.  -1. ]
     [ 1.5 -2.   0.5]]
    

    Tip

    If the matrix is singular or not square, a LinAlgError exception is raised. If you want, you can check the result manually. This is left as an exercise for the reader.

  3. Let...