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
51
52
53
54
55
56
57
58
| #include<stdio.h>
#include<conio.h>
#include<math.h>
double fun(double S, double R, int m, int n, int k, double q)
{double f,qmn,qmk;
double kp, fm, fk;
fm=m;
fk=k;
kp=fm/fk;
qmn=pow(q,m*n);
qmk=pow(q,kp);
f=qmn*qmk+qmk-S*k*(qmk-1.)/R;
return f;
}
double funp(double S, double R, int m, int n, int k, double q)
{double ff,qmn1,qmk1;
double kp,fm,fk;
fm=m;
fk=k;
kp=fm/fk;
qmn1=pow(q,m*n-1.);
qmk1=pow(q,kp-1.);
ff=kp*m*n*qmn1*qmk1+kp*qmk1-kp*S*k*qmk1/R;
return ff;
}
int main ()
{
double S,SS,q0,q,eps,R,t,y,kk;
double kp,fm,fk; int k,m,n;
S=1512,182;
m=2;
R=30;
n=12;
k=12;
eps=0.001;
kk=0;
fm=m;
fk=k;
kp=fm/fk;
q0=2.;
do
{
kk++;
q=q0-fun(S,R,m,n,k,q0)/funp(S,R,m,n,k,q0);
SS=(pow(q,m*n)-1.);
SS=SS/(pow(q,kp)-1.)*R/fk*pow(q,kp);
printf(" SS=%lf \n",SS);
printf("kk=%lf", kk);
printf("q =%lf", q);
}
// t=fabs(S-SS);
while(t>=eps);
q0=q;
y=(q-1.)*fm;
printf("y=%lf", y);
getch();
} |