//为什么只从61开始，不从0开始
#include<stdio.h> 
#include<math.h> 
#define PI 3.1415926 
void main(){ 
double a=0.125,b=0.600,c=0.150,d=0.275,e=0.575,f=1;
double B,C,E,F,G,I,L,M,O;/*B=θ3,C=θ4,E=Se，F =ω3，G=ω4，I=Ve，L=а3，M=а4，O=аe*/
double x=0; 
 printf(" @1     @3     @4      Se      W3     W4     Ve      A3     A4     Ae \n"); 
 for(x=0;x<6.3;x+=PI*1/180) { 
 B=atan((d+a*sin(x))/(a*cos(x))); /*求θ3*/ 
 if(B<0)
 B=PI+B; 
 C=PI-asin((e-b*sin(B))/c); /*求θ4*/ 
 if(C<0)
 C=PI+C; 
 E=b*cos(B)+c*cos(C); /*求 Se*/ 
 F=(a*f*(a+d*sin(x)))/(d*d+a*a+2*d*a*sin(x)); /*求ω3*/
 G=-(F*b*cos(B))/(c*cos(C)); /*求ω4*/ 
 I=-(F*b*sin(B-C))/cos(C); /*求Ve*/ 
 L=((d*d-a*a)*d*a*f*f*cos(x))/((d*d+a*a+2*d*a*sin(x))*(d*d+a*a+2*d*a*sin(x))); /*求а3*/ 
 M=(F*F*b*sin(B)+G*G*c*sin(C)-L*b*cos(B))/(c*cos(C)); /*求а4*/ 
 O=-(L*b*sin(B-C)+F*F*b*cos(B-C)-G*G*c)/cos(C);/*求аe*/
printf("%3.0f  %3.3f  %3.3f  %3.3f  %3.3f  %3.3f  %3.3f  %3.3f  %3.3f  %3.3f\n",x*180/PI,(B*180)/PI,(C*180)/PI,E,F,G,I,L,M,O);
} }