各位大侠，本人是一名机械工程师，是一个编程小白，但需要求解一个凸轮零件的曲线，我在网上找到了相关的技术资料，在电脑里下载了DEV-C++的编程软件，运行了这个文件，但不行；
源程序如下

 #include "stdio.h"
#include "math.h"
void main ()
{
 FILE*fp;
 double pi=3.1415926
 double L2,L1,L3,L4,L5;
 double L0;
 double xb,yb;
 double r1=20.0;//凸轮直径（可设定）
 double r2=23.5;//支持工件滚子直径（可设定）
 double Rmin=9.50,Rmax=69.0;//支持工件半径差值（可设定）
 double xita=165.0*pi/180.0;//两臂之间的夹角（可设定）
 double airfmin=0.0, airfmax=0.0, airf=0.0;
 double beita=0.0 omg=0.0;
 double t=0.0//加持工件半径差值
       xb=0.0;
    yb=0.0;
    L1=116.726, L2=116.726;//支撑臂长度（可设定）
    L3=80.0, L4=85.0;//工件回转中心与侧臂回转中心的距离（可设定）
    L5=104.0;//夹持最小工件时滚子到凸轮最小端面的距离（可设定）
    L0=sqrt(L3*L3+L4*L4);
    omg=atan(L3/L4);
    airfmin=acos((L1*L1+L0*L0-(Rmin+r2)*(Rmin+r2))/(2*L1*L0));//最小弧度
    airfmax=acos((L1*L1+L0*L0-(Rmax+r2)*(Rmax+r2))/(2*L1*L0));//最大弧度
       fp=fopen("c:\tulunquxian.dat","wb");
 for(t=0.0; t<=Rmax-Rmin; t=t+0.1)      
{
                 airf=acos((L1*L1+L0*L0-(Rmin+r2+t)*(Rmin+r2+1))/(2*L1*L0));
    if(airf<airfmin) break;
    if(airf>airfmax) break;
                beita=pi-xita-omg+airf;
    xb=L4+L2*cos(beita)-(L5+(Rmin+r2)/2+t);
    yb=L3-L2*sin(betia);
    fprintf(fp,"%1f%1f%1f\r\n",xb,yb,0.0);
    printf("xb=%1f,yb=%1f,\n",xb,yb);
}
     fclose(fp);
}

  请教各位大侠，我在电脑了安装什么软件来运行这个程序才能得到曲线的点的文件；我的电脑系统是WIN10;; 我的电话13589930219； QQ号3301141889；重谢