2017-04-04 15 views
2

누구도 /는 무엇을 알고 있습니까? 예 : Mathematica에서 V[n_][i_/;i<=imax,0];=V[n][i,0]=0?/는 무엇입니까? Mathematica에서 의미가 있습니까?

이것은 유한 차분 방법을 사용하여 직사각형 금속 골짜기의 잠재적 인 잠재력을 끌어내는 프로그램의 일부입니다.

개선 된 기능이지만 다음은 작동하지 않습니다.

V[0][i_, j_] := 0; 
tol = 0.0025; 
imax = 20; jmax = 20; 
V[n_][i_, 0] := V[n][i, 0] = 0; 
V[n_][i_, jmax] := V[n][i, jmax] = 100*Sin[0.05*Pi*i]; 
V[n_][0, j_] := V[n][0, j] = 0; 
V[n_][imax, j_] := V[n][imax, j] = 0; 
V[n_][i_, j_] := V[n][i, j] = (1/4)*(V[n - 1][i + 1, j] + 
       V[n - 1][i - 1, j] + V[n - 1][i, j + 1] + V[n - 1][i, j - 1]); 
For[m = 1, Max[Table[Abs[V[m][i, j] - V[m - 1][i, j]], {i, 1, imax - 1}, {j, 1, jmax - 1}]] >= tol, m++; 
    If[m > 20, Break[]]]; 
    Print[m]; 
    Show[ListPlot3D[Table[V[m][i, j], {j, 0, 20}, {i, 0, 20}]], BoxRatios -> {1, 1, 0.85}, PlotRange -> {0, 1}, Axes -> True, AxesLabel -> {"x", "y", "V"}] 
+1

https://reference.wolfram.com/language/ref/Condition.html – agentp

+0

@agentp 고마워, 도움이됩니다. 그리고'ListPlot3D'가 위에서 언급 한 코드에서 적절한 함수라고 생각합니까? –

+0

@Pengyy 귀하의 도움에 감사드립니다! –

답변

1

저는 이미 Wolfram의 참고 문헌을 통해 요점을 얻었습니다. 여기에 웹 https://reference.wolfram.com/language/ref/Condition.html이 있습니다. @agentp 도움에 감사드립니다.