400.0false600400600.01filenamerutherford$0.mmldate1362464770000url/Users/ddamelin/Documents/CC Stuff/Interactions/MWActiviites/rutherford/rutherford$0.mmlsize24221codebase/Users/ddamelin/Documents/CC Stuff/Interactions/MWActiviites/rutherforddefine %alphaCharge 0.5;
define %tbDx array(%number_of_textboxes);
define %tbDy array(%number_of_textboxes);
define %i 0;
define %flowOn 1;
while (%i < %number_of_textboxes);
set %tbDx[%i] %textbox[%i].x - %atom[0].rx;
set %tbDy[%i] %textbox[%i].y - %atom[0].ry;
%i++;
endwhile;
falseparticleINset %radius %ellipse[0].width/2;
set %particleIn nearest(0,%ellipse[0].x + %radius, %ellipse[0].y + %radius, %radius);
if (%particleIn != -1);
set atom[1].charge 0;
else;
set atom[1].charge %alphaCharge;
endif;stopParticlesset %i 1;
set %offset 1;
set %radius %ellipse[0].width/2;
select atom within (%offset, %offset, "%width-2*%offset", "%height-2*%offset");
charge %alphaCharge;
select atom within (%radius, 0);
charge 0;
select atom 0;
charge 1;
select atom within (0,0,%offset,%height);
set vx 0;
set vy 0;
charge 0;
select atom within (0,0,%width,%offset);
set vx 0;
set vy 0;
charge 0;
select atom within ("%width-%offset",0,%offset,%height);
set vx 0;
set vy 0;
charge 0;
select atom within (0,"%height-%offset",%width,%offset);
set vx 0;
set vy 0;
charge 0;
1000flowif (%flowOn);
set %y %height*rand();
add atom 0 (1.1, %y);
set %lastAtom %number_of_atoms - 1;
set atom[%lastAtom].vx 3000;
set atom[%lastAtom].charge %alphaCharge;
select atom %lastAtom;
traj on;
endif;1.0151000255None228223252255300.012552552552550300.0161.053.025502552550025525500025515.0110600.01.00.0398.011116.0586.00.014.0112599.01.0384.0397.011310.01.02.0falsetrue11Arial11411+292.0216.01112+203.0206.01113+363.0177.01114+321.0285.01115+230.0127.01116+369.0239.01117+334.0109.01118+236.0261.01119+287.0157.031.0false311.0203.0