write a program to find solution of linear equation using Power method in c++.

 write a program to find solution of linear equation using Power method in c++..


Source Code:




/********************************************
    language: C++ CONA

    Author : Tanveer Khan

    Rajasthan University, Jaipur

************************************************/
#include<iostream.h>
#include<conio.h>
void main()
{
clrscr();
float a[4][4],x[4],s[4],lamda,temp,iteration;
int i,j,n;
cout<<"enter number of equations: ";
cin>>n;
cout<<"enter number of iteratins: ";
cin>>iteration;
cout<<"enter coefficients equation wise:\n";
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
{
cout<<"position("<<i<<","<<j<<")";
cin>>a[i][j];
}
}
cout<<"set value of x:\n";
for(i=1;i<=n;i++)
{
x[i]=1.0;
cout<<x[i]<<endl;
}
for(int k=1;k<=iteration;k++)
{
for(i=1;i<=3;i++)
{
s[i]=0.0;
for(j=1;j<=n;j++)
{
s[i]+=a[i][j]*x[j];
}
}
lamda=s[1];
for(i=1;i<=n;i++)
{
 temp=s[i]/lamda;
x[i]=temp;
}
}
cout<<"\nlargest eigen value="<<lamda;
cout<<"\nassociated eigen vector=\n";
for(i=1;i<=n;i++)
{
cout<<x[i]<<endl;
}
getch();
}


OUTPUT: