d*z 发帖数: 150 | 1
try
k=3.584548 when t=0
k=-5.647731 when t=1
k=3.631065 when t=2
k=3.939163 when t=3
#include
#include
#define MIN_K_VALUE 1.989
double A(double k)
{
return sqrt(2.2 * k *k -3.95353);
}
double B(double k)
{
return sqrt(k*k-3.95353);
}
double Delta(double t,double k)
{
double a,b,da,db;
double v,w;
double t1,t2;
double c1,c2;
a=A(k);
b=B(k);
da=2.2*k/a;
db=k/b;
t1=tan((3.16-t)*b);
t2=tan(t*a);
c1=cos((3.16-t)*b);
c2=cos(t*a);
v=a*t1-b*t2;
w=da*t1-db*t2;
w+=a*(3.16-t)*db/(c1*c1 |
|