机械优化设计，曲柄摇杆机构解析法C程序

析法运动分析及力分析C语言源程序
#include "stdio.h"
#include "math.h"
#define PI 3.1415926
#define g 9.8
void main()
{double l1=0.13,l2=0.13,l3=0.13,l4=0.26,l5=0.11,l6=0.065,l7=0.13,G3=100,G4=240,
G5=220,Pr=2000,J3=0.012,J4=2.8;
double x=0.016,y=0.01,a1,a2,a3,a4,w1,w2,w3,w4,q,q1,A4x,A4y,A3x,A3y,F23x,F23y,
F34x,F34y,F45x,F45y,F63x,F63y,M;
double c,c4,s,v,a,A,B,C,A1,B1,C1;
int N;
float p[20]={0,30,60,90,120,150,180,210,240,270,300,330,360,239,60,229,70};
for (N=0;N<17;N++)