fy">=1.0/phi2;
( calculation of nz ... n );
/******************** nz ******************/( i=0; i lt; n; i ++)
{[i]=0.0; E1 [i]=0.0; E2 [i]=1.0;
} (j=0; j lt; step; j ++)
{printf ( step% d n , j);
(i=0; i lt; n; i ++)
{ee1=e; ee2=e; ee3=e; ( i =% d b b b b b b b , i); (i2= 0; i2 lt; n; i2 ++)
(r2 (i, i2,3)!=0)
{= (E0 [i] * r2 (i, i2,0) + E1 [i] * r2 (i, i2,1) + E2 [i] * r2 (i, i2, 2)) * a * a * a/(r2 (i, i2,3) * r2 (i, i
2,3) * r2 (i, i2,3) * r2 (i, i2,3)); +=S * r2 (i, i2,0);
+=S * r2 (i, i2,1);
+=S * r2 (i, i2,2);
} [i]=ee1/n; [i]=ee2/n; [i]=ee3/n;
}
}
r0=2.0 * ASIZE;=0.0;=0.0; (i=0; i lt; n; i ++)
{= sqrt ((z [i] -r0) * (z [i] -r0) + (y [i] -ASIZE/2.0) * (y [i] -ASIZE/2.0 ) + (x [i] -
ASIZE/2.0) * (x [i] -ASIZE/2.0)); (r!=0) +=(E2 [i]) * a * a * a/r0/r0/ r0;
}=(S);
=1.0/phi2; ("nx =% f
ny =% fnz =% f" , nx/(nx + ny + nz), ny/(nx + ny + nz), nz/(nx + ny + nz));
printf ( nDepolarization coefficients were successfully calculated! n );
0;
}
