cpp
#include<stdio.h>
#include<math.h>
double c(double T);
double c(double T)
{return ( (1.59e-11)*pow(T,4) - (1.58e-8)*pow(T,3) + (3.49e-6)*pow(T,2) + (4.72e-4)*T - (9.34e-3) );}
main()
{
float k,T1,T2,h;
double T,sum_odd,sum_even,ans;
int i,n;
printf("\nEnter the value of k:");
scanf("%f",&k);
n=2*k;
printf("n is the even number = %d\n",n);
printf("T1= 330K\n");
printf("T2 = 300K\n");
T1=330;
T2=300;
h=(T2-T1) / (n);
printf("the value of h is %1.2f\n",h);
sum_odd=0;
for(i=2;i<=n;i+=2){sum_odd+=c(T1+(i-1)*h);}
sum_even=0;
for(i=2;i<n;i+=2){sum_even+=c(T1+i*h);}
ans=(h/3)*(c(T1)+c(T2)+(4.0*sum_odd)+(2.0*sum_even));
printf("\nThe definite integral from T1 to T2");
printf("\nOf c(T)d(T) is approx.%.4f\n",ans);
system("pause");}
*edit: Please use code tags in the future, thanks!
This post has been edited by Martyr2: 2 Apr, 2008 - 09:39 PM
2 Apr, 2008 - 08:41 AM
| 