刚体转动惯量的测量

#include <stdio.h>
#include <stdlib.h>

int main()
{
    double j0,j1,j2,b1,b2,a,b3,b4,j3,j4,b5,b6,b7,b8,j5,j6,j7,j8,j9;
    double m1=0.05,g=9.8,r1=0.019,r2=0.1,m2=0.5829,r3=0.086,r4=0.1,m3=0.6358,m4=0.1726;
    scanf("%lf%lf%lf%lf",&b2,&b1,&b4,&b3);
    j0=m1*r1*(g-r1*b2)/(b2-b1);
    printf("空台%lf\n",j0);
    j2=m1*r1*(g-r1*b4)/(b4-b3);
    j3=j2-j0;
    printf("实际圆环：%lf\n",j3);
    a=m3*(r3*r3+r4*r4)/2;
    printf("理论圆环：%lf\n",a);
    scanf("%lf%lf",&b6,&b5);
    j4=m1*r1*(g-r1*b6)/(b6-b5);
    j5=j4-j0;
    printf("实际圆盘：%lf\n",j5);
    j6=m2*r4*r4/2;
    printf("理论圆盘：%lf\n",j6);
    scanf("%lf%lf",&b8,&b7);
    j6=m1*r1*(g-r1*b8)/(b8-b7);
    j6=j6-j0;
    j7=j6+m4*0.05*0.05;
    j8=j6+m4*0.075*0.075;
    printf("原位0：%lf\n",j6);
    printf("原位1：%lf\n",j7);
    printf("原位2：%lf\n",j8);




}