Построить график по численному решению системы диф.ур

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.DataVisualization.Charting;
using System.Windows.Forms;
 
namespace WindowsFormsApplication1
{
    public partial class Form1 : Form
    {
 
        public Form1()
        {
            InitializeComponent();
 
        }
 
        private double F(double a, double q, double t, double x, double y)
        {
            return ((2 * q * Math.Cos(2 * t) - a) * x);
        }
 
        
        
        
        private void button1_Click(object sender, EventArgs e)
        {
            double X1 = double.Parse(textBox1.Text);
            double X2 = double.Parse(textBox2.Text);
            double a = double.Parse(AtextBox.Text);
            double q = double.Parse(QtextBox.Text);
            double h = double.Parse(HtextBox.Text);
            double k1, k2, k3, k4, x0, y0, x1, y1, h1, l1, l2, l3, l4;
            h1 = (X2 - X1) / h;
            x0 = 0;
            y0 = 0;
            double t0 = 0;
 
            for (int i = 0; i < h1; i++)
            {
 
                l1 = y0;
                k1 = F(a, q, t0, x0, y0);
 
                l2 = y0 + k1 * h / 2;
                k2 = F(a, q, t0 + h / 2, x0 + l1 * h / 2, y0 + k1 * h / 2);
 
                l3 = 0 + k2 * h / 2;
                k3 = F(a, q, 0 + h / 2, 0 + l2 * h / 2, 0 + k2 * h / 2);
 
                l4 = 0 + h * l3;
                k4 = F(a, q, 0 + h, 0 + h * l3, 0 + h * k3);
 
                x1 = x0 + (l1 + 2 * l2 + 2 * l3 + l4) * (h / 6);
                y1 = y0 + (k1 + 2 * k2 + 2 * k3 + k4) * (h / 6);
 
                chart1.Series["Series1"].Points.AddXY(x1, y1);
 
                x0 = x1;
                y0 = y1;
 
            }
 
 
        }
 
 
 
        public void Form1_Load(object sender, EventArgs e)
        {
           
        }
    }
}