＜計算0＞　　u[0,T]=f[T], u[1,T]=-d/dT u[1,T]

Clrea[X, a, ramda] ; (*** ** ** ** ** **     計算パラメ - タ         ************)

L = 1 * 10^-6 ;             &nbs ... p;              (* 拡散時間 *)

S2 = 10^6 ;              &n ... sp;   係数an の計算         ************)

a[n_] := 2 * ramda[n]/(ramda[n] - Sin[ramda[n]] * Cos[ramda[n]]) * NIntegrate[(F[X] + N0 * ( ... nbsp;カルシウムイオン濃度 NCa[X, T] の計算         ************)

NCa[X_, T_] = N0 * (1 - L * S2/(1 + L * S2) * X) + Sum[a[n] * Sin[ramda[n] * X] * Exp[-ramda ... p; 右端でのカルシウムイオンの流速          ************)

flux[X_, T_] = D[NCa[X, T], X] ;

fluxCa[T_] = flux[1, T] ; (*** ** ** ** ** **     計算結果の表示         ************)

Plot[{NCa[X, 0], NCa[X, 0.001], NCa[X, 0.01], NCa[X, 0.1], NCa[X, 1], NCa[X, 10], NCa[X, 100], NCa[X, 1000]}, {X, 0, 1}, PlotRange -> All] ;

Plot[-fluxCa[T], {T, 0, 10}, PlotRange -> All] ;

Null

[Graphics:../HTMLFiles/index_11.gif]

[Graphics:../HTMLFiles/index_12.gif]

tau

1/1000

Plot3D[NCa[X, T], {X, 0, 1}, {T, 0, 10}, PlotRange -> All] ;

[Graphics:../HTMLFiles/index_16.gif]

Created by Mathematica  (July 18, 2007) Valid XHTML 1.1!