Решение дифференциальных уравнений методом Эйлера

#include <iostream>
#include <cmath>
#include <vector>
using namespace std;
double f1(double w, double x, double y) 
{
    return -0.804*((cos(w) + sin(w - y)) / (cos(x) + sin(x - y)));
}
double f2(double w, double x, double y) {
    return -(1.51*cos(w) + 1.88*cos(x)*f1(w,x,y)) / cos(y);
}
double f3(double w, double x, double y, double z)
{
    return -0.655*f2(w, x, y) / sin(z);
}
 
int main()
 
{
 
    double  h, b = 1.68;
    double pi = asin(1) * 2;
    
 
    int i, n = 24, a = 0;
    double t[24], w[24], x[24], y[24], z[24];
 
    t[0] = a;
    w[0] = 0;
    x[0] = 2.97;
    y[0] = 4.17;
    z[0] = 3.97;
 
    
 
    printf("t0=%2.2f\w0=%2.2f\x0=%2.2f\y0=%2.2f\z0=%2.2f\t\n", t[0], x[0], y[0], z[0]);
 
    h = (b - a) / n;
 
    cout << "\nh=" << h << endl;
 
    printf("\ni\t t\t w\t x\t y\t z\n");
    
    for (i = 0; i <= n; i++)
 
    {
 
        t[i + 1] = t[i] + h;
        w[i + 1] = w[i] + 2 * pi / 1.68;
        x[i + 1] = x[i] + h * f1(w, x, y);
        y[i + 1] = y[i] + h * f2(w, x, y);
        z[i + 1] = z[i] + h * f3(w, x, y, z);
 
 
 
        printf("%1d\t%8.7f\t%8.7f\t%8.7f\t%8.7f\t%8.7f\n", i, t[i], w[i], x[i], y[i], z[i]);
 
    }
 
    char st;
 
    cin >> st;
 
}

#include <iostream>
#include <cmath>
#include <vector>
using namespace std;
double f1(double w, double x, double y)
{
    return -0.804*((cos(w) + sin(w - y)) / (cos(x) + sin(x - y)));
}
double f2(double w, double x, double y) {
    return -(1.51*cos(w) + 1.88*cos(x)*f1(w, x, y)) / cos(y);
}
double f3(double w, double x, double y, double z)
{
    return -0.655*f2(w, x, y) / sin(z);
}
 
int main()
 
{
 
    double h, b = 1.68;
    double pi = asin(1) * 2;
 
 
    int i, n = 24, a = 0;
    double t[24], w[24], x[24], y[24], z[24];
 
    t[0] = a;
    w[0] = 0;
    x[0] = 2.97;
    y[0] = 4.17;
    z[0] = 3.97;
 
 
 
    printf("t0= %2.2f \tw0= %2.2f \tx0= %2.2f \ty0= %2.2f \tz0= %2.2f \t\n", t[0], w[0], x[0], y[0], z[0]);
 
    h = (b - a) / n;
 
    cout << "\nh=" << h << endl;
 
    printf("\ni\t t\t w\t x\t y\t z\n");
 
    for (i = 0; i <= n; i++)
 
    {
 
        t[i + 1] = t[i] + h;
        w[i + 1] = w[i] + 2 * pi / 1.68;
        x[i + 1] = x[i] + h * f1(w[i], x[i], y[i]);
        y[i + 1] = y[i] + h * f2(w[i], x[i], y[i]);
        z[i + 1] = z[i] + h * f3(w[i], x[i], y[i], z[i]);
 
 
 
        printf("%1d\t%8.7f\t%8.7f\t%8.7f\t%8.7f\t%8.7f\n", i, t[i], w[i], x[i], y[i], z[i]);
 
    }
 
    char st;
 
    cin >> st;
 
}