Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <iostream>
- #include <cmath>
- using namespace std;
- double f(double x, double y){
- double fi = 3 + x*x,
- psi = sqrt(1+x*x)/(1+0.4*x*x),
- L = 2*sin(1+y),
- F = y*fi + psi + L;
- return F;
- }
- double f1(double x, double y){
- double fi = 3 + x*x,
- L1 = 2*cos(1+y),
- F1 = fi + L1;
- return F1;
- }
- int main() {
- int n;
- cin >> n;
- double eps, x;
- cin >> eps;
- for(int j = 0; j< n; j++){
- double y,k;
- cin >> x;
- y = -(sqrt(1+x*x)/(1+0.4*x*x))/(3 + x*x);
- k = y - f(x,y)/f1(x,y);
- while(abs(k - y)>eps){
- y = k - f(x,k)/f1(x,k);
- k = y - f(x,y)/f1(x,y);
- }
- cout << x <<' '<< y << endl;
- cout<<f(x,k)<<endl;
- }
- system("pause");
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement