Решить уравнение при помощи методов Эйлера и Рунге Кутта

[Gladiator1995]

Здравсвуйте, дамы и господа! В институте дали такое задание:
Вот условие задачи
Решить методом Эйлера и Рунге-Кутта уравнение, описывающее смещение x электрона в атоме под действием электрического поля монохроматической световой волны

$$\ddot{x} + \frac{1}{\tau}\dot{x} + {\omega }^{2}x + \gamma {x}^{2} = eE_{0}\cos (pt)$$

\omega - Круговая частота собственных колебаний электрона в атоме
\tau - Их время затухания
\gamma - Коэффициент нелинейности возвращающей силы, действующей на электрон со стороны ядра
e, m - заряд и масса электрона
E_{0} - амплитудное значение напряженности электрического поля световой волны
p - круговая частота волны

Написал программу, касающуюся метода Эйлера, но ее фазовый портрет отвратителен! Помогите найти ошибку

Используется синтаксис C

<span class="dxdy-math dxdy-math-inline" data-math-render="interactive" data-display="0" data-tex-b64="I2lmbmRlZiBGVU5DVElPTgojZGVmaW5lIEZVTkNUSU9OCgpkb3VibGUgRjEoZG91YmxlIEUwLCBkb3VibGUgcCwgZG91YmxlIHQsIGRvdWJsZSB0YXUsIGRvdWJsZSBvbWVnYSwgZG91YmxlIGdhbW1hLCBkb3VibGUgeCwgZG91YmxlIHkpOwoKZG91YmxlIEYyKGRvdWJsZSBFMCwgZG91YmxlIHAsIGRvdWJsZSB0LCBkb3VibGUgdGF1LCBkb3VibGUgb21lZ2EsIGRvdWJsZSBnYW1tYSwgZG91YmxlIHgsIGRvdWJsZSB5KTsKCiNlbmRpZg==" data-fallback-src="https://dxdy.ru/math/77a90a0f3ce59a721411fa13cce0768782.png" data-fallback-hash="77a90a0f3ce59a721411fa13cce07687"><span class="dxdy-math-output"></span><img src='https://dxdy.ru/math/77a90a0f3ce59a721411fa13cce0768782.png' class="latex" alt="#ifndef FUNCTION
#define FUNCTION

double F1(double E0, double p, double t, double tau, double omega, double gamma, double x, double y);

double F2(double E0, double p, double t, double tau, double omega, double gamma, double x, double y);

#endif" title="#ifndef FUNCTION
#define FUNCTION

double F1(double E0, double p, double t, double tau, double omega, double gamma, double x, double y);

double F2(double E0, double p, double t, double tau, double omega, double gamma, double x, double y);

#endif" loading="lazy" decoding="async" /></span>

<span class="dxdy-math dxdy-math-inline" data-math-render="interactive" data-display="0" data-tex-b64="I2luY2x1ZGUgIkZ1bmN0aW9uLmgiCiNpbmNsdWRlIDxzdGRpby5oPgojaW5jbHVkZSA8bWF0aC5oPgoKI2RlZmluZSBxIDEsNkUtMTkKI2RlZmluZSBtIDksMUUtMjcKCmRvdWJsZSBGMShkb3VibGUgRTAsIGRvdWJsZSBwLCBkb3VibGUgdCwgZG91YmxlIHRhdSwgZG91YmxlIG9tZWdhLCBkb3VibGUgZ2FtbWEsIGRvdWJsZSB4LCBkb3VibGUgeSkKewoJLy8gdCAtINCy0YDQtdC80Y8KICAJLy8gRTAgLSDQsNC80L/Qu9C40YLRg9C00L3QvtC1INC30L3QsNGH0LXQvdC40LUg0L3QsNC/0YDRj9C20LXQvdC90L7RgdGC0Lgg0Y3Qu9C10LrRgtGA0LjRh9C10YHQutC+0LPQviDQv9C+0LvRjyDRgdCy0LXRgtC+0LLQvtC5INCy0L7Qu9C90YsKICAJLy8g0YAgLSDQutGA0YPQs9C+0LLQsNGPINGH0LDRgdGC0L7RgtCwINCy0L7Qu9C90YsKCS8vIHRhdSAtINCy0YDQtdC80Y8g0LfQsNGC0YPRhdCw0L3QuNGPINGB0LLQvtCx0L7QtNC90YvRhSDQutC+0LvQtdCx0LDQvdC40LkKCS8vIG9tZWdhIC0g0LrRgNGD0LPQvtCy0LDRjyDRh9Cw0YHRgtC+0YLQsCDRgdCy0L7QsdC+0LTQvdGL0YUg0LrQvtC70LXQsdCw0L3QuNC5INGN0LvQtdC60YLRgNC+0L3QsCDQsiDQsNGC0L7QvNC1CgkvLyBnYW1tYSAtINC60L7RjdGE0YTQuNGG0LjQtdC90YIg0L3QtdC70LjQvdC10LnQvdC+0YHRgtC4INCy0L7Qt9Cy0YDQsNGJ0LDRjtGJ0LXQuSDRgdC40LvRiwoJcmV0dXJuIChxICogRTAgKiBjb3MocCAqIHQpL20gLSB5Lyh0YXUpIC0gKG9tZWdhICogb21lZ2EgKiB4KSAtIGdhbW1hICogeCAqIHgpOyAvLyBGMSh4LCB5LCB0KQp9Cgpkb3VibGUgRjIoZG91YmxlIEUwLCBkb3VibGUgcCwgZG91YmxlIHQsIGRvdWJsZSB0YXUsIGRvdWJsZSBvbWVnYSwgZG91YmxlIGdhbW1hLCBkb3VibGUgeCwgZG91YmxlIHkpCnsKCS8vIHQgLSDQstGA0LXQvNGPCiAgCS8vIEUwIC0g0LDQvNC/0LvQuNGC0YPQtNC90L7QtSDQt9C90LDRh9C10L3QuNC1INC90LDQv9GA0Y/QttC10L3QvdC+0YHRgtC4INGN0LvQtdC60YLRgNC40YfQtdGB0LrQvtCz0L4g0L/QvtC70Y8g0YHQstC10YLQvtCy0L7QuSDQstC+0LvQvdGLCiAgCS8vINGAIC0g0LrRgNGD0LPQvtCy0LDRjyDRh9Cw0YHRgtC+0YLQsCDQstC+0LvQvdGLCgkvLyB0YXUgLSDQstGA0LXQvNGPINC30LDRgtGD0YXQsNC90LjRjyDRgdCy0L7QsdC+0LTQvdGL0YUg0LrQvtC70LXQsdCw0L3QuNC5CgkvLyBvbWVnYSAtINC60YDRg9Cz0L7QstCw0Y8g0YfQsNGB0YLQvtGC0LAg0YHQstC+0LHQvtC00L3Ri9GFINC60L7Qu9C10LHQsNC90LjQuSDRjdC70LXQutGC0YDQvtC90LAg0LIg0LDRgtC+0LzQtQoJLy8gZ2FtbWEgLSDQutC+0Y3RhNGE0LjRhtC40LXQvdGCINC90LXQu9C40L3QtdC50L3QvtGB0YLQuCDQstC+0LfQstGA0LDRidCw0Y7RidC10Lkg0YHQuNC70YsKCXJldHVybiB5OyAvLyBGMih4LCB5LCB0KQp9" data-fallback-src="https://dxdy.ru/math/7b0132310eef7b9bb7d934e13ec3269f82.png" data-fallback-hash="7b0132310eef7b9bb7d934e13ec3269f"><span class="dxdy-math-output"></span><img src='https://dxdy.ru/math/7b0132310eef7b9bb7d934e13ec3269f82.png' class="latex" alt="#include "Function.h"
#include <stdio.h>
#include <math.h>

#define q 1,6E-19
#define m 9,1E-27

double F1(double E0, double p, double t, double tau, double omega, double gamma, double x, double y)
{
        // t - время
        // E0 - амплитудное значение напряженности электрического поля световой волны
        // р - круговая частота волны
        // tau - время затухания свободных колебаний
        // omega - круговая частота свободных колебаний электрона в атоме
        // gamma - коэффициент нелинейности возвращающей силы
        return (q * E0 * cos(p * t)/m - y/(tau) - (omega * omega * x) - gamma * x * x); // F1(x, y, t)
}

double F2(double E0, double p, double t, double tau, double omega, double gamma, double x, double y)
{
        // t - время
        // E0 - амплитудное значение напряженности электрического поля световой волны
        // р - круговая частота волны
        // tau - время затухания свободных колебаний
        // omega - круговая частота свободных колебаний электрона в атоме
        // gamma - коэффициент нелинейности возвращающей силы
        return y; // F2(x, y, t)
}" title="#include "Function.h"
#include <stdio.h>
#include <math.h>

#define q 1,6E-19
#define m 9,1E-27

double F1(double E0, double p, double t, double tau, double omega, double gamma, double x, double y)
{
        // t - время
        // E0 - амплитудное значение напряженности электрического поля световой волны
        // р - круговая частота волны
        // tau - время затухания свободных колебаний
        // omega - круговая частота свободных колебаний электрона в атоме
        // gamma - коэффициент нелинейности возвращающей силы
        return (q * E0 * cos(p * t)/m - y/(tau) - (omega * omega * x) - gamma * x * x); // F1(x, y, t)
}

double F2(double E0, double p, double t, double tau, double omega, double gamma, double x, double y)
{
        // t - время
        // E0 - амплитудное значение напряженности электрического поля световой волны
        // р - круговая частота волны
        // tau - время затухания свободных колебаний
        // omega - круговая частота свободных колебаний электрона в атоме
        // gamma - коэффициент нелинейности возвращающей силы
        return y; // F2(x, y, t)
}" loading="lazy" decoding="async" /></span>

<span class="dxdy-math dxdy-math-inline" data-math-render="interactive" data-display="0" data-tex-b64="I2lmbmRlZiBFVUxFUgojZGVmaW5lIEVVTEVSCiNpbmNsdWRlIDxzdGRpby5oPgoKdm9pZCBFdWxlcihGSUxFICplZiwgZG91YmxlIEUwLCBkb3VibGUgcCwgZG91YmxlIHQsIGRvdWJsZSB0YXUsIGRvdWJsZSBvbWVnYSwgZG91YmxlIGdhbW1hLCBjb25zdCBkb3VibGUgYSwgY29uc3QgZG91YmxlIGIsIGNvbnN0IGRvdWJsZSBlcHMpOwoKI2VuZGlm" data-fallback-src="https://dxdy.ru/math/082bda10334f645d7c1390e23065a2af82.png" data-fallback-hash="082bda10334f645d7c1390e23065a2af"><span class="dxdy-math-output"></span><img src='https://dxdy.ru/math/082bda10334f645d7c1390e23065a2af82.png' class="latex" alt="#ifndef EULER
#define EULER
#include <stdio.h>

void Euler(FILE *ef, double E0, double p, double t, double tau, double omega, double gamma, const double a, const double b, const double eps);

#endif" title="#ifndef EULER
#define EULER
#include <stdio.h>

void Euler(FILE *ef, double E0, double p, double t, double tau, double omega, double gamma, const double a, const double b, const double eps);

#endif" loading="lazy" decoding="async" /></span>

<span class="dxdy-math dxdy-math-inline" data-math-render="interactive" data-display="0" data-tex-b64="I2luY2x1ZGUgIkZ1bmN0aW9uLmgiCiNpbmNsdWRlICJFdWxlci5oIgojaW5jbHVkZSA8bWF0aC5oPgoKdm9pZCBFdWxlcihGSUxFICplZiwgZG91YmxlIEUwLCBkb3VibGUgcCwgZG91YmxlIHQsIGRvdWJsZSB0YXUsIGRvdWJsZSBvbWVnYSwgZG91YmxlIGdhbW1hLCBjb25zdCBkb3VibGUgYSwgY29uc3QgZG91YmxlIGIsIGNvbnN0IGRvdWJsZSBlcHMpCnsKCWRvdWJsZSB4LCB4MiwgeSwgeTIsIGR0OwoJeCA9IGE7Cgl5ID0gYjsKCWR0ID0gMC4wMTsKCWZwcmludGYoZWYsICJJbml0aWFsIGRhdGEgeCgwKSA9IGEgPSAlZjsgeSgwKSA9IGIgPSAlZjsgRXBzaWxvbiA9IGVwcyA9ICVmO1xuIiwgeCwgeSwgZXBzKTsKCWZwcmludGYoZWYsICJFMCA9ICVmOyIsIEUwKTsKCWZwcmludGYoZWYsICJcbiIpOwoJZnByaW50ZihlZiwgInAgPSAlZjsiLCBwKTsKCWZwcmludGYoZWYsICJcbiIpOwoJZnByaW50ZihlZiwgInRhdSA9ICVmOyIsIHRhdSk7CglmcHJpbnRmKGVmLCAiXG4iKTsKCWZwcmludGYoZWYsICJvbWVnYSA9ICVmOyIsIG9tZWdhKTsKCWZwcmludGYoZWYsICJcbiIpOwoJZnByaW50ZihlZiwgImdhbW1hID0gJWY7IiwgZ2FtbWEpOwoJZm9yKGludCBpID0gMDtpIDwgMTAwMDA7ICsraSkKCQl7CQoJCQl3aGlsZSAoZmFicygoRjEoRTAsIHAsIHQsIHRhdSwgb21lZ2EsIGdhbW1hLCB4LCB5KSAqIGR0IC8gMikpID4gZXBzIHx8IGZhYnMoKEYyKEUwLCBwLCB0LCB0YXUsIG9tZWdhLCBnYW1tYSwgeCwgeSkgKiBkdCAvIDIpKSA+IGVwcykJZHQvPTIuMDsKCQkJeDIgPSB4ICsgRjEoRTAsIHAsIHQsIHRhdSwgb21lZ2EsIGdhbW1hLCB4LCB5KSAqIGR0OyAvLyA/Pz8/Pz8/PyB4ID8geSA/ID8/Pz8gPz8/Pz8/ID8/Pz8/PyA/Pz8/ID8/ID8/IC4uLgoJCQl5MiA9IHkgKyBGMihFMCwgcCwgdCwgdGF1LCBvbWVnYSwgZ2FtbWEsIHgsIHkpICogZHQ7IC8vIC4uLiA/Pz8gPyA/ID8/Pz8gPz8/Pz8/CgkJCXggPSB4MjsKCQkJeSA9IHkyOwoJCQl0ICs9IGR0OwoJCQlmcHJpbnRmKGVmLCAiJWYgICAlZiAgICVmXG4iLCB0ICwgeDIsIHkyKTsKCQkJZHQgPSAwLjAxOwoJCX0KfQ==" data-fallback-src="https://dxdy.ru/math/e72121ebdb2676ca5b8fa37ab1fcff8782.png" data-fallback-hash="e72121ebdb2676ca5b8fa37ab1fcff87"><span class="dxdy-math-output"></span><img src='https://dxdy.ru/math/e72121ebdb2676ca5b8fa37ab1fcff8782.png' class="latex" alt="#include "Function.h"
#include "Euler.h"
#include <math.h>

void Euler(FILE *ef, double E0, double p, double t, double tau, double omega, double gamma, const double a, const double b, const double eps)
{
        double x, x2, y, y2, dt;
        x = a;
        y = b;
        dt = 0.01;
        fprintf(ef, "Initial data x(0) = a = %f; y(0) = b = %f; Epsilon = eps = %f;\n", x, y, eps);
        fprintf(ef, "E0 = %f;", E0);
        fprintf(ef, "\n");
        fprintf(ef, "p = %f;", p);
        fprintf(ef, "\n");
        fprintf(ef, "tau = %f;", tau);
        fprintf(ef, "\n");
        fprintf(ef, "omega = %f;", omega);
        fprintf(ef, "\n");
        fprintf(ef, "gamma = %f;", gamma);
        for(int i = 0;i < 10000; ++i)
                {       
                        while (fabs((F1(E0, p, t, tau, omega, gamma, x, y) * dt / 2)) > eps || fabs((F2(E0, p, t, tau, omega, gamma, x, y) * dt / 2)) > eps)    dt/=2.0;
                        x2 = x + F1(E0, p, t, tau, omega, gamma, x, y) * dt; // ???????? x ? y ? ???? ?????? ?????? ???? ?? ?? ...
                        y2 = y + F2(E0, p, t, tau, omega, gamma, x, y) * dt; // ... ??? ? ? ???? ??????
                        x = x2;
                        y = y2;
                        t += dt;
                        fprintf(ef, "%f   %f   %f\n", t , x2, y2);
                        dt = 0.01;
                }
}" title="#include "Function.h"
#include "Euler.h"
#include <math.h>

void Euler(FILE *ef, double E0, double p, double t, double tau, double omega, double gamma, const double a, const double b, const double eps)
{
        double x, x2, y, y2, dt;
        x = a;
        y = b;
        dt = 0.01;
        fprintf(ef, "Initial data x(0) = a = %f; y(0) = b = %f; Epsilon = eps = %f;\n", x, y, eps);
        fprintf(ef, "E0 = %f;", E0);
        fprintf(ef, "\n");
        fprintf(ef, "p = %f;", p);
        fprintf(ef, "\n");
        fprintf(ef, "tau = %f;", tau);
        fprintf(ef, "\n");
        fprintf(ef, "omega = %f;", omega);
        fprintf(ef, "\n");
        fprintf(ef, "gamma = %f;", gamma);
        for(int i = 0;i < 10000; ++i)
                {       
                        while (fabs((F1(E0, p, t, tau, omega, gamma, x, y) * dt / 2)) > eps || fabs((F2(E0, p, t, tau, omega, gamma, x, y) * dt / 2)) > eps)    dt/=2.0;
                        x2 = x + F1(E0, p, t, tau, omega, gamma, x, y) * dt; // ???????? x ? y ? ???? ?????? ?????? ???? ?? ?? ...
                        y2 = y + F2(E0, p, t, tau, omega, gamma, x, y) * dt; // ... ??? ? ? ???? ??????
                        x = x2;
                        y = y2;
                        t += dt;
                        fprintf(ef, "%f   %f   %f\n", t , x2, y2);
                        dt = 0.01;
                }
}" loading="lazy" decoding="async" /></span>

<span class="dxdy-math dxdy-math-inline" data-math-render="interactive" data-display="0" data-tex-b64="I2luY2x1ZGUgIkV1bGVyLmgiCiNpbmNsdWRlIDxzdGRpby5oPgoKI2RlZmluZSBxIDEsNkUtMTkKI2RlZmluZSBtIDksMUUgLTI3CgppbnQgbWFpbigpCnsKCWRvdWJsZSB4LCB5LCBFMCwgcCwgdCwgdGF1LCBvbWVnYSwgZ2FtbWEsIGVwczsKCS8vIHQgLSDDosOww6XDrMO/CiAgCS8vIEUwIC0gw6DDrMOvw6vDqMOyw7PDpMOtw67DpSDDp8Otw6DDt8Olw63DqMOlIMOtw6DDr8Oww7/DpsOlw63DrcOuw7HDssOoIMO9w6vDpcOqw7LDsMOow7fDpcOxw6rDrsOjw64gw6/DrsOrw78gw7HDosOlw7LDrsOiw67DqSDDosOuw6vDrcO7CiAgCS8vIMOwIC0gw6rDsMOzw6PDrsOiw6DDvyDDt8Ogw7HDssOuw7LDoCDDosOuw6vDrcO7CgkvLyB0YXUgLSDDosOww6XDrMO/IMOnw6DDssOzw7XDoMOtw6jDvyDDscOiw67DocOuw6TDrcO7w7Ugw6rDrsOrw6XDocOgw63DqMOpCgkvLyBvbWVnYSAtIMOqw7DDs8Ojw67DosOgw78gw7fDoMOxw7LDrsOyw6Agw7HDosOuw6HDrsOkw63Du8O1IMOqw67Dq8Olw6HDoMOtw6jDqSDDvcOrw6XDqsOyw7DDrsOtw6Agw6Igw6DDssOuw6zDpQoJLy8gZ2FtbWEgLSDDqsOuw73DtMO0w6jDtsOow6XDrcOyIMOtw6XDq8Oow63DpcOpw63DrsOxw7LDqCDDosOuw6fDosOww6DDucOgw77DucOlw6kgw7HDqMOrw7sKCS8vIGVwc2lsb24gPSBlcHMgLSAKCXByaW50ZigiRW50ZXIgZXBzaWxvblxuIik7CglzY2FuZigiJWxmIiwgJmVwcyk7CglwcmludGYoIkVudGVyIHgsIHlcbiIpOwoJc2NhbmYoIiVsZiAlbGYiLCAmeCwgJnkpOwoJcHJpbnRmKCJFbnRlciBFMDpcbiIpOwoJc2NhbmYoIiVsZiIsICZFMCk7CglwcmludGYoIkVudGVyIHA6XG4iKTsKCXNjYW5mKCIlbGYiLCAmcCk7CglwcmludGYoIkVudGVyIHQ6XG4iKTsKCXNjYW5mKCIlbGYiLCAmdCk7CglwcmludGYoIkVudGVyIHRhdTpcbiIpOwoJc2NhbmYoIiVsZiIsICZ0YXUpOwoJcHJpbnRmKCJFbnRlciBvbWVnYTpcbiIpOwoJc2NhbmYoIiVsZiIsICZvbWVnYSk7CglwcmludGYoIkVudGVyIGdhbW1hOlxuIik7CglzY2FuZigiJWxmIiwgJmdhbW1hKTsKCXByaW50ZigiXG4iKTsKCUZJTEUgKmVmOwoJZWYgPSBmb3BlbigiRXVsZXIudHh0IiwgInciKTsKCUV1bGVyKGVmLCBFMCwgcCwgdCwgdGF1LCBvbWVnYSwgZ2FtbWEsIHgsIHksIGVwcyk7CglmY2xvc2UoZWYpOwoJcmV0dXJuIDA7Cn0=" data-fallback-src="https://dxdy.ru/math/d41d8cd98f00b204e9800998ecf8427e82.png" data-fallback-hash="d41d8cd98f00b204e9800998ecf8427e"><span class="dxdy-math-output"></span><img src='https://dxdy.ru/math/d41d8cd98f00b204e9800998ecf8427e82.png' class="latex" alt="#include "Euler.h"
#include <stdio.h>

#define q 1,6E-19
#define m 9,1E -27

int main()
{
        double x, y, E0, p, t, tau, omega, gamma, eps;
        // t - âðåìÿ
        // E0 - àìïëèòóäíîå çíà÷åíèå íàïðÿæåííîñòè ýëåêòðè÷åñêîãî ïîëÿ ñâåòîâîé âîëíû
        // ð - êðóãîâàÿ ÷àñòîòà âîëíû
        // tau - âðåìÿ çàòóõàíèÿ ñâîáîäíûõ êîëåáàíèé
        // omega - êðóãîâàÿ ÷àñòîòà ñâîáîäíûõ êîëåáàíèé ýëåêòðîíà â àòîìå
        // gamma - êîýôôèöèåíò íåëèíåéíîñòè âîçâðàùàþùåé ñèëû
        // epsilon = eps - 
        printf("Enter epsilon\n");
        scanf("%lf", &eps);
        printf("Enter x, y\n");
        scanf("%lf %lf", &x, &y);
        printf("Enter E0:\n");
        scanf("%lf", &E0);
        printf("Enter p:\n");
        scanf("%lf", &p);
        printf("Enter t:\n");
        scanf("%lf", &t);
        printf("Enter tau:\n");
        scanf("%lf", &tau);
        printf("Enter omega:\n");
        scanf("%lf", &omega);
        printf("Enter gamma:\n");
        scanf("%lf", &gamma);
        printf("\n");
        FILE *ef;
        ef = fopen("Euler.txt", "w");
        Euler(ef, E0, p, t, tau, omega, gamma, x, y, eps);
        fclose(ef);
        return 0;
}" title="#include "Euler.h"
#include <stdio.h>

#define q 1,6E-19
#define m 9,1E -27

int main()
{
        double x, y, E0, p, t, tau, omega, gamma, eps;
        // t - âðåìÿ
        // E0 - àìïëèòóäíîå çíà÷åíèå íàïðÿæåííîñòè ýëåêòðè÷åñêîãî ïîëÿ ñâåòîâîé âîëíû
        // ð - êðóãîâàÿ ÷àñòîòà âîëíû
        // tau - âðåìÿ çàòóõàíèÿ ñâîáîäíûõ êîëåáàíèé
        // omega - êðóãîâàÿ ÷àñòîòà ñâîáîäíûõ êîëåáàíèé ýëåêòðîíà â àòîìå
        // gamma - êîýôôèöèåíò íåëèíåéíîñòè âîçâðàùàþùåé ñèëû
        // epsilon = eps - 
        printf("Enter epsilon\n");
        scanf("%lf", &eps);
        printf("Enter x, y\n");
        scanf("%lf %lf", &x, &y);
        printf("Enter E0:\n");
        scanf("%lf", &E0);
        printf("Enter p:\n");
        scanf("%lf", &p);
        printf("Enter t:\n");
        scanf("%lf", &t);
        printf("Enter tau:\n");
        scanf("%lf", &tau);
        printf("Enter omega:\n");
        scanf("%lf", &omega);
        printf("Enter gamma:\n");
        scanf("%lf", &gamma);
        printf("\n");
        FILE *ef;
        ef = fopen("Euler.txt", "w");
        Euler(ef, E0, p, t, tau, omega, gamma, x, y, eps);
        fclose(ef);
        return 0;
}" loading="lazy" decoding="async" /></span>

[Pphantom]

Тема перемещена из форума «Программирование» в форум «Карантин»
по следующим причинам:

- неправильно набраны формулы (краткие инструкции: «Краткий FAQ по тегу [math]» и видеоролик Как записывать формулы);
- неправильно оформлен код.

Исправьте все Ваши ошибки и сообщите об этом в теме Сообщение в карантине исправлено.
Настоятельно рекомендуется ознакомиться с темами Что такое карантин и что нужно делать, чтобы там оказаться и Правила научного форума.