BL-Physics | Python Fiddle

def position2Mass3Springs (mass1InitialPosition, mass2InitialPosition, mass1, mass2, springConstant1, springConstant2, springConstant3, time):
	DataSet2Mass3Springs = open("DataSet2Mass3Springs.txt", "w")
	mass1Equilibrium = 10
	mass2Equilibrium = 20
	mass1Position = mass1Equilibrium + mass1InitialPosition
	mass2Position = mass2Equilibrium + mass2InitialPosition
	mass1Velocity = 0
	mass2Velocity = 0
	timeStamp = 0
	while timeStamp <= time:
		if timeStamp >= 0:
			mass1NetForce = -1 * springConstant1 * (mass1Position-mass1Equilibrium) + -1 * springConstant2 * (10-(mass2Position-mass1Position))
			mass2NetForce = -1 * springConstant3 * (mass2Position-mass2Equilibrium) + -1 * springConstant2 * ((mass2Position-mass1Position)-10)
		mass1Acceleration = mass1NetForce/mass1
		mass2Acceleration = mass2NetForce/mass2
		mass1Position = mass1Position + mass1Velocity * 0.001 + 0.5 * mass1Acceleration * 0.001 * 0.001
		mass2Position = mass2Position + mass2Velocity * 0.001 + 0.5 * mass2Acceleration * 0.001 * 0.001
		mass1Velocity = mass1Velocity + mass1Acceleration * 0.001
		mass2Velocity = mass2Velocity + mass2Acceleration * 0.001
		mass1PotentialEnergy = 0.5 * springConstant1 * (mass1Position-mass1Equilibrium) * (mass1Position-mass1Equilibrium) + 0.5 * springConstant2 * (10-(mass2Position-mass1Position)) * (10-(mass2Position-mass1Position))
		mass2PotentialEnergy = 0.5 * springConstant3 * (mass2Position-mass2Equilibrium) * (mass2Position-mass2Equilibrium) + 0.5 * springConstant2 * ((mass2Position-mass1Position)-10) * ((mass2Position-mass1Position)-10)
		mass1KineticEnergy = 0.5 * mass1 * mass1Velocity * mass1Velocity
		mass2KineticEnergy = 0.5 * mass2 * mass2Velocity * mass2Velocity
		mass1TotalEnergy = mass1PotentialEnergy + mass1KineticEnergy
		mass2TotalEnergy = mass2PotentialEnergy + mass2KineticEnergy
		systemTotalEnergy = mass1TotalEnergy + mass2TotalEnergy
		timeStamp = timeStamp + 0.001
		DataSet2Mass3Springs.write(str(timeStamp) + " " + str(mass1Position) + " " + str(mass2Position) + " " + str(mass1KineticEnergy) + " " + str(mass1PotentialEnergy) + " " + str(mass2KineticEnergy) + " " + str(mass2PotentialEnergy)+ " " + str(systemTotalEnergy) + "\n")
position2Mass3Springs(-3,0,50,50,100,100,100,10)

Python Fiddle

Python Cloud IDE