include
include
include
include
define pi 3.14159265359
lo*lo*p-2*mu,freq=2.25*1e6,wavelength=(long double)lo/freq,dh=(long double)wavelength/
30.0,dt=(long double)dh/(lo*1.5);
(1000*dh));
(p*dh),lambdaplus2mudtbydh=(lambda+2*mu)*dt/dh,lambdadtbydh=lambda*dt/dh,dtmubydh=dt*mu/
dh;
double**U,long double**V){
for(int k=0,l=0;k<=yno-1 && l<=yno;k++,l++){
U[i+1][l]+=dtbyrhodh*(X[i+1][l+1]-X[i+1][l]+Z[i+1][l]-
Z[i][l]);
[k+1]-Y[j][k+1]);
}
double**U,long double**V){
for(int k=0,l=0;k<=yno-1 && l<=yno;k++,l++){
U[i+1][k])+lambdadtbydh*(V[i+1][k+1]-V[i][k+1]);
V[i][k+1])+lambdadtbydh*(U[i+1][k+1]-U[i+1][k]);
U[j][l]);
int main(){
clock_t start,end;
long double time_taken;
start=clock();
long double **X,**Y,**U,**V,**Z;int n=1;
X=Make2DDoubleArray(xno+2,yno+2);
Y=Make2DDoubleArray(xno+2,yno+2);
Z=Make2DDoubleArray(xno+1,yno+1);
U=Make2DDoubleArray(xno+2,yno+2);
V=Make2DDoubleArray(xno+2,yno+2);
for (n=1;n<=timesteps;n++){
}
end=clock();
time_taken=(long double)(end-start)/CLOCKS_PER_SEC;
printf("Time elapsed is %Lf\nGRID Size:%Lf*%Lf\nTime Steps
Taken:%d\n",time_taken,(xno),floor(yno),n);
return 0;
}
