Scribd
Upload
9 views

Eigen Value Using Power Method

Uploaded by

ad599066
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
9 views

Eigen Value Using Power Method

Uploaded by

ad599066
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

#include<stdio.

h>
#include<conio.h>
#include<math.h>
void main()
{

int i,j,n;
float A[40][40],x[40],z[40],e[40],zmax,emax;
printf("\nEnter the order of matrix:");
scanf("%d",&n);
printf("\nEnter matrix elements row-wise\n");
for(i=1; i<=n; i++){
for(j=1; j<=n; j++){
printf("A[%d][%d]=", i,j);
scanf("%f",&A[i][j]);
}
}
printf("\nEnter the column vector\n");
for(i=1; i<=n; i++){
printf("X[%d]=",i);
scanf("%f",&x[i]);
}
do
{
for(i=1; i<=n; i++){
z[i]=0;
for(j=1; j<=n; j++){
z[i]=z[i]+A[i][j]*x[j];
}
}
zmax=fabs(z[1]);
for(i=2; i<=n; i++){
if((fabs(z[i]))>zmax)
zmax=fabs(z[i]);
}
for(i=1; i<=n; i++){
z[i]=z[i]/zmax;
}
for(i=1; i<=n; i++){
e[i]=0;
e[i]=fabs((fabs(z[i]))-(fabs(x[i])));
}
emax=e[1];
for(i=2; i<=n; i++){
if(e[i]>emax)
emax=e[i];
}
for(i=1; i<=n; i++){
x[i]=z[i];
}
}
while(emax>0.001);
printf("\n The required eigen value is %f",zmax);
printf("\n\nThe required eigen vector is :\n");
for(i=1; i<=n; i++){
printf("%f\t",z[i]);
}
getch();
}
OUTPUT:
Enter the order of matrix:3 Enter the order of matrix:2

Enter matrix elements row-wise Enter matrix elements row-wise


A[1][1]=2 A[1][1]=2
A[1][2]=1 A[1][2]=1
A[1][3]=3 A[2][1]=3
A[2][1]=-2 A[2][2]=-1
A[2][2]=1
A[2][3]=0 Enter the column vector
A[3][1]=4 X[1]=2
A[3][2]=2 X[2]=3
A[3][3]=1
The required eigen value is 2.791796
Enter the column vector
X[1]=2 The required eigen vector is :
X[2]=1 1.000000 0.790962
X[3]=3
...Program finished with exit code 0
The required eigen value is 4.003649 Press ENTER to exit console.

The required eigen vector is :


1.000000 -0.665705 0.888973

...Program finished with exit code 0


Press ENTER to exit console.

You might also like