如何用c语言编写这个微积分程序(第二题)求代码

如何用c语言编写这个微积分程序(第二题)求代码,第1张

#include <stdio.h>

#include <math.h>

double integ(double sub,double sup,double f(double))

double f(double x)

int main(void) 

{

printf("%E\n",integ(5,7,f))

return 0

}

double integ(double sub,double sup,double f(double))

{

double i,sum=0

for(i=subi<=supi+=1e-6)

{

sum+=f(i)*1e-6

}

return sum

}

double f(double x)

{

return exp(5*x)*pow(x,4)-3*x

}

辛普森法

#include <stdio.h>

#include <math.h>

#define F(X) (4.0/(1+X*X))

static float EPS = 1.0E-14

static int COUNT=1

static double a=0.0, b=1.0

static double M[32],T[32],S[32],C[32],R[32],E[32]

double GETM(int K)

{ unsigned long j,n=1

double x=0,y=0,step=0

for(j=0j<Kj++) n*=2

step = (b-a)/n

x = a+step/2

for(j=0j<nj++){y+=F(x)x+=step}

return (y*step)

}

Simpson(float EPS)

{ int k=0

T[0]= (b-a)*(F(a)+F(b))/2.0

S[0]=T[0]

E[0]=1.0E10

for(k=0k<20k++)

{ M[k] = GETM(k)

T[k+1]=(M[k]+T[k])/2.0

S[k+1]=(4.0*T[k+1]-T[k])/3.0

COUNT++

E[k+1]=fabs(S[k+1]-S[k])

if(k<3)continue

if(E[k+1]<EPS)break

}

return

}

ShawResult()

{ int k

system("cls")

printf("\n K M[K] T[K] S[K]")

printf(" E[K]")

printf("\n-------------------------------------------")

printf("-----------------------------------")

for(k=0k<COUNTk++)

printf("\n%2d %20.15lf%20.15lf%20.15lf%12.4e",k,M[k],T[k],S[k],E[k])

printf("\n--------------------------------------------")

printf("-----------------------------------")

getch()

system("cls")

return

}

SaveResult()

{ int k

FILE * fp

fp=fopen("Simpson.htm","w")

if(!fp) return

fprintf(fp,"<html><head></head>")

fprintf(fp,"<body bgcolor = #006699 text = #FFff00>")

fprintf(fp,"<pre><font size=\"6\">")

fprintf(fp,"\n K M[K] T[K] S[K]")

fprintf(fp," E[K]")

fprintf(fp,"\n----------------------------------------------")

fprintf(fp,"----------------------------------")

for(k=0k<COUNTk++)

fprintf(fp,"\n%2d %20.15f%20.15f%20.15f%12.4le",k,M[k],T[k],S[k],E[k])

fprintf(fp,"\n----------------------------------------------")

fprintf(fp,"----------------------------------")

fprintf(fp,"</pre></font></body>")

fclose(fp)

return

}

main()

{ Simpson(EPS)

ShawResult()

SaveResult()

}

1、积分工具(App)

英文名:CalculusTools数学软件,用来计算导数、定积分、弧长、泰勒级数,绘制函数对应的图,包括极坐标和参数以及斜率场。此外,还包括积分表和自定义键盘。

2、matlab7.0

是MathWorks出品一款的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。

是一款积分,微分,函数等都可以进行计算的软件。

扩展资料:

关于matlab7.0的主要功能

1、交互式工具可以按迭代的方式探查、设计及求解问题

2、此高级语言可用于技术计算

3、此开发环境可对代码、文件和数据进行管理

4、各种工具可用于构建自定义的图形用户界面

5、各种函数可将基于MATLAB的算法与外部应用程序和语言(如C、C++、Fortran、Java、COM以及MicrosoftExcel)集成

6、数学函数可用于线性代数、微积分、统计、傅立叶分析、筛选、优化以及数值积分等

7、二维和三维图形函数可用于可视化数据。

参考资料来源:百度百科——微积分工具

参考资料来源:百度百科——matlab7.0


欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/yw/11693353.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-18
下一篇 2023-05-18

发表评论

登录后才能评论

评论列表(0条)

保存