密立根油滴实验科大奥瑞计算代码

静态法：

#include <iostream>
#include <cmath>
using namespace std;
int main(){
    double q, t, r, u, x, x1, x2,n1,e1,e0=1.602;
    double p1 = 981,g = 9.80,l = 2.00e-3,b = 8.23e-3,p = 0.101e6,d = 5.00e-3,n = 1.83e-5;
    cout<<"先输入时间,再输入电压:"<<endl;
    cin>>t>>u;
    r = sqrt(9 * n*l / (2 * p1*g*t));
    x1 = ((n*l) / (t*(1 + b / (p*r))));
    x2 = x1*x1*x1;
    x = sqrt(x2);
    q = 18 * 3.14*x*d / (sqrt(2 * p1*g)*u);
    n1=q/(1.6e-19);
    if(n1-int(n1)>=0.5)    n1=int(n1)+1;
    else if(n1-int(n1)<0.5||n1-int(n1)>=0)    n1=int(n1);
    else cout<<"error";
    e1=q/n1;
    cout<<"油滴所带电量q="<<q<<endl;
    cout<<"实验数据所得电荷量e="<<e1<<endl;