Autor
|
Tema: Problema codigo 0x4463c4 (Leído 1,531 veces)
|
oscarariza
Desconectado
Mensajes: 3
|
#include<iostream.h> #include<fstream.h> #include<math.h> int main() { double h,f,pi,bm,ac,A,B,Ma,densidad,Vmolar,R,T,P,P2,Q,D,th,Z1,Z2,Z3,v,u,xr1,xr2,xr3,xi2,xi3,g,y,mayor,sumacomposiciones; int i,j,n,k; R=10.737; pi=3.141592; cout<<"Presion:"; cin>>P; cout<<endl; cout<<"Temperatura:"; cin>>T; cout<<endl; fstream entrada,salida,salida2; entrada.open("datos.txt",ios::in); if(!entrada) { cout<<"no encontrado"<<endl; system("pause"); return 0; } entrada>>n; cout<<"filas: "<<n<<endl<<endl; entrada>>k; cout<<"columnas: "<<k<<endl<<endl; double Dat[n][k]; cout<<endl<<"***************Datos******************"<<endl<<endl; for(i=0;i<n;i=i+1) { for(j=0;j<k;j=j+1) { entrada>>Dat [j]; cout<<Dat[j]<<" "; } cout<<endl; } entrada.close();
sumacomposiciones=0.0; for(i=0;i<n;i=i+1) { sumacomposiciones=sumacomposiciones+Dat}
salida.open("calculos.txt",ios::out); salida2.open("resultados.txt",ios::out);
if(sumacomposiciones==1.0) { cout<<endl<<"***************Vector a***************"<<endl<<endl; salida<<endl<<"***************Vector a***************"<<endl<<endl;
double a;
for(i=0;i<n;i=i+1) { a=0.45724*((pow((R*Dat[2]),2))/Dat[3]); cout<<a<<" "; salida<<a<<" "; }
cout<<endl; salida<<endl;
cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl;
cout<<endl<<"***************Vector b***************"<<endl<<endl; salida<<endl<<"***************Vector b***************"<<endl<<endl;
double b;
for(i=0;i<n;i=i+1) { b=0.07780*((R*Dat[2])/Dat[3]); cout<<b<<" "; salida<<b<<" "; }
cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl;
cout<<endl<<"***************Vector m***************"<<endl<<endl; salida<<endl<<"***************Vector m***************"<<endl<<endl; double m;
for(i=0;i<n;i=i+1) { m=0.3746+1.5423*Dat[4]-0.2699*pow(Dat[4],2); cout<<m<<" "; salida<<m<<" "; }
cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl<<"***************Vector alfa***************"<<endl<<endl; salida<<endl<<"***************Vector alfa***************"<<endl<<endl; double alfa;
for(i=0;i<n;i=i+1) { alfa=pow((1+m*(1-pow((T/Dat[2]),0.5))),2); cout<<alfa<<" "; salida<<alfa<<" "; }
cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl; salida<<endl;
bm=0.0; for(i=0;i<n;i=i+1) { bm=bm+Dat}
cout<<"bm=";cout<<bm<<endl<<endl; salida<<"bm=";salida<<bm<<endl<<endl; ac=0.0; for(i=0;i<n;i=i+1) { for(j=0;j<n;j=j+1) { ac=ac+Dat- *Dat[j][0]*pow(a*a[j]*alfa*alfa[j],0.5);
} }
cout<<"ac=";cout<<ac<<endl<<endl; salida<<"ac=";salida<<ac<<endl<<endl;
A=(ac*P)/(pow(R*T,2)); cout<<"A=";cout<<A<<endl<<endl;
B=(bm*P)/(R*T); cout<<"B=";cout<<B<<endl<<endl;
P2=(3*(A-3*pow(B,2)-2*B)-pow((B-1),2))/3; cout<<"P2=";cout<<P2<<endl<<endl;
Q=(2*pow(B-1,3)-9*(B-1)*(A-3*pow(B,2)-2*B)+27*-1*(A*B-pow(B,2)-pow(B,3)))/27; cout<<"Q=";cout<<Q<<endl<<endl;
D=pow(P2/3,3)+pow(Q/2,2); cout<<"D=";cout<<D<<endl<<endl;
f=-Q/2; g=pow(D,0.5); h=fabs(f+g); u=pow(h,1.0/3.0); y=fabs(f-g); v=pow(y,1.0/3.0);
Ma=0.0; for(i=0;i<n;i=i+1) { Ma=Ma+(Dat} cout<<"Ma="<<cout<<Ma<<endl<<endl;
salida2<<P<<endl<<endl; salida2<<T<<endl<<endl;
if(D<=0.0) { th=acos(-Q/(2.0*pow(pow(fabs(P2),3.0)/27.0,0.5))); Z1=(2.0*pow(fabs(P2)/3.0,0.5)*cos(th/3.0))-((B-1.0)/3.0); Z2=(-2.0*pow(fabs(P2)/3.0,0.5)*cos(th+pi/3.0))-((B-1.0)/3.0); Z3=(-2.0*pow(fabs(P2)/3.0,0.5)*cos(th-pi/3.0))-((B-1.0)/3.0); cout<<"Z1=";cout<<Z1<<endl<<endl; cout<<"Z2=";cout<<Z2<<endl<<endl; cout<<"Z3=";cout<<Z3<<endl<<endl; if(Z3>Z1 and Z3>Z2) { mayor=Z3; } else { if(Z1>Z2 and Z1>Z3) { mayor=Z1; } } if(Z2>Z3 and Z2>Z1) { mayor=Z2; } cout<<"Zg="<<cout<<mayor<<endl<<endl; salida2<<"Zg="<<cout<<mayor<<endl<<endl; densidad=(P*Ma)/(mayor*R*T); Vmolar=(mayor*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<endl<<endl; salida2<<"densidad="<<salida2<<densidad<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<endl<<endl; salida2<<"Volumen molar="<<salida2<<Vmolar<<endl<<endl; } else { xr1=(u+v)-((B-1.0)/3.0); xr2=-1.0*((u+v)/2.0)-((B-1.0)/3.0); xr3=-1.0*((u+v)/2.0)-((B-1.0)/3.0); xi2=(pow(3.0,0.5)*(u-v))/2.0; xi3=-1*(pow(3.0,0.5)*(u-v))/2.0; cout<<"Z1="<<cout<<xr1<<endl<<endl; cout<<"Z2="<<cout<<xr2<<" + "<<xi2<<" i "<<endl<<endl; cout<<"Z3="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; if(xr1>xr2 and xr1>xr3) { mayor=xr1; cout<<"Zg="<<cout<<xr1<<endl<<endl; salida2<<"Zg="<<salida2<<xr1<<endl<<endl; densidad=(P*Ma)/(xr1*R*T); Vmolar=(xr1*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<endl<<endl; salida2<<"densidad="<<salida2<<densidad<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<endl<<endl; salida2<<"Volumen molar="<<salida2<<Vmolar<<endl<<endl; } else { if(xr2>xr1 and xr2>xr3) { mayor=xr2; cout<<"Zg="<<cout<<xr2<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"Zg="<<salida2<<xr2<<" + "<<xi2<<" i "<<endl<<endl; densidad=(P*Ma)/(xr2*R*T); Vmolar=(xr2*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"densidad="<<cout<<densidad<<" + "<<xi2<<" i "<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi2<<" i "<<endl<<endl; } } if(xr3>xr2 and xr3>xr1) { mayor=xr3; cout<<"Zg="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"Zg="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; densidad=(P*Ma)/(xr3*R*T); Vmolar=(xr3*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"densidad="<<cout<<densidad<<" + "<<xi3<<" i "<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi3<<" i "<<endl<<endl; } } } else { cout<<endl; cout<<"Suma de las composiciones es diferente a 1, no es valido:"<<endl; } salida.close(); salida2.close(); cout<<endl; system ("pause"); return 0; }
Lo que esta en rojo es lo que presenta problemas, al dar la respuesta en el programa sale con 0x4463c4
por ejemplo
Ma=0x4463c420.313 lo que esta en azul es lo que debe salir
|
|
|
En línea
|
|
|
|
|
oscarariza
Desconectado
Mensajes: 3
|
Disculpe es que soy nuevo en esto
|
|
« Última modificación: 17 Abril 2012, 19:16 pm por oscarariza »
|
En línea
|
|
|
|
xiruko
Desconectado
Mensajes: 438
|
#include<iostream.h> #include<fstream.h> #include<math.h> int main() { double h,f,pi,bm,ac,A,B,Ma,densidad,Vmolar,R,T,P,P2,Q,D,th,Z1,Z2,Z3,v,u,xr1,xr2,xr3,xi2,xi3,g,y,mayor,sumacomposiciones; int i,j,n,k; R=10.737; pi=3.141592; cout<<"Presion:"; cin>>P; cout<<endl; cout<<"Temperatura:"; cin>>T; cout<<endl; fstream entrada,salida,salida2; entrada.open("datos.txt",ios::in); entrada>>n; cout<<"filas: "<<n<<endl<<endl; entrada>>k; cout<<"columnas: "<<k<<endl<<endl; double Dat[n][k]; cout<<endl<<"***************Datos******************"<<endl<<endl; for(i=0;i<n;i=i+1) { for(j=0;j<k;j=j+1) { entrada>>Dat[j]; cout<<Dat[j]<<" "; } cout<<endl; } entrada.close(); sumacomposiciones=0.0; for(i=0;i<n;i=i+1) { sumacomposiciones=sumacomposiciones+Dat ; } salida.open("calculos.txt",ios::out); salida2.open("resultados.txt",ios::out); if(sumacomposiciones==1.0) { cout<<endl<<"***************Vector a***************"<<endl<<endl; salida<<endl<<"***************Vector a***************"<<endl<<endl; double a[i]; for(i=0;i<n;i=i+1) { a[i]=0.45724*((pow((R*Dat[2]),2))/Dat[3]); cout<<a[i]<<" "; salida<<a[i]<<" "; } cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl<<"***************Vector b***************"<<endl<<endl; salida<<endl<<"***************Vector b***************"<<endl<<endl; double b[i]; for(i=0;i<n;i=i+1) { b[i]=0.07780*((R*Dat[2])/Dat[3]); cout<<b[i]<<" "; salida<<b[i]<<" "; } cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl<<"***************Vector m***************"<<endl<<endl; salida<<endl<<"***************Vector m***************"<<endl<<endl; double m[i]; for(i=0;i<n;i=i+1) { m[i]=0.3746+1.5423*Dat[4]-0.2699*pow(Dat[4],2); cout<<m[i]<<" "; salida<<m[i]<<" "; } cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl<<"***************Vector alfa***************"<<endl<<endl; salida<<endl<<"***************Vector alfa***************"<<endl<<endl; double alfa[i]; for(i=0;i<n;i=i+1) { alfa=pow((1+m*(1-pow((T/Dat[2]),0.5))),2); cout<<alfa[i]<<" "; salida<<alfa[i]<<" "; } cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl; salida<<endl; bm=0.0; for(i=0;i<n;i=i+1) { bm=bm+Dat[i][0]*b[i]; } cout<<"bm=";cout<<bm<<endl<<endl; salida<<"bm=";salida<<bm<<endl<<endl; ac=0.0; for(i=0;i<n;i=i+1) { for(j=0;j<n;j=j+1) { ac=ac+Dat[i][0]*Dat[j][0]*pow(a[i]*a[j]*alfa[i]*alfa[j],0.5); } } cout<<"ac=";cout<<ac<<endl<<endl; salida<<"ac=";salida<<ac<<endl<<endl; A=(ac*P)/(pow(R*T,2)); cout<<"A=";cout<<A<<endl<<endl; B=(bm*P)/(R*T); cout<<"B=";cout<<B<<endl<<endl; P2=(3*(A-3*pow(B,2)-2*B)-pow((B-1),2))/3; cout<<"P2=";cout<<P2<<endl<<endl; Q=(2*pow(B-1,3)-9*(B-1)*(A-3*pow(B,2)-2*B)+27*-1*(A*B-pow(B,2)-pow(B,3)))/27; cout<<"Q=";cout<<Q<<endl<<endl; D=pow(P2/3,3)+pow(Q/2,2); cout<<"D=";cout<<D<<endl<<endl; f=-Q/2; g=pow(D,0.5); h=fabs(f+g); u=pow(h,1.0/3.0); y=fabs(f-g); v=pow(y,1.0/3.0); Ma=0.0; for(i=0;i<n;i=i+1) { Ma=Ma+(Dat[i][0]*Dat[i][1]); } cout<<"Ma="<<cout<<Ma<<endl<<endl; salida2<<P<<endl<<endl; salida2<<T<<endl<<endl; if(D<=0.0) { th=acos(-Q/(2.0*pow(pow(fabs(P2),3.0)/27.0,0.5))); Z1=(2.0*pow(fabs(P2)/3.0,0.5)*cos(th/3.0))-((B-1.0)/3.0); Z2=(-2.0*pow(fabs(P2)/3.0,0.5)*cos(th+pi/3.0))-((B-1.0)/3.0); Z3=(-2.0*pow(fabs(P2)/3.0,0.5)*cos(th-pi/3.0))-((B-1.0)/3.0); cout<<"Z1=";cout<<Z1<<endl<<endl; cout<<"Z2=";cout<<Z2<<endl<<endl; cout<<"Z3=";cout<<Z3<<endl<<endl; if(Z3>Z1 and Z3>Z2) { mayor=Z3; } else { if(Z1>Z2 and Z1>Z3) { mayor=Z1; } } if(Z2>Z3 and Z2>Z1) { mayor=Z2; } cout<<"Zg="<<cout<<mayor<<endl<<endl; salida2<<"Zg="<<cout<<mayor<<endl<<endl; densidad=(P*Ma)/(mayor*R*T); Vmolar=(mayor*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<endl<<endl; salida2<<"densidad="<<salida2<<densidad<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<endl<<endl; salida2<<"Volumen molar="<<salida2<<Vmolar<<endl<<endl; } else { xr1=(u+v)-((B-1.0)/3.0); xr2=-1.0*((u+v)/2.0)-((B-1.0)/3.0); xr3=-1.0*((u+v)/2.0)-((B-1.0)/3.0); xi2=(pow(3.0,0.5)*(u-v))/2.0; xi3=-1*(pow(3.0,0.5)*(u-v))/2.0; cout<<"Z1="<<cout<<xr1<<endl<<endl; cout<<"Z2="<<cout<<xr2<<" + "<<xi2<<" i "<<endl<<endl; cout<<"Z3="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; if(xr1>xr2 and xr1>xr3) { mayor=xr1; cout<<"Zg="<<cout<<xr1<<endl<<endl; salida2<<"Zg="<<salida2<<xr1<<endl<<endl; densidad=(P*Ma)/(xr1*R*T); Vmolar=(xr1*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<endl<<endl; salida2<<"densidad="<<salida2<<densidad<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<endl<<endl; salida2<<"Volumen molar="<<salida2<<Vmolar<<endl<<endl; } else { if(xr2>xr1 and xr2>xr3) { mayor=xr2; cout<<"Zg="<<cout<<xr2<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"Zg="<<salida2<<xr2<<" + "<<xi2<<" i "<<endl<<endl; densidad=(P*Ma)/(xr2*R*T); Vmolar=(xr2*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"densidad="<<cout<<densidad<<" + "<<xi2<<" i "<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi2<<" i "<<endl<<endl; } } if(xr3>xr2 and xr3>xr1) { mayor=xr3; cout<<"Zg="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"Zg="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; densidad=(P*Ma)/(xr3*R*T); Vmolar=(xr3*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"densidad="<<cout<<densidad<<" + "<<xi3<<" i "<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi3<<" i "<<endl<<endl; } } } else { cout<<endl; cout<<"Suma de las composiciones es diferente a 1, no es valido:"<<endl; } salida.close(); salida2.close(); cout<<endl; system ("pause"); return 0; }
Desde Ma para abajo los resultados me salen con 0x4463c4 por ejemplo Ma=0x4463c420.313 lo que esta en azul es lo que debe salir pd. tan dificil es usar las etiquetas de codigo?
|
|
|
En línea
|
|
|
|
oscarariza
Desconectado
Mensajes: 3
|
#include<iostream.h> #include<fstream.h> #include<math.h> int main() { double h,f,pi,bm,ac,A,B,Ma,densidad,Vmolar,R,T,P,P2,Q,D,th,Z1,Z2,Z3,v,u,xr1,xr2,xr3,xi2,xi3,g,y,mayor,sumacomposiciones; int i,j,n,k; R=10.737; pi=3.141592; cout<<"Presion:"; cin>>P; cout<<endl; cout<<"Temperatura:"; cin>>T; cout<<endl; fstream entrada,salida,salida2; entrada.open("datos.txt",ios::in); entrada>>n; cout<<"filas: "<<n<<endl<<endl; entrada>>k; cout<<"columnas: "<<k<<endl<<endl; double Dat[n][k]; cout<<endl<<"***************Datos******************"<<endl<<endl; for(i=0;i<n;i=i+1) { for(j=0;j<k;j=j+1) { entrada>>Dat[i][j]; cout<<Dat[i][j]<<" "; } cout<<endl; } entrada.close(); sumacomposiciones=0.0; for(i=0;i<n;i=i+1) { sumacomposiciones=sumacomposiciones+Dat[i][0]; } salida.open("calculos.txt",ios::out); salida2.open("resultados.txt",ios::out); if(sumacomposiciones==1.0) { cout<<endl<<"***************Vector a***************"<<endl<<endl; salida<<endl<<"***************Vector a***************"<<endl<<endl; double a[i]; for(i=0;i<n;i=i+1) { a[i]=0.45724*((pow((R*Dat[i][2]),2))/Dat[i][3]); cout<<a[i]<<" "; salida<<a[i]<<" "; } cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl<<"***************Vector b***************"<<endl<<endl; salida<<endl<<"***************Vector b***************"<<endl<<endl; double b[i]; for(i=0;i<n;i=i+1) { b[i]=0.07780*((R*Dat[i][2])/Dat[i][3]); cout<<b[i]<<" "; salida<<b[i]<<" "; } cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl<<"***************Vector m***************"<<endl<<endl; salida<<endl<<"***************Vector m***************"<<endl<<endl; double m[i]; for(i=0;i<n;i=i+1) { m[i]=0.3746+1.5423*Dat[i][4]-0.2699*pow(Dat[i][4],2); cout<<m[i]<<" "; salida<<m[i]<<" "; } cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl<<"***************Vector alfa***************"<<endl<<endl; salida<<endl<<"***************Vector alfa***************"<<endl<<endl; double alfa[i]; for(i=0;i<n;i=i+1) { alfa=pow((1+m[i]*(1-pow((T/Dat[i][2]),0.5))),2); cout<<alfa[i]<<" "; salida<<alfa[i]<<" "; } cout<<endl; salida<<endl; cout<<endl<<"***************************************"<<endl; salida<<endl<<"***************************************"<<endl; cout<<endl; salida<<endl; bm=0.0; for(i=0;i<n;i=i+1) { bm=bm+Dat[i][0]*b[i]; } cout<<"bm=";cout<<bm<<endl<<endl; salida<<"bm=";salida<<bm<<endl<<endl; ac=0.0; for(i=0;i<n;i=i+1) { for(j=0;j<n;j=j+1) { ac=ac+Dat[i][0]*Dat[j][0]*pow(a[i]*a[j]*alfa[i]*alfa[j],0.5); } } cout<<"ac=";cout<<ac<<endl<<endl; salida<<"ac=";salida<<ac<<endl<<endl; A=(ac*P)/(pow(R*T,2)); cout<<"A=";cout<<A<<endl<<endl; B=(bm*P)/(R*T); cout<<"B=";cout<<B<<endl<<endl; P2=(3*(A-3*pow(B,2)-2*B)-pow((B-1),2))/3; cout<<"P2=";cout<<P2<<endl<<endl; Q=(2*pow(B-1,3)-9*(B-1)*(A-3*pow(B,2)-2*B)+27*-1*(A*B-pow(B,2)-pow(B,3)))/27; cout<<"Q=";cout<<Q<<endl<<endl; D=pow(P2/3,3)+pow(Q/2,2); cout<<"D=";cout<<D<<endl<<endl; f=-Q/2; g=pow(D,0.5); h=fabs(f+g); u=pow(h,1.0/3.0); y=fabs(f-g); v=pow(y,1.0/3.0); Ma=0.0; for(i=0;i<n;i=i+1) { Ma=Ma+(Dat[i][0]*Dat[i][1]); } cout<<"Ma="<<cout<<Ma<<endl<<endl; salida2<<P<<endl<<endl; salida2<<T<<endl<<endl; if(D<=0.0) { th=acos(-Q/(2.0*pow(pow(fabs(P2),3.0)/27.0,0.5))); Z1=(2.0*pow(fabs(P2)/3.0,0.5)*cos(th/3.0))-((B-1.0)/3.0); Z2=(-2.0*pow(fabs(P2)/3.0,0.5)*cos(th+pi/3.0))-((B-1.0)/3.0); Z3=(-2.0*pow(fabs(P2)/3.0,0.5)*cos(th-pi/3.0))-((B-1.0)/3.0); cout<<"Z1="<<cout<<Z1<<endl<<endl; cout<<"Z2="<<cout<<Z2<<endl<<endl; cout<<"Z3="<<cout<<Z3<<endl<<endl; if(Z3>Z1 and Z3>Z2) { mayor=Z3; } else { if(Z1>Z2 and Z1>Z3) { mayor=Z1; } } if(Z2>Z3 and Z2>Z1) { mayor=Z2; } cout<<"Zg="<<cout<<mayor<<endl<<endl; salida2<<"Zg="<<cout<<mayor<<endl<<endl; densidad=(P*Ma)/(mayor*R*T); Vmolar=(mayor*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<endl<<endl; salida2<<"densidad="<<salida2<<densidad<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<endl<<endl; salida2<<"Volumen molar="<<salida2<<Vmolar<<endl<<endl; } else { xr1=(u+v)-((B-1.0)/3.0); xr2=-1.0*((u+v)/2.0)-((B-1.0)/3.0); xr3=-1.0*((u+v)/2.0)-((B-1.0)/3.0); xi2=(pow(3.0,0.5)*(u-v))/2.0; xi3=-1*(pow(3.0,0.5)*(u-v))/2.0; cout<<"Z1="<<cout<<xr1<<endl<<endl; cout<<"Z2="<<cout<<xr2<<" + "<<xi2<<" i "<<endl<<endl; cout<<"Z3="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; if(xr1>xr2 and xr1>xr3) { mayor=xr1; cout<<"Zg="<<cout<<xr1<<endl<<endl; salida2<<"Zg="<<salida2<<xr1<<endl<<endl; densidad=(P*Ma)/(xr1*R*T); Vmolar=(xr1*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<endl<<endl; salida2<<"densidad="<<salida2<<densidad<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<endl<<endl; salida2<<"Volumen molar="<<salida2<<Vmolar<<endl<<endl; } else { if(xr2>xr1 and xr2>xr3) { mayor=xr2; cout<<"Zg="<<cout<<xr2<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"Zg="<<salida2<<xr2<<" + "<<xi2<<" i "<<endl<<endl; densidad=(P*Ma)/(xr2*R*T); Vmolar=(xr2*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"densidad="<<cout<<densidad<<" + "<<xi2<<" i "<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi2<<" i "<<endl<<endl; salida2<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi2<<" i "<<endl<<endl; } } if(xr3>xr2 and xr3>xr1) { mayor=xr3; cout<<"Zg="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"Zg="<<cout<<xr3<<" + "<<xi3<<" i "<<endl<<endl; densidad=(P*Ma)/(xr3*R*T); Vmolar=(xr3*R*T)/P; cout<<endl<<"***************Resultados***************"<<endl<<endl; salida2<<endl<<"***************Resultados***************"<<endl<<endl; cout<<"densidad="<<cout<<densidad<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"densidad="<<cout<<densidad<<" + "<<xi3<<" i "<<endl<<endl; cout<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi3<<" i "<<endl<<endl; salida2<<"Volumen molar="<<cout<<Vmolar<<" + "<<xi3<<" i "<<endl<<endl; } } } else { cout<<endl; cout<<"Suma de las composiciones es diferente a 1, no es valido:"<<endl; } salida.close(); salida2.close(); cout<<endl; system ("pause"); return 0; }
Desde Ma para abajo en el programa las respuestas me salen acompañadas del 0x4463c4
|
|
« Última modificación: 17 Abril 2012, 19:21 pm por oscarariza »
|
En línea
|
|
|
|
|
|