Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include "pch.h"
- #include <iostream>
- #include "math.h"
- using namespace std;
- double psi(double x)
- {
- double s;
- s = (1 - x - +pow(x, 2)) / sqrt(1 + (2 * pow(x, 2)));
- return s;
- }
- double fi(double x)
- {
- double s;
- s = (5 * pow(x, 2) - x + 7.9);
- return s;
- }
- double lam(double x)
- {
- const double e = 2.71;
- double s;
- s = 3.1*pow(e, -pow(x, 2));
- return s;
- }
- double lamh(double x)
- {
- const double e = 2.71;
- double s;
- s = -6.2*pow(e, -pow(x, 2))*x;
- return s;
- }
- const double eps = 0.0000000001;
- int main()
- {
- int n, i, j;
- double x, y, y1;
- cout << "Enter n = ";
- cin >> n;
- for (i = 0; i < n; i++)
- {
- y1 = 0;
- y = 0;
- cout << "Enter x = ";
- cin >> x;
- y1 = -(psi(x)/fi(x));
- for (j = 0; true; j++)
- {
- y = y1 - ((y1*fi(x) + psi(x) + lam(y1)) / (fi(x) + lamh(y1)));
- if (abs(y - y1) <= eps)
- {
- break;
- }
- else
- {
- y1 = y;
- y = 0;
- }
- }
- cout << '(' << x << ',' << y << ')' << endl;
- cout << "F(" << x << ',' << y << ")= " << y * fi(x) + psi(x) + lam(y) << endl;
- }
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement