The question is not clear to me, why not simply do:
sage: Y = [1,2,4,3,7,2,1]
sage: h = 0.1
sage: [(Y[k+1]-Y[k-1])/h for k in range(1,len(Y)-1)]
[30.0000000000000,
10.0000000000000,
30.0000000000000,
-10.0000000000000,
-60.0000000000000]
$y_{k+1} $ and $y_k$ needs to remain unknown till the formation of equation, so that they can act as variables in the linear system. My interest is to automate the finite difference method in sage. Should I need to define some new data structure. Thanks for responding though!!!
Asked: 2016-01-25 12:09:28 -0500
Seen: 95 times
Last updated: Jan 26 '16
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.
