use imsl
PARAMETER TS=3000, XS=100
REAL C(0:XS),CTR(0:XS)
REAL T,X, DT,DX, C0, CS,T0,TMAX,X0,XMAX ! C(T,X)
WRITE(*,*)"PLEASE INPUT TIME STEP(0-3000) "
READ(*,*)TIME
T0=0.0
TMAX=20*3600
DT=(TMAX-T0)/TS
X0=0
XMAX=5
DX=(XMAX-X0)/XS
D=2.0E-5
C0=0.0
CS=0.8
OPEN (1,FILE="C:\DIF.DAT")
DO I=1,TS
T=T0+DT*I
DO J=1,XS-1
X=X0+DX*J
TEMP1=X/(2.0*SQRT(1.0*D*T))
C(J)=CS*(1-ERF(TEMP1))
IF(I==TIME)WRITE(1,*) X,C(J)
END DO
END DO
CLOSE(1)
END
1) Start the Microsoft Developer Studio.
2) Tool File pull-down menu.
3) options
4) Click on the "Directories" tab.
5) From the Category pull-down menu, select "Lib files".
6) Add the directory of the IMSL Fortran 90 MP Library
7) Click on "OK" to exit the dialog and save the settings.