f=open('Results.dat', 'w')
v=0.000005
h=0.001
dx=h/10000
dt=dx/(2*v)
c=[0]*100
c[0]=1
C=[0]*100
k=0
while k<=100:
for i in range(1, len(c)):
C[i]=c[i]+(dt*v*(c[i-1]-c[i]))/dx
res=repr(c)
f.write(res+'\n\n')
for j in range(1, len(C)):
c[j]=C[j]
k+=1
print c
f.close()