#include
#include
#include
using namespace std;
int main()
{
cout.precision(2);
cout.setf(ios::fixed);
int n,i,j;
cout<<"\nEnter the order of the matrix:\n"; //Get the order of the square matrix
cin>>n;
double a[n][n],b[n],c[n],k,eps,y; //declare arrays and other variables to be used
cout<<"\nEnter the elements of matrix row-wise:\n";
for (i=0;i>a[i][j];
cout<<"\nEnter the initial values of the eigen-vector:\n";
for (i=0;i>b[i];
cout<<"\nEnter the accuracy desired:\n";
cin>>eps;
k=b[0]; //Assign some initial value to the eigen value, 'k'
do
{
y=k;
for (i=0;ik?abs(c[i]):k; //Find the largest element of c and assign it to k where k is the eigen value
for (i=0;ieps); //Check if the error in eigen value is within the tolerance limit
cout<<"The eigen-value is: "<