Нарисовать молекулу ДНК

using System;
using System.Collections.Generic;
using System.Drawing;
using System.Drawing.Drawing2D;
using System.Windows.Forms;
using System.Linq;
 
namespace WindowsFormsApplication324
{
    public partial class Form1 : Form
    {
        private DNA dna;
 
        public Form1()
        {
            InitializeComponent();
 
            SetStyle(ControlStyles.AllPaintingInWmPaint | ControlStyles.OptimizedDoubleBuffer | ControlStyles.UserPaint | ControlStyles.ResizeRedraw, true);
 
            //new TrackBar { Parent = this }.ValueChanged += (o, e) => { dna.Angle = (o as TrackBar).Value / 10f; Invalidate(); };
            new Timer { Interval = 150, Enabled = true }.Tick += delegate { dna.Angle += 0.03f; Invalidate(); };
 
            dna = new DNA();
            dna.Build();
        }
 
        protected override void OnPaint(PaintEventArgs e)
        {
            e.Graphics.SmoothingMode = SmoothingMode.HighQuality;
            new DNADrawer().Draw(e.Graphics, dna);
        }
    }
 
    class Ball
    {
        public Color Color { get; set; }
        public float Radius { get; set; }
    }
 
    /// <summary>
    /// Molecular chain item
    /// </summary>
    class DNAItem
    {
        /// <summary>
        /// Coordinate along molecular axis
        /// </summary>
        public float Index { get; set; }
        /// <summary>
        /// Rotation angle around molecular axis
        /// </summary>
        public float Angle { get; set; }
        
        public Ball Ball1 { get; set; }
        public Ball Ball2 { get; set; }
 
        public DNAItem()
        {
            Ball1 = new Ball {Color = Color.Red};
            Ball2 = new Ball { Color = Color.Blue };
        }
    }
 
    /// <summary>
    /// DNA molecule
    /// </summary>
    class DNA
    {
        public PointF StartPoint { get; set; }
        /// <summary>
        /// Direction and scale of molecular axis
        /// </summary>
        public PointF Direction { get; set; }
        /// <summary>
        /// Rotation angle around molecular axis
        /// </summary>
        public float Angle { get; set; }
        /// <summary>
        /// Radius 
        /// </summary>
        public float Radius { get; set; }
        /// <summary>
        /// Ball radius
        /// </summary>
        public float BallRadius { get; set; }
        /// <summary>
        /// Chain items
        /// </summary>
        public List<DNAItem> Items { get; private set; }
 
        public DNA()
        {
            Items = new List<DNAItem>();
            Direction = new PointF(14, -14);
            StartPoint = new PointF(50, 500);
            Radius = 50;
            BallRadius = 10;
        }
 
        public void Build(int itemCount = 50, float angleStep = 0.5f)
        {
            Items.Clear();
            var a = 0f;
 
            for(int i=0; i<itemCount; i++)
            {
                Items.Add(new DNAItem(){Angle = a, Index = i});
                a += angleStep;
            }
        }
    }
 
    /// <summary>
    /// Drawer
    /// </summary>
    class DNADrawer
    {
        public void Draw(Graphics gr, DNA dna)
        {
            gr.TranslateTransform(dna.StartPoint.X, dna.StartPoint.Y);
            gr.RotateTransform((float)(180 * Math.Atan2(dna.Direction.Y, dna.Direction.X) / Math.PI));
            
            var scale = (float)Math.Sqrt(dna.Direction.X*dna.Direction.X + dna.Direction.Y*dna.Direction.Y);
 
            foreach(var item in dna.Items)
            {
                Draw(gr, dna, item, scale);
            }
        }
 
        private void Draw(Graphics gr, DNA dna, DNAItem item, float scale)
        {
            var a = item.Angle + dna.Angle;
            var z = (float)Math.Sin(a);
            var y = (float)Math.Cos(a);
            var x = scale*item.Index;
            var dr = 0.7f * z;
 
            if (z > 0)
            {
                Draw(gr, item.Ball1, new PointF(x, y - dna.Radius * y), dna.BallRadius - dr);
                Draw(gr, item.Ball2, new PointF(x, y + dna.Radius * y), dna.BallRadius + dr);
            }else
            {
                Draw(gr, item.Ball2, new PointF(x, y + dna.Radius * y), dna.BallRadius + dr);
                Draw(gr, item.Ball1, new PointF(x, y - dna.Radius * y), dna.BallRadius - dr);
            }
        }
 
        private void Draw(Graphics gr, Ball ball, PointF p, float r)
        {
            using (var brush = new LinearGradientBrush(new PointF(p.X - r, p.Y - r), new PointF(p.X + 2* r, p.Y + 2* r), ball.Color, Color.White))
            using (var penBrush = new LinearGradientBrush(new PointF(p.X - 3, 0), new PointF(p.X + 2, 0), ball.Color, Color.White))
            using (var pen = new Pen(penBrush, 3))
            {
                gr.DrawLine(pen, p.X, p.Y, p.X, 0);
                gr.FillEllipse(brush, p.X - r, p.Y - r, 2*r, 2*r);
            }
        }
    }
}