home *** CD-ROM | disk | FTP | other *** search
/ Piper's Pit BBS/FTP: ibm 0020 - 0029 / ibm0020-0029 / ibm0028.tar / ibm0028 / MCAD254.ZIP / APPL / DIFFN.MCD < prev    next >
Encoding:
Text File  |  1988-07-14  |  1.3 KB  |  73 lines
  1. .MCD 20000 0
  2. .CMD PLOTFORMAT logs=0,0 subdivs=1,1 size=5,15 type=l
  3. .CMD FORMAT  rd=d ct=10 im=i et=3 zt=15 pr=3 mass length time charge
  4. .CMD SET ORIGIN 0
  5. .CMD SET TOL 0.001000
  6. .CMD MARGIN 0
  7. .CMD LINELENGTH 78
  8. .CMD SET PRNCOLWIDTH 8
  9. .CMD SET PRNPRECISION 4
  10. .TXT 0 42 1 39 
  11. a1,38,49,37
  12. Copyright (c) 1988 by MathSoft, Inc.
  13. .TXT 1 -41 1 12 
  14. a1,11,78,10
  15. DIFFUSION
  16. .TXT 0 80 1 13 
  17. a1,12,78,11
  18. /equations
  19. .TXT 2 -80 1 25 
  20. a1,24,77,23
  21. diffusion coefficient:
  22. .EQN 0 25 1 9 
  23. c~.01
  24. .TXT 2 -25 1 29 
  25. a1,28,77,27
  26. distribution at time zero:
  27. .EQN 0 30 1 32 
  28. f(x)~if(|(x-.5)<.125,1,0)
  29. .TXT 2 -30 1 10 
  30. a1,9,77,8
  31. x step:
  32. .EQN 0 9 1 9 
  33. k~.05
  34. .TXT 0 14 1 13 
  35. a1,12,77,11
  36. time step:
  37. .EQN 0 12 1 8 
  38. h~.1
  39. .TXT 0 11 1 24 
  40. a1,23,52,22
  41. number of time steps:
  42. .EQN 0 24 1 8 
  43. T~30
  44. .EQN 0 10 3 14 
  45. X~floor(1/k)
  46. .EQN 0 18 1 11 
  47. j~0;X
  48. .EQN 0 17 1 11 
  49. i~0;T
  50. .EQN 0 16 2 10 
  51. x[j~j*k
  52. .EQN 1 -33 1 15 
  53. J~1;X-1
  54. .EQN 0 17 1 11 
  55. I~1;T
  56. .EQN 1 -63 4 11 
  57. k^2/(2*c)=?
  58. .EQN 1 46 2 10 
  59. u[(i,0)~0
  60. .EQN 0 17 2 10 
  61. u[(i,X)~0
  62. .EQN 0 16 2 14 
  63. u[(0,j)~f(x[j)
  64. .TXT 1 -131 1 51 
  65. a1,50,52,49
  66. For stability, h should be less than or equal to
  67. .EQN 1 80 4 9 
  68. ▐~(h*c)/k^2
  69. .EQN 1 18 2 52 
  70. u[(I,J)~u[(I-1,J)+▐*(u[(I-1,J+1)-2*u[(I-1,J)+u[(I-1,J-1))
  71. .EQN 1 -88 10 50 
  72. 1.1&0&u[(0,j),u[(10,j),u[(30,j){1,1,8,25,l}@1&0&x[j
  73.