Euler's method for a differential
... i am only interested in writting the information if its at one of these
values... problem is, i keep getting a syntax error with the if/then statements
--------------------------------
fout = open('/selftest.txt','w')
var('t y')
f(t,y) = 5-3*y^(1/2)
t = 0.0
y = 2.0
steps = 60
h = 0.05
for i in range(steps +1):
if t == 0:
print 'i=', i, 't=', t, 'y=', y
s ='i='+repr(i)+' t='+repr(t)+' y='+repr(y) +'\n'
fout.write(s)
y = (f(t,y)*h)+y
t = t+h
elif (t == 1.0):
print 'i=', i, 't=', t, 'y=', y
s ='i='+repr(i)+' t='+repr(t)+' y='+repr(y) +'\n'
fout.write(s)
y = (f(t,y)*h)+y
t = t+h
elif (t == 1.5):
print 'i=', i, 't=', t, 'y=', y
s ='i='+repr(i)+' t='+repr(t)+' y='+repr(y) +'\n'
fout.write(s)
y = (f(t,y)*h)+y
t = t+h
else:
y = (f(t,y)*h)+y
t = t+h
fout.close()
--
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
