印刷用メモ

#include<stdio.h>
#include<stdlib.h>
#include<math.h>

int main(void)
{
FILE *Z;
double a,b,c,d,e,n,x,u,sum;

printf("start");
scanf("%lf",&a);
printf("end");
scanf("%lf",&b);
printf("division");
scanf("%lf",&n);

d=(b-a)/n;


Z=fopen("integral2.txt","w");
for(x=-3.0;x<=0;x=x+d){
sum=0;
for(u=a;u<=b;u=u+d){

e=0.5*(pow(u,0.5)/(exp(-x)*exp(u)-1.0)+pow(u+d,0.5)/(exp(-x)*exp(u+d)-1.0))*d;

sum=sum+e;

}

fprintf(Z,"%lf %lf\n",x,sum);


}
fclose(Z);

return 0;

}
~