fork download
  1. #include <iostream>
  2. #include <cmath>
  3.  
  4. using std::cin;
  5.  
  6. int main() {
  7. double A[2][2], B[2][2];
  8. int I[2][2] = {1,0,0,1};
  9. for (int i=0; i<2; i++) {
  10. for (int j=0; j<2; j++) {
  11. cin>>A[i][j];
  12. }
  13. }
  14. double tr_A = A[0][0]+ A[1][1];
  15. double sqrt_det_A = sqrt(A[0][0]*A[1][1] - A[0][1]*A[1][0]);
  16. double divisor = sqrt(tr_A + 2*sqrt_det_A);
  17. for (int i=0; i<2; i++) {
  18. for (int j=0; j<2; j++) {
  19. B[i][j] = (A[i][j] + sqrt_det_A*I[i][j]) / divisor;
  20. }
  21. }
  22. printf("Matrix A:\n%f %f\n%f %f\n", A[0][0], A[0][1], A[1][0], A[1][1]);
  23. printf("Square root of matrix A:\n%f %f\n%f %f", B[0][0], B[0][1], B[1][0], B[1][1]);
  24. return 0;
  25. }
Success #stdin #stdout 0.01s 5272KB
stdin
2
1
2
3
stdout
Matrix A:
2.000000 1.000000
2.000000 3.000000
Square root of matrix A:
1.333333 0.333333
0.666667 1.666667