解決しました。

解決しました。

sqr(z)関数で変換する前の、電荷位置として仮の位置(da,db)を想定してsqr(z)変換
し、更に全体を45度回転させた後の電荷位置が(a,b)になる、と言う条件式からを
(da,db)を求めて解決しました。
また、電気力線がダブって描かれたので、-pi<=v<=piと修正しました。
修正後のプログラムは下記のとおりです。

---直角平板の場合の等電位線と電気力線の作図----
option arithmetic complex
set Window -1,10,-1,10
draw grid
LET j=sqr(-1)
LET a=2 !電荷位置x座標の設定
LET b=2 !電荷位置y座標の設定
LET theta=2*(ATN(b/a)-PI/4)
LET dl=a*a/COS(theta/2+PI/4)/COS(theta/2+PI/4)
LET da=dl*COS(theta)
LET db=dl*SIN(theta)
for u=0 to 10 step 0.1 !等電位値（パラメータ）の設定
for v=-pi to pi step pi/300
LET w=u+j*v
LET z=(1+EXP(w))/(1-EXP(w))*(-da)+j*db
LET q=SQR(z)*EXP(j*PI/4)
plot lines:Re(q),Im(q);
next V
plot lines
next U

set line color "red"
for v=-pi to pi step pi/30 !等電気力線値（パラメータ）の設定
for u=0 to 10 step 0.02
LET w=u+j*v
LET z=(1+EXP(w))/(1-EXP(w))*(-da)+j*db
LET q=SQR(z)*EXP(j*PI/4)
plot lines:Re(q),Im(q);
next U
plot lines
next V
end

Re: 続き	返事を書くノートメニュー
島村1243 <bjllmpcujp> 2007/11/11 10:18:12 ** この記事は3回修正されてます
解決しました。 sqr(z)関数で変換する前の、電荷位置として仮の位置(da,db)を想定してsqr(z)変換し、更に全体を45度回転させた後の電荷位置が(a,b)になる、と言う条件式からを (da,db)を求めて解決しました。また、電気力線がダブって描かれたので、-pi<=v<=piと修正しました。修正後のプログラムは下記のとおりです。 ---直角平板の場合の等電位線と電気力線の作図---- option arithmetic complex set Window -1,10,-1,10 draw grid LET j=sqr(-1) LET a=2 !電荷位置x座標の設定 LET b=2 !電荷位置y座標の設定 LET theta=2(ATN(b/a)-PI/4) LET dl=aa/COS(theta/2+PI/4)/COS(theta/2+PI/4) LET da=dlCOS(theta) LET db=dlSIN(theta) for u=0 to 10 step 0.1 !等電位値（パラメータ）の設定 for v=-pi to pi step pi/300 LET w=u+jv LET z=(1+EXP(w))/(1-EXP(w))(-da)+jdb LET q=SQR(z)EXP(jPI/4) plot lines:Re(q),Im(q); next V plot lines next U set line color "red" for v=-pi to pi step pi/30 !等電気力線値（パラメータ）の設定 for u=0 to 10 step 0.02 LET w=u+jv LET z=(1+EXP(w))/(1-EXP(w))(-da)+jdb LET q=SQR(z)EXP(jPI/4) plot lines:Re(q),Im(q); next U plot lines next V end