#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <math.h>

const int D=300;

int main (){

	srand(time(NULL));
	int energy[D+1],x,i;
	double R;
	double P=.1;
	//P is the probability to get  a random kick
	int left=0,right=0;

	for(i=0;i<=D;i++)
		energy[i]=i%2;
	//this defines the energy landscape

	double T=2.0;
	//T is the temperature

	for(i=0;i<10000;i++){
		x=150;//starting position
		do{
			R=(double)rand()/RAND_MAX;
			if(R<.5){
				if( (double)rand()/RAND_MAX < P)
					x=x+1;
				else if(energy[x+1]<energy[x])
					x=x+1;
				else if( (double)rand()/RAND_MAX < exp(-(energy[x+1]-energy[x])/T) )
					x=x+1;
			}
			else{
				if( (double)rand()/RAND_MAX < P)
					x=x-1;
				else if(energy[x-1]<energy[x])
					x=x-1;
				else if( (double)rand()/RAND_MAX < exp(-(energy[x-1]-energy[x])/T) )
					x=x-1;
			}
		}while( x!=0 && x!=D );//exits when you reach 0 or D

		if( x==0 )
			left++;
		if ( x==D )
			right++;
	}//for
	printf("Left %d Right %d\n",left,right);

}//main