fork download
copy 
  1. from math import sin, cos
  2.  
  3.  
  4. def f1(x, y):
  5.     return ((y**2)*(sin(x+y)-1.6*x)-(cos(x+y))*(x**2+y**2-1)) / ((cos(x+y)-1.6)*(y**2)-(cos(x+y)*(x**2)))
  6.  
  7.  
  8. def f2(x, y):
  9.     return (-(x**2)*(sin(x+y)-1.6*x)+(cos(x+y)-1.6)*(x**2+y**2-1)) / ((cos(x+y)-1.6)*(y**2)-(cos(x+y)*(x**2)))
  10.  
  11.  
  12. x1 = 0.6
  13. x2 = 0.6
  14. for i in range(100):
  15.     print(i, x1, x2)
  16.     temp1 = x1 - f1(x1, x2)
  17.     temp2 = x2 - f2(x1, x2)
  18.     if abs(x1 - temp1) < 0.000001 and abs(x2 - temp2) < 0.000001:
  19.         print(temp1, temp2)
  20.         break
  21.     x1 = temp1
  22.     x2 = temp2
Success #stdin #stdout 0.02s 9764KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
0 0.6 0.6
1 0.7586705593778994 1.2191072183998783
2 0.7517489928885273 0.5073567440076151
3 0.7316327250712791 1.240934258271905
4 0.7439210856479251 0.5384408867695838
5 0.7020198831697975 1.1587934303282954
6 0.7093006688168685 0.5338150640256161
7 0.7370739701272804 1.2285164625085057
8 0.7435655114431248 0.5287943756365269
9 0.7129436985817128 1.188314508816299
10 0.7221675822904888 0.5332092488717408
11 0.7259602115170797 1.2091703603004107
12 0.7335545406027876 0.5299293230941687
13 0.7199752524514789 1.2007372479772576
14 0.7286248282768614 0.5316860291825689
15 0.7221251425958245 1.2033279282488454
16 0.7302847578908701 0.5308700196106699
17 0.7216657563105909 1.2031283505798647
18 0.7300203797657827 0.5311708982213574
19 0.7215335797439598 1.202642139675719
20 0.7298348781587565 0.5311021716497688
21 0.7217767614017496 1.2031484249794049
22 0.7300770673206792 0.5310879605943656
23 0.7215840170712614 1.2028021452811424
24 0.729899159748093 0.5311190831959329
25 0.7217007418968482 1.2029938705296226
26 0.7300023147947774 0.5310936731599739
27 0.7216418895985848 1.202904992026475
28 0.7299522935056304 0.5311092941808193
29 0.7216664635848856 1.2029380968430954
30 0.7299721575383225 0.531101316068375
31 0.7216588349430305 1.2029303237740887
32 0.7299666301424647 0.5311047000582945
33 0.7216595603133278 1.2029289621683343
34 0.7299666195141974 0.531103616751063
35 0.7216608699920827 1.2029322594068597
36 0.729968069442663 0.5311037473034228
37 0.7216594573636033 1.2029295334650094
38 0.7299667175185733 0.5311039075479094
39 0.7216604360096401 1.2029312213661876
40 0.7299676029106587 0.5311037238151075
41 0.7216598896617479 1.2029303544967218
42 0.7299671278779013 0.5311038537320951
43 0.7216601451665785 1.202930724696138
44 0.7299673410467443 0.5311037802142671
45 0.7216600487787024 1.2029306046499013
46 0.7299672656351912 0.5311038150571693
47 0.7216600722713007 1.2029306203666779
48 0.7299672805586603 0.531103801655731
49 0.721660075432782 1.2029306369567503
50 0.7299672862616685 0.5311038050996365
51 0.7216600663586441 1.2029306173618013
52 0.7299672770453722 0.5311038053732863
53 0.7216600740131427 1.2029306313476567
54 0.7299672841705716 0.5311038042204017
55 0.7216600692254079 1.2029306233853112
56 0.729967279914406 0.5311038052262711
57 0.721660071704519 1.2029306271810682
58 0.7299672820347559 0.5311038045868195
59 0.7216600706355452 1.2029306257089805
60 0.7299672811624943 0.5311038049222209
61 0.7216600709884466 1.2029306260955046
62 0.7299672814250686 0.5311038047754396
63 0.7216600709372146 1.2029306261178783
64 0.7299672814069833 0.5311038048252008
65 0.7216600708931433 1.2029306259962764
66 0.7299672813554746 0.5311038048169476
67 0.721660070947667 1.2029306261040627
68 0.7299672814083118 0.5311038048116948
69 0.7216600709082062 1.2029306260350359
70 0.7299672813723641 0.5311038048187529
71 0.7216600709308705 1.2029306260714476
72 0.7299672813921853 0.5311038048135268
73 0.7216600709199754 1.202930626055411
74 0.7299672813830315 0.5311038048165708
75 0.7216600709242501 1.2029306260609078
76 0.72996728138642 0.5311038048150883
77 0.7216600709230947 1.20293062605995
78 0.7299672813856388 0.5311038048156802
79 0.7216600709230608 1.2029306260594135
80 0.7299672813854862 0.5311038048155133
81 0.7216600709233933 1.202930626060164
82 0.7299672813858322 0.5311038048155152
83 0.7216600709230913 1.2029306260596022
84 0.7299672813855484 0.5311038048155572
85 0.7216600709232868 1.2029306260599313
86 0.7299672813857233 0.5311038048155172
87 0.7216600709231833 1.2029306260597714
88 0.7299672813856344 0.5311038048155434
89 0.7216600709232288 1.2029306260598347
90 0.7299672813856717 0.5311038048155295
91 0.7216600709232132 1.2029306260598167
92 0.7299672813856598 0.5311038048155358
93 0.7216600709232159 1.2029306260598174
94 0.7299672813856612 0.5311038048155337
95 0.7216600709232173 1.2029306260598212
96 0.7299672813856629 0.531103804815534
97 0.7216600709232155 1.2029306260598176
98 0.7299672813856611 0.5311038048155342
99 0.7216600709232167 1.2029306260598196
https://ideone.com/qpWJV1
language:
Python 3 nbc (python 3.7.3)
created:
4 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!