I think Derintegral is okay, I'm focusing on FunctionGrapher5.py--but you can comment on both.
(To Johan Nilsson: I haven't had time to implement Simpson's rule instead of Reimann's sum yet... but I have gotten it to work faster anyway.)
I will probably expand the code yet... Add extra things like find area between curve and x axis in interval [a,b]...
Any suggestions are greatly appreciated. If you want the code to work on your machine, you'll need to have my Derintegral.py sys.path, obviously.
#### Start of FunctionGrapher5.py ################### from __future__ import division from visual import * import Derintegral import os, random from math import * ja = 0
scenexdefault = 375
def start():
for objects in scene.objects:
objects.visible = 0
scene.title = "Function Grapher by Jacob, Inc."
scene.x = scenexdefault
scene.visible=1
scene.exit=0
scene.userspin = 0
scene.range=(10,10,1)
scene.background=(1,1,1)
global xaxis
global yaxis
xaxis = curve(pos=[(100,0,0),(-100,0,0)],color=color.black)
yaxis = curve(pos=[(0,100,0),(0,-100,0)],color=color.black)
global radiusaxis
global radiusaxis2
radiusaxis = curve(pos=[(-100,-100),(100,100)],color=color.black)
radiusaxis2 = curve(pos=[(-100,100),(100,-100)],color=color.black)
radiusaxis.visible = 0
radiusaxis2.visible = 0start() d = 1 print """\ List of Commands: clear quit remove lines return lines gotl # (Graph One Tangent Line) w/optional x coordinate gatl # (Graph All Tangent Lines) w/optional step gonl # (Graph One Normal Line) w/optional x coordinate ganl # (Graph All Normal Lines) w/optional step
Function Syntax: [y,x,or r] = function [range] ["s" followed by float for step] Brackets mean that it's not required. So, in effect, you don't have to type a function """
def graphit(type,f,range2,step):
li = curve(color=color.blue)
if type in ['y','x']:
x = -15
radiusaxis.visible = 0
radiusaxis2.visible = 0
while -15 <= x <= 15:
if eval(range2):
try:
a = f(x)
if -15 <= a <= 15:
if type == 'y':
li.append(pos=(x,a,0))
elif type == 'x':
li.append(pos=(a,x,0))
else:
li = curve(color=color.blue)
except:
pass
else:
li = curve(color=color.blue)
x = x+step
elif type == 'r':
exec "def m(t): return f(t)*cos(t)"
exec "def n(t): return f(t)*sin(t)"
t = 0
while 0 <= t <= 10*pi:
if eval(range2):
try:
li.append(pos=(m(t),n(t),0))
except:
li = curve(color=color.blue)
t = t+stepprint 'Please type in functions in below. '
while 1:
lists=[]
y = raw_input(">")
if y == 'clear':
scene.visible=0
start()
print "-"*36
continue
elif y == 'quit':
scene.visible = 0
del scene
break
elif y == 'remove lines':
a = [radiusaxis,radiusaxis2,xaxis,yaxis]
for x in a:
x.visible = 0
d = 0
continue
elif y == 'return lines':
a = [radiusaxis,radiusaxis2,xaxis,yaxis]
for x in a:
x.visible = 1
d = 1
continue
elif y.startswith('gatl'):
try:
step2 = float(y.lstrip("gatl "))
except:
step2 = 0.5
x = -20
while -20 <=x<= 20:
try:
der = Derintegral.diffatpt(f,x)
if abs(der)<=25:
if type == 'y':
curve(pos=[(-15,eval("der*(-15-x)+f(x)")),(15,eval("der*(15-x)+f(x)"))],color=color.red)
else:
curve(pos=[(eval("der*(-15-x)+f(x)"),-15),(eval("der*(15-x)+f(x)"),15)],color=color.red)
except:
pass
x = x + step2
continue
elif y.startswith('ganl'):
try:
step2 = float(y.lstrip("ganl "))
except:
step2 = 0.5
x = -20
while -20 <=x<= 20:
try:
der = Derintegral.diffatpt(f,x)
if abs(der)<=25:
if type == 'y':
curve(pos=[(-15,eval("(-1/der)*(-15-x)+f(x)")),(15,eval("(-1/der)*(15-x)+f(x)"))],color = color.red)
else:
curve(pos=[(eval("(-1/der)*(-15-x)+f(x)"),-15),(eval("(-1/der)*(15-x)+f(x)"),15)],color = color.red)
except:
pass
x = x + step2
continue
elif y.startswith('gotl'):
try:
x = float(y.lstrip('gotl '))
except:
x = float(raw_input('What x coordinate do you want the tangent line at? '))
der = Derintegral.diffatpt(f,x)
try:
if abs(der)<= 25:
if type == 'y':
curve(pos=[(-15,eval("der*(-15-x)+f(x)")),(15,eval("der*(15-x)+f(x)"))],color=color.red)
else:
curve(pos=[(eval("der*(-15-x)+f(x)"),-15),(eval("der*(15-x)+f(x)"),15)],color=color.red)
except:
pass
continue
elif y.startswith('gonl'):
try:
x = float(y.lstrip('gonl '))
except:
x = float(raw_input('What x coordinate do you want the tangent line at? '))
der = Derintegral.diffatpt(f,x)
try:
if abs(der)<= 25:
if type == 'y':
curve(pos=[(-15,eval("(-1/der)*(-15-x)+f(x)")),(15,eval("(-1/der)*(15-x)+f(x)"))],color=color.red)
else:
curve(pos=[(eval("(-1/der)*(-15-x)+f(x)"),-15),(eval("(-1/der)*(15-x)+f(x)"),15)],color=color.red)
except:
pass
continue
elif y.startswith('getpt'):
y = y.lstrip("getpt ")
if y:
m = float(eval(y))
else:
m = float(eval(raw_input("f(?) ")))
try:
print f(m)
except ValueError:
print "Math Domain Error"
continue
elif y == 'regraph':
graphit(type,f,range2,step)
continue
if y.count(" = ") == 1:
y = y.split(" = ")
type = y[0].lower()
y = y[1]
y = y.replace("y","x")
if type == 'r':
defaultrange = '0<=t<=5*pi'
y = y.replace('x','t')
if d == 1:
radiusaxis.visible = 1
radiusaxis2.visible = 1
else:
defaultrange = '-15<=x<=15'
else:
type = 'y'
defaultrange = '-15<=x<=15'
y = y.split(" ",1)
tempfunct = y.pop(0)
if y:
y = y[0]
if y.count('st'):
if y.startswith('st'):
y = y.lstrip('st')
step = float(y)
else:
y = y.split(" ")
step = float(y.pop().lstrip('st'))
y = " ".join(y)
range2 = defaultrange
else:
range2 = y
step = 0.005
else:
range2 = defaultrange
step = 0.005
y = tempfunct
range2 = range2.replace('and','or')
range2 = range2.replace(',','<=x<=')
try:
exec "def f(x): return %s" % y
except:
continue
graphit(type,f,range2,step)
###### End of FunctionGrapher5.py ###########
#### Start of Derintegral.py ########## from __future__ import division import psyco psyco.full() from math import *
def diffatpt(funct,pt): """Return the derivative of a function at a specific point. """ sminc = 1e-10 x = pt+sminc y2 = funct(x) x2 = x x = pt-sminc y1 = funct(x) x1 = x slope = (y2-y1)/(x2-x1) return slope
def defintegral(fofx,x,max1):
total = 0
step = 1e-5
exec "def f(x): return %s" % fofx
while x <= max1:
try:
total = total+f(x)
except:
pass
x = x+step
return abs(total*step)
######End of Derintegral.py ##########Like I said, any comments appreciated.
Jacob Schmidt
_______________________________________________ Tutor maillist - [email protected] http://mail.python.org/mailman/listinfo/tutor
