双曲正弦函数计算编程

#include

#include

main()

{

double a=0,a1,a2,a3,a4,M1,Q1;

a1=cosh(a)*cos(a);

a2=0.5*(sin(a)*cosh(a)+cos(a)*sinh(a));

a3=sinh(a)*sin(a)*0.5;

a4=0.25*(sin(a)*cosh(a)-cos(a)*sinh(a));

printf("a=%.2f a1=%.4f a2=%.4f a3=%.4f a4=%.4f\n",a,a1,a2,a3,a4);

M1=38016*a1+29236*a2-86345*a3+100600*a4;

printf("M1=%.0f\n",M1);

Q1=9648*a1-28494*a2+40239.9*a3-50181*a4;

printf("Q1=%.0f\n",Q1);

double b=0.1,b1,b2,b3,b4,M2,Q2;

b1=cosh(b)*cos(b);

b2=0.5*(sin(b)*cosh(b)+cos(b)*sinh(b));

b3=sinh(b)*sin(b)*0.5;

b4=0.25*(sin(b)*cosh(b)-cos(b)*sinh(b));

printf("b=%.2f b1=%.4f b2=%.4f b3=%.4f b4=%.4f\n",b,b1,b2,b3,b4);

M2=38016*b1+29236*b2-86344.6*b3+100600*b4;

printf("M2=%.0f\n",M2);

Q2=9648*b1-28494*b2+40239.9*b3-50181*b4;

printf("Q2=%.0f\n",Q2);

double c=0.15,c1,c2,c3,c4,M3,Q3;

c1=cosh(c)*cos(c);

c2=0.5*(sin(c)*cosh(c)+cos(c)*sinh(c));

c3=sinh(c)*sin(c)*0.5;

c4=0.25*(sin(c)*cosh(c)-cos(c)*sinh(c));

printf("c=%.2f c1=%.4f c2=%.4f c3=%.4f c4=%.4f\n",c,c1,c2,c3,c4);

M3=38016*c1+29236*c2-86344.6*c3+100600*c4;

printf("M3=%.0f\n",M3);

Q3=9648*c1-28494*c2+40239.9*c3-50181*c4;

printf("Q3=%.0f\n",Q3);

double d=0.17,d1,d2,d3,d4,M4,Q4;

d1=cosh(d)*cos(d);

d2=0.5*(sin(d)*cosh(d)+cos(d)*sinh(d));

d3=sinh(d)*sin(d)*0.5;

d4=0.25*(sin(d)*cosh(d)-cos(d)*sinh(d));

printf("d=%.2f d1=%.4f d2=%.4f d3=%.4f d4=%.4f\n",d,d1,d2,d3,d4); M4=38016*d1+29236*d2-86344.6*d3+100600*d4;

printf("M4=%.0f\n",M4);

Q4=9648*d1-28494*d2+40239.9*d3-50181*d4;

printf("Q4=%.0f\n",Q4);

double e=0.2,e1,e2,e3,e4,M5,Q5;

e1=cosh(e)*cos(e);

e2=0.5*(sin(e)*cosh(e)+cos(e)*sinh(e));

e3=sinh(e)*sin(e)*0.5;

e4=0.25*(sin(e)*cosh(e)-cos(e)*sinh(e));

printf("e=%.2f e1=%.4f e2=%.4f e3=%.4f e4=%.4f\n",e,e1,e2,e3,e4); M5=38016*e1+29236*e2-86344.6*e3+100600*e4;

printf("M5=%.0f\n",M5);

Q5=9648*e1-28494*e2+40239.9*e3-50181*e4;

printf("Q5=%.0f\n",Q5);

double f=0.23,f1,f2,f3,f4,M6,Q6;

f1=cosh(f)*cos(f);

f2=0.5*(sin(f)*cosh(f)+cos(f)*sinh(f));

f3=sinh(f)*sin(f)*0.5;

f4=0.25*(sin(f)*cosh(f)-cos(f)*sinh(f));

printf("f=%.2f f1=%.4f f2=%.4f f3=%.4f f4=%.4f\n",f,f1,f2,f3,f4);

M6=38016*f1+29236*f2-86344.6*f3+100600*f4;

printf("M6=%.0f\n",M6);

Q6=9648*f1-28494*f2+40239.9*f3-50181*f4;

printf("Q6=%.0f\n",Q6);

double g=0.3,g1,g2,g3,g4,M7,Q7;

g1=cosh(g)*cos(g);

g2=0.5*(sin(g)*cosh(g)+cos(g)*sinh(g));

g3=sinh(g)*sin(g)*0.5;

g4=0.25*(sin(g)*cosh(g)-cos(g)*sinh(g));

printf("g=%.2f g1=%.4f g2=%.4f g3=%.4f g4=%.4f\n",g,g1,g2,g3,g4); M7=38016*g1+29236*g2-86344.6*g3+100600*g4;

printf("M7=%.0f\n",M7);

Q7=9648*g1-28494*g2+40239.9*g3-50181*g4;

printf("Q7=%.0f\n",Q7);

double h=0.33,h1,h2,h3,h4,M8,Q8;

h1=cosh(h)*cos(h);

h2=0.5*(sin(h)*cosh(h)+cos(h)*sinh(h));