- 締切済み
分数の分母ゼロ問題
分数の計算で分母がゼロになると当然発散するわけですが、その回避の方法として以下のようなものを考えています。 x=1.0 y=0.0 a=x/y ! a=∞ ゼロ割計算が一旦発生する b=y/x ! b=0.0 if(y.eq.0.) then c=a ! しかしここでaを採用しないのでプログラムには関係なし。 else c=b ! こちらが選択される。 endif というようなものです。一旦ゼロ割が生じるけれどもそれを採用しないのでトラブル回避という格好です。このような方法は問題ないと考えられるでしょうか。 もちろん、このような簡単なプログラムではありません。複雑なアルゴリズムの中の１つの局面でこのような処理をするということです。それとも絶対このような処理は行うべきではないのでしょうか。 実験すべきかと思いますが、本当の方は複雑なので実験結果があてはまるかどうかも分からないものですので。 よろしくお願いします。
- skmsk19410
- お礼率54% (276/507)
- その他（プログラミング・開発）
- 回答数6
- ありがとう数4
- みんなの回答 （6）
- 専門家の回答
みんなの回答
- notnot
- ベストアンサー率47% (4854/10269)
- Tacosan
- ベストアンサー率23% (3656/15482)
- ultraCS
- ベストアンサー率44% (3956/8947)
- notnot
- ベストアンサー率47% (4854/10269)
数学の話じゃなくて、プログラミングの話なら、実数のゼロ除算でどうなるかは言語によって違いますので、特定言語に限定して話さないと意味がないです。 ・例外が発生。例外処理を記述できる ・例外が発生。例外処理を記述できない ・結果が浮動小数点の無限大になり、処理続行 などなど 例外が発生するケースでは、a=x/y の時点で発生するので、その後のif文の実行がどうなるかは関係ないです。
お礼
回答、有り難うございます。 ご指摘の内容はコンパイラ、処理系、言語など諸々の条件によっていろいろ変わるということかと思います。 このプログラムではゼロ割しても後の処理でその問題を回避する処置がなされるため体勢に影響はないという解釈が成り立ちます。ゼロ割した変数の内容をその後参照しないからですね。コンパイラがそのように判断するというようなことではないかと思います。 一方、その逆の解釈、すなわち、ゼロ割が生じるようなプログラムはそれがどのような局面でさえ、絶対に動作しない、どんな言語、どんな処理系でもこれが生じたらそのシステムは異常終了するのである、というものです。この場合、このアルゴリズムは何が何でも排除しなければならないわけです。 数千行程度のプログラムのバグフィックスをおこなっており、最初の質問に示したような趣旨のアルゴリズムを見つけました。この箇所こそその元凶であると特定したかったのですが、必ずしもそういうことではないということのようですね。もちろん、私がこのアルゴリズムを積極的に採用したいと思っている訳ではありません。私が期待した回答は後者、すなわち絶対回避せよ、なのです。
- gyrocompas
- ベストアンサー率23% (24/104)
- buriburi3
- ベストアンサー率44% (353/792)
関連するQ&A
- 初期条件を代入すると分母が0になる初期値問題
2階線形常微分方程式 (x^2 y')' = -x^2 を解くと、両辺を x で積分して x^2 y' = -x^3/3 + C 両辺を x^2 で割って y' = -x/3 + C/x^2 ... (b) 両辺を x で積分して y = -x^2/6 - C/x + D ... (a) となります。 ここで初期条件 y(0) = 1, y'(0) = 0 だと y = -x^2/6 + 1 になるらしいのですが、(a) で y(0) = 1 とおいて D について解こうと思っても、 y(0) = -C/0 + D = 1 と分母が0になってしまい計算できません。 (b) y' = -x/3 + C/x^2 の場合も同様で、 y'(0) = C/0 = 0 となってしまいます。 このような場合はどうやって積分定数 C, D を求めればよいのでしょうか？ 解き方が間違っていますか？ ご教示ください。よろしくお願いします。
- ベストアンサー
- 数学・算数
- BASIC言語の問題（初心者）
こんばんは(>_<) 学校からの課題で非常に困っています!!!! BASIC言語の仕組みは読めても、どうしても組み立てが出来ないのです。。。 *１* ２の平方根を出力しなさい（二分法による平方根の計算） ただしδ＝0.0001とする *２* ２の平方根を出力しなさい（ニュートン・ラプソン法による平方根の計算） ただしδ＝0.0001とする ヒントとして、それぞれのアルゴリズム（？）があるので紹介します(>_<) *１*（xの平方根を精度δで求める） a←0 b←x while b-a < δ do c←(a+b)/2 if C*C > x then b←c else a←c endif done return a *２*（xの平方根を精度δで求める） y←x while │x-y*y│ < 2yδ do y←（y*y+x）/2y done return y 自分は文系で、パソコンの知識まったくないので、 どうかお助け下さいo(;△;)o涙
- ベストアンサー
- その他（プログラミング・開発）
- 分数の足し算をさせるプログラムが分かりません。どなたか分かりませんか？
分数の足し算をさせるプログラムが分かりません。 Ｃ言語の問題で分数の足し算までは一応できるんですが、答えがでたときに整数で出すやり方と約分して表す方法が分かりません。 どなたか知恵を貸してくれませんか？ ユーザから4つの整数を入力し、はじめに入力された2個の整数と後に入力された2個の整数を分数と考え、その分数の和を表示するプログラムを作成せよ。 例えば、「3」「4」「5」「6」と入力されたときは、3/4 + 5/6を計算する。 そのプログラム内では分数の和を計算する関数を作成する。 さらに、 約分を行う関数を 再帰呼び出しを利用して作成する。 void yakubun(int *a1, int *a2) 例えば、以下の場合1/2と表示される。 int i=10,j=20; yakubun(&i,&j); printf(“%d / %d”, i, j); ちなみにここまでできました↓ #include<stdio.h> void bunsu_tasizan(int a1,int a2,int b1,int b2, int *c1,int *c2 ) { *c1=(a1*b2)+(b1*a2); *c2=(a2*b2); } int main() { int x1,x2,y1,y2,z1,z2; printf("整数を入力してください"); scanf("%d",&x1); scanf("%d",&x2); scanf("%d",&y1); scanf("%d",&y2); if(x2==0||y2==0||x2==0&&y2==0) printf("0以外を入力してください"); else{ bunsu_tasizan(x1,x2,y1,y2,&z1,&z2); printf("%d/%d",z1,z2);} return (0); }
- ベストアンサー
- C・C++・C#
- 中学校で習う比の問題で分からないところがあります。
（X－２Y）：（３X－４Y）＝５：１のときX：Yを求めなさい。との問題なのですが、この場合は比の性質「A:B＝C：Dのとき、AD＝BC」で解けるのでしょうか。 それとも「A:B＝C：DのときAC＝BD」又「A:B＝C：DのときBC：AD」で解くのでしょうか？どれで解けばよいでしょうか？ 一応手元に回答はあるのですが、答えのみなのであまり詳しく載っていません。 回答はX：Y＝９：７です。 念のため自分で計算したのですが、数字が逆になってしまいます。 ※※Ａ※※※Ｂ※※※Ｃ※Ｄ※※ （X－２Y）：（３X－４Y）＝５：１ ↓の式は「A:B＝C：DのときBC：AD」を使いました。 ５（３X－４Y）＝１（X－２Y） １５X－２０Y＝X－２Y １５X－X＝－２Y＋２０Y １４X＝１８Y 両方を２で割る ７X＝９Y X：Y＝７：９ になりましたが、誤りだそうです。 回答は反対のX：Y＝９：７になっています。 質問なんですがX：Y＝７：９が→X：Y＝９：７にするには 何か計算したり、動かしたり、公式に入れたりしているの でしょうか？そのあたりを※詳しく詳しく教えてください。 （※は投稿するとどうも端に寄ってしまうので、わざと※を付けて寄らないようにしています。計算とはまったく関係ありませんので） ご指導よろしくお願いします。
- ベストアンサー
- 数学・算数
- UWSCのスクリプトについて
画像Dを認識したら画面上に画像Bがあっても 処理Bを実行しない為にはどうしたら良いでしょうか。 -------- while True BTN(LEFT,click,150,10,110) //処理A Sleep(2) ifb CHKIMG("A.bmp") x=G_IMG_X y=G_IMG_Y BTN(left,click,x,y,110) Sleep(1.5) endif Ifb ChkImg("B.bmp") //処理B BTN(LEFT,click,1120,90,110) Sleep(1) KBD(VK_return,DOWN,30) endif ifb CHKIMG("C.bmp") //処理C x=G_IMG_X y=G_IMG_Y BTN(left,click,x,y,110) Sleep(1.5) wend --------
- ベストアンサー
- その他（プログラミング・開発）
- 分数式での疑問
x^2-xy-cy+cx / x^2+bx-xy-by =x(x-y)+c(x-y) / x(x-y)+b(x-y) =(x-y)(x+c) / (x-y)(x+b) =x+c / x+b ↑の式が納得できませんでした。 2行目の x(x-y)+c(x-y) / x(x-y)+b(x-y) を、 x(x-y)で約分して ⇒c(x-y) / b(x-y) 更に(x-y)で約分して c / b と答えたのですが、どこで間違ったのでしょう。 x(x-y)+c(x-y) / x(x-y)+b(x-y) ↑ここで分母と分子を (x-y) で括らなければ答えが間違ってしまうのは何故ですか？
- 締切済み
- 数学・算数
- 数学（ベクトル）の問題
http://okwave.jp/qa/q8022847.html のNO.7の回答より、 さらに途中式を書いたのですが、 L^2 = m(t-n(s))^2-(a1^2+b1^2+c1^2 ) {(c1(z2-z1)+b1(y2-y1)+a1(x2-x1)+(a1a2+b1b2+c1c2)s)/(a1^2+b1^2+c1^2 )}^2+(a2^2+b2^2+c2^2 ){(s^2 )+2s{c2(z2-z1)+b2(y2-y1)+a2(x2-x1)}/{(a2^2+b2^2+c2^2 )} }+{(x2-x1)^2 }+{(y2-y1)^2 }+{(z2-z1)^2 } n(s)={c1(z2-z1)+b1(y2-y1)+a1(x2-x1)+(a1a2+b1b2+c1c2)s}/(a1^2+b1^2+c1^2 ) = m(t-n(s))^2-{c1(z2-z1)+b1(y2-y1)+a1(x2-x1)+(a1a2+b1b2+c1c2)s}^2/(a1^2+b1^2+c1^2 )+(a2^2+b2^2+c2^2 ){(s^2 )+2s{c2(z2-z1)+b2(y2-y1)+a2(x2-x1)}/{(a2^2+b2^2+c2^2 )} }+{(x2-x1)^2 }+{(y2-y1)^2 }+{(z2-z1)^2 } = m(t-n(s))^2-{c1(z2-z1)+b1(y2-y1)+a1(x2-x1)+(a1a2+b1b2+c1c2)s}^2/(a1^2+b1^2+c1^2 )+(a2^2+b2^2+c2^2 )(s^2 )+2s{c2(z2-z1)+b2(y2-y1)+a2(x2-x1)}+{(x2-x1)^2 }+{(y2-y1)^2 }+{(z2-z1)^2 } ところで {c1(z2-z1)+b1(y2-y1)+a1(x2-x1)+(a1a2+b1b2+c1c2)s}^2 = {c1(z2-z1)+b1(y2-y1)+a1(x2-x1)}^2 +2{c1(z2-z1)+b1(y2-y1)+a1(x2-x1)}(a1a2+b1b2+c1c2)s +{(a1a2+b1b2+c1c2)s}^2 = {c1(z2-z1)+b1(y2-y1)+a1(x2-x1)}^2 +2{c1(z2-z1)+b1(y2-y1)+a1(x2-x1)}(a1a2+b1b2+c1c2)s +(a1a2+b1b2+c1c2)^2 (*s)^2 これより、 L^2= m(t-n(s))^2+s^2 {(a2^2+b2^2+c2^2 )-(a1a2+b1b2+c1c2)^2/(a1^2+b1^2+c1^2 )} +s[2{(c2(z2-z1)+b2(y2-y1)+a2(x2-x1))-{c1(z2-z1)+b1(y2-y1)+a1(x2-x1)}(a1a2+b1b2+c1c2)s/(a1^2+b1^2+c1^2 )}] +{(x2-x1)^2 }+{(y2-y1)^2 }+{(z2-z1)^2 }-{c1(z2-z1)+b1(y2-y1)+a1(x2-x1)}^2/(a1^2+b1^2+c1^2 ) 簡単、 L^2 =m(t-n(s))^2+ps^2+p1s+p2 =m(t-n(s))^2+p(s^2+p1s/p)+p2 =m(t-n(s))^2+p(s^2+p1s/p+(p1/p)^2-(p1/p)^2 )+p2 =m(t-n(s))^2+p(s^2+p1s/p+(p1/p)^2 )-(p1)^2/p+p2 まで、計算したのですが（間違っていたら申し訳ありません）、 ここから、どのように q=-p1/2p が導出できるのかがわからないです。 （rは導出できました。） 数式だらけで分かりづらいと思いますが、計算ミスを指摘しつつ、導出過程も分かりやすくお願いします。
- 締切済み
- 数学・算数
- エジプトの分数問題。どうしても解きたい。
4abc-b-4c∈K , 4abc-b-c∈H , x*y∈T , x≠1、y≠1 Tは素数以外の数で可解集合とする。 ここで、P+k=4abc-b-c∈H とすると 1≦c=4m-k≦4 とする。 P+3=4ab(4m-3)-b-4m+3 , P∈K P+6=4ab(4m-6)-b-4m+6 , P∈K P+7=4ab(4m-7)-b-4m+7 , P∈K P+ 12=4ab(4m-12)-b-4m+12 , P∈K ここで、15ｎ+aの形の関数を考える。 すると １５ｎ＋１∈K、P+12 １５ｎ＋２∈K、P+12 １５ｎ＋３∈T、３ １５ｎ＋４∈K、P＋7 １５ｎ＋５∈T、５ １５ｎ＋６∈T、３ １５ｎ＋７∈K、P+7 １５ｎ＋８∈K、P+3 １５ｎ＋９∈T、３ １５ｎ＋１０∈T、５ １５ｎ＋１１∈H １５ｎ＋１２∈T、３ １５ｎ＋１３∈H １５ｎ＋１４∈H １５ｎ＋１５∈T、１５ ここで具体的に１５ｎ＋１＝ｐを計算してみましょう。 P+12=4abc-b-c=(4ac-1)*b-c b=n+1,a=2,c=2 P+12=15n+15-2=15n+13∈H P+12=4ab*4-b-4∈H P=4ab*4-b-4*4∈K このように計算していくと15n+aは解けることになります。 自信はないのですが。 また、わけわからないことを言って迷惑をかけています。 先に謝っておきます。どうもすみません。 本人はこんな簡単に解けるはずはないと思っています。
- ベストアンサー
- 数学・算数
お礼
回答有り難うございます。 私も同感です。ゼロ割問題のフォローが計算の後だったら絶対にダメなのでしょうか。私の希望はこのアルゴリズムの完全否定です。処理系や言語にかかわらず、可能性としての否定ではなく、完全否定です。