#include <iostream> #include <iomanip> #include <cmath> using namespace std; double f1(double x1, double x2, double t) // первая функция { return -52 * x1 - 100 * x2 + exp(-t); } double f2(double x1, double t) // вторая функция { return x1 + sin(t); } int main() { //setlocale(LC_ALL,"rus"); double h = 0.01; // шаг double x1 = 1.0; double x2 = 0.0; double t = 0.0; cout << setw(10) << t << setw(12) << x1 << setw(12) << x2 << endl; while (t <= 2.0) { double k1 = h*f1(x1,x2,t); double m1 = h*f2(x1,t); double k2 = h*f1(x1+k1/2,x2+m1/2,t+h/2); double m2 = h*f2(x1+k1/2,t+h/2); double k3 = h*f1(x1+k2/2,x2+m2/2,t+h/2); double m3 = h*f2(x1+k2/2,t+h/2); double k4 = h*f1(x1+k3,x2+m3,t+h); double m4 = h*f2(x1+k2,t+h); x1 += (k1+2*k2+2*k3+k4)/6; x2 += (m1+2*m2+2*m3+m4)/6; t += h; cout << setw(10) << t << setw(12) << x1 << setw(12) << x2 << endl; } //system("pause"); return 0; }
Standard input is empty
0 1 0 0.01 0.598934 0.00792517 0.02 0.355535 0.012775 0.03 0.207742 0.0157975 0.04 0.117897 0.017749 0.05 0.0631435 0.019088 0.06 0.0296148 0.0200918 0.07 0.00889575 0.0209281 0.08 -0.00411902 0.0216978 0.09 -0.0125264 0.0224615 0.1 -0.0182044 0.0232552 0.11 -0.0222913 0.0241001 0.12 -0.0254767 0.0250084 0.13 -0.0281777 0.0259867 0.14 -0.030646 0.0270384 0.15 -0.0330331 0.028165 0.16 -0.0354297 0.0293665 0.17 -0.0378895 0.0306425 0.18 -0.0404439 0.031992 0.19 -0.0431106 0.0334138 0.2 -0.0458993 0.0349065 0.21 -0.0488146 0.0364687 0.22 -0.051858 0.0380989 0.23 -0.0550294 0.0397956 0.24 -0.0583275 0.0415574 0.25 -0.0617503 0.0433826 0.26 -0.0652956 0.04527 0.27 -0.0689609 0.0472179 0.28 -0.0727435 0.0492249 0.29 -0.0766408 0.0512896 0.3 -0.0806501 0.0534106 0.31 -0.0847687 0.0555866 0.32 -0.088994 0.057816 0.33 -0.0933231 0.0600976 0.34 -0.0977537 0.0624299 0.35 -0.102283 0.0648118 0.36 -0.106908 0.0672418 0.37 -0.111627 0.0697187 0.38 -0.116438 0.0722411 0.39 -0.121337 0.0748079 0.4 -0.126322 0.0774178 0.41 -0.13139 0.0800694 0.42 -0.136541 0.0827617 0.43 -0.14177 0.0854934 0.44 -0.147076 0.0882634 0.45 -0.152456 0.0910703 0.46 -0.157909 0.0939132 0.47 -0.163431 0.0967907 0.48 -0.169021 0.0997019 0.49 -0.174677 0.102646 0.5 -0.180395 0.10562 0.51 -0.186175 0.108626 0.52 -0.192014 0.11166 0.53 -0.19791 0.114723 0.54 -0.20386 0.117812 0.55 -0.209864 0.120928 0.56 -0.215918 0.124068 0.57 -0.222021 0.127233 0.58 -0.228171 0.13042 0.59 -0.234365 0.13363 0.6 -0.240603 0.13686 0.61 -0.246881 0.14011 0.62 -0.253199 0.143379 0.63 -0.259554 0.146667 0.64 -0.265944 0.149971 0.65 -0.272367 0.153291 0.66 -0.278823 0.156627 0.67 -0.285308 0.159977 0.68 -0.291821 0.16334 0.69 -0.298361 0.166716 0.7 -0.304926 0.170103 0.71 -0.311513 0.173501 0.72 -0.318122 0.176909 0.73 -0.32475 0.180326 0.74 -0.331396 0.183752 0.75 -0.338058 0.187184 0.76 -0.344735 0.190623 0.77 -0.351426 0.194067 0.78 -0.358127 0.197517 0.79 -0.364838 0.20097 0.8 -0.371558 0.204427 0.81 -0.378285 0.207886 0.82 -0.385016 0.211347 0.83 -0.391752 0.214808 0.84 -0.39849 0.21827 0.85 -0.405228 0.221731 0.86 -0.411966 0.225191 0.87 -0.418701 0.228648 0.88 -0.425434 0.232103 0.89 -0.432161 0.235554 0.9 -0.438881 0.239001 0.91 -0.445594 0.242443 0.92 -0.452298 0.245879 0.93 -0.458992 0.249309 0.94 -0.465673 0.252731 0.95 -0.472341 0.256146 0.96 -0.478995 0.259553 0.97 -0.485633 0.26295 0.98 -0.492254 0.266337 0.99 -0.498857 0.269715 1 -0.50544 0.273081 1.01 -0.512003 0.276435 1.02 -0.518543 0.279777 1.03 -0.52506 0.283106 1.04 -0.531552 0.286422 1.05 -0.538019 0.289723 1.06 -0.544459 0.293009 1.07 -0.550871 0.296281 1.08 -0.557254 0.299536 1.09 -0.563606 0.302775 1.1 -0.569928 0.305996 1.11 -0.576216 0.3092 1.12 -0.582471 0.312386 1.13 -0.588691 0.315552 1.14 -0.594876 0.3187 1.15 -0.601023 0.321827 1.16 -0.607133 0.324934 1.17 -0.613204 0.328021 1.18 -0.619235 0.331085 1.19 -0.625224 0.334128 1.2 -0.631172 0.337148 1.21 -0.637077 0.340145 1.22 -0.642938 0.343119 1.23 -0.648754 0.346068 1.24 -0.654524 0.348993 1.25 -0.660247 0.351893 1.26 -0.665923 0.354768 1.27 -0.671549 0.357616 1.28 -0.677127 0.360439 1.29 -0.682654 0.363234 1.3 -0.688129 0.366002 1.31 -0.693552 0.368742 1.32 -0.698923 0.371455 1.33 -0.704239 0.374138 1.34 -0.7095 0.376793 1.35 -0.714706 0.379418 1.36 -0.719856 0.382013 1.37 -0.724948 0.384578 1.38 -0.729983 0.387112 1.39 -0.734958 0.389615 1.4 -0.739874 0.392087 1.41 -0.74473 0.394526 1.42 -0.749525 0.396934 1.43 -0.754257 0.399309 1.44 -0.758928 0.401651 1.45 -0.763535 0.403959 1.46 -0.768078 0.406234 1.47 -0.772557 0.408475 1.48 -0.77697 0.410681 1.49 -0.781318 0.412853 1.5 -0.785598 0.41499 1.51 -0.789812 0.417091 1.52 -0.793957 0.419156 1.53 -0.798034 0.421186 1.54 -0.802042 0.423179 1.55 -0.80598 0.425135 1.56 -0.809848 0.427055 1.57 -0.813645 0.428937 1.58 -0.81737 0.430782 1.59 -0.821024 0.432589 1.6 -0.824605 0.434358 1.61 -0.828112 0.436088 1.62 -0.831546 0.43778 1.63 -0.834906 0.439433 1.64 -0.838191 0.441047 1.65 -0.841401 0.442621 1.66 -0.844535 0.444156 1.67 -0.847593 0.445651 1.68 -0.850575 0.447106 1.69 -0.853479 0.44852 1.7 -0.856306 0.449894 1.71 -0.859055 0.451227 1.72 -0.861726 0.452519 1.73 -0.864318 0.45377 1.74 -0.86683 0.45498 1.75 -0.869264 0.456148 1.76 -0.871617 0.457274 1.77 -0.87389 0.458359 1.78 -0.876082 0.459401 1.79 -0.878194 0.460401 1.8 -0.880224 0.461359 1.81 -0.882172 0.462273 1.82 -0.884038 0.463146 1.83 -0.885823 0.463975 1.84 -0.887524 0.464761 1.85 -0.889143 0.465504 1.86 -0.890679 0.466204 1.87 -0.892131 0.46686 1.88 -0.893499 0.467472 1.89 -0.894784 0.468041 1.9 -0.895985 0.468566 1.91 -0.897101 0.469048 1.92 -0.898133 0.469485 1.93 -0.89908 0.469878 1.94 -0.899942 0.470227 1.95 -0.900719 0.470531 1.96 -0.901411 0.470791 1.97 -0.902017 0.471007 1.98 -0.902538 0.471178 1.99 -0.902973 0.471305 2 -0.903322 0.471387