Hola que tal !!.
Estoy haciendo una simulación a partículas en un gas. La verdad no sé mucho de programación y en laboratorio nos dieron el código del programa. Después de mucho ya lo logré compilar sin embargo no logro ejecutarlo, me sale un error de "Segmentation fault"; segun me dicen que es error del programa pero no encuentro nada.
Ahi va el código
-------------------------------------
//simulacion de gases
#include<stdio.h>
#include<math.h>
int N;
float lambda=1.5, energia=0.0,energiaold,n,sig;
float T,A,rx[26],ry[26],L;
void inicio(void);
void montecarlo(void);
float ran(void);
void salida(void);
void energia1(void);
void periodicas(int);
int traslape(int);
main ()
{
inicio ();
montecarlo();
salida();
energia1();
}
void inicio(void)
{
int i,m,j;
float rho,sig=1.0;
N=25;
rho=0.1;
A=(float)N*sig*sig/rho;
L=sqrt(A);
//generar las coordenadas
m=1;
for(i=1;j<=5;++i){
for(j=1;j<=5;++j){
rx[m]=(0.01+(L-0.02/4)*(j-1));
ry[m]=(0.01+(L-0.02/4)*(i-1));
++m;
}}
}
void montecarlo(void)
{
int i,j=1;
float desp=0.4;
energia1();
energiaold=energia;
for(j=1;j<1000;j++)
i=random()%25+1;
do{
rx[i]=rx[i]+desp*(ran()-0.5);
ry[i]=ry[i]+desp*(ran()-0.5);
periodicas(i);
traslape(i);
energia1();
}
while (energia>=energiaold);
}
float ran(void)
{
float r;
r=(float)(random(1000))/1000.0;
return(r);
}
void salida(void){
}
int traslape(int i){
int t=0,k,n,sig;
float dx,dy,d2;
for(k=1; k<=n; ++k)
if (k!=i){
dx=rx[i]-rx[k];
dy=ry[i]-ry[k];
d2=dx*dx+dy*dy;
if(d2<sig*sig) t=1;}
return(t);}
void configuracion(void){
int i;
for(i=1;i<=25;++i) printf("%d %f\t %f\n",i,rx[i],ry[i]);
}
void energia1(void){
int i,j,n;
float rr;
for(i=1;i<=n-1;++i);
for(j=i+1;j<=n;++j){
//distancia entre moleculas.
rr=(rx[i]-rx[j])*(rx[i]-rx[j])+(ry[i]-ry[j])*(ry[i]-ry[j]);
if(rr<=lambda*lambda) energia+=-1;
}}
//subrutinas periodicas.
void periodicas(int i){
if(rx[i]>L) rx[i]-=L;
else if(rx[i]<0.0) rx[i]+=L;
if(ry[i]>L) ry[i]-=L;
else if(ry[i]<0.0) rx[i]+=L;
}
--------------------------
Si alguien sabe como lo puedo arreglar sería muy bueno.
Gracias!!
