问题及代码:
根据sin的泰勒展开式编写程序,求出sin(π/2)与sin(56°)的值。
(x取弧度值,而非角度值)
#include <iostream>
#include <iomanip>
using namespace std;
const double pi=3.1415926;
double mysin(double x);
int main( )
{
cout<<"sin(π/2)的值为"<<setiosflags(ios::fixed)<<setprecision(6)<<mysin(pi/2.0)<<endl;
cout<<"sin(56°)的值为"<<setiosflags(ios::fixed)<<setprecision(6)<<mysin(pi*56.0/180)<<endl;
return 0;
}
//下面定义mysin函数,求sin值
double mysin(double x)
{
double myabs(double y);
double Xn,sum=x,a=x;
int f=1,b=1;
for (;myabs(Xn)>0.00001;)
{
a=-a*x*x;
f=f*(b+1)*(b+2);
b+=2;
Xn=a/f;
sum=Xn+sum;
}
return sum;
}
//下面定义myabs函数,求绝对值
double myabs(double y)
{
if (y<0)
y=-y;
return y;
}
运行结果:
学习心得:
写得有点艰难…总是想依赖贺老 囧~
a,f的值不必每次都累乘 b+=2也是之前的遗忘点