分享几个自己对编程成果（C语言版）

shumabaolong · 发表于 2010-9-11 22:00:32

马上注册，结识高手，享用更多资源，轻松玩转三维网社区。

您需要登录才可以下载或查看，没有帐号？注册

x

好久没来论坛发言了，来分享几个我编写对程序，对工科学生也许会有点启发。

用迭代法求解方程超越方程的解
#include "math.h"
main()
{
double x1,x2;
x1=1;
x2=2;
while(x2-x1>0.000001||x1-x2>0.000001)
{
x1=x2;
x2=log10(x1+2);
printf("x2=%lf\n",x2);
}
printf("result is %lf\n",x2);
}

对分法解方程
double f(double x)
{
double y;
y=x*x*x+2*x-5;
return y;
}
main()
{
double a=1,b=2;
double c,d,r;
int i;
for(i=0;i<20;i++)
{
c=(a+b)/2;
d=f(c);
if(d==0)
  {
  r=c;
  break;
  }
else
  {
  if(d>0)
  b=c;
  else
  a=c;
  }
}
r=c;
printf("The result of ""x^3+2*x-5=0"" is %lf\n",r);
}

一下程序段将以0为起点0.05为步长以2*PI为终点计算各个自变量的正弦值，并将计算结果以点的坐标的形式（x,y）打印到文件中，方便绘图软件处理。使用这种方法，即使不懂C语言的图形操作也可以绘制出高质量的曲线来。

#define PI 3.1415926
#include<stdlib.h>
#include<stdio.h>
#include<math.h>
main()
{
FILE * pf;
double x,y;
char xx,yy;
x=0.0;
if((pf=fopen("d:\\file1.txt","w"))==NULL)
{
printf("cannot open this file\n");
exit(1);
}
while(x<=(PI*2.0))
{
y=sin(x);
fprintf(pf,"%lf,%lf\n",x,y);
x+=0.05;
}
fprintf(pf,"%lf,%lf\n",2*PI,sin(2*PI));
fclose(pf);
}

		自动登录	找回密码
密码			注册

通知

[原创] 分享几个自己对编程成果（C语言版）

马上注册，结识高手，享用更多资源，轻松玩转三维网社区。

浏览过的版块