Математический маятник С++ Windows Form (нарисовать график в панели)

private: System::Void button1_Click(System::Object^  sender, System::EventArgs^  e) {
                   l = System::Convert::ToDouble(textBox1->Text);//длина
                    m = System::Convert::ToDouble(textBox2->Text);//масса
                     k = System::Convert::ToDouble(textBox3->Text);//коэффициент сопротивления среды
                      phi0 = 3.14*(System::Convert::ToDouble(textBox4->Text))/180;// частота вынуждающей силы
                       F0 = System::Convert::ToDouble(textBox5->Text);
                       //   W = System::Convert::ToDouble(textBox6->Text);//амплитуда вынуждающей силы
                 int width = panel1->Width;
                 int height = panel1->Height;
                 x0 = width/2;
                 y0 = height/2;
                
                 R = 5;
                // phi0 = Pi/6;//начальная фаза откуда шарик начинает колебаться
                 g = 10;
             W=1;
                 beta = 0.5*k/m;
                 f0 = -F0/(m*l);
                // W=sqrt(abs(w0*w0-beta*beta));
                 
                 t = 0;
                 //рисунок с шариком рисуется только для частоты W
                 //tetta и B нужны чтобы нарисовать рисунок, они и вычисляются специально
                 //для этой частоты
                 w0 = sqrt(g/l);//собств.частота своб.колебаний                
                 teta = TETTA(W);//фаза колебаний, на которую различаются колебания 
                    //вынудж.силы и собственные
                 B = AmpVS(W);//Амплитуда колебаний вынуждающей силы в уравнении для 
                 //угла отклонения phi
                 A = phi0 - B*cos(teta);//часть амплитуды члена, отвечающего за своб.колебания
                 
                 col = gcnew Color();
                 mypen = gcnew Pen(col->Red);
                 image = panel1->CreateGraphics();
                     mypen1 = gcnew Pen(col->Blue);
                 image1= panel1->CreateGraphics();
                 image->Clear(col->FromArgb(255, 225, 200));
                image1->DrawLine(mypen1, 0, y0, width, y0);
                 x = l*sin(phi(0));//координаты X и Y в полярных координатах
                 y = l*cos(phi(0));
 
                 Draw(x,y-3,255,0,0);
 
                 chart1->Series["Коор-та маятника"]->Points->Clear();
                 chart1->Series["амплитуда своб кол"]->Points->Clear();
                 chart1->Series["амплитуда своб колеб"]->Points->Clear();
            //   chart1->Series["Под действием вын силы"]->Points->Clear();
            //   chart1->Series["не действ вын сила"]->Points->Clear();
                 chart2->Series["Фазовая диаграмма"]->Points->Clear();
                 chart3->Series["Series1"]->Points->Clear();
                 chart4->Series["Series1"]->Points->Clear();
                 //графики амплитуды и фазы 
                 //частота меняется от нуля до выбранного w с шагом 0.001
                 double w=1.;
                 double dW = 0.0;
                 while (dW < w)
                 {
                     chart3->Series["Series1"]->Points->AddXY(dW,abs(AmpVS(dW)));
                     if (TETTA(dW)/Pi > 0) { 
                         chart4->Series["Series1"]->Points->AddXY(dW,TETTA(dW)/Pi-1); }
                     else  chart4->Series["Series1"]->Points->AddXY(dW,TETTA(dW)/Pi);
                     dW+= 0.001;
                 }
                 
 
                 timer1->Start();
             }
private: System::Void timer1_Tick(System::Object^  sender, System::EventArgs^  e) {
              Draw(x,y-3,255,225,200);
                 t+= 0.1;
 
                 x = l*sin(phi(t));
                // y = l*cos(phi(t));
 
                 Draw(x,y-3,255,0,0);
 
                 chart1->Series["Коор-та маятника"]->Points->AddXY(t, x);
            //   chart1->Series["Series4"]->Points->AddXY(t,l*sin(Ampfree(t) + phiforced(t)));
            //  chart1->Series["Под действием вын силы"]->Points->AddXY(t,l*sin(phiforced(t)));
            //  chart1->Series["не действ вын сила"]->Points->AddXY(t,l*sin(phifree(t)));
 
             chart1->Series["амплитуда своб кол"]->Points->AddXY(t,l*sin(Ampfree(t)));
             chart1->Series["амплитуда своб колеб"]->Points->AddXY(t,l*sin(-Ampfree(t)));
 
                 chart2->Series["Фазовая диаграмма"]->Points->AddXY(x,Xderivative(t));
                 
             }
 
            double TETTA(double w)
            {return -atan(beta*w/(sqr(w0) - sqr(w)));}
 
            double AmpVS(double w)
            {return f0/sqrt(sqr(sqr(w0) - sqr(w)) + sqr(beta*w));}
 
            double phi(double s)
            {return phifree(s) + phiforced(s);}
 
            double Ampfree(double s)
            {return A*exp(-beta*s);}
 
            double phifree(double s)
            {return Ampfree(s)*cos(w0*s);}
 
            double phiforced(double s)
            {return B*cos(W*s + teta);}
 
            double sqr(double s)
            {return s*s;}
 
            double Phiderivative(double s)
            {return -beta*phifree(s) - Ampfree(s)*w0*sin(w0*s) -
                    B*W*sin(W*t + teta);}
 
            double Xderivative(double s)
            {return l*cos(phi(s))*Phiderivative(s);}
            //а я ленивая панда! бу-га-га!
            void Draw(double X, double Y,int a1,int a2,int a3)
            { 
              mypen->Color = col->FromArgb(a1, a2, a3);
              image->DrawEllipse(mypen, x0 + (X) - R,
                                 y0 + (Y) - R, 2*R, 2*R);
              
              image->DrawLine(mypen, x0, y0, x0 + (X), y0 + (Y));
            }
    private: System::Void button2_Click(System::Object^  sender, System::EventArgs^  e) {
                 timer1->Stop();
             }