Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include "pch.h"
- #include <iostream>
- #include <cmath>
- #include <locale>
- using namespace std;
- double F(double x, double y) {
- return ((-0.5)*(y / x)) + x * x;
- }
- double tocnoe_rechenie(double x) {
- return ((2.0/7)*x*x*x) + ((5.0/7) / sqrt(x));
- }
- int main() {
- setlocale(LC_ALL, "Russian");
- double a = 1, b = 2, h = 0.1;
- double n = (b - a) / h;
- double X[12], Y[12], Y1[12], raz[12];
- X[0] = a, Y[0]=1;
- for (int i = 1; i < n+1; i++) {
- X[i] = a + i * h;
- Y[i] = Y[i - 1] + h * F(X[i - 1], Y[i - 1]);
- //Y[1] = 1;
- }
- for (int i = 1; i < n+1; i++) {
- cout << "X[" << i << "]=" << X[i] << " "<<"\n";
- }
- cout << endl;
- for (int i = 1; i < n+1; i++) {
- cout << "Y[" << i << "]=" << Y[i] << " "<<"\n";
- }
- cout << "\nТочное решение\n" << endl;
- for (int i = 1; i < n+1; i++) {
- Y1[i] = tocnoe_rechenie(X[i]);
- cout << "Y1[" << i << "]=" << Y1[i] << " " << "\n";
- }
- cout << "\nПогрешность\n" << endl;
- for (int i = 1; i < n+1; i++) {
- raz[i] = abs(Y1[i] - Y[i]);
- cout << raz[i] << "\n";
- }
- system("pause");
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement