Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include<stdio.h>
- #include<math.h>
- void kiir (double tomb[], int db){
- for (int i=0; i<db; i++)
- printf("%lf\n", tomb[i]);
- }
- double tavolsag(double PR[], double PRV[],int n){
- double osszeg=0.0;
- for (int i=0;i<n;++i){
- osszeg += (PRV[i] - PR[i]) * (PRV[i] - PR[i]);
- return sqrt(osszeg);
- }
- }
- int main(void){
- double L[4][4]={ //J JP JPL M
- {0.0, 0.0, 1.0/3.0, 0.0},
- {1.0, 1.0/2.0, 1.0/3.0, 1.0},
- {0.0, 1.0/2.0, 0.0, 0.0},
- {0.0, 0.0, 1.0/3.0, 0.0}
- };
- double PR[4]={0.0, 0.0, 0.0 ,0.0};
- double PRV[4]={1.0/4.0, 1.0/4.0, 1.0/4.0, 1.0/4.0};
- int i,j;
- for(;;){
- for(i=0;i<4;i++)
- PR[i] = PRV[i];
- for (i=0;i<4;i++){
- double ideg=0;
- for (j=0;j<4;j++)
- ideg+=L[i][j]*PR[j];
- PRV[i]=ideg;
- }
- if(tavolsag(PR,PRV,4)<0.000000001)
- break;
- }
- kiir(PR,4);
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement