计算机作业2008级.doc

在化学反应CO(g)+1/2O2(g)=CO2(g)中1克分子CO2分解百分比依赖于方程 (P/K2 – 1)x3+3x – 2=0, 其中P是CO2在大气中的压强, K是一依赖于温度的平衡常数(在2800K时, K=1.648), P=1标准大气压时, 求x (0.7586) #include math.h #include stdio.h main( ) { float a=0.0, b=1.0; float x; x=half(a,b); printf(x=%f\n,x); } float half(float a,float b) { float p=1,k=1.648,x,f,fa; fa=(p/(k*k)-1.0)*a*a*a+3.0*a-2.0; do { x=(a+b)/2.0; f=(p/(k*k)-1.0)*x*x*x+3.0*x-2.0; if(f*fa0) a=x; else b=x; } while(fabs(f)1.0e-5); return(x); } 程序运行结果为：x=0.758601 假定氨气符合范德华方程(P+n2a/v2)(v-nb)=nRT, 已知a=4.19大气压.升2/摩尔2, b=0.0373升/摩尔2, n=1, P=1大气压, R=0.08206升.大气压/摩尔.K, T=423.2K, 求氨气的1克分子体积. (34.64升) #include math.h #include stdio.h main() { printf(v=%f,newt()); } float f(float x ) { float p=1.0,n=1.0,a=4.19,b=0.0373,R=0.08206,T=423.2; return (p+n*n*a/(x*x))*(x-n*b)-n*R*T; } float df(float x) { float p=1.0,n=1.0,a=4.19,b=0.0373,R=0.08206,T=423.2; return p-n*n*a/(x*x)+2*n*n*n*a*b/(x*x*x); } float newt() { float eps=1.0e-5; float p=1.0,n=1.0,a=4.19,b=0.0373,R=0.08206,T=423.2; float v=n*R*T/p; float v0; do{ v0=v; v=v0-f(v0)/df(v0); }while(fabs(v0-v)eps); return(v); } 程序运行结果为：v= 34.644276 3. 将氯甲烷（CH3Cl）、氯乙烯（C2H5Cl）、氰化氢（HCN）与氨(NH3)混合，然后对混合物进行元素分析，各个元素所占质量分数如下：C：18.71%，H：11.38%，Cl：35.68%。请计算各化合物占混合物总量的质量分数。已知各元素的相对原子质量为：C：12.01，H：1.008，N：14.01，Cl：35.45 (编程求解和使用Excel规划求解均可，建议使用Excel方法做) #include stdlib.h #include math.h #include stdio.h int sgauss(n,a,b) int n; double a[],b[]; { int *js,key,k,i,j,is,p,q; double d,t; js=malloc(n * sizeof(int)); /*开辟用于记忆列交换信息的动态空间*/ key=1; /*置非奇异标志*/ for(k=0;k=n-2,k++) { d=0.0; for(i=k;i=n-1;i++) /*全选主元*/ for(j=k;j=n-1;j++) { t=fabs(a[i * n+j]; if(td) {d=t;js[k]=j;is=i;} /*记忆行、列交换信息*/ } if(d+1.0==1.0) key=0; /*置奇异标志*/ else { if(js[k]!=k) f