# HG changeset patch
# User Madhusudan.C.S <madhusudancs@gmail.com>
# Date 1257843831 -19800
# Node ID a835affb14471200a122ef063928daae5aa11bb5
# Parent  2d08c45681a1aba2407e8e318372255e2f03bec9
Added final changes to all day 2 slides for Calicut workshop.

diff -r 2d08c45681a1 -r a835affb1447 day2/exercises.tex
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/day2/exercises.tex	Tue Nov 10 14:33:51 2009 +0530
@@ -0,0 +1,201 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%Tutorial slides on Python.
+%
+% Author: Prabhu Ramachandran <prabhu at aero.iitb.ac.in>
+% Copyright (c) 2005-2009, Prabhu Ramachandran
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+\documentclass[14pt,compress]{beamer}
+%\documentclass[draft]{beamer}
+%\documentclass[compress,handout]{beamer}
+%\usepackage{pgfpages} 
+%\pgfpagesuselayout{2 on 1}[a4paper,border shrink=5mm]
+
+% Modified from: generic-ornate-15min-45min.de.tex
+\mode<presentation>
+{
+  \usetheme{Warsaw}
+  \useoutertheme{infolines}
+  \setbeamercovered{transparent}
+}
+
+\usepackage[english]{babel}
+\usepackage[latin1]{inputenc}
+%\usepackage{times}
+\usepackage[T1]{fontenc}
+
+% Taken from Fernando's slides.
+\usepackage{ae,aecompl}
+\usepackage{mathpazo,courier,euler}
+\usepackage[scaled=.95]{helvet}
+
+\definecolor{darkgreen}{rgb}{0,0.5,0}
+
+\usepackage{listings}
+\lstset{language=Python,
+    basicstyle=\ttfamily\bfseries,
+    commentstyle=\color{red}\itshape,
+  stringstyle=\color{darkgreen},
+  showstringspaces=false,
+  keywordstyle=\color{blue}\bfseries}
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% Macros
+\setbeamercolor{emphbar}{bg=blue!20, fg=black}
+\newcommand{\emphbar}[1]
+{\begin{beamercolorbox}[rounded=true]{emphbar} 
+      {#1}
+ \end{beamercolorbox}
+}
+\newcounter{time}
+\setcounter{time}{0}
+\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}}
+
+\newcommand{\typ}[1]{\texttt{#1}}
+
+\newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}}  }
+
+%%% This is from Fernando's setup.
+% \usepackage{color}
+% \definecolor{orange}{cmyk}{0,0.4,0.8,0.2}
+% % Use and configure listings package for nicely formatted code
+% \usepackage{listings}
+% \lstset{
+%    language=Python,
+%    basicstyle=\small\ttfamily,
+%    commentstyle=\ttfamily\color{blue},
+%    stringstyle=\ttfamily\color{orange},
+%    showstringspaces=false,
+%    breaklines=true,
+%    postbreak = \space\dots
+% }
+
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% Title page
+\title[Exercises]{Exercises}
+
+\author[FOSSEE] {FOSSEE}
+
+\institute[IIT Bombay] {Department of Aerospace Engineering\\IIT Bombay}
+\date[] {8 November, 2009\\Day 2}
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+%\pgfdeclareimage[height=0.75cm]{iitmlogo}{iitmlogo}
+%\logo{\pgfuseimage{iitmlogo}}
+
+
+%% Delete this, if you do not want the table of contents to pop up at
+%% the beginning of each subsection:
+\AtBeginSubsection[]
+{
+  \begin{frame}<beamer>
+    \frametitle{Outline}
+    \tableofcontents[currentsection,currentsubsection]
+  \end{frame}
+}
+
+
+% If you wish to uncover everything in a step-wise fashion, uncomment
+% the following command: 
+%\beamerdefaultoverlayspecification{<+->}
+
+%\includeonlyframes{current,current1,current2,current3,current4,current5,current6}
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% DOCUMENT STARTS
+\begin{document}
+
+\begin{frame}
+  \titlepage
+\end{frame}
+
+\begin{frame}{Problem 1.1}
+  The aliquot of a number is defined as: the sum of the \emph{proper} divisors of the number. For example:\\
+aliquot(12) = 1 + 2 + 3 + 4 + 6 = 16.\\
+  Write a function that returns the aliquot number of a given number. 
+\end{frame}
+
+\begin{frame}{Problem 1.2}
+  Pair of numbers (a, b) is said to be \alert{amicable} if aliquot number of a is b and aliquot number of b is a.\\
+  Example: \texttt{220, 284}\\
+  Write a program that prints all four digit amicable pairs.
+  
+\inctime{20}
+\end{frame}
+
+%% \begin{frame}{Problem 2}
+%%   Given an empty chessboard and one Bishop placed in any s%quare, say (r, c), generate the list of all squares the Bi%shop could move to.
+
+%% \end{frame}
+
+\begin{frame}[fragile]
+  \frametitle{Problem 2}
+  Given a string of numbers like, ``1, 3-7, 12, 15, 18-21'', produce the following list \\
+  \begin{lstlisting}
+    [1,3,4,5,6,7,12,15,18,19,20,21]
+  \end{lstlisting}
+\inctime{10}
+\end{frame}
+
+\begin{frame} 
+  \frametitle{Problem 3}
+  Count frequencies of words in a file.
+\inctime{25}
+\end{frame}
+
+\begin{frame}[fragile]
+  \frametitle{Problem set 4}
+  Central difference
+  \begin{equation*}
+  \frac{sin(x+h)-sin(x-h)}{2h}
+  \end{equation*}
+  \begin{lstlisting}
+  In []: x = linspace(0, 2*pi, 100)
+  In []: y = sin(x)
+  In []: deltax = x[1] - x[0]
+  \end{lstlisting}
+  \pause
+    \begin{enumerate}
+      \item Given this, get the finite difference of sin in the range 0 to 2*pi
+    \end{enumerate}
+\end{frame}
+
+\begin{frame}
+  \frametitle{Problem Set 5}
+  \begin{itemize}
+      \item[5.1] Write a function that plots any regular n-gon given \typ{n}.
+      \item[5.2] Consider the logistic map, $f(x) = kx(1-x)$, plot it for
+          $k=2.5, 3.5$ and $4$ in the same plot.
+\end{itemize}
+\end{frame}
+
+\begin{frame}[fragile] 
+\frametitle{Problem Set 5}
+  \begin{columns}
+    \column{0.6\textwidth}
+    \small{
+    \begin{itemize}
+      \item[3] Consider the iteration $x_{n+1} = f(x_n)$ where $f(x) = kx(1-x)$.  Plot the successive iterates of this process as explained below. 
+    \end{itemize}}
+    \column{0.35\textwidth}
+    \hspace*{-0.5in}
+  \includegraphics[height=1.6in, interpolate=true]{data/cobweb}  
+\end{columns}
+\end{frame}
+
+\begin{frame}
+  \frametitle{Problem Set 5.3}
+  Plot the cobweb plot as follows:
+  \begin{enumerate}
+    \item Start at $(x_0, 0)$ ($\implies$ i=0)
+    \item Draw a line to $(x_i, f(x_i))$
+    \item Set $x_{i+1} = f(x_i)$
+    \item Draw a line to $(x_{i+1}, x_{i+1})$
+    \item $(i\implies i+1)$ 
+    \item Repeat from 2 for as long as you want 
+  \end{enumerate}
+\inctime{20}
+\end{frame}
+
+\end{document}
diff -r 2d08c45681a1 -r a835affb1447 day2/session1.tex
--- a/day2/session1.tex	Tue Nov 10 14:33:06 2009 +0530
+++ b/day2/session1.tex	Tue Nov 10 14:33:51 2009 +0530
@@ -51,7 +51,8 @@
 \setcounter{time}{0}
 \newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}}
 
-\newcommand{\typ}[1]{\texttt{#1}}
+\newcommand{\typ}[1]{\lstinline{#1}}
+
 
 \newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}}  }
 
@@ -130,7 +131,7 @@
   \begin{itemize}
     \item Numbers: float, int, complex
     \item Strings
-    \item Boolean
+    \item Booleans
   \end{itemize}
 \end{frame}
 
@@ -138,7 +139,7 @@
 \begin{frame}[fragile]
   \frametitle{Numbers}
   \begin{itemize}
-    \item \kwrd{int}\\ \kwrd{int} = whole number, no matter what the size!
+    \item \kwrd{int}\\ whole number, no matter what the size!
   \begin{lstlisting}
 In []: a = 13
 
@@ -167,9 +168,9 @@
   \end{lstlisting}
 \end{frame}
 
-\subsection{Boolean}
+\subsection{Booleans}
 \begin{frame}[fragile]
-  \frametitle{Boolean}
+  \frametitle{Booleans}
   \begin{lstlisting}
 In []: t = True
 
@@ -211,15 +212,19 @@
 In []: print w[0] + w[2] + w[-1]
 Out[]: hlo
 
-In []: len(w) # guess what
+In []: len(w)
 Out[]: 5
   \end{lstlisting}
 \end{frame}
 
 \begin{frame}[fragile]
   \frametitle{Strings \ldots}
+  \emphbar{Strings are immutable}
   \begin{lstlisting}
-In []: w[0] = 'H' # Can't do that!
+In []: w[0] = 'H' 
+  \end{lstlisting}
+  \pause
+  \begin{lstlisting}
 --------------------------------------------
 TypeError  Traceback (most recent call last)
 
@@ -248,7 +253,7 @@
   \end{lstlisting}
 \end{frame}
 
-\begin{frame}
+\begin{frame}[fragile]
   \frametitle{A bit about IPython}
   \begin{itemize}
     \item IPython provides better help
@@ -265,7 +270,7 @@
 \end{frame}
 
 \begin{frame}[fragile]
-\frametitle{Still with strings}
+  \frametitle{Still with strings}
   \begin{itemize}
     \item We saw split() yesterday
     \item join() is the opposite of split()
@@ -294,6 +299,7 @@
 \section{Operators}
 \begin{frame}[fragile]
   \frametitle{Arithmetic operators}
+  \small
   \begin{lstlisting}
 In []: 1786 % 12
 Out[]: 10
@@ -351,15 +357,15 @@
 \begin{frame}[fragile]
   \frametitle{String operations}
   \begin{lstlisting}
-In []: s = 'Hello '
+In []: s = 'Hello'
 
 In []: p = 'World'
 
 In []: s + p 
-Out[]: 'Hello World'
+Out[]: 'HelloWorld'
 
 In []: s * 4
-Out[]: 'Hello Hello Hello Hello'
+Out[]: 'HelloHelloHelloHello'
   \end{lstlisting}
 \end{frame}
 
@@ -367,13 +373,16 @@
   \frametitle{String operations \ldots}
   \begin{lstlisting}
 In []: s * s
+  \end{lstlisting}
+  \pause
+  \begin{lstlisting}
 --------------------------------------------
 TypeError  Traceback (most recent call last)
 
 /<ipython console> in <module>()
 
-TypeError: can't multiply sequence by
-                non-int of type 'str'
+TypeError: can`t multiply sequence by
+                non-int of type `str`
   \end{lstlisting}
 \end{frame}
 
@@ -401,7 +410,7 @@
 In []: int(17 / 2.0)
 Out[]: 8
 
-In []: float(17 / 2)  # Recall
+In []: float(17 / 2)
 Out[]: 8.0
 
 In []: str(17 / 2.0)
@@ -455,7 +464,7 @@
   \frametitle{Simple IO: Console output}
   \begin{itemize}
     \item \texttt{print} is straight forward
-    \item Put the following code snippet in a file
+    \item Put the following code snippet in a file \emph{hello1.py}
   \end{itemize}
   \begin{lstlisting}
 print "Hello"
@@ -470,7 +479,7 @@
 
 \begin{frame}[fragile]
   \frametitle{Simple IO: Console output \ldots}
-Put the following code snippet in a file
+Put the following code snippet in a file \emph{hello2.py}
   \begin{lstlisting}
 print "Hello",
 print "World"
@@ -522,7 +531,7 @@
   \begin{lstlisting}
 ...
 a = raw_input('Enter number(Q to quit):')
-num = int(a) if a != 'Q' else break
+num = int(a) if a != 'Q' else 0
 ...
   \end{lstlisting}
 \end{frame}
diff -r 2d08c45681a1 -r a835affb1447 day2/session2.tex
--- a/day2/session2.tex	Tue Nov 10 14:33:06 2009 +0530
+++ b/day2/session2.tex	Tue Nov 10 14:33:51 2009 +0530
@@ -51,7 +51,7 @@
 \setcounter{time}{0}
 \newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}}
 
-\newcommand{\typ}[1]{\texttt{#1}}
+\newcommand{\typ}[1]{\lstinline{#1}}
 
 \newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}}  }
 
@@ -127,10 +127,10 @@
 \subsection{Basic Looping}
 \begin{frame}[fragile]
   \frametitle{\typ{while}}
-Example: Fibonacci series
+\begin{block}{Example: Fibonacci series}
+  Sum of previous two elements defines the next
+\end{block}
   \begin{lstlisting}
-# the sum of two elements
-# defines the next
 In []: a, b = 0, 1
 In []: while b < 10:
   ...:     print b,
@@ -145,23 +145,24 @@
 \frametitle{\typ{range()}}
 \kwrd{range([start,] stop[, step])}\\
 \begin{itemize}
-  \item range() returns a list of integers
-  \item The \emph{start} and the \emph{step} arguments are optional
-  \item \emph{stop} argument is not included in the list
+  \item \typ{range()} returns a list of integers
+  \item The \typ{start} and the \typ{step} arguments are optional
+  \item \typ{stop} is not included in the list
 \end{itemize}
 \vspace*{.5in}
-\begin{itemize}
-  \item \alert{Anything within \typ{[]} is optional}
+\begin{block}{Documentation convention}
   \begin{itemize}
-    \item Nothing to do with Python.
+    \item \alert{Anything within \typ{[]} is optional}
+    \begin{itemize}
+      \item Nothing to do with Python.
+    \end{itemize}
   \end{itemize}
-\end{itemize}
-
+\end{block}
 \end{frame}
 
 \begin{frame}[fragile]
-  \frametitle{\typ{for} \ldots \typ{range()}}
-Example: print squares of first \typ{n} numbers
+  \frametitle{\texttt{for} \ldots \typ{range()}}
+Example: print squares of first \typ{5} numbers
   \begin{lstlisting}
 In []: for i in range(5):
  ....:     print i, i * i
@@ -176,7 +177,7 @@
 \end{frame}
 
 \begin{frame}[fragile]
-  \frametitle{\typ{for} \ldots \typ{range()}}
+  \frametitle{\texttt{for} \ldots \typ{range()}}
 Example: print squares of odd numbers from 3 to 9
   \begin{lstlisting}
 In []: for i in range(3, 10, 2):
@@ -240,6 +241,9 @@
   \begin{lstlisting}
 In []: num = [1, 2, 3, 4]
 
+In []: num + [9, 10, 11]
+Out[]: [1, 2, 3, 4, 9, 10, 11]
+
 In []: num.append([9, 10, 11])
 
 In []: num
@@ -325,7 +329,6 @@
 In []: t[0] + t[3] + t[-1]
 Out[]: 13
 
-# Try the following!
 In []: t[4] = 7 
 \end{lstlisting}
 \pause
@@ -361,8 +364,8 @@
 Out[]: 52.530000000000001
   \end{lstlisting}
   \begin{block}{Note!}
-    Duplicate keys are not allowed!\\
-    Dictionaries are iterable through keys.
+    Duplicate keys $\Rightarrow$ overwritten!\\
+    You can iterate through a dictionary using keys.
   \end{block}
 \end{frame}
 
@@ -375,6 +378,12 @@
 In []: 'Econ' in player
 Out[]: False
   \end{lstlisting}
+  \begin{block}{Note}
+    \begin{itemize}
+      \item We can check for the containership of keys only
+      \item Not values
+    \end{itemize}
+  \end{block}
 \end{frame}
 
 \begin{frame}[fragile]
@@ -384,17 +393,18 @@
 Out[]: ['Runs', 'Inn', 'Avg', 'Mat']
 
 In []: player.values()
-Out[]: [10823, 233, 52.530000000000001, 134]
+Out[]: [10823, 233, 
+        52.530000000000001, 134]
   \end{lstlisting}
 \end{frame}
 
 \begin{frame} {Problem Set 2.1: Problem 2.1.1}
-You are given date strings of the form ``29, Jul 2009'', or ``4 January 2008''. In other words a number, a string and another number, with a comma sometimes separating the items.\\Write a function that takes such a string and returns a tuple (yyyy, mm, dd) where all three elements are ints.
+You are given date strings of the form ``29 Jul, 2009'', or ``4 January 2008''. In other words a number, a string and another number, with a comma sometimes separating the items.\\Write a function that takes such a string and returns a tuple (yyyy, mm, dd) where all three elements are ints.
 \end{frame}
 
-\subsection{Set}
+\subsection{Sets}
 \begin{frame}[fragile]
-  \frametitle{Set}
+  \frametitle{Sets}
     \begin{itemize}
       \item Simplest container, mutable
       \item No ordering, no duplicates
@@ -447,18 +457,14 @@
 \end{frame}
 
 \begin{frame}
-  \frametitle{Problem set 2.2}
-  \begin{description}
-    \item[2.2.1] Given a dictionary of the names of students and their marks, identify how many duplicate marks are there? and what are these?
-\end{description}
-\inctime{15}
+  \frametitle{Problem set 2.2: Problem 2.2.1}
+Given a dictionary of the names of students and their marks, identify how many duplicate marks are there? and what are these?
 \end{frame}
 
 \begin{frame}
-  \frametitle{Problem set 2.2}
-  \begin{description}
-    \item[2.2.2] Given a list of words, find all the anagrams in the list
-\end{description}
+  \frametitle{Problem 2.2.2}
+Given a list of words, find all the anagrams in the list.
+
 \inctime{15}
 \end{frame}
 
@@ -468,7 +474,6 @@
   \begin{itemize}
     \item \kwrd{def} - keyword to define a function
     \item Arguments are local to a function
-    \item Docstrings are important!
     \item Functions can return multiple values
   \end{itemize}
 \end{frame}
@@ -487,6 +492,7 @@
     else:
         return 0
   \end{lstlisting}
+  \emphbar{Note docstrings}
 \end{frame}
 
 \begin{frame}[fragile]
@@ -517,8 +523,15 @@
   \frametitle{What did we learn?}
   \begin{itemize}
     \item Loops: \kwrd{while}, \kwrd{for}
-    \item Advanced Data structures
+    \item Advanced Data structures:
+    \begin{itemize}
+      \item Lists
+      \item Tuples
+      \item Dictionaries
+      \item Sets
+    \end{itemize}
     \item Functions
+    \item Docstrings
   \end{itemize}
 \end{frame}
 
diff -r 2d08c45681a1 -r a835affb1447 day2/session3.tex
--- a/day2/session3.tex	Tue Nov 10 14:33:06 2009 +0530
+++ b/day2/session3.tex	Tue Nov 10 14:33:51 2009 +0530
@@ -51,7 +51,7 @@
 \setcounter{time}{0}
 \newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}}
 
-\newcommand{\typ}[1]{\texttt{#1}}
+\newcommand{\typ}[1]{\lstinline{#1}}
 
 \newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}}  }
 
@@ -73,7 +73,7 @@
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % Title page
-\title[Basic Python]{Python language: Data structures and functions}
+\title[Basic Python]{Python language: Functions, modules and objects}
 
 \author[FOSSEE Team] {The FOSSEE Group}
 
@@ -143,20 +143,15 @@
 
 \begin{frame}[fragile]
   \frametitle{Functions: default arguments \ldots}
-  \small
   \begin{lstlisting}
-def ask_ok(prompt, complaint='Yes or no!'):
-    while True:
-        ok = raw_input(prompt)
-        if ok in ('y', 'ye', 'yes'):
-            return True
-        if ok in ('n', 'no', 'nop',
-                            'nope'):
-            return False
-        print complaint
+In []: def welcome(greet, name="World"):
+  ....     print greet, name
 
-ask_ok('?')
-ask_ok('?', '[Y/N]')
+In []: welcome("Hello")
+Hello World
+
+In []: welcome("Hi", "Guido")
+Hi Guido
   \end{lstlisting}
 \end{frame} 
 
@@ -164,43 +159,43 @@
 \begin{frame}[fragile]
   \frametitle{Functions: Keyword arguments}
 We have seen the following
-  \begin{lstlisting}
-In []: legend(['sin(2y)'],
-              loc='center')
+\begin{lstlisting}
+In []: legend(['sin(2y)'], 
+                 loc = 'center')
+
 In []: plot(y, sin(y), 'g',
-                 linewidth=2)
+                 linewidth = 2)
+
 In []: annotate('local max',
-                 xy=(1.5, 1))
+                 xy = (1.5, 1))
+
 In []: pie(science.values(),
-            labels=science.keys())
+            labels = science.keys())
   \end{lstlisting}
 \end{frame}
 
 \begin{frame}[fragile]
   \frametitle{Functions: keyword arguments \ldots}
-  \small
   \begin{lstlisting}
-def ask_ok(prompt, complaint='Yes or no!'):
-    while True:
-        ok = raw_input(prompt)
-        if ok in ('y', 'ye', 'yes'):
-            return True
-        if ok in ('n', 'no', 'nop',
-                            'nope'):
-            return False
-        print complaint
+In []: def welcome(greet, name="World"):
+  ....     print greet, name
+
+In []: welcome("Hello", "James")
+Hello James
 
-ask_ok(prompt='?')
-ask_ok(prompt='?', complaint='[y/n]')
-ask_ok(complaint='[y/n]', prompt='?')
-\end{lstlisting}
+In []: welcome("Hi", name="Guido")
+Hi Guido
+
+In []: welcome(name="Guido", greet="Hey")
+Hey Guido
+  \end{lstlisting}
 \end{frame}
 
 \subsection{Built-in functions}
 \begin{frame}
   {Before writing a function}
   \begin{itemize}
-      \item Variety of builtin functions are available
+      \item Variety of built-in functions are available
       \item \typ{abs, any, all, len, max, min}
       \item \typ{pow, range, sum, type}
       \item Refer here:
@@ -221,36 +216,32 @@
 
 \begin{frame}{Problem 3.3}
   Write a program that generates a list of all four digit numbers that have all their digits even and are perfect squares.\newline\\\emph{For example, the output should include 6400 but not 8100 (one digit is odd) or 4248 (not a perfect square).}
+
 \inctime{15}
 \end{frame}
 
 \section{Modules}
 \begin{frame}[fragile]
-  \frametitle{\typ{from} \ldots \typ{import} magic}
+  \frametitle{\texttt{from} \ldots \texttt{import} magic}
   \begin{lstlisting}
-from scipy.interpolate import splrep
-from scipy.interpolate import splev
-
-from scipy.integrate import quad
 from scipy.integrate import odeint
 
 from scipy.optimize import fsolve
   \end{lstlisting}
-\emphbar{All the above statements import one function into your namespace}
+\emphbar{Above statements import a function to our namespace}
 \end{frame}
 
 \begin{frame}[fragile]
   \frametitle{Running scripts from command line}
   \small
   \begin{itemize}
-    \item Start cmd
-    \item cd to Desktop
-    \item python sine\_plot.py
+    \item Fire up a terminal
+    \item python four\_plot.py
   \end{itemize}
   \pause
   \begin{lstlisting}
 Traceback (most recent call last):
-  File "sine_plot.py", line 1, in <module>
+  File "four_plot.py", line 1, in <module>
     x = linspace(-5*pi, 5*pi, 500)
 NameError: name 'linspace' is not defined
   \end{lstlisting}
@@ -258,15 +249,14 @@
 
 \begin{frame}[fragile]
   \frametitle{Remedy}
-  \emphbar{Adding what lines to sine\_plot.py makes this program work?}
   \begin{lstlisting}
 from scipy import *
   \end{lstlisting}
-\alert{Now run python sine\_plot.py again!}
+\alert{Now run python four\_plot.py again!}
   \pause
   \begin{lstlisting}
 Traceback (most recent call last):
-  File "sine_plot.py", line 4, in <module>
+  File "four_plot.py", line 4, in <module>
     plot(x, x, 'b')
 NameError: name 'plot' is not defined
   \end{lstlisting}
@@ -274,11 +264,10 @@
 
 \begin{frame}[fragile]
   \frametitle{Remedy \ldots}
-  \emphbar{What should we add now?}
   \begin{lstlisting}
 from pylab import *
   \end{lstlisting}
-\alert{Now run python sine\_plot.py again!!}
+\alert{Now run python four\_plot.py again!!}
 \end{frame}
 
 \begin{frame}[fragile]
@@ -298,16 +287,14 @@
 \begin{frame}[fragile]
   \frametitle{Package hierarchies}
   \begin{lstlisting}
-from scipy.interpolate import splev
-
-from scipy.integrate import quad
+from scipy.integrate import odeint
 
 from scipy.optimize import fsolve
   \end{lstlisting}
 \end{frame}
 
 \begin{frame}[fragile]
-  \frametitle{from \ldots import in a conventional way!}
+  \frametitle{\texttt{from} \ldots \texttt{import} - conventional way!}
   \small
   \begin{lstlisting}
 from scipy import linspace, pi, sin
@@ -327,7 +314,7 @@
 \end{frame}
 
 \begin{frame}[fragile]
-  \frametitle{from \ldots import in a conventional way!}
+  \frametitle{\texttt{from} \ldots \texttt{import} - conventional way!}
   \small
   \begin{lstlisting}
 import scipy
@@ -368,13 +355,13 @@
 \begin{frame}[fragile]
   \frametitle{Modules of special interest}
   \begin{description}[matplotlibfor2d]
-    \item[\typ{pylab}] Easy, interactive, 2D plotting
+    \item[\texttt{pylab}] Easy, interactive, 2D plotting
 
-    \item[\typ{scipy}] arrays, statistics, optimization, integration, linear
+    \item[\texttt{scipy}] arrays, statistics, optimization, integration, linear
             algebra, Fourier transforms, signal and image processing,
             genetic algorithms, ODE solvers, special functions, and more
 
-    \item[\typ{Mayavi}] Easy, interactive, 3D plotting
+    \item[\texttt{Mayavi}] Easy, interactive, 3D plotting
   \end{description}
 \end{frame}
 
diff -r 2d08c45681a1 -r a835affb1447 day2/session4.tex
--- a/day2/session4.tex	Tue Nov 10 14:33:06 2009 +0530
+++ b/day2/session4.tex	Tue Nov 10 14:33:51 2009 +0530
@@ -144,7 +144,7 @@
 \begin{frame}[fragile] 
   \frametitle{gcd revisited!}
   \begin{itemize}
-  \item Open gcd.py
+  \item Open \texttt{gcd.py}
   \end{itemize}  
 \begin{lstlisting}
     def gcd(a, b):
@@ -156,15 +156,15 @@
     print gcd(16, 76)
 \end{lstlisting}
   \begin{itemize}
-  \item python gcd.py
+  \item \texttt{python gcd.py}
   \end{itemize}
 \end{frame}
 
 \begin{frame}[fragile] 
   \frametitle{Find lcm using our gcd module}
   \begin{itemize}
-  \item Open lcm.py  
-  \item $lcm = \frac{a*b}{gcd(a,b)}$
+  \item Open \texttt{lcm.py}  
+  \item $lcm = \frac{a * b}{gcd(a,b)}$
   \end{itemize}  
 \begin{lstlisting}
     from gcd import gcd    
@@ -174,7 +174,7 @@
     print lcm(14, 56)
 \end{lstlisting}
   \begin{itemize}
-  \item python lcm.py
+  \item \texttt{python lcm.py}
   \end{itemize}
   \begin{lstlisting}
 5
@@ -185,7 +185,7 @@
 
 \begin{frame}[fragile] 
   \frametitle{Writing stand-alone module}  
-Edit gcd.py file to:
+Edit \texttt{gcd.py} file to:
 \begin{lstlisting}
     def gcd(a, b):
         if a % b == 0: 
@@ -197,22 +197,23 @@
         print gcd(16, 76)
 \end{lstlisting}
   \begin{itemize}
-  \item python gcd.py
-  \item python lcm.py
+  \item \texttt{python gcd.py}
+  \item \texttt{python lcm.py}
   \end{itemize}
 \end{frame}
 
 \begin{frame}[fragile]
-  \frametitle{More use of main}
-  For automating tests.
+  \frametitle{Automating tests}
   \begin{lstlisting}
 if __name__ == '__main__':
     for line in open('numbers.txt'):
         numbers = line.split()
         x = int(numbers[0])
         y = int(numbers[1])
-        result = (int(numbers[2]))
-        assert gcd(x, y) == result
+        result = int(numbers[2])
+        if gcd(x, y) != result:
+            print "Failed gcd test
+                          for", x, y
   \end{lstlisting}  
 \end{frame}
 
@@ -242,7 +243,7 @@
   \frametitle{Code Layout}
   \begin{itemize}
         \item Indentation
-        \item Tabs or Spaces??
+        \item Tabs or Spaces?
         \item Maximum Line Length
         \item Blank Lines
         \item Encodings
@@ -251,8 +252,8 @@
 
 \begin{frame}{Whitespaces in Expressions}
   \begin{itemize}
-        \item When to use extraneous whitespaces??
-        \item When to avoid extra whitespaces??
+        \item When to use extraneous whitespaces?
+        \item When to avoid extra whitespaces?
         \item Use one statement per line
    \end{itemize}
 \end{frame}
@@ -319,43 +320,36 @@
 \end{frame}
 
 \begin{frame}[fragile]
+  \frametitle{Processing user input}
+  \begin{lstlisting}
+prompt = 'Enter a number(Q to quit): '
+
+a = raw_input(prompt)
+
+num = int(a) if a != 'Q' else 0
+  \end{lstlisting}
+  \emphbar{What if the user enters some other alphabet?}
+\end{frame}
+
+
+\begin{frame}[fragile]
   \frametitle{Handling Exceptions}
-  Python uses \typ{try} and \typ{except} clause.
-  %%Revisiting the raw\_input
+  Python provides \typ{try} and \typ{except} clause.
   \begin{lstlisting}
-a = raw_input('Enter number(Q to quit):')
+prompt = 'Enter a number(Q to quit): '
+
+a = raw_input(prompt)
 try:
     num = int(a)
     print num
 except:
     if a == 'Q':
-        print 'Exiting...'
+        print "Exiting ..."
     else:
-        print 'Wrong input!'      
-  \end{lstlisting}
-  
-  
+        print "Wrong input ..."
+  \end{lstlisting}  
 \end{frame}
 
-%% \begin{frame}[fragile]
-%%   \frametitle{Solving it with \typ{try} and \typ{except}}
-%% \vspace{-0.2in}
-%%   \begin{lstlisting}
-%% highest = 0
-%% for record in open('sslc1.txt'):
-%%     fields = record.split(';')
-%%     try:
-%%         total = 0
-%%         for score_str in fields[3:8]:
-%%             score = int(score_str)
-%%             total += score
-%%         if total > highest:
-%%             highest = total
-%%     except:        
-%%         pass
-%% print highest
-%%   \end{lstlisting}
-%% \end{frame}
 \subsection{Strategy}
 \begin{frame}[fragile]
     \frametitle{Debugging effectively}
@@ -421,6 +415,20 @@
 \inctime{10}
 \end{frame}
 
+\begin{frame}
+  \frametitle{Summary}
+We have covered:
+  \begin{itemize}
+  \item Following and Resolving Error Messages.
+  \item Exceptions.
+  \item Handling exceptions
+  \item Approach for Debugging.
+%  \item Writting and running tests.
+  \end{itemize}
+\end{frame}
+
+\end{document}
+
 %% \begin{frame}
 %%     \frametitle{Testing}
    
@@ -433,130 +441,116 @@
 %%     \end{itemize}
 %% \end{frame}
 
-\section{Test Driven Approach}
-\begin{frame}
-    \frametitle{Need for Testing!}
-   
-    \begin{itemize}
-        \item Quality
-        \item Regression
-        \item Documentation
-    \end{itemize}
-    %% \vspace*{0.25in}
-    %% \emphbar{It is to assure that section of code is working as it is supposed to work}
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{Example}
-    \begin{block}{Problem Statement}
-      Write a function to check whether a given input
-      string is a palindrome.
-    \end{block}
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{Function: palindrome.py}
-\begin{lstlisting}    
-def is_palindrome(input_str):
-  return input_str == input_str[::-1]
-\end{lstlisting}    
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{Test for the palindrome: palindrome.py}
-\begin{lstlisting}    
-def test_function_normal_words():
-  input = "noon"
-  assert is_palindrome(input) == True
-
-if __name__ == "main'':
-  test_function_normal_words()
-\end{lstlisting}    
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{Running the tests.}
-\begin{lstlisting}    
-$ nosetests palindrome.py 
-.
-----------------------------------------------
-Ran 1 test in 0.001s
-
-OK
-\end{lstlisting}    
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{Exercise: Including new tests.}
-\begin{lstlisting}    
-def test_function_ignore_cases_words():
-  input = "Noon"
-  assert is_palindrome(input) == True
-\end{lstlisting}
-     \vspace*{0.25in}
-     Check\\
-     \PythonCode{$ nosetests palindrome.py} \\
-     \begin{block}{Task}
-     Tweak the code to pass this test.
-     \end{block}
-\end{frame}
-
-%\begin{frame}[fragile]
-%    \frametitle{Lets write some test!}
-%\begin{lstlisting}    
-%#for form of equation y=mx+c
-%#given m and c for two equation,
-%#finding the intersection point.
-%def intersect(m1,c1,m2,c2):
-%    x = (c2-c1)/(m1-m2)
-%    y = m1*x+c1
-%    return (x,y)
-%\end{lstlisting}
-%
-%Create a simple test for this
-%
-%function which will make it fail.
-%
-%\inctime{15} 
-%\end{frame}
-%
-
-%% \begin{frame}[fragile]
-%%     \frametitle{Exercise}
-%%     Based on Euclid's algorithm:
-%%     \begin{center}
-%%     $gcd(a,b)=gcd(b,b\%a)$
-%%     \end{center}
-%%     gcd function can be written as:
-%%     \begin{lstlisting}
-%%     def gcd(a, b):
-%%       if a%b == 0: return b
-%%       return gcd(b, a%b)
-%%     \end{lstlisting}
-%%     \vspace*{-0.15in}
-%%     \begin{block}{Task}
-%%       \begin{itemize}
-%%       \item Write at least 
-%%         two tests for above mentioned function.
-%%       \item Write a non recursive implementation
-%%       of gcd(), and test it using already 
-%%       written tests.
-%%       \end{itemize}
-%%     \end{block}
-    
-%% \inctime{15} 
-%% \end{frame}
-
-\begin{frame}
-  \frametitle{Summary}
-We have coverd:
-  \begin{itemize}
-  \item Following and Resolving Error Messages.
-  \item Exceptions.
-  \item Handling exceptions
-  \item Approach for Debugging.
-  \item Writting and running tests.
-  \end{itemize}
-\end{frame}
-
-\end{document}
+% \section{Test Driven Approach}
+% \begin{frame}
+%     \frametitle{Need for Testing!}
+%    
+%     \begin{itemize}
+%         \item Quality
+%         \item Regression
+%         \item Documentation
+%     \end{itemize}
+%     %% \vspace*{0.25in}
+%     %% \emphbar{It is to assure that section of code is working as it is supposed to work}
+% \end{frame}
+% 
+% \begin{frame}[fragile]
+%     \frametitle{Example}
+%     \begin{block}{Problem Statement}
+%       Write a function to check whether a given input
+%       string is a palindrome.
+%     \end{block}
+% \end{frame}
+% 
+% \begin{frame}[fragile]
+%     \frametitle{Function: palindrome.py}
+% \begin{lstlisting}    
+% def is_palindrome(input_str):
+%   return input_str == input_str[::-1]
+% \end{lstlisting}    
+% \end{frame}
+% 
+% \begin{frame}[fragile]
+%     \frametitle{Test for the palindrome: palindrome.py}
+% \begin{lstlisting}    
+% def test_function_normal_words():
+%   input = "noon"
+%   assert is_palindrome(input) == True
+% 
+% if __name__ == "main'':
+%   test_function_normal_words()
+% \end{lstlisting}    
+% \end{frame}
+% 
+% \begin{frame}[fragile]
+%     \frametitle{Running the tests.}
+% \begin{lstlisting}    
+% $ nosetests palindrome.py 
+% .
+% ----------------------------------------------
+% Ran 1 test in 0.001s
+% 
+% OK
+% \end{lstlisting}    
+% \end{frame}
+% 
+% \begin{frame}[fragile]
+%     \frametitle{Exercise: Including new tests.}
+% \begin{lstlisting}    
+% def test_function_ignore_cases_words():
+%   input = "Noon"
+%   assert is_palindrome(input) == True
+% \end{lstlisting}
+%      \vspace*{0.25in}
+%      Check\\
+%      \PythonCode{$ nosetests palindrome.py} \\
+%      \begin{block}{Task}
+%      Tweak the code to pass this test.
+%      \end{block}
+% \end{frame}
+% 
+% %\begin{frame}[fragile]
+% %    \frametitle{Lets write some test!}
+% %\begin{lstlisting}    
+% %#for form of equation y=mx+c
+% %#given m and c for two equation,
+% %#finding the intersection point.
+% %def intersect(m1,c1,m2,c2):
+% %    x = (c2-c1)/(m1-m2)
+% %    y = m1*x+c1
+% %    return (x,y)
+% %\end{lstlisting}
+% %
+% %Create a simple test for this
+% %
+% %function which will make it fail.
+% %
+% %\inctime{15} 
+% %\end{frame}
+% %
+% 
+% %% \begin{frame}[fragile]
+% %%     \frametitle{Exercise}
+% %%     Based on Euclid's algorithm:
+% %%     \begin{center}
+% %%     $gcd(a,b)=gcd(b,b\%a)$
+% %%     \end{center}
+% %%     gcd function can be written as:
+% %%     \begin{lstlisting}
+% %%     def gcd(a, b):
+% %%       if a%b == 0: return b
+% %%       return gcd(b, a%b)
+% %%     \end{lstlisting}
+% %%     \vspace*{-0.15in}
+% %%     \begin{block}{Task}
+% %%       \begin{itemize}
+% %%       \item Write at least 
+% %%         two tests for above mentioned function.
+% %%       \item Write a non recursive implementation
+% %%       of gcd(), and test it using already 
+% %%       written tests.
+% %%       \end{itemize}
+% %%     \end{block}
+%     
+% %% \inctime{15} 
+% %% \end{frame}
\ No newline at end of file
diff -r 2d08c45681a1 -r a835affb1447 day2/session5.tex
--- a/day2/session5.tex	Tue Nov 10 14:33:06 2009 +0530
+++ b/day2/session5.tex	Tue Nov 10 14:33:51 2009 +0530
@@ -1,33 +1,48 @@
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%Tutorial slides on Python.
+% Tutorial slides on Python.
 %
 % Author: Prabhu Ramachandran <prabhu at aero.iitb.ac.in>
 % Copyright (c) 2005-2009, Prabhu Ramachandran
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-\documentclass[14pt,compress]{beamer}
-%\documentclass[draft]{beamer}
-%\documentclass[compress,handout]{beamer}
-%\usepackage{pgfpages} 
-%\pgfpagesuselayout{2 on 1}[a4paper,border shrink=5mm]
+\documentclass[compress,14pt]{beamer}
+% \documentclass[handout]{beamer}
+% \usepackage{pgfpages}
+% \pgfpagesuselayout{4 on 1}[a4paper,border, shrink=5mm,landscape]
+\usepackage{tikz}
+\newcommand{\hyperlinkmovie}{}
+%\usepackage{movie15}
 
-% Modified from: generic-ornate-15min-45min.de.tex
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% Note that in presentation mode 
+% \paperwidth  364.19536pt
+% \paperheight 273.14662pt
+% h/w = 0.888
+
+
 \mode<presentation>
 {
   \usetheme{Warsaw}
+  %\usetheme{Boadilla}
+  %\usetheme{default}
   \useoutertheme{infolines}
   \setbeamercovered{transparent}
 }
 
+% To remove navigation symbols
+\setbeamertemplate{navigation symbols}{}
+
+\usepackage{amsmath}
 \usepackage[english]{babel}
 \usepackage[latin1]{inputenc}
-%\usepackage{times}
+\usepackage{times}
 \usepackage[T1]{fontenc}
 
 % Taken from Fernando's slides.
 \usepackage{ae,aecompl}
 \usepackage{mathpazo,courier,euler}
 \usepackage[scaled=.95]{helvet}
+\usepackage{pgf}
 
 \definecolor{darkgreen}{rgb}{0,0.5,0}
 
@@ -40,50 +55,65 @@
   keywordstyle=\color{blue}\bfseries}
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% Macros
+% My Macros
+\setbeamercolor{postit}{bg=yellow,fg=black}
 \setbeamercolor{emphbar}{bg=blue!20, fg=black}
 \newcommand{\emphbar}[1]
 {\begin{beamercolorbox}[rounded=true]{emphbar} 
       {#1}
  \end{beamercolorbox}
 }
+%{\centerline{\fcolorbox{gray!50} {blue!10}{
+%\begin{minipage}{0.9\linewidth}
+%    {#1} 
+%\end{minipage}
+%    }}}
+
+\newcommand{\myemph}[1]{\structure{\emph{#1}}}
+\newcommand{\PythonCode}[1]{\lstinline{#1}}
+
+\newcommand{\tvtk}{\texttt{tvtk}}
+\newcommand{\mlab}{\texttt{mlab}}
+
 \newcounter{time}
 \setcounter{time}{0}
-\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\tiny \thetime\ m}}
-
-\newcommand{\typ}[1]{\texttt{#1}}
-
-\newcommand{\kwrd}[1]{ \texttt{\textbf{\color{blue}{#1}}}  }
+\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\vspace*{0.1in}\tiny \thetime\ m}}
 
-%%% This is from Fernando's setup.
-% \usepackage{color}
-% \definecolor{orange}{cmyk}{0,0.4,0.8,0.2}
-% % Use and configure listings package for nicely formatted code
-% \usepackage{listings}
-% \lstset{
-%    language=Python,
-%    basicstyle=\small\ttfamily,
-%    commentstyle=\ttfamily\color{blue},
-%    stringstyle=\ttfamily\color{orange},
-%    showstringspaces=false,
-%    breaklines=true,
-%    postbreak = \space\dots
-% }
+\newcommand\BackgroundPicture[1]{%
+  \setbeamertemplate{background}{%
+      \parbox[c][\paperheight]{\paperwidth}{%
+      \vfill \hfill
+ \hfill \vfill
+}}}
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% Configuring the theme
+%\setbeamercolor{normal text}{fg=white}
+%\setbeamercolor{background canvas}{bg=black}
+
 
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % Title page
-\title[Exercises]{Exercises}
+\title[3D Plotting]{3D data Visualization}
 
 \author[FOSSEE] {FOSSEE}
 
 \institute[IIT Bombay] {Department of Aerospace Engineering\\IIT Bombay}
 \date[] {8 November, 2009\\Day 2, Session 5}
+
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-%\pgfdeclareimage[height=0.75cm]{iitmlogo}{iitmlogo}
-%\logo{\pgfuseimage{iitmlogo}}
+%\pgfdeclareimage[height=0.75cm]{iitblogo}{iitblogo}
+%\logo{\pgfuseimage{iitblogo}}
 
+\AtBeginSection[]
+{
+  \begin{frame}<beamer>
+    \frametitle{Outline}      
+    \tableofcontents[currentsection,currentsubsection]
+  \end{frame}
+}
 
 %% Delete this, if you do not want the table of contents to pop up at
 %% the beginning of each subsection:
@@ -95,109 +125,434 @@
   \end{frame}
 }
 
-
-% If you wish to uncover everything in a step-wise fashion, uncomment
-% the following command: 
-%\beamerdefaultoverlayspecification{<+->}
-
-%\includeonlyframes{current,current1,current2,current3,current4,current5,current6}
-
+\AtBeginSection[]
+{
+  \begin{frame}<beamer>
+    \frametitle{Outline}
+    \tableofcontents[currentsection,currentsubsection]
+  \end{frame}
+}
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 % DOCUMENT STARTS
 \begin{document}
 
 \begin{frame}
-  \titlepage
+  \maketitle
+\end{frame}
+
+\begin{frame}
+  \frametitle{Outline}
+  \tableofcontents
+  % You might wish to add the option [pausesections]
 \end{frame}
 
-\begin{frame}{Problem 1.1}
-  The aliquot of a number is defined as: the sum of the \emph{proper} divisors of the number. \\For example: 
-\center{aliquot(12) = 1 + 2 + 3 + 4 + 6 = 16.}\\
-  Write a function that returns the aliquot number of a given number. 
+\section{3D Data Visualization}
+
+\begin{frame}
+    \frametitle{What is visualization?}
+    \Large
+    \begin{center}
+    Visual representation of data
+    \end{center}
 \end{frame}
 
-\begin{frame}{Problem 1.2}
-  Pair of numbers (a, b) is said to be \alert{amicable} if aliquot number of a is b and aliquot number of b is a.\\
-  Example: \texttt{220, 284}\\
-  Write a program that prints all four digit amicable pairs.
-  
-\inctime{20}
+
+%% \begin{frame}
+%%     \frametitle{Is this new?}    
+%%     \begin{center}
+%%     We have moved from:
+%%     \end{center}
+%%     \begin{columns}
+%%     \column{}
+%%     \hspace*{-1in}    
+%%     \includegraphics[width=1.75in,height=1.75in, interpolate=true]{data/3832}      
+%%     \column{}\hspace*{-0.25in}
+%%     To
+%%     \column{}
+%%     \hspace*{-1in}
+%%     \includegraphics[width=1.75in, height=1.75in, interpolate=true]{data/torus}  
+%%     \end{columns}
+%% \end{frame}
+
+\begin{frame}
+    \frametitle{3D visualization}
+    \Large
+    \begin{center}
+        Harder but important
+    \end{center}
 \end{frame}
 
-%% \begin{frame}{Problem 2}
-%%   Given an empty chessboard and one Bishop placed in any s%quare, say (r, c), generate the list of all squares the Bi%shop could move to.
-
-%% \end{frame}
+\begin{frame}
+    \frametitle{Is this Graphics?}
+    \Large
+    \begin{center}
+        Visualization is about data!
+    \end{center}
+\end{frame}
 
-\begin{frame}[fragile]
-  \frametitle{Problem Set 2}
-  Given a string like, ``1, 3-7, 12, 15, 18-21'', produce the list \\
-  \begin{lstlisting}
-    [1,3,4,5,6,7,12,15,18,19,20,21]
-  \end{lstlisting}
+\begin{frame}
+    \frametitle{Examples: trajectory in space}
+    \Large
+    \begin{center}
+        \pgfimage[width=2.5in]{MEDIA/m2/mlab/plot3d_ex}
+    \end{center}
+\end{frame}
+
+\begin{frame}
+    \frametitle{Examples: Fire in a room}
+    \Large
+    \begin{center}
+        Demo of data
+    \end{center}
 \inctime{10}
 \end{frame}
 
-\begin{frame} 
-  \frametitle{Problem Set 3}
-  \begin{description}
-    \item[3.1] Count word frequencies in a file.
-\end{description}
-\inctime{5}
+\section{Tools available}
+
+\subsection{mlab}
+
+\begin{frame}
+    {Overview}
+    \Large
+    \begin{itemize}
+        \item Simple
+        \item Convenient
+        \item Full-featured
+    \end{itemize}
+\end{frame}
+
+\begin{frame}[fragile]
+
+    \frametitle{Getting started}
+    \myemph{\Large Vanilla:}
+    \begin{lstlisting}[language=bash]
+        $ ipython -wthread
+    \end{lstlisting}
+    \myemph{\Large with Pylab:}
+    \begin{lstlisting}[language=bash]
+        $ ipython -pylab -wthread
+    \end{lstlisting}
+\end{frame}
+
+\begin{frame}[fragile]
+    \frametitle{Using mlab}
+
+    \begin{lstlisting}
+In []:from enthought.mayavi import mlab
+    \end{lstlisting}
+
+    \vspace*{0.5in}
+
+    \myemph{\Large Try these}
+
+    \vspace*{0.25in}
+
+    \begin{lstlisting}
+In []: mlab.test_<TAB>
+In []: mlab.test_contour3d()
+In []: mlab.test_contour3d??
+    \end{lstlisting}
+\end{frame}
+
+\begin{frame}
+    {Exploring the view}
+    \begin{columns}
+        \column{0.6\textwidth}
+    \pgfimage[width=3in]{MEDIA/m2/contour3d}
+        \column{0.4\textwidth}
+        \begin{itemize}
+            \item Mouse
+            \item Keyboard
+            \item Toolbar
+            \item Mayavi icon\pgfimage[width=0.2in]{MEDIA/m2/m2_icon}
+        \end{itemize}
+    \end{columns}
+\end{frame}
+
+\begin{frame}[fragile]
+    \frametitle{\mlab\ plotting functions}
+    \begin{columns}
+        \column{0.25\textwidth}
+        \myemph{\Large 0D data}
+        \column{0.5\textwidth}
+    \pgfimage[width=2in]{MEDIA/m2/mlab/points3d_ex}
+    \end{columns}
+
+    \begin{lstlisting}
+In []: t = linspace(0, 2*pi, 50)
+In []: u = cos(t) * pi
+In []: x, y, z = sin(u), cos(u), sin(t)
+    \end{lstlisting}
+    \emphbar{\PythonCode{In []: mlab.points3d(x, y, z)}}
+\end{frame}
+
+\begin{frame}
+  \begin{columns}
+        \column{0.25\textwidth}
+        \myemph{\Large 1D data}
+        \column{0.5\textwidth}
+        \pgfimage[width=2.5in]{MEDIA/m2/mlab/plot3d_ex}
+  \end{columns}
+  \emphbar{\PythonCode{In []: mlab.plot3d(x, y, z, t)}}
+
+    Plots lines between the points
+    
+\end{frame}
+
+\begin{frame}[fragile]
+    \begin{columns}
+        \column{0.25\textwidth}
+        \myemph{\Large 2D data}
+        \column{0.5\textwidth}
+        \pgfimage[width=2in]{MEDIA/m2/mlab/surf_ex}
+    \end{columns}            
+    \begin{lstlisting}
+In []: x, y = mgrid[-3:3:100j,-3:3:100j]
+In []: z = sin(x*x + y*y)
+    \end{lstlisting}
+
+    \emphbar{\PythonCode{In []: mlab.surf(x, y, z)}}
+
+    \alert{Assumes the points are rectilinear}
+
+\end{frame}
+
+\begin{frame}[fragile]
+  \frametitle{mgrid}
+  \begin{lstlisting}
+In []: mgrid[0:3,0:3]
+Out[]: 
+array([[[0, 0, 0],
+        [1, 1, 1],
+        [2, 2, 2]],
+
+       [[0, 1, 2],
+        [0, 1, 2],
+        [0, 1, 2]]])
+
+In []: mgrid[-1:1:5j]
+Out[]: array([-1., -0.5,  0.,  0.5,  1.])
+\end{lstlisting}
 \end{frame}
 
 \begin{frame}[fragile]
-  \frametitle{Problem set 4}
-  Central difference
-  \begin{equation*}
-  \frac{sin(x+h)-sin(x-h)}{2h}
-  \end{equation*}
+  \frametitle{Example}
   \begin{lstlisting}
-  In []: x = linspace(0, 2*pi, 100)
-  In []: y = sin(x)
-  In []: deltax = x[1] - x[0]
-  \end{lstlisting}
+In []: x, y = mgrid[-1:1:5j, -1:1:5j]
+In []: z = x*x + y*y
+
+In []: z
+Out[]: 
+array([[ 2.  , 1.25, 1.  , 1.25, 2.  ],
+       [ 1.25, 0.5 , 0.25, 0.5 , 1.25],
+       [ 1.  , 0.25, 0.  , 0.25, 1.  ],
+       [ 1.25, 0.5 , 0.25, 0.5 , 1.25],
+       [ 2.  , 1.25, 1.  , 1.25, 2.  ]])
+\end{lstlisting}
+\end{frame}
+
+\begin{frame}[fragile]
+    \myemph{\Large 2D data: \texttt{mlab.mesh}}
+    \vspace*{0.25in}
+
+    \emphbar{\PythonCode{In []: mlab.mesh(x, y, z)}}
+
+    \alert{Points needn't be regular}
+
+    \vspace*{0.25in}
+\begin{lstlisting}
+In []: phi, theta = mgrid[0:pi:20j, 
+...                         0:2*pi:20j]
+In []: x = sin(phi)*cos(theta)
+In []: y = sin(phi)*sin(theta)
+In []: z = cos(phi)
+In []: mlab.mesh(x, y, z, 
+...           representation=
+...           'wireframe')
+\end{lstlisting}
+
+\end{frame}
+
+\begin{frame}[fragile]
+
+  \begin{columns}
+        \column{0.25\textwidth}
+        \myemph{\Large 3D data}
+        \column{0.5\textwidth}
+        \pgfimage[width=1.5in]{MEDIA/m2/mlab/contour3d}\\        
+    \end{columns}
+\begin{lstlisting}
+In []: x, y, z = mgrid[-5:5:64j, 
+...                -5:5:64j, 
+...                -5:5:64j]
+In []: mlab.contour3d(x*x*0.5 + y*y + 
+                   z*z*2)
+\end{lstlisting}
+\end{frame}
+
+\begin{frame}[fragile]
+
+    \myemph{\Large 3D vector data: \PythonCode{mlab.quiver3d}}
+    \vspace*{0.25in}
+
+    \pgfimage[width=2in]{MEDIA/m2/mlab/quiver3d_ex}\\
+    
+\begin{lstlisting}
+In []: mlab.test_quiver3d()
+\end{lstlisting}
+
+\emphbar{\PythonCode{obj = mlab.quiver3d(x, y, z, u, v, w)}}
+\inctime{20}
+\end{frame}
+
+
+\subsection{Mayavi2}
+
+\begin{frame}
+  \frametitle{Introduction to Mayavi}
+  \begin{itemize}
+  \item Most scientists not interested in details of visualization
+  \item Visualization of data files with a nice UI
+  \item Interactive visualization of data (think Matlab)
+  \item Embedding visualizations in applications
+  \item Customization
+  \end{itemize}
   \pause
-    \begin{enumerate}
-      \item Given this, get the finite difference of sin in the range 0 to 2*pi
-    \end{enumerate}
+  \begin{block}{The Goal}
+      Provide a \alert{flexible} library/app for all of these needs!
+  \end{block}
+\end{frame}
+
+\begin{frame}
+    {Overview of features}
+      \vspace*{-0.3in}
+  \begin{center}    
+    \hspace*{-0.2in}\pgfimage[width=5in]{MEDIA/m2/m2_app3_3}
+  \end{center}    
+\end{frame}
+
+
+\begin{frame}
+    \frametitle{Mayavi in applications}
+      \vspace*{-0.3in}
+  \begin{center}    
+    \hspace*{-0.2in}\pgfimage[width=4.5in]{MEDIA/m2/m2_envisage}
+  \end{center}
 \end{frame}
 
 \begin{frame}
-  \frametitle{Problem Set 5}
-  \begin{itemize}
-      \item[5.1] Write a function that plots any regular n-gon given \typ{n}.
-      \item[5.2] Consider the logistic map, $f(x) = kx(1-x)$, plot it for
-          $k=2.5, 3.5$ and $4$ in the same plot.
-\end{itemize}
+    \frametitle{Live in your dialogs}
+      \vspace*{0.1in}
+  \begin{center}    
+    \hspace*{-0.2in}\pgfimage[width=2.5in]{MEDIA/m2/mlab_tui}
+  \end{center}
+\end{frame}
+
+\begin{frame}
+    {Exploring the documentation}
+    \begin{center}
+    \pgfimage[width=4in]{MEDIA/m2/m2_ug_doc}
+    \end{center}
 \end{frame}
 
-\begin{frame}[fragile] 
-\frametitle{Problem Set 5}
-  \begin{columns}
-    \column{0.6\textwidth}
-    \small{
-    \begin{itemize}
-      \item[3] Consider the iteration $x_{n+1} = f(x_n)$ where $f(x) = kx(1-x)$.  Plot the successive iterates of this process as explained below. 
-    \end{itemize}}
-    \column{0.35\textwidth}
-    \hspace*{-0.5in}
-  \includegraphics[height=1.6in, interpolate=true]{data/cobweb}  
-\end{columns}
+
+\begin{frame}
+  \frametitle{Summary}
+      \begin{itemize}
+          \item \url{http://code.enthought.com/projects/mayavi}
+          \item Uses VTK (\url{www.vtk.org})
+          \item BSD license
+          \item Linux, win32 and Mac OS X
+          \item Highly scriptable
+          \item Embed in Traits UIs (wxPython and PyQt4)
+          \item Envisage Plugins
+          \item Debian/Ubuntu/Fedora
+          \item \alert{Pythonic}
+      \end{itemize}
+    
+      \inctime{10}
+
 \end{frame}
 
 \begin{frame}
-  \frametitle{Problem Set 5.3}
-  Plot the cobweb plot as follows:
-  \begin{enumerate}
-    \item Start at $(x_0, 0)$ ($\implies$ i=0)
-    \item Draw a line to $(x_i, f(x_i))$
-    \item Set $x_{i+1} = f(x_i)$
-    \item Draw a line to $(x_{i+1}, x_{i+1})$
-    \item $(i\implies i+1)$ 
-    \item Repeat from 2 for as long as you want 
-  \end{enumerate}
+    {Getting hands dirty!}
+
+        \begin{block}{Motivational problem}
+        Atmospheric data of temperature over the surface of the earth.
+        Let temperature ($T$) vary linearly with height ($z$):
+        \begin{center}            
+        $T = 288.15 - 6.5z$
+        \end{center}
+        \end{block}
+\end{frame}
+
+\begin{frame}[fragile]
+    \frametitle{Simple solution}
+
+    \begin{lstlisting}
+lat = linspace(-89, 89, 37)
+lon = linspace(0, 360, 37)
+z = linspace(0, 100, 11)
+    \end{lstlisting}
+\pause
+    \begin{lstlisting}
+x, y, z = mgrid[0:360:37j,-89:89:37j,
+                0:100:11j]
+t = 288.15 - 6.5*z
+mlab.contour3d(x, y, z, t)
+mlab.outline()
+mlab.colorbar()
+    \end{lstlisting}
+\end{frame}
+
+\begin{frame}[fragile]
+    \frametitle{Exercise: Lorenz equation}
+    \begin{columns}
+        \column{0.25\textwidth}
+        \begin{eqnarray*}
+        \frac{d x}{dt} &=& s (y-x)\\
+        \frac{d y}{d t} &=& rx -y -xz\\
+        \frac{d z}{d t} &=& xy - bz\\
+        \end{eqnarray*}
+        \column{0.25\textwidth}
+        Let $s=10,$
+        $r=28,$ 
+        $b=8./3.$
+    \end{columns}
+    \structure{\Large Region of interest}
+  \begin{lstlisting}
+x, y, z = mgrid[-50:50:20j,-50:50:20j,
+                -10:60:20j]
+  \end{lstlisting}
 \inctime{20}
+
+\end{frame}
+\begin{frame}[fragile]
+    \frametitle{Solution}
+  \begin{lstlisting}
+def lorenz(x,y,z,s=10.,r=28.,b=8./3.):
+    u = s*(y-x)
+    v = r*x-y-x*z
+    w = x*y-b*z
+    return u,v,w
+x,y,z = mgrid [-50:50:20j,-50:50:20j,
+                    -10:60:20j ]
+u,v,w = lorenz( x , y , z )
+# Your plot here
+#
+mlab.show()
+
+  \end{lstlisting}
+\end{frame}
+
+\begin{frame}
+  \frametitle{We have covered:}
+  \begin{itemize}
+  \item Need of visualization.
+  \item Using mlab to create 3 D plots.
+  \item Mayavi Toolkit.
+  \end{itemize}
 \end{frame}
 
 \end{document}
+
diff -r 2d08c45681a1 -r a835affb1447 day2/session6.tex
--- a/day2/session6.tex	Tue Nov 10 14:33:06 2009 +0530
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,558 +0,0 @@
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% Tutorial slides on Python.
-%
-% Author: Prabhu Ramachandran <prabhu at aero.iitb.ac.in>
-% Copyright (c) 2005-2009, Prabhu Ramachandran
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-\documentclass[compress,14pt]{beamer}
-% \documentclass[handout]{beamer}
-% \usepackage{pgfpages}
-% \pgfpagesuselayout{4 on 1}[a4paper,border, shrink=5mm,landscape]
-\usepackage{tikz}
-\newcommand{\hyperlinkmovie}{}
-%\usepackage{movie15}
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% Note that in presentation mode 
-% \paperwidth  364.19536pt
-% \paperheight 273.14662pt
-% h/w = 0.888
-
-
-\mode<presentation>
-{
-  \usetheme{Warsaw}
-  %\usetheme{Boadilla}
-  %\usetheme{default}
-  \useoutertheme{infolines}
-  \setbeamercovered{transparent}
-}
-
-% To remove navigation symbols
-\setbeamertemplate{navigation symbols}{}
-
-\usepackage{amsmath}
-\usepackage[english]{babel}
-\usepackage[latin1]{inputenc}
-\usepackage{times}
-\usepackage[T1]{fontenc}
-
-% Taken from Fernando's slides.
-\usepackage{ae,aecompl}
-\usepackage{mathpazo,courier,euler}
-\usepackage[scaled=.95]{helvet}
-\usepackage{pgf}
-
-\definecolor{darkgreen}{rgb}{0,0.5,0}
-
-\usepackage{listings}
-\lstset{language=Python,
-    basicstyle=\ttfamily\bfseries,
-    commentstyle=\color{red}\itshape,
-  stringstyle=\color{darkgreen},
-  showstringspaces=false,
-  keywordstyle=\color{blue}\bfseries}
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% My Macros
-\setbeamercolor{postit}{bg=yellow,fg=black}
-\setbeamercolor{emphbar}{bg=blue!20, fg=black}
-\newcommand{\emphbar}[1]
-{\begin{beamercolorbox}[rounded=true]{emphbar} 
-      {#1}
- \end{beamercolorbox}
-}
-%{\centerline{\fcolorbox{gray!50} {blue!10}{
-%\begin{minipage}{0.9\linewidth}
-%    {#1} 
-%\end{minipage}
-%    }}}
-
-\newcommand{\myemph}[1]{\structure{\emph{#1}}}
-\newcommand{\PythonCode}[1]{\lstinline{#1}}
-
-\newcommand{\tvtk}{\texttt{tvtk}}
-\newcommand{\mlab}{\texttt{mlab}}
-
-\newcounter{time}
-\setcounter{time}{0}
-\newcommand{\inctime}[1]{\addtocounter{time}{#1}{\vspace*{0.1in}\tiny \thetime\ m}}
-
-\newcommand\BackgroundPicture[1]{%
-  \setbeamertemplate{background}{%
-      \parbox[c][\paperheight]{\paperwidth}{%
-      \vfill \hfill
- \hfill \vfill
-}}}
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% Configuring the theme
-%\setbeamercolor{normal text}{fg=white}
-%\setbeamercolor{background canvas}{bg=black}
-
-
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% Title page
-\title[3D Plotting]{3D data Visualization}
-
-\author[FOSSEE] {FOSSEE}
-
-\institute[IIT Bombay] {Department of Aerospace Engineering\\IIT Bombay}
-\date[] {8 November, 2009\\Day 2, Session 6}
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-%\pgfdeclareimage[height=0.75cm]{iitblogo}{iitblogo}
-%\logo{\pgfuseimage{iitblogo}}
-
-\AtBeginSection[]
-{
-  \begin{frame}<beamer>
-    \frametitle{Outline}      
-    \tableofcontents[currentsection,currentsubsection]
-  \end{frame}
-}
-
-%% Delete this, if you do not want the table of contents to pop up at
-%% the beginning of each subsection:
-\AtBeginSubsection[]
-{
-  \begin{frame}<beamer>
-    \frametitle{Outline}
-    \tableofcontents[currentsection,currentsubsection]
-  \end{frame}
-}
-
-\AtBeginSection[]
-{
-  \begin{frame}<beamer>
-    \frametitle{Outline}
-    \tableofcontents[currentsection,currentsubsection]
-  \end{frame}
-}
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% DOCUMENT STARTS
-\begin{document}
-
-\begin{frame}
-  \maketitle
-\end{frame}
-
-\begin{frame}
-  \frametitle{Outline}
-  \tableofcontents
-  % You might wish to add the option [pausesections]
-\end{frame}
-
-\section{3D Data Visualization}
-
-\begin{frame}
-    \frametitle{What is visualization?}
-    \Large
-    \begin{center}
-    Visual representation of data
-    \end{center}
-\end{frame}
-
-
-%% \begin{frame}
-%%     \frametitle{Is this new?}    
-%%     \begin{center}
-%%     We have moved from:
-%%     \end{center}
-%%     \begin{columns}
-%%     \column{}
-%%     \hspace*{-1in}    
-%%     \includegraphics[width=1.75in,height=1.75in, interpolate=true]{data/3832}      
-%%     \column{}\hspace*{-0.25in}
-%%     To
-%%     \column{}
-%%     \hspace*{-1in}
-%%     \includegraphics[width=1.75in, height=1.75in, interpolate=true]{data/torus}  
-%%     \end{columns}
-%% \end{frame}
-
-\begin{frame}
-    \frametitle{3D visualization}
-    \Large
-    \begin{center}
-        Harder but important
-    \end{center}
-\end{frame}
-
-\begin{frame}
-    \frametitle{Is this Graphics?}
-    \Large
-    \begin{center}
-        Visualization is about data!
-    \end{center}
-\end{frame}
-
-\begin{frame}
-    \frametitle{Examples: trajectory in space}
-    \Large
-    \begin{center}
-        \pgfimage[width=2.5in]{MEDIA/m2/mlab/plot3d_ex}
-    \end{center}
-\end{frame}
-
-\begin{frame}
-    \frametitle{Examples: Fire in a room}
-    \Large
-    \begin{center}
-        Demo of data
-    \end{center}
-\inctime{10}
-\end{frame}
-
-\section{Tools available}
-
-\subsection{mlab}
-
-\begin{frame}
-    {Overview}
-    \Large
-    \begin{itemize}
-        \item Simple
-        \item Convenient
-        \item Full-featured
-    \end{itemize}
-\end{frame}
-
-\begin{frame}[fragile]
-
-    \frametitle{Getting started}
-    \myemph{\Large Vanilla:}
-    \begin{lstlisting}[language=bash]
-        $ ipython -wthread
-    \end{lstlisting}
-    \myemph{\Large with Pylab:}
-    \begin{lstlisting}[language=bash]
-        $ ipython -pylab -wthread
-    \end{lstlisting}
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{Using mlab}
-
-    \begin{lstlisting}
-In []:from enthought.mayavi import mlab
-    \end{lstlisting}
-
-    \vspace*{0.5in}
-
-    \myemph{\Large Try these}
-
-    \vspace*{0.25in}
-
-    \begin{lstlisting}
-In []: mlab.test_<TAB>
-In []: mlab.test_contour3d()
-In []: mlab.test_contour3d??
-    \end{lstlisting}
-\end{frame}
-
-\begin{frame}
-    {Exploring the view}
-    \begin{columns}
-        \column{0.6\textwidth}
-    \pgfimage[width=3in]{MEDIA/m2/contour3d}
-        \column{0.4\textwidth}
-        \begin{itemize}
-            \item Mouse
-            \item Keyboard
-            \item Toolbar
-            \item Mayavi icon\pgfimage[width=0.2in]{MEDIA/m2/m2_icon}
-        \end{itemize}
-    \end{columns}
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{\mlab\ plotting functions}
-    \begin{columns}
-        \column{0.25\textwidth}
-        \myemph{\Large 0D data}
-        \column{0.5\textwidth}
-    \pgfimage[width=2in]{MEDIA/m2/mlab/points3d_ex}
-    \end{columns}
-
-    \begin{lstlisting}
-In []: t = linspace(0, 2*pi, 50)
-In []: u = cos(t) * pi
-In []: x, y, z = sin(u), cos(u), sin(t)
-    \end{lstlisting}
-    \emphbar{\PythonCode{In []: mlab.points3d(x, y, z)}}
-\end{frame}
-
-\begin{frame}
-  \begin{columns}
-        \column{0.25\textwidth}
-        \myemph{\Large 1D data}
-        \column{0.5\textwidth}
-        \pgfimage[width=2.5in]{MEDIA/m2/mlab/plot3d_ex}
-  \end{columns}
-  \emphbar{\PythonCode{In []: mlab.plot3d(x, y, z, t)}}
-
-    Plots lines between the points
-    
-\end{frame}
-
-\begin{frame}[fragile]
-    \begin{columns}
-        \column{0.25\textwidth}
-        \myemph{\Large 2D data}
-        \column{0.5\textwidth}
-        \pgfimage[width=2in]{MEDIA/m2/mlab/surf_ex}
-    \end{columns}            
-    \begin{lstlisting}
-In []: x, y = mgrid[-3:3:100j,-3:3:100j]
-In []: z = sin(x*x + y*y)
-    \end{lstlisting}
-
-    \emphbar{\PythonCode{In []: mlab.surf(x, y, z)}}
-
-    \alert{Assumes the points are rectilinear}
-
-\end{frame}
-
-\begin{frame}[fragile]
-  \frametitle{mgrid}
-  \begin{lstlisting}
-In []: mgrid[0:3,0:3]
-Out[]: 
-array([[[0, 0, 0],
-        [1, 1, 1],
-        [2, 2, 2]],
-
-       [[0, 1, 2],
-        [0, 1, 2],
-        [0, 1, 2]]])
-
-In []: mgrid[-1:1:5j]
-Out[]: array([-1., -0.5,  0.,  0.5,  1.])
-\end{lstlisting}
-\end{frame}
-
-\begin{frame}[fragile]
-  \frametitle{Example}
-  \begin{lstlisting}
-In []: x, y = mgrid[-1:1:5j, -1:1:5j]
-In []: z = x*x + y*y
-
-In []: z
-Out[]: 
-array([[ 2.  , 1.25, 1.  , 1.25, 2.  ],
-       [ 1.25, 0.5 , 0.25, 0.5 , 1.25],
-       [ 1.  , 0.25, 0.  , 0.25, 1.  ],
-       [ 1.25, 0.5 , 0.25, 0.5 , 1.25],
-       [ 2.  , 1.25, 1.  , 1.25, 2.  ]])
-\end{lstlisting}
-\end{frame}
-
-\begin{frame}[fragile]
-    \myemph{\Large 2D data: \texttt{mlab.mesh}}
-    \vspace*{0.25in}
-
-    \emphbar{\PythonCode{In []: mlab.mesh(x, y, z)}}
-
-    \alert{Points needn't be regular}
-
-    \vspace*{0.25in}
-\begin{lstlisting}
-In []: phi, theta = mgrid[0:pi:20j, 
-...                         0:2*pi:20j]
-In []: x = sin(phi)*cos(theta)
-In []: y = sin(phi)*sin(theta)
-In []: z = cos(phi)
-In []: mlab.mesh(x, y, z, 
-...           representation=
-...           'wireframe')
-\end{lstlisting}
-
-\end{frame}
-
-\begin{frame}[fragile]
-
-  \begin{columns}
-        \column{0.25\textwidth}
-        \myemph{\Large 3D data}
-        \column{0.5\textwidth}
-        \pgfimage[width=1.5in]{MEDIA/m2/mlab/contour3d}\\        
-    \end{columns}
-\begin{lstlisting}
-In []: x, y, z = mgrid[-5:5:64j, 
-...                -5:5:64j, 
-...                -5:5:64j]
-In []: mlab.contour3d(x*x*0.5 + y*y + 
-                   z*z*2)
-\end{lstlisting}
-\end{frame}
-
-\begin{frame}[fragile]
-
-    \myemph{\Large 3D vector data: \PythonCode{mlab.quiver3d}}
-    \vspace*{0.25in}
-
-    \pgfimage[width=2in]{MEDIA/m2/mlab/quiver3d_ex}\\
-    
-\begin{lstlisting}
-In []: mlab.test_quiver3d()
-\end{lstlisting}
-
-\emphbar{\PythonCode{obj = mlab.quiver3d(x, y, z, u, v, w)}}
-\inctime{20}
-\end{frame}
-
-
-\subsection{Mayavi2}
-
-\begin{frame}
-  \frametitle{Introduction to Mayavi}
-  \begin{itemize}
-  \item Most scientists not interested in details of visualization
-  \item Visualization of data files with a nice UI
-  \item Interactive visualization of data (think Matlab)
-  \item Embedding visualizations in applications
-  \item Customization
-  \end{itemize}
-  \pause
-  \begin{block}{The Goal}
-      Provide a \alert{flexible} library/app for all of these needs!
-  \end{block}
-\end{frame}
-
-\begin{frame}
-    {Overview of features}
-      \vspace*{-0.3in}
-  \begin{center}    
-    \hspace*{-0.2in}\pgfimage[width=5in]{MEDIA/m2/m2_app3_3}
-  \end{center}    
-\end{frame}
-
-
-\begin{frame}
-    \frametitle{Mayavi in applications}
-      \vspace*{-0.3in}
-  \begin{center}    
-    \hspace*{-0.2in}\pgfimage[width=4.5in]{MEDIA/m2/m2_envisage}
-  \end{center}
-\end{frame}
-
-\begin{frame}
-    \frametitle{Live in your dialogs}
-      \vspace*{0.1in}
-  \begin{center}    
-    \hspace*{-0.2in}\pgfimage[width=2.5in]{MEDIA/m2/mlab_tui}
-  \end{center}
-\end{frame}
-
-\begin{frame}
-    {Exploring the documentation}
-    \begin{center}
-    \pgfimage[width=4in]{MEDIA/m2/m2_ug_doc}
-    \end{center}
-\end{frame}
-
-
-\begin{frame}
-  \frametitle{Summary}
-      \begin{itemize}
-          \item \url{http://code.enthought.com/projects/mayavi}
-          \item Uses VTK (\url{www.vtk.org})
-          \item BSD license
-          \item Linux, win32 and Mac OS X
-          \item Highly scriptable
-          \item Embed in Traits UIs (wxPython and PyQt4)
-          \item Envisage Plugins
-          \item Debian/Ubuntu/Fedora
-          \item \alert{Pythonic}
-      \end{itemize}
-    
-      \inctime{10}
-
-\end{frame}
-
-\begin{frame}
-    {Getting hands dirty!}
-
-        \begin{block}{Motivational problem}
-        Atmospheric data of temperature over the surface of the earth.
-        Let temperature ($T$) vary linearly with height ($z$):
-        \begin{center}            
-        $T = 288.15 - 6.5z$
-        \end{center}
-        \end{block}
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{Simple solution}
-
-    \begin{lstlisting}
-lat = linspace(-89, 89, 37)
-lon = linspace(0, 360, 37)
-z = linspace(0, 100, 11)
-    \end{lstlisting}
-\pause
-    \begin{lstlisting}
-x, y, z = mgrid[0:360:37j,-89:89:37j,
-                0:100:11j]
-t = 288.15 - 6.5*z
-mlab.contour3d(x, y, z, t)
-mlab.outline()
-mlab.colorbar()
-    \end{lstlisting}
-\end{frame}
-
-\begin{frame}[fragile]
-    \frametitle{Exercise: Lorenz equation}
-    \begin{columns}
-        \column{0.25\textwidth}
-        \begin{eqnarray*}
-        \frac{d x}{dt} &=& s (y-x)\\
-        \frac{d y}{d t} &=& rx -y -xz\\
-        \frac{d z}{d t} &=& xy - bz\\
-        \end{eqnarray*}
-        \column{0.25\textwidth}
-        Let $s=10,$
-        $r=28,$ 
-        $b=8./3.$
-    \end{columns}
-    \structure{\Large Region of interest}
-  \begin{lstlisting}
-x, y, z = mgrid[-50:50:20j,-50:50:20j,
-                -10:60:20j]
-  \end{lstlisting}
-\inctime{20}
-
-\end{frame}
-\begin{frame}[fragile]
-    \frametitle{Solution}
-  \begin{lstlisting}
-def lorenz(x,y,z,s=10.,r=28.,b=8./3.):
-    u = s*(y-x)
-    v = r*x-y-x*z
-    w = x*y-b*z
-    return u,v,w
-x,y,z = mgrid [-50:50:20j,-50:50:20j,
-                    -10:60:20j ]
-u,v,w = lorenz( x , y , z )
-# Your plot here
-#
-mlab.show()
-
-  \end{lstlisting}
-\end{frame}
-
-\begin{frame}
-  \frametitle{We have covered:}
-  \begin{itemize}
-  \item Need of visualization.
-  \item Using mlab to create 3 D plots.
-  \item Mayavi Toolkit.
-  \end{itemize}
-\end{frame}
-
-\end{document}
-