using System;

namespace Science.Physics.GeneralPhysics
{
	/// <summary>
	/// Nucleus
	/// </summary>
	public class Nucleus
	{
		public Nucleus()
		{
		}
		private double C1 = 15.7, C2 = 17.8, C3 = 0.71, C4 = 23.6;
		private int Z, N, A;
		private double mass, u;
		public int AtomicNumber
		{
			get{return Z;}
			set{Z = value;}
		}
		public int NeutronNumber
		{
			get{return N;}
			set{N = value;}
		}
		public int MassNumber
		{
			get{return A;}
			set{A = value;}
		}
		public double Radius
		{
			get{return 1.2E-15*Math.Pow(A,1.0/3.0);}
		}
		public double BindingEnergyMeV
		{
			get
			{
				return C1*A 
					 - C2*Math.Pow(A,2.0/3.0)
    				 - C3*Z*(Z-1)/Math.Pow(A,1.0/3.0)
	    			 - C4*(N-Z)*(N-Z)/A;
			}
		}
		public double MassMeV
		{
			get{return mass;}
			set{mass = value;}
		}
		public double MassINu
		{
			get{return u;}
			set{u = value;}
		}
	}
}