#include <stdio.h>
#include <math.h>
#include <mem.h>
#include <conio.h>
#include <stdlib.h>
#include <string.h>

#define N 100
#define Q 225
#define PI 3.14159
int xx[N][N];
float Qb=150;
static float n=0,j0=0,j2=0,Fcr=0,F=0,H=0,u=0,x=0;
float Fp=200,f=0.11, H1=210,E=2.06*10*10*10*10*10;
float D,D0,D1,D2,D3,D4,D5,T1,T2,Lp,L,d1,d;
float

i=0,j,k,F,h1,H,H1,Fp,d,I,l,d0,d1,d2,p,j0,j1=PI/6,j2,S,t,b,u,qb,dp;

main()
{ int i=0;
loop: i++;
{if(luoganxiaojing())
goto loop;

if(luoganshengjiaodangliangmocajiao())
goto loop;
H=2.5*d2;
if(luoganwendingxing())
goto loop;

if(luowenyaqiangdu())
goto loop;
S=n*p;
h1=2*d;
D1=0.8*d;
D=1.8*d;
D2=1.5*d;
D3=1.4*D2;
D5=D3+(H1+H+h1)*tan(75/180.0*PI);
D4=1.4*D5;
j0=atan(S/(PI*d2));
j2=atan(f/cos(j1/2));
T1=F*d2*tan(j0+j2)/2.0;
d0=D1+1;
D0=D-1;
T2=1/3.0*f*0.15*F*(D0*D0*D0-d0*d0*d0)/(D0*D0-d0*d0);
Lp=(T1+T2)/Fp;
L=Lp+D/2.0;
xiaolv();
}
printf("%d",L);
}

int luoganxiaojing()/*螺纹小径*/
{
if(i<1.25*1000*sqrt(F*1.0/(Q*10*10*10)))
return(1);
}

int luoganshengjiaodangliangmocajiao()
/*螺纹升角和当量摩擦角*/
{
printf("Please input the value of d1,d,D,d2,p,where

d1=%d\n",i);
scanf("%d%d%d%d%d",&d1,&d,&D,&d2,&p);
j0=atan(S/(PI*d2));
j2=atan(f/cos(j1/2));
if(j0>j2||d2<0.8*sqrt(F/(2.5*7.5)))
return(1);
}

int luoganwendingxing() /*螺纹杆稳定型*/
{
I=PI*d1*d1*d1*d1/64.0;
Fcr=PI*PI*E*I/(4*(H1+h1+H/2.0)*(H1+h1+H/2.0));

if(Fcr/F<3.5)
return(1);
}

int luowenyaqiangdu()/*螺纹牙强度*/
{
t=F/(PI*D*b*u);
qb=6*F*l/(PI*D*b*b*u);
if(t>(Q*0.6)||qb>Q)
return(1);
}

int xiaolv()/*效率*/
{dp=pow((Fp*Lp/(0.1*Qb)),-3);
x=Fp/(2*PI*(T1+T2));
}