Advertisement
Guest User

biraj bro

a guest
Feb 26th, 2020
130
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.07 KB | None | 0 0
  1. #include<stdio.h>
  2. #include<conio.h>
  3. #include<math.h>
  4. main()
  5. {
  6. int i,j,n;
  7. float A[40][40],x[40],z[40],e[40],zmax,emax;
  8. printf("\nEnter the order of matrix:");
  9. scanf("%d",&n);
  10. printf("\nEnter matrix elements row-wise\n");
  11. for(i=1; i<=n; i++)
  12. {
  13. for(j=1; j<=n; j++)
  14. {
  15. printf("A[%d][%d]=", i,j);
  16. scanf("%f",&A[i][j]);
  17. }
  18. }
  19. printf("\nEnter the column vector\n");
  20. for(i=1; i<=n; i++)
  21. {
  22. printf("X[%d]=",i);
  23. scanf("%f",&x[i]);
  24. }
  25. do
  26. {
  27. for(i=1; i<=n; i++)
  28. {
  29. z[i]=0;
  30. for(j=1; j<=n; j++)
  31. {
  32. z[i]=z[i]+A[i][j]*x[j];
  33. }
  34. }
  35. zmax=fabs(z[1]);
  36. for(i=2; i<=n; i++)
  37. {
  38. if((fabs(z[i]))>zmax)
  39. zmax=fabs(z[i]);
  40. }
  41. for(i=1; i<=n; i++)
  42. {
  43. z[i]=z[i]/zmax;
  44. }
  45. for(i=1; i<=n; i++)
  46. {
  47. e[i]=0;
  48. e[i]=fabs((fabs(z[i]))-(fabs(x[i])));
  49. }
  50. emax=e[1];
  51. for(i=2; i<=n; i++)
  52. {
  53. if(e[i]>emax)
  54. emax=e[i];
  55. }
  56. for(i=1; i<=n; i++)
  57. {
  58. x[i]=z[i];
  59. }
  60. }
  61. while(emax>0.001);
  62. printf("\n The required eigen value is %f",zmax);
  63. printf("\n\nThe required eigen vector is :\n");
  64. for(i=1; i<=n; i++)
  65. {
  66. printf("%f\t",z[i]);
  67. }
  68. getch();
  69. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement