#include stdafx.h

#include "stdafx.h"

#include <cmath>

#include <iostream>

using namespace std;

#define PI 3.14

int _tmain(int argc, _TCHAR* argv[])

{

double phi1 = 0.00174532925, phi2 = 0.00174532925;

double Amp = 1.5, l(40.0);

int w1(3), w2(3);

double secDer1, secDer2;

double fder1(0.1 * PI/180), fder2(0.1 * PI /180);

double step = 0.05;

for (int j = 1; j <= 189; ++j)

{

double t = j * step;

secDer1 = -(((9.8 / l) + ((Amp * w1 * w1) / l) * cos(w1 * t))*phi1);

fder1 += secDer1 * step;

phi1 += fder1 * step;

if (Amp*w2 > sqrt(2 * 2.8*l))

phi2 = 0;

else

{

secDer2 = -(((9.8 / l) + ((Amp * w2 * w2) / l) * cos(w2 * t))*phi2);

fder2 += secDer2 * step;

phi2 += fder2 * step;

}

cout << t << '\t' << phi1 << '\t' << phi2 <<

endl;

Amp += step;

}

return 0;

}

Дата добавления: 2015-09-29