#!/usr/bin/python
from numpy import *
from pylab import *
from ising import *
h=open('data.txt','w')
hs=open('stru.txt','w')
h1=open('data1.txt','w')
hs1=open('stru1.txt','w')
h2=open('data2.txt','w')
hs2=open('stru2.txt','w')
h3=open('data3.txt','w')
hs3=open('stru3.txt','w')
h4=open('data4.txt','w')
hs4=open('stru4.txt','w')
h5=open('data5.txt','w')
hs5=open('stru5.txt','w')
h6=open('data6.txt','w')
hs6=open('stru6.txt','w')
h7=open('data7.txt','w')
hs7=open('stru7.txt','w')
h8=open('data8.txt','w')
hs8=open('stru8.txt','w')
h9=open('data9.txt','w')
hs9=open('stru9.txt','w')
h10=open('data10.txt','w')
hs10=open('stru10.txt','w')
h11=open('data11.txt','w')
hs11=open('stru11.txt','w')


value,S=ising(3000,100,-4.,0.5,0.)
h.write(str(value))
hs.write(str(S))
value1,S1=ising(3000,100,-4.,1.0,0.)
value2,S2=ising(3000,100,-4.,1.5,0.)
value3,S3=ising(3000,100,-4.,2.0,0.)
value4,S4=ising(3000,100,-4.,2.1,0.)
value5,S5=ising(3000,100,-4.,2.2,0.)
value6,S6=ising(3000,100,-4.,2.269,0.)
value7,S7=ising(3000,100,-4.,2.4,0.)
value8,S8=ising(3000,100,-4.,2.5,0.)
value9,S9=ising(3000,100,-4.,2.6,0.)
value10,S10=ising(3000,100,-4.,2.8,0.)
value11,S11=ising(3000,100,-4.,3.0,0.)
#pcolor(S)
#show()
#t=arange(1,101,1)
#figure()
#plot(t,value[:,-1],'b')
#figure()
#plot(t,value[:,1],'g')