- #36
xts
- 885
- 0
Don't try to be vicious! It took me about 3 mins to write and run the following C program:Misr said:I hope you spend some time checking these calculations
Code:
#define DEGTORAD (2*3.141592654/360)
#define alpha (10*DEGTORAD)
#define n (1.5)
void misr(void) {
for (int i=-40; i<=40; i+=10) {
double beta = DEGTORAD*i;
double phi = asin(sin(beta)/n);
double psi = phi + alpha;
double theta = asin(n*sin(psi));
double delta = theta - alpha - beta;
printf ("%3.0f %2.1f\n", beta/DEGTORAD, delta/DEGTORAD);
}
}
-40 6.6
-30 5.7
-20 5.2
-10 5.0
0 5.1
10 5.5
20 6.2
30 7.6
40 10.3
Great! So now make remaining plots (for other apical angles) and compare all of them.I tried to be more accurate this time