import sys

if len(sys.argv) == 1:
	print "Usage: \"python rungekutta.py \'$ODE\' '$IC'\" where $ODE is the ODE with y as dependent variable and x as independent variable and IC the initial condition"
	sys.exit(1)

def step(F,y0, h, p):
	k1 = 
	k2 = 
	k3 = 
	y1 = y0 + 

h = 0.25
N = 2
A = [0.,2.]	 #example Moore1965interval p. 97
FA = (float(F.subs(x==A[0]))*h/p,float(F.subs(x==A[1]))*h/p)
print "FA", FA
y0 = [IC,IC]
step(F,y0, h, p)