Программа, рассчитывающая задачу по молекулярной физике

using System;
 
class Program
{
    static void Main()
    {
        Console.Write("Введите площадь: ");
        double s = double.Parse(Console.ReadLine());
        Console.Write("Введите расстояние: ");
        double d = double.Parse(Console.ReadLine());
        double c = s / (4 * Math.PI * d);
        Console.Write("Ёмкость конденсатора: {0}", c);
        Console.ReadKey();
    }
}

using System;
using System.Collections.Generic;
using System.Drawing;
using System.Windows.Forms;
 
namespace WindowsFormsApplication231
{
    public partial class Form1 : Form
    {
        private Solver solver;
        public float Scale = 10000f;//шкала (пикселей на метр)
        public float dT = 0.0001f;//шаг моделирования (секунд на шаг)
 
        public Form1()
        {
            InitializeComponent();
 
            SetStyle(ControlStyles.AllPaintingInWmPaint | ControlStyles.OptimizedDoubleBuffer | ControlStyles.UserPaint | ControlStyles.ResizeRedraw, true);
            WindowState = FormWindowState.Maximized;
 
            solver = new Solver();
 
            var bt = new Button() {Parent = this, Text = "Start", Left = 100};
            bt.Click += (o, e) => Start();
 
            Application.Idle += new EventHandler(Application_Idle);
        }
 
        private void Start()
        {
            solver.Capacitor = new Capacitor(0.003, 0.03);//задаем площадь и расст между пластин конденсатора
            solver.E = 0.001;//задаем напряженность эл поля
 
            solver.Particles.Clear();
            solver.Particles.Add(new Particle(1, 0.001, 5));//создаем частицу
            solver.Particles.Add(new Particle(3, 0.002, 3));//создаем частицу
            solver.Particles.Add(new Particle(1.5, 0.0005, -1));//создаем частицу
        }
 
        void Application_Idle(object sender, EventArgs e)
        {
            solver.DoStep(dT);
            Invalidate();
        }
 
        protected override void OnPaint(PaintEventArgs e)
        {
            if (solver == null)
                return;
 
            var center = new Point(ClientRectangle.Left + 20, (ClientRectangle.Top + ClientRectangle.Bottom)/2);
            var gr = e.Graphics;
            gr.SmoothingMode = System.Drawing.Drawing2D.SmoothingMode.HighQuality;
 
            using(var pen = new Pen(Color.Gray))
            {
                pen.EndCap = System.Drawing.Drawing2D.LineCap.ArrowAnchor;
                gr.DrawLine(pen, ClientRectangle.Left, center.Y, ClientRectangle.Right, center.Y);
                gr.DrawLine(pen, center.X, ClientRectangle.Bottom, center.X, ClientRectangle.Top);
            }
 
            if (solver.Particles.Count == 0)
                return;
 
            using (var pen = new Pen(Color.Blue, 3))
                gr.DrawLine(pen, (float)center.X, (float)(center.Y - Scale * solver.Capacitor.D / 2), (float)(center.X + Scale * solver.Capacitor.Width), (float)(center.Y - Scale * solver.Capacitor.D / 2));
 
            using (var pen = new Pen(Color.Red, 3))
                gr.DrawLine(pen, (float)center.X, (float)(center.Y + Scale * solver.Capacitor.D / 2), (float)(center.X + Scale * solver.Capacitor.Width), (float)(center.Y + Scale * solver.Capacitor.D / 2));
 
            foreach (var p in solver.Particles)
            {
                var x = (float) (center.X + Scale*p.X - 1);
                var y = (float) (center.Y - Scale*p.Y - 1);
                gr.FillEllipse(p.Q > 0 ? Brushes.Magenta : Brushes.Blue, x, y, 3, 3);
                if (p.Stop)
                {
                    gr.DrawString(string.Format("tg(a) = {0:0.0000} (теор. {1:0.0000})", p.Vy / p.Vx, solver.GetTheorE(p)), Font, Brushes.Red, x + 5, y + 5);
                }
            }
 
        }
    }
 
    /// Решатель
    class Solver
    {
        public List<Particle> Particles = new List<Particle>();//частицы
        public Capacitor Capacitor;//конденсатор
        public double E;//Напряженность поля
        public const double g = 9.8;//ускорение свободного падения
 
        //Щаг моделирования
        public void DoStep(double dt)
        {
            foreach(var p in Particles)
            if(!p.Stop)
            {
                var fY = p.Q * E;
                var aY = fY / p.M;//без учета силы тяжести, если нужно учитывать - нужно отнять g
                p.Vy += aY * dt;
                p.X += p.Vx * dt;
                p.Y += p.Vy * dt;
 
                if (p.X > Capacitor.Width)
                    p.Stop = true;
            }
        }
 
        /// Теоретический угол выхода частицы из конденсатора
        public double GetTheorE(Particle p)
        {
            return (p.Q/p.M)*Capacitor.Width*E/(p.Vx0*p.Vx0);
        }
    }
 
    /// <summary>
    /// Частица
    /// </summary>
    class Particle
    {
        public double X;//координата X (м)
        public double Y;//координата Y (м)
        public double Vx0;//начальная скорость  по X (м/с)
        public double Vx;//текущая скорость по X (м/с)
        public double Vy;//текущая скорость по Y (м/с)
        public double M;//масса (кг)
        public double Q;//заряд (кулон)
 
        public bool Stop = false;
 
        public Particle(double vx0, double m, double q)
        {
            this.X = 0;
            this.Y = 0;
            this.Vx0 = this.Vx = vx0;
            this.Vy = 0;
            this.M = m;
            this.Q = q;
        }
    }
 
    /// <summary>
    /// Конденсатор
    /// </summary>
    class Capacitor
    {
        public double Square;//площадь (м^2)
        public double D;//расст между пластинами (м)
 
        // Ширина пластин, считая, что конденсатор квадратный
        public double Width
        {
            get { return Math.Sqrt(Square); }
        }
 
        public Capacitor(double square, double d)
        {
            this.Square = square;
            this.D = d;
        }
    }
}