Помогите с программой

[anna.menshutina]

Условие:Разработать модель движения практически невесомой заряженной частицы в электрическом поле, созданном системой нескольких фиксированных в пространстве заряженных тел, в случае, когда заряженные тела находятся в одной плоскости и в ней же находится движущаяся частица.


Есть весь код программы, НО нужно добавить 3 фиксированные заряженные частицы на графике( т.е добавть жирную точку или закрашенный круг на график) + заряженную частицу, которую позьзаватель может изменять, чтобы траектория этой не фиксированной частицы изменялась и это было видно на графике.

НА ЯЗЫКЕ C#
И ЭТО ООП

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;

namespace WindowsFormsApplication1
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
            Calc();
            chart1.DataBind();
        }

        public void Calc()
        {
            double[] R = new double[120];
            double[] F = new double[120];
            double[] V = new double[120];
            int j;
            double p, e0, m, q1, q2, q3, q4;
            p = 3.1415;
            e0 = 8.85 * Math.Pow(10, -12);
            m = Math.Pow(9, -10);
            q1 = Math.Pow(9, -9);
            q2 = (2 * Math.Pow(9, -9));
            q3 = (2 * Math.Pow(9, -9));
            q4 = (2 * Math.Pow(9, -9));
            double i;
            for (i = 0.2, j = 0; i <= 24; i += 0.2, j++)
                R[j] = i;
            for (int k = 0; k < 120; k++)
                F[k] = (q1 / (4 * p * e0 * Math.Pow(R[k], 2))) * (q3 + q2 * 0.958 + q4 * 0.958);
            for (int k = 0; k < 120; k++)
                V[k] = Math.Sqrt(F[k] / m * 2 * R[k]);
            chart1.Series["Series1"].ChartType = System.Windows.Forms.DataVisualization.Charting.SeriesChartType.Line;
            chart1.Series["Series1"].Points.DataBindXY(R, V);
        }

        private void chart1_Click(object sender, EventArgs e)
        {

        }
    }
}