 
 	//Set integration parameters (end time, number of steps, etc.)
 	const int blockSize = 1024;
-	unsigned int steps = 8*blockSize;
+	unsigned int steps = 32*blockSize;
 	float dt = 0;
 
 	dim3 BIGgridSize(32,32);

			cudaFree(d_beta);
			cudaFree(d_radius);
			cudaFree(d_lengths);
			cudaFree(d_origins);
			cudaFree(d_lines);

//Print some data to screen


	cudaFree(d_radiimin);
	cudaFree(d_radiimax);
	
	//Declaring arrays to save the winding data
	float *d_alpha, *d_beta;
	checkCudaErrors(cudaMalloc(&d_alpha, dataCount*sizeof(float)));
	checkCudaErrors(cudaMalloc(&d_beta, dataCount*sizeof(float)));

	//Set up streams for independent execution
	cudaStream_t RK4, windings, windings2, lengths;
