# Python skript for Taylor-rekke eksempel,
# funksjonen f(x,y) = (1-y^2)*cos(x) Taylor-utvikles til andre ordens ledd
# 1 - 1/2*x^2 - y^2

from scitools.std import *
# Det ser ut til at matplotlib backend ikke klarer flere flater samtidig,
# derfor bruker vi heller gnuplot backend
from scitools.easyviz.gnuplot_ import *

# rutenett med -5 <= x <= 5 og -4 <= y <= 4 og med henholdsvis 26 og 21 punkt
x,y = ndgrid (linspace (-5, 5, 26), linspace (-4, 4, 21))
surf (x, y, (1 - y**2)*cos(x))
xlabel ('x')
ylabel ('y')
hold('on')
surf (x, y, 1 - 1./2*x**2 - y**2)
hold('off')
raw_input('Press Return key to quit: ')

# Dessverre ser rotasjon av de to flatene ikke riktig ut med gnuplot backend,
# og med matplotlib backend er kun en flate synlig om gangen.

# Fra kommandolinje: scitools file2interactive taylor2.py