没改多少,就是一些类型转换
#include
#include
#include
#define PI 3。1415926
#define t PI/180
void main()
{
float a0=85。
0f,b0=23。0f,m=0。588792f,n=1。813319f,l=2。398730f;
float A,B,C,s,j,k1,k2,k;
float x0=1。0f, y0=0。0f,u1=1。
0f/45,u2=0。301f/85;
float x[100],y1[100],y2[100],a1[100],f1[100],f2[100],r[100];
int i;
FILE *p;
if((p=fopen("d:\yj。
txt","w"))==NULL)
{
printf("can't open the file!");
exit(0);
}
printf(" i a1[i] f1[i] r[i] k x[i] y1[i] y2[i]
");
fprintf(p," i a1[i] f1[i] r[i] k x[i] y1[i] y2[i]
");
for(i=0;a1[i]<=45;i )
{
a1[0]=0;
A=(float)sin((a1[i] a0)*(float)t);
B=(float)cos((a1[i] a0)*(float)t)-n;
C=(1 m*m n*n-l*l)/(2*m)-n*(float)cos((a1[i] a0)*(float)t)/m;
j=x0 u1*a1[i];
s=(float)sqrt(A*A B*B-C*C);
f1[i]=2*((float)atan((A s)/(B C)))/((float)t)-b0;
f2[i]=2*((float)atan((A-s)/(B C)))/((float)t)-b0;
r[i]=((float)log(j)-y0)/u2;
k1=f1[i]-r[i];
k2=f2[i]-r[i];
x[i]=a1[i]*u1 x0;
y2[i]=(float)log(x[i]);
if(fabs(k1) {
k=k1;
y1[i]=f1[i]*u2 y0;
printf(" %-4d %-5。
1f %-10。4f %-8。4f %-8。4f %-7。4f %-8。4f %0。4f
",
i,a1[i],f1[i],r[i],k,x[i],y1[i],y2[i]);
fprintf(p," %-4d %-5。
1f %-10。4f %-8。4f %-8。4f %-7。4f %-8。4f %0。4f
",
i,a1[i],f1[i],r[i],k,x[i],y1[i],y2[i]);
}
else
{
k=k2;
y1[i]=f2[i]*u2 y0;
printf(" %-6d %-7。
1f %-12。4f %-10。4f %-10。4f %-9。4f %-10。4f %2。4f
",
i,a1[i],f2[i],r[i],k,x[i],y1[i],y2[i]);
fprintf(p," %-6d %-7。
1f %-12。4f %-10。4f %-10。4f %-9。4f %-10。4f %2。4f
",
i,a1[i],f2[i],r[i],k,x[i],y1[i],y2[i]);
}
a1[i 1]=a1[i] 0。
5f;
}
fclose(p);
}
结果:
i a1[i] f1[i] r[i] k x[i] y1[i] y2[i]
0 0。
0 -0。6893 0。0000 -0。6893
1。0000 -0。0024 0。0000
1 0。5 0。5527 3。1204 -2。5676
1。0111 0。0020 0。
0110
2 1。0 1。7708 6。2067 -4。4359
1。0222 0。0063 0。0220
3 1。5 2。9664 9。2596 -6。2931
1。
0333 0。0105 0。0328
4 2。0 4。1410 12。2798 -8。1389
1。0444 0。0147 0。0435
5 2。5 5。2956 15。2682 -9。
9726
1。0556 0。0188 0。0541
6 3。0 6。4315 18。2252 -11。7937
1。0667 0。0228 0。0645
7 3。
5 7。5496 21。1515 -13。6019
1。0778 0。0267 0。0749
8 4。0 8。6509 24。0479 -15。3970
1。0889 0。0306 0。
0852
9 4。5 9。7363 26。9148 -17。1786
1。1000 0。0345 0。0953
10 5。0 10。8064 29。7530 -18。9466
1。
1111 0。0383 0。1054
11 5。5 11。8620 32。5629 -20。7008
1。1222 0。0420 0。1153
12 6。0 12。9039 35。
3451 -22。4412
1。1333 0。0457 0。1252
13 6。5 13。9325 38。1001 -24。1676
1。1444 0。0493 0。1349
14 7。
0 14。9485 40。8286 -25。8801
1。1556 0。0529 0。1446
15 7。5 15。9525 43。5309 -27。5784
1。1667 0。
0565 0。1542
16 8。0 16。9448 46。2076 -29。2628
1。1778 0。0600 0。1636
17 8。5 17。9261 48。8592 -30。
9332
1。1889 0。0635 0。1730
18 9。0 18。8966 51。4862 -32。5896
1。2000 0。0669 0。1823
19 9。
5 19。8568 54。0889 -34。2320
1。2111 0。0703 0。1915
20 10。0 20。8072 56。6678 -35。8606
1。2222 0。
0737 0。2007
21 10。5 21。7480 59。2234 -37。4755
1。2333 0。0770 0。2097
22 11。0 22。6795 61。7561 -39。
0766
1。2444 0。0803 0。2187
23 11。5 23。6022 64。2662 -40。6640
1。2556 0。0836 0。2276
24 12。
0 24。5163 66。7543 -42。2380
1。2667 0。0868 0。2364
25 12。5 25。4221 69。2206 -43。7985
1。2778 0。
0900 0。2451
26 13。0 26。3198 71。6656 -45。3458
1。2889 0。0932 0。2538
27 13。5 27。2098 74。0896 -46。
8798
1。3000 0。0964 0。2624
28 14。0 28。0922 76。4929 -48。4007
1。3111 0。0995 0。2709
29 14。
5 28。9673 78。8760 -49。9087
1。3222 0。1026 0。2793
30 15。0 29。8352 81。2391 -51。4039
1。3333 0。
1057 0。2877
31 15。5 30。6963 83。5826 -52。8863
1。3444 0。1087 0。2960
32 16。0 31。5506 85。9069 -54。
3562
1。3556 0。1117 0。3042
33 16。5 32。3984 88。2121 -55。8137
1。3667 0。1147 0。3124
34 17。
0 33。2399 90。4987 -57。2588
1。3778 0。1177 0。3205
35 17。5 34。0752 92。7669 -58。6918
1。3889 0。
1207 0。3285
36 18。0 34。9044 95。0171 -60。1127
1。4000 0。1236 0。3365
37 18。5 35。7278 97。2494 -61。
5217
1。4111 0。1265 0。3444
38 19。0 36。5454 99。4643 -62。9189
1。4222 0。1294 0。3522
39 19。
5 37。3574 101。6619 -64。3045
1。4333 0。1323 0。3600
40 20。0 38。1640 103。8425 -65。6786
1。4444 0。
1351 0。3677
41 20。5 38。9652 106。0065 -67。0413
1。4556 0。1380 0。3754
42 21。0 39。7612 108。1540 -68。
3928
1。4667 0。1408 0。3830
43 21。5 40。5521 110。2852 -69。7331
1。4778 0。1436 0。3905
44 22。
0 41。3380 112。4005 -71。0625
1。4889 0。1464 0。3980
45 22。5 42。1190 114。5001 -72。3811
1。5000 0。
1492 0。4055
46 23。0 42。8952 116。5842 -73。6890
1。5111 0。1519 0。4128
47 23。5 43。6668 118。6530 -74。
9863
1。5222 0。1546 0。4202
48 24。0 44。4337 120。7068 -76。2731
1。5333 0。1573 0。4274
49 24。
5 45。1961 122。7457 -77。5497
1。5444 0。1600 0。4347
50 25。0 45。9540 124。7701 -78。8160
1。5556 0。
1627 0。4418
51 25。5 46。7076 126。7800 -80。0723
1。5667 0。1654 0。4490
52 26。0 47。4570 128。7757 -81。
3187
1。5778 0。1681 0。4560
53 26。5 48。2021 130。7574 -82。5553
1。5889 0。1707 0。4630
54 27。
0 48。9431 132。7253 -83。7822
1。6000 0。1733 0。4700
55 27。5 49。6800 134。6796 -84。9995
1。6111 0。
1759 0。4769
56 28。0 50。4130 136。6204 -86。2074
1。6222 0。1785 0。4838
57 28。5 51。1419 138。5480 -87。
4061
1。6333 0。1811 0。4906
58 29。0 51。8670 140。4625 -88。5955
1。6444 0。1837 0。4974
59 29。
5 52。5883 142。3642 -89。7759
1。6556 0。1862 0。5041
60 30。0 53。3058 144。2531 -90。9473
1。6667 0。
1888 0。5108
61 30。5 54。0196 146。1295 -92。1098
1。6778 0。1913 0。5175
62 31。0 54。7297 147。9934 -93。
2637
1。6889 0。1938 0。5241
63 31。5 55。4362 149。8452 -94。4090
1。7000 0。1963 0。5306
64 32。
0 56。1391 151。6849 -95。5458
1。7111 0。1988 0。5371
65 32。5 56。8385 153。5127 -96。6742
1。7222 0。
2013 0。5436
66 33。0 57。5343 155。3287 -97。7944
1。7333 0。2037 0。5500
67 33。5 58。2267 157。1331 -98。
9064
1。7444 0。2062 0。5564
68 34。0 58。9157 158。9261 -100。0104
1。7556 0。2086 0。5628
69 34。
5 59。6013 160。7078 -101。1065
1。7667 0。2111 0。5691
70 35。0 60。2835 162。4783 -102。1947
1。7778 0。
2135 0。5754
71 35。5 60。9624 164。2377 -103。2753
1。7889 0。2159 0。5816
72 36。0 61。6380 165。
9863 -104。3483
1。8000 0。2183 0。5878
73 36。5 62。3104 167。7241 -105。4137
1。8111 0。2207 0。5939
74 37。
0 62。9795 169。4512 -106。4718
1。8222 0。2230 0。6001
75 37。5 63。6454 171。1679 -107。5226
1。8333 0。
2254 0。6061
76 38。0 64。3080 172。8742 -108。5662
1。8444 0。2277 0。6122
77 38。5 64。9675 174。
5703 -109。6027
1。8556 0。2301 0。6182
78 39。0 65。6239 176。2562 -110。6323
1。8667 0。2324 0。6242
79 39。
5 66。2771 177。9321 -111。6550
1。8778 0。2347 0。6301
80 40。0 66。9272 179。5982 -112。6709
1。8889 0。
2370 0。6360
81 40。5 67。5742 181。2544 -113。6802
1。9000 0。2393 0。6419
82 41。0 68。2181 182。
9010 -114。6829
1。9111 0。2416 0。6477
83 41。5 68。8590 184。5381 -115。6791
1。9222 0。2438 0。6535
84 42。
0 69。4967 186。1657 -116。6690
1。9333 0。2461 0。6592
85 42。5 70。1315 187。7840 -117。6525
1。9444 0。
2483 0。6650
86 43。0 70。7631 189。3931 -118。6300
1。9556 0。2506 0。6707
87 43。5 71。3918 190。
9930 -119。6013
1。9667 0。2528 0。6763
88 44。0 72。0174 192。5840 -120。5666
1。9778 0。2550 0。6820
89 44。
5 72。6399 194。1660 -121。5261
1。9889 0。2572 0。6876
90 45。0 73。2595 195。7392 -122。4797
2。0000 0。
2594 0。6931
Press any key to continue。
全部