工数上机作业

发布 2022-09-05 09:25:28 阅读 9773

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 ...