1.题目:计算的函数值。计算结果取7位有效数字。
算法:数值积分方法(复化牛顿-科特斯公式):
程序:#include <>
#include <>
float func(float x)
int n,i;
float t,h,t1;
t=x/2*(1+exp(-x*x/2));
for (n=1;n<=500;n++)
for (i=0,h=0;i<=n-1;i++)t1=(2*i+1)*x/(2*n);
h=h+exp(-t1*t1/2)*x/n;
t=(t+h)/2;
return t;
int main()
int k;
float f;
for (k=1;k<=30;k++)
f=func(0.1*k)/sqrt(2*3.14)+0.5;
printf("x=%f,f(x)=%f",0.1*k,f);
printf("");
结果:分析:根据结果可知,该函数是一个单调增加的函数。
2.题目:人造地球卫星的轨道可视为平面上的椭圆,地心位于椭圆的一个焦点处。
已知一颗人造地球卫星近地点距地球表面439 km,远地点距地球表面 2384 km,地球半径为6371 km。求该卫星的轨道长度。
算法:利用平面曲线的弧长公式, ds=
程序:#include <>
#include <>
#define n 1000000
double a,b;
int main()
double func(double x);
double x=0,dx,g=0,a,c,c,dy;
double d1=6371+2384,d2=6371+439;
int i;
a=(d1+d2)/2;
c=(d1-d2)/2;
a=a*a;
c=c*c;
b=a-c;
dx=a/n;
for(i=0;i
g=4*g;
printf("%1f",l);
return 0;
double func(double x)double y;
y=x*x/a;
y=b*(1-y);
y=sqrt(y);
return y;
结果:48707.3
通过题目给出的条件求出半长轴、半短轴长,再利用公式ds=3. 利用切线法求方程的近似解(取n=2)算法:过(,f())作为曲线y=f(x)的切线,交于x轴于点。
重复以上过程n次,得到切线法求近似解的公式。
程序:#include <>
int main()
float x0=2,f,f;
int i;
for(i=0;i<3;i++)
printf("%4.3f",x0);
结果:2.095
分析:取,利用公式得到方程的近似解x=2.095。
高数上机作业
高等数学实验作业。1 两个重要极限 计算下列函数的函数值并画出图形,观察两个重要极限值 利用matlab进行函数作图 输入命令 x 10 0.01 10 y1 sin x x figure,plot x,y1 x 0 0.01 100 y2 1 x 1.x figure,plot x,y2 利用计算...
高数上机作业
数学实验报告 实验人员 院 系学号 姓名 实验地点 计算机中心机房。实验一。实验名称 观察f x cos x的各阶泰勒展开的图形。实验目的 利用mathematica计算函数的各阶泰勒多项式,并通过绘制曲线图形,来进一步掌握泰勒展开与函数逼近的思想。实验程序 1 固定x0 0,观察阶数n的影响。输入...
高数上机作业
1 两个重要极限 计算下列函数的函数值并画出图形,观察两个重要极限值 1 程序 x 10 0.01 10 y sin x x plot x,y x 0 0.01 1 y sin x x plot x,y x属于 10到10图像为 x属于0到1的图像为 故x趋近于0时,f x 趋近于1.2 程序 x ...