physics 6a final

using System;
class HelloWorld {
    static double pp(double x, double pow) {
        return x*x;
    }
    static double M = 0.198;
    static double l = 0.12;
    static double m = 0.00323;
    //static double teta1 = 0.13126;
    //static double teta2 = 0.1078;

    static double K(double R2, double R1, double T2, double T1) {
        return (( 8 * Math.PI * Math.PI * M * (R2 * R2 - R1 * R1))/(T2 * T2 - T1 * T1));
    }
    /*static double gffgv1(double R2, double R1, double T2, double T1) {
        return ((K(R2, R1, T2, T1) * T1 )/ (2.0 * Math.PI * m * l)) * teta1;
    }
    static double fgfgv2(double R2, double R1, double T2, double T1) {
        return ((K(R2, R1, T2, T1) * T2 )/ (2.0 * Math.PI * m * l)) * teta2;
    }*/

    static double v1(double kk, double T2, double T1, double teta1, double teta2) {
        return ((kk * T1 )/ (2.0 * Math.PI * m * l)) * teta1;
    }
    static double v2(double kk, double T2, double T1, double teta1, double teta2) {
        return ((kk * T2 )/ (2.0 * Math.PI * m * l)) * teta2;
    }
  static void Main() {

    double T1 = 0.954;
    double T2 = 1.574;
    double M = 0.198;
    double R1 = 0.041;
    double R2 = 0.091;
    double ddR1 = 0.001;
    double ddR2 = 0.001;
    double ddT1 = 0.073;
    double ddT2 = 0.166;

    double ddK = Math.Sqrt(
          pp((K(R2, R1, T2, T1) - K(R2 + ddR2, R1, T2, T1)),2) +
          pp((K(R2, R1, T2, T1) - K(R2, R1 + ddR1, T2, T1)),2) +
          pp((K(R2, R1, T2, T1) - K(R2, R1, T2 + ddT2, T1)),2) +
          pp((K(R2, R1, T2, T1) - K(R2, R1, T2, T1 + ddT1)),2)
    );
    /*
    Console.WriteLine(K(R2, R1, T2, T1));
    Console.WriteLine();
    Console.WriteLine(K(R2 + ddR2, R1, T2, T1));
    Console.WriteLine(K(R2, R1 + ddR1, T2, T1));
    Console.WriteLine(K(R2, R1, T2 + ddT2, T1));
    Console.WriteLine(K(R2, R1, T2, T1 + ddT1));
    Console.WriteLine();
    Console.WriteLine(bobble);
    Console.WriteLine();
    Console.WriteLine();
    Console.WriteLine();
    Console.WriteLine();*/

    double teta1 = 0.13126;
    double teta2 = 0.1078;
    double ddteta1 = 0.00451;
    double ddteta2 = 0.01057;
    double kk = K(R2, R1, T2, T1);
    double ddV2 = Math.Sqrt(
          pp((v2(kk, T2, T1, teta1, teta2) - v2(kk + ddK, T2, T1, teta1, teta2)),2) +
          pp((v2(kk, T2, T1, teta1, teta2) - v2(kk, T2 + ddT2, T1, teta1, teta2)),2) +
          pp((v2(kk, T2, T1, teta1, teta2) - v2(kk, T2, T1, teta1, teta2 + ddteta2)),2)
    );

    Console.WriteLine(v2(kk, T2, T1, teta1, teta2));
    Console.WriteLine();
    Console.WriteLine(v2(kk + ddK, T2, T1, teta1, teta2));
    Console.WriteLine(v2(kk, T2 + ddT2, T1, teta1, teta2));
    Console.WriteLine(v2(kk, T2, T1, teta1, teta2 + ddteta2));
    Console.WriteLine();
    Console.WriteLine(ddV2);
    /*
    double ddV2 = Math.Sqrt(
          pp((v2(kk, T2, T1) - v2(kk + ddK, T2, T1)),2) +
          pp((v2(kk, T2, T1) - v2(kk, T2 + ddT2, T1)),2) +
          pp((v2(kk, T2, T1) - v2(kk, T2, T1 + ddT1)),2)
    );

    Console.WriteLine(v2(kk, T2, T1));
    Console.WriteLine();
    Console.WriteLine(v2(kk + ddK, T2, T1));
    Console.WriteLine(v2(kk, T2 + ddT2, T1));
    Console.WriteLine(v2(kk, T2, T1 + ddT1));
    Console.WriteLine();
    Console.WriteLine(ddV2);*/
  }
}