using System;
using GP = Science.Physics.GeneralPhysics;

namespace ScienceTest.PhysicsTest.GeneralPhysicsTest
{
	/// <summary>
	/// ForceTest
	/// </summary>
	public class ForceTest
	{
		public ForceTest()
		{
		}
		private string result;
		public string Result
		{
			get{return result;}
		}

		public void Compute()
		{
            GP.Scalar.FunctionOfPosition func = new GP.Scalar.FunctionOfPosition(xtoU);
			GP.PotentialEnergy U = new GP.PotentialEnergy();
            U.ScalarFunctionOfPosition = func;
            
			GP.Position r = new GP.Position();
			r.X = 1.0;
			r.Y = 3.0;
			r.Z = 2.0;
		
			GP.Force f = new GP.Force(U, r);
			
			result += f.ToString();	
		}
		private GP.Scalar xtoU(GP.Position x)
		{
            GP.Scalar s = new GP.Scalar();
			s.Magnitude = 0.5*10.0*x.X*x.X+13.0*x.Y*x.Y-12.0*x.Z*x.Z;
            return s;
		}	
	}
}
// -10.000000000006 +/- 0 i    -77.9999999999844 +/- 0 j     +48.0000000000018 +/- 0 k    (N)