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){
- int i;
- for (i = 0; i < db; ++i){
- printf("%lf\n", tomb[i]);
- }
- }
- double tavolsag (double PR[], double PRv[], int n){
- int i;
- double osszeg = 0.0;
- for (i = 0; i < n; ++i)
- osszeg += (PRv[i] - PR[i]) * (PRv[i] - PR[i]);
- return sqrt(osszeg);
- }
- int main (){
- double L[4][4] = {
- {0.0, 0.0, 1.0 / 3.0, 0.0},
- {1.0, 1.0 / 2.0, 1.0 / 3.0, 0.0},
- {0.0, 1.0 / 2.0, 0.0, 0.0},
- {0.0, 0.0, 1.0 / 3.0, 1.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;
- double d = 0.5;
- for (;;){
- for (i = 0; i < 4; ++i)
- {
- PR[i]=0.0;
- for (j = 0; j < 4 ; ++j){
- PR[i] += (PRv[j] * L[i][j] + (1-d)/4);
- }
- PR[i] *= d;
- }
- if(tavolsag(PR,PRv, 4) < 0.000000001)
- break;
- for (i = 0; i < 4; ++i){
- PRv[i]=PR[i];
- }
- }
- kiir (PR, 4);
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement