import numpy as np
import pandas as pd 
import math
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt



def test_run():
	d1=np.random.randint(500,size=(1000,2))
	pd1=pd.DataFrame(d1, columns=['x1','x2'])
	
	for i in range(len(pd1)):
		pd1.ix[i,'x2']=np.random.randint(500)
		pd1['y']=1000*pd1['x1']+(-1000)*pd1['x2']+np.random.uniform()
	'''
	threedee = plt.figure().gca(projection='3d')
	threedee.scatter(pd1['x1'], pd1['x2'], pd1['y'])
	threedee.set_xlabel('x1')
	threedee.set_ylabel('x2')
	threedee.set_zlabel('y')
	plt.show()
	'''
	
	pd1.to_csv('Data/best4linreg.csv',header=False,index=False)
	
	
	
	

if __name__=="__main__":
    test_run()
