1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50
|
long double F(long double a,long double ap,long double b,long double bp,long double c,long double cp,long double d,long double dp)
{ long double r;
long double R0=1;
long double Ho=65000;
long double cv=3*powl(10,8);
/* r=d^2(ct)/ds^2 */
r=-((2./3)*R0*R0*powl((3*Ho)/(2*cv),4./3)*powl((a),1./3))*(powl(bp,2)+powl(b,2)*powl(cp,2)+powl(b,2)*powl(sinl(c),2)*powl(dp,2));
return r;
}
long double G(long double a,long double ap,long double b,long double bp,long double c,long double cp,long double d,long double dp)
{ long double r;
long double Ho=65000;
long double cv=3*powl(10,8);
/* r=d^2(r)/ds^2 */
r=-((4./(3*(a)))*ap*bp)+b*(powl(cp,2)+powl(sinl(c),2)*powl(dp,2));
return r;
}
long double H(long double a,long double ap,long double b,long double bp,long double c,long double cp,long double d,long double dp)
{ long double r;
long double Ho=65000;
long double cv=3*powl(10,8);
/* r=d^2(theta)/ds^2 */
r=-(2./b)*bp*cp-(4./(3*(a)))*ap*cp+(sinl(2*c)*powl(dp,2))/2;
return r;
}
long double I(long double a,long double ap,long double b,long double bp,long double c,long double cp,long double d,long double dp)
{ long double r;
long double Ho=65000;
long double cv=3*powl(10,8);
/* r=d^2(phi)/ds^2 */
r=-2*dp*((2./(3*(a)))*ap+(1./tanl(c))*cp+(1./b)*bp);
return r;
} |
Partager