Home

<b>Executing the code will leave the page unresponsive for a few minutes.</b> Spin Description // VERY INPORTANT // you have to let the program run for a long time // then plot particle size vs the average in an excel // spread sheet and do scatter plot followed by trend line // to get the linear equation, solve for x when y=2.0023193 // and you should get something near 1820 // to get more accuracy do the following // make m instead of 200 to 100 and set do=d1=1500 // then do the line fitting again . This willl take 5 times longer (1 hour) // run several times and average x to get more accurate results //javascript program must compute all points before displaying all results

This is just a replacement in case Javascript is not available or used for SEO purposes

// VERY INPORTANT

// you have to let the program run for a long time 
// then plot particle size vs the average in an excel 
// spread sheet and do scatter plot followed by trend line
// to get the linear equation, solve for x when y=2.0023193
// and you should get something near 1820
// to get more accuracy do the following
// make m instead of 200 to 100 and set do=d1=1500
// then do the line fitting again . This willl take 5 times longer (1 hour)
// run several times and average x to get more accurate results

//javascript program must compute all points before displaying all results
 
//insert code here
 function GraphIt() {
 
	var newElement = document.createElement('p');
    var L = 1000000;
    var w=1;
    var f1  =0;
    var f = 0;
    var ii = 0;
    var q = 0;
    var g = 0;
    var en = 0;
    var en1 = 0;
    var edx = 0;
    var edx1 = 0;
    var kj = 100000000; // increase for accuracy
    var m = 0;
    var km = 500000;
 
    var d0 = 1700;
    var d1 = 1700;
    var intr = .50 ;
 
 
 
    var rand = new Random();
 
    // create an array 's' and 'l' and initialize all elements to 0
 
//    var fr = new Array();

// for (var i = 0; i <= km;i++) {
//	fr.push([0,0]);
 //  }

for ( var m = 1; m <20; m++) {

d0=d0+10;
				d1=d1+10;

var S = new Array();
    var Sy = new Array();

for (var i = 0; i <= d1;i++) {
		Sy.push(0);
		
 }

for (var i = 0; i <= d0;i++) {
		
		S.push(0);
 }
   
    var dist = 2*m*intr + d0;
  //  var st1 =(L/2)  - m*intr;
  //  var st0 =(L/2) + d0 + m*intr ;

var st1 = (L/2);

var st0 = (L/2) ;
 
	f = 0;	f1 = 0;	edx = 0;i = 0;	edx1 = 0;	en = 0.0;  en1=0;  ent=0;

ii=0; g=0;
	//document.lf.log.value += m+" "+dist+" "+en+"\n";

for ( var i = 1; i <kj; i++) {
        // throw random lines
        var p = d0 * rand.random();
        var li =  ( d0) *rand.random();
        var p1 =  d1 *rand.random();
        var li1 = (d1)* rand.random();

p=Math.floor(p)
          li=Math.floor(li)
           p1=Math.floor(p1)
       li1=Math.floor(li1)

//document.lf.log.value += f+" "+dist+" "+en+"\n";
        // put random lines through conditions
if (st0+p>st1+p1  )

{//e
//document.lf.log.value += f+" "+dist+" "+en+"\n";
if ( st1+p1 + li1 > st0+ p - li ) // this one is em

{//a

g = g + 1;

}//a1

else

{//d

if ((d0-p-li)/d0 > rand.random()  && (p-li)/d0 > rand.random() )

{//b
en = en + li1;
////{//c
Sy[Math.floor(p1)] = Sy[Math.floor(p1)] + 1;
S[Math.floor(p)] = S[Math.floor(p)] +1 ;
f1 = f1 + 1;
f = f + 1;
/////}//c1

}//b1

}//d1
                                                                ///}//e1
}//if p>p1
{//hh
if (st1+p1>st0+p )
{//gg
//document.lf.log.value += f+" "+dist+" "+en+"\n";

if ( st1+p1 - li1 < st0+ p + li) // this one is gravity 
{
g = g + 1;
}
else
{//ff
if ((d0-p-li)/d0 > rand.random()  && (p-li)/d0 > rand.random() )
{
en = en + li1;

Sy[Math.floor(p1)] = Sy[Math.floor(p1)] + 1;
S[Math.floor(p)] = S[Math.floor(p)] +1 ;
f1 = f1 + 1;
f = f + 1;
ii = ii + 1;
     	 // document.lf.log.value += f+" "+S[p]+" "+en+"\n";                                                               //}

}
}//ff
 } //gg                                                                   ///}
}//if p1>p hh
//} //if outside
} // Next j
  //     en = f/en;
       en1 = f/en1;

for ( var ind = 0; ind<=d0; ind++) {
			// CALCULATE EXPECTATION VALUE AFTER INTERACTION HAS TAKEN PLACE
			edx = edx + (ind) * S[ind];
			edx1 = edx1 + (ind) * Sy[ind];

}

var sp= f/ii;
var se = d0/(edx1/f);

//  g  (probability)      spin(expectation)      Average
	document.lf.log.value +=d0+"   "+sp+" "+se+" "+(se+sp)/2+"\n";

}