約 17 件
https://w.atwiki.jp/resistance3/pages/39.html
91KrGR edquylyiucsx, [url=http //mgsdxqlpmclo.com/]mgsdxqlpmclo[/url], [link=http //skgxyfkezonv.com/]skgxyfkezonv[/link], http //deztmwsoyioe.com/
https://w.atwiki.jp/ecovip/pages/607.html
http //buy-benicar.science/#2912 how much does benicar cost,
https://w.atwiki.jp/editguide/pages/76.html
HiWD9Z a href="http //becakiktbvtf.com/" becakiktbvtf /a , [url=http //ftmxbuxodrgl.com/]ftmxbuxodrgl[/url], [link=http //uuvurbsledpl.com/]uuvurbsledpl[/link], http //rkdzmsaxgphq.com/
https://w.atwiki.jp/broadcast/pages/52.html
MIYUKI TORII 美人育成プロジェクト20080308 1/3 MIYUKI TORII 美人育成プロジェクト20080308 2/3 MIYUKI TORII 美人育成プロジェクト20080308 3/3 MIYUKI TORII 美人育成プロジェクト20080308 1/3 MIYUKI TORII ç¾äººè²æããã¸ã§ã¯ã20080308 1/3 Uploaded by toritoritorii MIYUKI TORII 美人育成プロジェクト20080308 2/3 MIYUKI TORII ç¾äººè²æããã¸ã§ã¯ã20080308 2/3 Uploaded by toritoritorii MIYUKI TORII 美人育成プロジェクト20080308 3/3 MIYUKI TORII ç¾äººè²æããã¸ã§ã¯ã20080308 3/3 Uploaded by toritoritorii
https://w.atwiki.jp/evan/pages/74.html
äºåãã±ã´ã¡ã³ã²ãªã²ã³æ°åå ´ç ç ´Uploaded by UDF 【コメント】 名前
https://w.atwiki.jp/broadcast/pages/34.html
鳥居みゆき(エンタの神様/捨てられた女) MIYUKI TORII エンタの神様 20080301 鳥居みゆき(エンタの神様/捨てられた女) MIYUKI TORII エンタの神様 20080301 MIYUKI TORII ã¨ã³ã¿ã®ç¥æ§ 20080301 Uploaded by toritoritorii
https://w.atwiki.jp/tohohopeacewalk/pages/193.html
思01_6522t "HN変更とチャンネル桜" http //www21.atwiki.jp/tohohopeacewalk/pages/168.html http //pipponan.fc2web.com/shikosakugo_01/shi-6522.htm -------------------------------------------------------------------------------- [6522]HN変更とチャンネル桜 dozugawa 09/7/9(木) 2 01 [6523]Re(1) HN変更とチャンネル桜 平社員 09/7/9(木) 20 18 [6524]Re(2) HN変更とチャンネル桜 dozugawa 09/7/9(木) 20 37 [6526]南京 (戦線後方記録映画) ピッポ 09/7/12(日) 20 24 [6530]Re(1) 南京 (戦線後方記録映画) 平社員 09/7/13(月) 22 26 [6534]Re(2) 南京 (戦線後方記録映画) ピッポ 09/7/14(火) 23 51 [6537]Re(3) 南京 (戦線後方記録映画) 平社員 09/7/15(水) 17 36 [6538]Re(4) 南京 (戦線後方記録映画) ピッポ 09/7/15(水) 19 35 [6539]Re(5) 南京 (戦線後方記録映画) 平社員 09/7/18(土) 18 46 [6536]Re(2) 南京 (戦線後方記録映画) 熊猫 09/7/15(水) 8 59 [6531]Re(1) HN変更とチャンネル桜 熊猫 09/7/14(火) 3 18 [6532]Re(2) HN変更とチャンネル桜 dozugawa 09/7/14(火) 21 49 [6535]Re(3) HN変更とチャンネル桜 熊猫 09/7/15(水) 8 32 -------------------------------------------------------------------------------- [6522]HN変更とチャンネル桜 ←back ↑menu ↑top forward→ dozugawa - 09/7/9(木) 2 01 - -------------------------------------------------------------------------------- 引用なし パスワード bykokudo改めdozugawaですよろしく。 チャンネル桜の歴史捏造家、水島総氏が本来のチャンネさえピンチなのに又エセ映画「南京の真実」四部作を作ろうとしている様です。 資金が無いなら止めればいいのに支援金を集めています。この中でプロパガンダとして有名な「南京~戦線後方記録映画」が使用される様です。 なんとかならんものですか?。 ttp //www.nankinnoshinjitsu.com/kaisetsu.html#interview [管理人修正回:文内URL] 441 hits -------------------------------------------------------------------------------- [6523]Re(1) HN変更とチャンネル桜 ←back ↑menu ↑top forward→ 平社員 - 09/7/9(木) 20 18 - -------------------------------------------------------------------------------- 引用なし パスワード ▼dozugawaさん: bykokudo改めdozugawaですよろしく。 チャンネル桜の歴史捏造家、水島総氏が本来のチャンネさえピンチなのに又エセ映画「南京の真実」四部作を作ろうとしている様です。 資金が無いなら止めればいいのに支援金を集めています。この中でプロパガンダとして有名な「南京~戦線後方記録映画」が使用される様です。 なんとかならんものですか?。 http //www.nankinnoshinjitsu.com/kaisetsu.html#interview 直接関係はない話なのですが「戦線後方記録映画」はどのような内容なのでしょうか? どのへんがウソ(というより事実に触れていない)のでしょう? 313 hits -------------------------------------------------------------------------------- [6524]Re(2) HN変更とチャンネル桜 ←back ↑menu ↑top forward→ dozugawa - 09/7/9(木) 20 37 - -------------------------------------------------------------------------------- 引用なし パスワード ▼平社員さん: ▼dozugawaさん: bykokudo改めdozugawaですよろしく。 チャンネル桜の歴史捏造家、水島総氏が本来のチャンネさえピンチなのに又エセ映画「南京の真実」四部作を作ろうとしている様です。 資金が無いなら止めればいいのに支援金を集めています。この中でプロパガンダとして有名な「南京~戦線後方記録映画」が使用される様です。 なんとかならんものですか?。 ttp //www.nankinnoshinjitsu.com/kaisetsu.html#interview 直接関係はない話なのですが「戦線後方記録映画」はどのような内容なのでしょうか? どのへんがウソ(というより事実に触れていない)のでしょう? >http //ja.wikipedia.org/wiki/å—京_(æˆ¦ç·šå¾Œæ–¹è¨˜éŒ²æ˜ ç”») 松井石根大将の南京入城から始まる国策(?)映画のようです、負傷兵を看護する看護婦、日本軍による戦没者のための合同慰霊祭等死体は全く出てきません。 >http //ilovenippon.jugem.jp/?eid=559こんなサイトでもオススメの様です。 [管理人修正1回:文内URL] 309 hits -------------------------------------------------------------------------------- [6526]南京 (戦線後方記録映画) ←back ↑menu ↑top forward→ ピッポ - 09/7/12(日) 20 24 - -------------------------------------------------------------------------------- 引用なし パスワード ▼dozugawaさん: >http //ja.wikipedia.org/wiki/å—京_(æˆ¦ç·šå¾Œæ–¹è¨˜éŒ²æ˜ ç”») 「売国奴」さんは、この映画をオススメのようですね(笑) >http //ilovenippon.jugem.jp/?eid=559こんなサイトでもオススメの様です。 そしてこのサイトも(笑々) 松井石根大将の南京入城から始まる国策(?)映画のようです、負傷兵を看護する看護婦、日本軍による戦没者のための合同慰霊祭等死体は全く出てきません。 それにしては、見てもいないようですね。こんな嘘8を書いているところをみると。 ~~~~~ この映画の中で愛国者である私がお勧めするシーンは、『査問会』のシーンです。 カメラマンの白井茂は、査問を受ける中国人の心理を他のシーンにはないクローズアップショットで捉えています。写された中国人の目が語るものは、軍部の検閲を受けたナレーションとは違うものです。 軍部の宣伝映画としての大きな制約の中で、「味方」と「敵」との双方の人間を撮えようと努めたスタッフの姿勢が垣間見えます。 282 hits -------------------------------------------------------------------------------- [6530]Re(1) 南京 (戦線後方記録映画) ←back ↑menu ↑top forward→ 平社員 - 09/7/13(月) 22 26 - -------------------------------------------------------------------------------- 引用なし パスワード この映画の中で愛国者である私がお勧めするシーンは、『査問会』のシーンです。 カメラマンの白井茂は、査問を受ける中国人の心理を他のシーンにはないクローズアップショットで捉えています。写された中国人の目が語るものは、軍部の検閲を受けたナレーションとは違うものです。 軍部の宣伝映画としての大きな制約の中で、「味方」と「敵」との双方の人間を撮えようと努めたスタッフの姿勢が垣間見えます。 うーん。私には行間を読む能力が欠けているので、確信はできないのですが「やらせ」があるように思えるということでしょうか? 当然、このような映画には「都合の悪い部分」を写さないと言うことがあるとは思いますけれど。 265 hits -------------------------------------------------------------------------------- [6531]Re(1) HN変更とチャンネル桜 ←back ↑menu ↑top forward→ 熊猫 - 09/7/14(火) 3 18 - -------------------------------------------------------------------------------- 引用なし パスワード ▼dozugawaさん: オフの生活が忙しすぎて、ネット生活を休止中の熊猫です。今日は久々に酒量が少なかったせいか、夜中に目が覚めたついでにレスさせていただきます。 水島総氏の映画作りそのものは支持されるべきものであり、氏の姿勢は尊重すべきではないかと思います。 民主主義とは自分の考えとは違う意見を持った人がその主張をすることを保証することにあるのではないかと思います。その主張の手段として映画とは実に平和的な方法です。映画の内容はともかく言論や表現の自由を守るという意味でも、水島総氏の行為は尊重されるべきでしょう。 「完成した作品を観てから批評する」という立場をとるべきでしょう。映画を評価できるのは、その作品を観た人のみであり、ここはひとつ映画の完成を楽しみにしましょう。 自分が観た映画を批評するのも言論の自由なのですから、楽しみに待っては如何でしょうか。 271 hits -------------------------------------------------------------------------------- [6532]Re(2) HN変更とチャンネル桜 ←back ↑menu ↑top forward→ dozugawa - 09/7/14(火) 21 49 - -------------------------------------------------------------------------------- 引用なし パスワード ▼熊猫さん: ▼dozugawaさん: オフの生活が忙しすぎて、ネット生活を休止中の熊猫です。今日は久々に酒量が少なかったせいか、夜中に目が覚めたついでにレスさせていただきます。 水島総氏の映画作りそのものは支持されるべきものであり、氏の姿勢は尊重すべきではないかと思います。 民主主義とは自分の考えとは違う意見を持った人がその主張をすることを保証することにあるのではないかと思います。その主張の手段として映画とは実に平和的な方法です。映画の内容はともかく言論や表現の自由を守るという意味でも、水島総氏の行為は尊重されるべきでしょう。 「完成した作品を観てから批評する」という立場をとるべきでしょう。映画を評価できるのは、その作品を観た人のみであり、ここはひとつ映画の完成を楽しみにしましょう。 自分が観た映画を批評するのも言論の自由なのですから、楽しみに待っては如何でしょうか。 >水島氏はチャンネル桜の放映資金もないのに、募金を募って「南京の真実」を4部作で作ろうとしています。 >実際のところその様な映画が必要か?、支援者は東中野某や櫻井某といった怱々(?)たる面々です。 >私は右翼と聞いただけでアレルギーが....の体質なので。関係ないですが拉致被害者横田さんのご両親の公演の際ブルーリボンを買おうとしたのですが「新しい教科書をつくる会」主催の為、公演だけ聞いて帰ってきました。 249 hits -------------------------------------------------------------------------------- [6534]Re(2) 南京 (戦線後方記録映画) ←back ↑menu ↑top forward→ ピッポ - 09/7/14(火) 23 51 - -------------------------------------------------------------------------------- 引用なし パスワード ▼平社員さん: うーん。私には行間を読む能力が欠けているので、確信はできないのですが「やらせ」があるように思えるということでしょうか? 東宝映画社、そして後の日本映画社の撮影クルーは、三脚を立てた35ミリの大型カメラで大人数のスタッフを擁して撮影していたようです。 ですから、大砲の発射、歩兵の突撃、城門突入など戦闘シーンなの殆どが、カメラの前で「よーい、すたーと」、号令によって進められていたでしょう。スタッフは、小隊長さんや分隊長さんと仲良くなり、その号令で「よいショット」を稼いだに違いありません。 しかし、これをもって「ヤラセ」というのは酷です。完成した映画は劇場で上映されたわけですから、アクションが画面の中央で起こらなくては迫力のある映画にはなりません。現実を反映した「再現」を「ヤラセ」といったら、ビデオカメラ以前のフィルム時代の動画は、すべてヤラセということになってしまいます(笑)。 当然、このような映画には「都合の悪い部分」を写さないと言うことがあるとは思いますけれど。 おそらく、東京の陸軍省の段階で「台本」または「梗概」を提出して許可を貰い、戦場では現地軍の指示で撮影項目と日程が決められていたでしょう。そうして更に、明日の許可を確認する毎日だったことでしょう。 「都合の悪い部分」を写してしまうということは、それでかなり制約できますが、なにしろ戦場ですから、「写しちゃった」という残酷シーンはあります。 検閲の中心は、編集と、ナレーション録音です。何度もなんども軍の情報部将校に試写し、一字一句細かくチェックされたでしょう。完成した後も上層部の裁可が必要だったようです。もし、出来上がってから不許可になると、ビデオと違って映画の場合は物凄い費用がかかりますから、そうなった場合は製作会社は大損です。 したがって検閲は絶対です。 加えて映画の場合は、軍がスポンサーのようなものです。戦争PR映画といってもよく、報道とは大違いでした。しかも、映画、新聞、ラジオ、当時の報道はみな「報導」と書き、ジャーナリズムではなく国民を"報せ導く"ものでした。戦争を支える世論をつくる「宣伝」機関でした。 そのような中では、映画のストーリーやナレーションは、なかなかスタッフの思うとおりにはなりません。しかしスタッフは、映像の中で人間の真実を伝えようとしました。 当時のスタッフは、大正デモクラシーの影響をうけて、社会主義リアリズムやジャーナリズムの理想に燃えた世代でした。生きた人間のダイナミックな心を写し取る映画の魅力にあこがれた人たちでした。ですから戦場での兵士の葛藤、敵国民の感情を、画面に捉えようと最大限の努力をしたのです。 そうした努力の結果が、亀井文夫監督の「戦ふ兵隊」「上海」に現れており、戦争体験世代の中には、あれは戦争映画ではなくヒューマニズム映画だった、と言う人もいました。白井茂さんの「南京」も、私は『査問』シーンの中にヒューマニズムを感じます。徳川夢声のナレーションを消して映像だけをサイレンとで見ると、中国人の目が語るものは何か、想像力を掻き立てられ、ゾクゾクするものがあります。 映画「南京」などを見るだけでなく、スタッフだった人の手記なども大いに参考にしてください。 260 hits -------------------------------------------------------------------------------- [6535]Re(3) HN変更とチャンネル桜 ←back ↑menu ↑top forward→ 熊猫 - 09/7/15(水) 8 32 - -------------------------------------------------------------------------------- 引用なし パスワード ▼dozugawaさん: 水島氏はチャンネル桜の放映資金もないのに、募金を募って「南京の真実」を4部作で作ろうとしています。 募金を募って映画製作するなんて非常に素晴らしいじゃないですか。右翼の皆さんには、どんどんお金を出してもらいましょうよ。これこそ草の根運動の基本です! 実際のところその様な映画が必要か?、支援者は東中野某や櫻井某といった怱々(?)たる面々です。 必要か不必要かではなく、表現の自由または言論の自由の問題です。更に同じ考えを持った人同士が力を合わせることは悪いことではありません。ただ東中野先生は支援したくても一身上の都合でお金が無いと思います。 必要無いと思った人は、協力しなければ良いだけのことです。完成した作品を観て素晴らしければ認めるし、駄作であれば批判しましょう。 259 hits -------------------------------------------------------------------------------- [6536]Re(2) 南京 (戦線後方記録映画) ←back ↑menu ↑top forward→ 熊猫 - 09/7/15(水) 8 59 - -------------------------------------------------------------------------------- 引用なし パスワード ▼平社員さん: うーん。私には行間を読む能力が欠けているので、確信はできないのですが「やらせ」があるように思えるということでしょうか? 南京が陥落したのは何時? 映画の撮影スタッフが南京に到着したのは何時? 1937年にタイムマシンが実用化されていたという記録は存在しません。過去の出来事を撮影することはできません。つまり「やらせ」ということです。 当然、このような映画には「都合の悪い部分」を写さないと言うことがあるとは思いますけれど。 【白井茂 証言】 虐殺の現場は二度見た。一度はサクがあったように思う。はるか離れているところで、銃殺していた。数は憶えていない。揚子江でない川のところで、機関銃で撃っているところも見た。私なら抵抗すると思ったが、彼等は従順に死を待っていたようだ。川にとび込んで、向うに泳きついた者もいた。二百人ぐらいいたと思う。場所は憶えていない。 【藤井慎一 証言】 小さな川の傍の門の中で捕虜らしき者を撃っているのを見た。白井氏と一緒だった。いくつかの死体に石油をまき、火をつけた。中に生きている兵隊がいて"早く射て”と胸を指さし、"蒋介石万歳”といったので、大変驚いたのを憶えている。それ以外にも、銀行の裏で百人以上が殺されているのを見た。胸のあたりを銃剣で突いていたように思う。虐殺の噂はきいたように思うが、見たのはこの時だけである。 残念ながらこれらのシーンは撮影されていません。 369 hits -------------------------------------------------------------------------------- [6537]Re(3) 南京 (戦線後方記録映画) ←back ↑menu ↑top forward→ 平社員 - 09/7/15(水) 17 36 - -------------------------------------------------------------------------------- 引用なし パスワード なるほど・・・一概にウソ映画と言い切ることはだめですね。 もう中古品しかなさそうですが、興味があるので今度見てみることにします。 あともう一つ、日本のように「写させない」ではなくて「写して削除する」という方法は珍しいのでしょうか? 358 hits -------------------------------------------------------------------------------- [6538]Re(4) 南京 (戦線後方記録映画) ←back ↑menu ↑top forward→ ピッポ - 09/7/15(水) 19 35 - -------------------------------------------------------------------------------- 引用なし パスワード ▼平社員さん: なるほど・・・一概にウソ映画と言い切ることはだめですね。 私の場合は、ナレーションは「宣伝」であって「事実」ではないと思っています。 しかし、そうした宣伝物の中にも、あるいは軍の命令書の中にも、透かしてみえる真実が隠されていると思います。 あともう一つ、日本のように「写させない」ではなくて「写して削除する」という方法は珍しいのでしょうか? スイマセン、意味がとれません。 358 hits -------------------------------------------------------------------------------- [6539]Re(5) 南京 (戦線後方記録映画) ←back ↑menu ↑top forward→ 平社員 - 09/7/18(土) 18 46 - -------------------------------------------------------------------------------- 引用なし パスワード いや、日本語がおかしかったですね。 スルーしてください 337 hits --------------------------------------------------------------------------------
https://w.atwiki.jp/suchikeisan/pages/7.html
1.Newton@Ì f(x), df(x) ÌÖW𦷮ð±¯D restart; y-f(x1) = (df(x1))(x-x1); y - f(x1) = df(x1)(x - x1) y = (df(x1))(x-x1)+f(x1); y = df(x1)(x - x1) + f(x1) y=0; df(x1)(x-x1)+f(x1)=0; df(x1)(x-x1)=-f(x1); x-x1=-f(x1)/df(x1); x=-f(x1)/df(x1)+x1; y = 0 df(x1)(x - x1) + f(x1) = 0 df(x1)(x - x1) = -f(x1) f(x1) x - x1 = - ------ df(x1) f(x1) x = - ------ + x1 df(x1) ñª@Ej [g@ restart; Digits = 10; f = unapply(exp(-x)-2*exp(-2*x), x); plot(f(x), x = -1 .. 2); x - exp(-x) - 2 exp(-2 x) ñª@ x1 =-1.0 x2 =2.0 f1 =f(x1) f2 =f(x2) epsilon =0.001 delta =0.001 for i from 1 to 100 do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.0) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; printf("%10.5f, %10.5f, %10.5f, %10.5f\n", x1, x2, f1, f2); if (abs(x2-x1) delta) then break; end if; if (abs(f2-f1) epsilon) then break; end if; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 0.50000, 2.00000, -0.12923, 0.09870 0.50000, 1.25000, -0.12923, 0.12233 0.50000, 0.87500, -0.12923, 0.06931 0.68750, 0.87500, -0.00285, 0.06931 0.68750, 0.78125, -0.00285, 0.03861 0.68750, 0.73438, -0.00285, 0.01938 0.68750, 0.71094, -0.00285, 0.00866 0.68750, 0.69922, -0.00285, 0.00301 0.68750, 0.69336, -0.00285, 0.00011 0.69043, 0.69336, -0.00136, 0.00011 0.69189, 0.69336, -0.00063, 0.00011 11åç®, x = 0.69189, f(x) = -0.00063 Newton@ x4 = -1.0; df = unapply(diff(f(x), x), x); epsilon = 0.1e-2; delta = 0.1e-2; for i to 100 do x5 = x4-f(x4)/df(x4); printf("%10.5f,~%10.5f\n", x5, f(x5)); if abs(x5-x4) delta then break end if; if abs(f(x5)-f(x4)) epsilon then break end if; x4 = x5 end do; printf("%d aring; ccedil; reg;,~x~=~%.5f,~f(x)~=~%.5f\n", i, x4, f(x4)); x - -exp(-x) + 4 exp(-2 x) -0.55064, -4.28169 -0.13485, -1.47479 0.22538, -0.47607 0.49737, -0.13151 0.64833, -0.02397 0.69032, -0.00142 0.69314, -0.00001 0.69315, 0.00000 8åç®, x = 0.69314, f(x) = -0.00001 absÍâÎl 3(2004ú) 3(a) restart; with(plots) with(plottools) Digits =10 f =unapply(exp(-x)-x^2,x); plot(f(x),x=-1..1); 2 x - exp(-x) - x x1 =-1.0 x2 =1.0 f1 =f(x1) f2 =f(x2) epsilon =0.001 delta =0.001 for i from 1 to 100 do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; printf("%10.5f, %10.5f, %10.5f, %10.5f\n", x1, x2, f1, f2); if (abs(x2-x1) delta) then break; end if; if (abs(f2-f1) epsilon) then break; end if; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 0.00000, 1.00000, 1.00000, -0.63212 0.50000, 1.00000, 0.35653, -0.63212 0.50000, 0.75000, 0.35653, -0.09013 0.62500, 0.75000, 0.14464, -0.09013 0.68750, 0.75000, 0.03018, -0.09013 0.68750, 0.71875, 0.03018, -0.02924 0.70312, 0.71875, 0.00065, -0.02924 0.70312, 0.71094, 0.00065, -0.01425 0.70312, 0.70703, 0.00065, -0.00679 0.70312, 0.70508, 0.00065, -0.00307 0.70312, 0.70410, 0.00065, -0.00121 11åç®, x = 0.70410, f(x) = -0.00121 x4 =-1.0 df =unapply(diff(f(x),x),x); epsilon =0.001 delta =0.001 for i from 1 to 100 do x5 =x4-f(x4)/df(x4); printf("%10.5f, %10.5f\n", x5, f(x5)); if (abs(x5-x4) delta) then break; end if; if (abs(f(x5)-f(x4)) epsilon) then break; end if; x4 =x5; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x4, f(x4)); x - -exp(-x) - 2 x 1.39221, -1.68973 0.83509, -0.26353 0.70983, -0.01214 0.70348, -0.00003 0.70347, -0.00000 5åç®, x = 0.70348, f(x) = -0.00003 3(b) äºåæ³ x1 =-1.0 x2 =1.0 f1 =f(x1) f2 =f(x2) x_f =0.7034674225; lst1 =[] for i from 1 to 100 do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst1 =[op(lst1), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 0.7034674225 101åç®, x = 0.70347, f(x) = 0.00000 Newtonæ³ x4 =-1.0 df =unapply(diff(f(x),x),x); x_f =0.7034674225; lst2 =[] for i from 1 to 100 do x5 =x4-f(x4)/df(x4); lst2 =[op(lst2), log(abs(x_f-x5))]; x4 =x5; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x4, f(x4)); x - -exp(-x) - 2 x 0.7034674225 101åç®, x = 0.70347, f(x) = 0.00000 l1 =listplot(lst1,linestyle=solid,legend="äºåæ³"") l2 =listplot(lst2,linestyle=dashdot,legend="Newtonæ³"") display([l1,l2]); 4(2008Nx) restart; Digits =10 f =unapply(x^4-x-0.12,x); plot(f(x),x=0..2); 4 x - x - x - 0.12 x1 =0.0 x2 =2.0 f1 =f(x1) f2 =f(x2) epsilon =1.0*10^(-8) delta =1.0*10^(-8) for i from 1 to 1000 do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; if (abs(x2-x1) delta) then break; end if; if (abs(f2-f1) epsilon) then break; end if; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 28åç®, x = 1.03717, f(x) = 0.00000 6(ü@) restart; Digits =10 func =x- x^2-4*x+1; x1 =0 x2 =2 f1 =func(x1) f2 =func(x2) plot({(z-x1)*(f1-f2)/(x1-x2)+f1,func(z)},z=0..2); x_f =fsolve(func(x)=0,x)[1]; 2 x - x - 4 x + 1 0.2679491924 with(plots) with(plottools) f =func N =40 å²ç·æ³ x1 =0 x2 =2 f1 =f(x1) f2 =f(x2) lst1 =[] for i from 1 to N do y1 =((f2-f1)/(x2-x1))*(x-x1)+f1 x3 =fsolve(y1=0,x); f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst1 =[op(lst1), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 41åç®, x = 0.26795, f(x) = 0.00000 äºåæ³ x1 =0.0 x2 =2.0 f1 =f(x1) f2 =f(x2) lst2 =[] for i from 1 to N do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst2 =[op(lst2), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 41åç®, x = 0.26795, f(x) = -0.00000 Newtonæ³ x4 =0.0 df =unapply(diff(f(x),x),x) lst3 =[] for i from 1 to N do x5 =x4-f(x4)/df(x4); lst3 =[op(lst3), log(abs(x_f-x5))]; x4 =x5; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x4, f(x4)); 41åç®, x = 0.26795, f(x) = 0.00000 l1 =listplot(lst1,linestyle=solid,legend="ü@") l2 =listplot(lst2,linestyle=dashdot,legend="ñª@") l3 =listplot(lst3,linestyle=spacedash,legend="Newton") display([l1,l2,l3]); 7(2009Nx) restart; Digits =10 with(plots) with(plottools) f =unapply(cos(x)-x^2,x); x_f =fsolve(f(x)=0,x); N =50 plot(f(x),x=0..1); 2 f = x - cos(x) - x x_f = 0.8241323123 äºåæ³ x1 =0.0 x2 =1.0 f1 =f(x1) f2 =f(x2) lst1 =[] for i from 1 to N do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst1 =[op(lst1), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 51åç®, x = 0.82413, f(x) = -0.00000 å²ç·æ³ x1 =0.0 x2 =1.0 f1 =f(x1) f2 =f(x2) lst2 =[] for i from 1 to N do y1 =((f2-f1)/(x2-x1))*(x-x1)+f1 x3 =fsolve(y1=0,x); f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst2 =[op(lst2), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 51åç®, x = 0.82413, f(x) = -0.00000 l1 =listplot(lst1,linestyle=solid,legend="ñª@") l2 =listplot(lst2,linestyle=dashdot,legend="ü@") display([l1,l2]); 8(2010Nx) restart; Digits =10 f =unapply(x^3-3*x+3,x); plot(f(x),x=-5..5); 3 x - x - 3 x + 3 x1 =-3; df =unapply(diff(f(x),x),x); epsilon =0.001 delta =0.001 for i from 1 to 100 do x2 =x1-f(x1)/df(x1); printf("%10.5f, %10.5f\n", x2, f(x2)); if (abs(x2-x1) delta) then break; end if; if (abs(f(x2)-f(x1)) epsilon) then break; end if; x1 =x2; end do printf("%dñÚ, x = %.5f, f(x) = %.5f\n", i, x1, f(x1)); -3 2 x - 3 x - 3 -2.37500, -3.27148 -2.14001, -0.38046 -2.10458, -0.00801 -2.10380, -0.00000 4ñÚ, x = -2.10458, f(x) = -0.00801 x3 =2; df =unapply(diff(f(x),x),x); epsilon =0.001 delta =0.001 for i from 1 to 100 do x4 =x3-f(x3)/df(x3); printf("%10.5f, %10.5f\n", x3, f(x3)); if (abs(x4-x3) delta) then break; end if; if (abs(f(x4)-f(x3)) epsilon) then break; end if; x3 =x4; end do printf("%dñÚ, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 2 2 x - 3 x - 3 2.00000, 5.00000 1.44444, 1.68038 0.92887, 1.01482 3.39447, 31.92927 2.38301, 9.38346 1.71449, 2.89625 1.21672, 1.15109 0.41805, 1.81891 1.15275, 1.07357 0.06453, 2.80668 1.00400, 1.00005 -40.56380, -66619.86560 -27.05959, -19729.42266 úlðx=-3Ƶ½êÍ4ñÚÅðªoܵ½ªúlðx=2Ƶ½êÍvZʪoܹñŵ½B
https://w.atwiki.jp/bc5656/pages/1650.html
夏凛さんが入室しました 夏凛 (同刻。キアシスの街を歩きスマホするガール 夏凛 っはぁ? 容量5GB? このペースで増やして大丈夫かっての。 夏凛 (なんかスマホ見て独り言言いながら歩いている 輝空さんが入室しました 輝空 おおっと危ない!!(夏凛の前に空中からロングコートはためかせて現れる青年 夏凛 っはぁ?(スマホスクロールしながら立ち止まる。輝空の方は見ない 輝空 こんな街中で歩きスマホだなんて都会育ちは全くマナーがなってないものだね! 輝空 ボクがすんでの所で気づいて声をかけたからよかったものの! 輝空 君が全く気づかず一般人にぶつかっていたら事だったよ! 夏凛 っはぁ。じゃあ今度は「歩道での浮遊魔術及び着陸禁止」ってのもマナーに加えとけっての。未来の議長さん。 輝空 ボクを議長と認めて法案を提案してくるのは立派な精神だが、内容は褒められないねえ! 輝空 古来より箒にまたがり空をかける魔女ぐらいは都会者の君でも絵本で見たことがあるだろう? 輝空 浮遊を禁ずる魔法使いの街が何処にあるって話さ! 現代人の身勝手な歩きスマホとはワケが違う! 夏凛 あーっそ。じゃあその古来より伝わる浮遊方で華麗にアタシを避けろって話よね。 夏凛 どいつもこいつも歩きスマホしてる都会でもどんくさい奴以外はぶつからないっての。 夏凛 自由自在に空飛び回るぐらいの元気さがあるってのに歩行人程度の体当たり食らうわけがないわよね? 輝空 それは歩きスマホ側のどんくささ次第としか言いようがないがね。 輝空 そもそも論で歩きスマホはするべきじゃあない。君の都会では勿論。ボクのこのキアシスでもね。 夏凛 っはぁー・・・で。この無限に広がらない話をするためだけにわざわざ地上まで降りてきたってワケ? 輝空 っはぁー・・・思い上がるのもいい加減にしてくれないかい? 輝空 新魔術好きの中年に毎日チヤホヤされてるからって誰でも接待対応してくれるとは思わない事だよ。 輝空 ボクが降りてきたのは何でも無い。ただ降りる場所だったというだけさ。 輝空 君に用なんて全く無い。 夏凛 あーっそ。 輝空 でもまあどうしてもというならキアシス新米の君に少し昔話を聞かせてあげないこともない。 輝空 ああ、そうか、聞きたいか。なら仕方ない。少しばかり話してあげよう。 夏凛 ぜーんぜん興味ないけど。 輝空 「魔女が目覚めていた」 輝空 「魔女は封印されていたわけじゃない。遠くに距離を置いていただけなんだ。」 輝空 「来ようと思えば来ることができた。でもしなかった。でも、今回、来た。」 輝空 「歴史の長い魔術都市キアシスでは、このような昔話はいくつもある。」 輝空 「いや、忘れ去られた話もあるかもしれない。それこそが次かも。」 輝空 「せいぜいアンテナを高くしておくことだ。」 輝空 ボクの街を守る気があるんならね。 夏凛 っはぁー・・・? 輝空 やっぱり無駄な話だったみたいだね! 輝空 全く、何故君がチヤホヤされるのか理解に苦しむよ! 輝空 そう。チヤホヤされ続けて今の地位を騙し取ったのなら、この話に興味が沸かざる負えないはずなのにね! 輝空 嗚呼、嘆かわしい! 輝空 新入りがコレじゃあ、この街の治安も危ぶまれるってものさ! 夏凛 っはぁ・・・ 夏凛 最近物騒だから八首としての自覚があるんなら街を守れとでも言いたいワケ? 輝空 ばぁーーーっかじゃないのか?! だーーーれが君なんか頼るものかよ! 夏凛 ふーーーーーん。意外なもんね。 輝空 勝手に勘違いして納得かい? 自意識過剰で妄想癖ときた君をおだてて使い回すのは中年たちの絶好の娯楽だろうね! 夏凛 っはぁーーーーー。 夏凛 あのね。アタシだってそんな我慢強い方じゃあないんだけど? 輝空 そりゃあ良かった。この街に合わないと思ったら一刻一秒でも早く故郷に帰ってもらって構わないよ? 夏凛 ふーーーん。考えておくわ。 夏凛 この街で暮らしててもアンタに会わずに済む方法をね。 夏凛 (踵返して来た道を戻っていく 輝空 それは明暗だ。君が牢屋にでも入ってくれればボクも会わずに済むんだけどね! 夏凛 牢屋で会ってみたい人は居るけど、一生行くのはごめんだっての。 輝空 それはそれは残念だ。 夏凛 (会話も途中に勝手に歩いて見えなくなっていく 夏凛さんが退室しました 輝空 ふん。たったあれだけのやり取りで道を変えて帰ってしまうとは。 輝空 案外やりがいのない子供だよ。全く。喫茶にでも向かう途中だったのかな? 輝空 それならそれで丁度いいとも言えるがね。ああ。だが残念。 キアシスのお空さんが入室しました 輝空 おそらく今この時間は彼がいるだろうからね。鉢合わせておくのも一興だったか。 キアシスのお空 (時刻は同刻、つまり朝方。 星は無いはずなのに、朝空にキラリと何かが光る 輝空 全く。悪運の強いラッキーガールだこと。 キアシスのお空 (再びキラリと何かが光る 先程より少し大きい キアシスのお空 (代わりに光ったのではなく、最初に光った朝星はそのまま―― キアシスのお空 (キラ、キラキラ――と次々と大きく――近付いている? 輝空 輝く 空 か。(天を見上げて 輝空 この僕に今度は何の祝福の光かな? キアシスのお空 (輝きが”視認”できるように なんと 星型 キアシスのお空 (☆ ← これである → ☆ キアシスのお空 (スター型のリング……? けっこうでかい! 直径10tdぐらいあるんじゃないか キアシスのお空 (スター次々と連鎖し、空から地に続くトンネルを描いているようにも見える 輝空 (両足ついて天を仰ぎ両手を広げる 輝空 歓迎しよう。 輝く 空の 星よ キアシスのお空 (輝空の上空に最後のスターが出現する 輝空 輝かしい未来と血統を持つ僕に何の用かな? キアシスのお空 (スターが輝き、虹色のベールで包まれ 文字通りトンネルを描く! キアシスのお空 (ものっそい収束音 何かチャージしてるような、ワープしてるような キアシスのお空 (瞬間 何かが せり出て キアシスのお空さんが退室しました マンガみたいなロケットさんが入室しました マンガみたいなロケット (意外! それはロケット!宇宙までいくアレ! 輝空 ・・・はぁ。(表情曇らせて溜息 輝空 宇宙まで行こう、あるいは宇宙から来ようという気概は認めるが、 輝空 なんだねそのデザインは。まるで宇宙ロケットそのものじゃないか。 マンガみたいなロケット (クレヨンとかで描くような色合い アポロチョコみたいなアレ 輝空 もっとこう夢と希望と魔術と神秘で満ちた幻想的なデザインにできなかったものかね? マンガみたいなロケット (そうこう言ってるうちに落下してきます>輝空 輝空 デザインは大事さ。性能だけでは民の心までは動かせない。 マンガみたいなロケット (直撃すると10D6みたいなダメージ判定が出ます 回避どうぞ 輝空 嘆かわしい。夢半ばで墜落するロケットだなんて。(一歩だけ、前へ歩む マンガみたいなロケット (ずっどーーんむ!! 輝空 ( その場 には、長い間地面に踵をつけ、深く刻まれた装飾魔法陣 マンガみたいなロケット (ロケットのトサカから地面に突き刺さ――突き刺さるれ――何やそれ 輝空 その夢、ボクが再び叶えてあげよう。 輝空 星流大魔法-(詠唱は既に唱えていた 輝く 空の 星よ ) 輝空 スター・ダスト・リヴァース。(魔法陣から煌めく星屑が溢れ出てロケットを空へと押し返す マンガみたいなロケット (オラオラ言いかねない勢いで星が溢れかえ、ロケットが押し返される! マンガみたいなロケット (が、妙に重い――いや魔法弾いてるっぽいぞ? 微妙に マンガみたいなロケット (20tdほど浮き上がったが、つるっと滑って マンガみたいなロケット (ぐるぐる回転しながら落下してくる!>輝空 輝空 ボクの好意が受けれないなんてねっっ、(両足で強く大地を踏んで大きく距離を離すダッシュ 輝空 (ボクは見逃さなかった。あのふざけたデザインのロケットはボクの大魔法を弾いていた。 マンガみたいなロケット (ずずーん!と地面に堕ち転がる巨塊 マンガみたいなロケット .......(砂煙 輝空 (それならば無理に意地を張って攻撃や防御の魔法を重ねるのは二流のやることさ 輝空 (スマートに回避する。戦闘中に機敏に動ける魔術師こそ現代魔術師だからね。 マンガみたいなロケット (砂煙が晴れる頃……横向きの扉がウィイイーンと開く 輝空 (空中に魔法陣を貼って立っている。砂煙を観察。 マンガみたいなロケット (扉に手がかかる 何かが 出てくる出てくる マンガみたいなロケットさんが退室しました 謎の美少女さんが入室しました 謎の美少女 (色テス 色テス 謎の美少女 (こんぐらい 謎の美少女 (やっぱこっち… 謎の美少女 (こんぐらい 謎の美少女 ……、 。。(出にくい横向きの扉から出てくるのは…美少女! 輝空 ごきげんよう! 魔術学園都市キアシスへようこそ! 輝空 旅行者かな? それならば観光の手配をしよう。 謎の美少女 (翡翠色のミディアムボブだが毛先はピンクで、所々にイエローが交じる 輝空 それとも漂流者かな? それならば手厚く保護をしなくちゃね。 謎の美少女 (何か……触覚? アクセ? 頭から二本ほど生えて、先っぽに丸いのついてる 輝空 あるいは…侵略者かな? そうでない事を願っているよ? 謎の美少女 …………(目が虚ろ 謎の美少女 ―――――― !!!!(ゔ!と口を抑える 謎の美少女 (いわゆるロケット酔いのようだ 多分 謎の美少女 ~~~~~~!!!(虹色の[自主規制]が[自主規制]する 輝空 宇宙人幼女の嘔吐シーンか。 筋金入りのフェチズム持ちには刺さる内容かもしれないが、ボクには全く響かないね。 輝空 悪いがボクは幼女と見れば誰でも助けるフェミニストでもヒーローでもない。 謎の美少女 (転がり落ちるようにロケットから降りてくる 謎の美少女 (幼女っていうか女性ぐらいかな 10代後半 輝空 (輝空くんは見下してあえて低めにいったり、高めに言ったりする子なのだ!) 謎の美少女 (ならOKOK 輝空 でも救急隊ぐらいは手配してあげようか。ボクの愛するこの街を得体の知れない吐瀉物で汚すわけにはいかないからね。 輝空 (内ポケットからスマホを取り出し、すばやく画面をタップ 謎の美少女 (フラフラと輝空に近付いてくる 輝空 (ボクはキアシス好きだが、懐古主義者じゃあない。知っていると思うがスマホもサーチライトも使いこなすよ。 輝空 (空中の魔法陣からふわりと降りて、両つま先立ちで着地 輝空 魔導警察の匿名無言回線にこの場所の位置座標を送った。 輝空 早ければすぐにでも迎えが来るだろう。 モブ導警察さんが入室しました 輝空 君が漂流者ならそこで手厚く保護してもらう事だね。 モブ導警察 デアエー!デアエー!(モブモブと現れるモブ警察 輝空 これはこれは。お早いお付きで。 謎の美少女 æ–‡å—åŒ–ã ……? あなたは誰……?(キョトンと輝空を見る モブ導警察 「信号は此処からかー!?」「うわ!なんか落ちてる!」「でっかいなー!」(モブモブとモブめいた感想でアワアワするモブ 輝空 何を言っているかさっぱりわからないが、君はヴァースに感謝するんだね。 輝空 創世神柳の威光により、君はすぐにでもこの世界の言語を手に入れる事だろうさ。 謎の美少女 ãã ã«ã¡ã¯ ……? あれ この反応……? 輝空 その力さえも弾く変わり者だっていうんなら、運がなかったね。 謎の美少女 ããªãã¯ãç §ã®ãéå½ã®ä………ººã§ããï¼? 運命の人 もしかして……あなたが? 輝空 見ての通りだ。ロケットが墜落し、宇宙人が出てきた。 輝空 戦闘能力、敵意があるかは不明だが、街に放って置くわけにもいくまいよ。 謎の美少女 (輝空と警察を交互に見る 何か悪いことしたかも顔 輝空 餅は餅屋。不審者は魔導警察へ。 モブ導警察 「ほうほう、なるりょです!」「あの子かわいいでござる」「街の平和にご協力ありがとうございます!リーズベルト氏!」 輝空 礼を言われる程のことはない。街に住まうものとして当然の行いをしたまでだよ。 モブ導警察 「それでは、詳しいお話は署で伺いましょう」「ねえねえ君歳はいくつ?天パは好み?」(モブモブと美少女に寄っていくモブ警察 謎の美少女 ……ããã ã⦅!(首をふるふるしながら後ずさる い いや…! 輝空 先程も言ったがボクはフェミストではないし、どちらかといえばリアリストだ。 輝空 得体の知れない者は公的機関で保護する。その方針に疑念はないし、君を見過ごすつもりもない。 謎の美少女 ãã®äººéã¯èª°ã§ãã…!? あの人達は誰なんですか…!? (モブ警察達を指さしながら輝空を見る モブ導警察 「こら、軽乃助!被疑者をナンパするな!」「そうだぞ!リーズベルト氏の前だぞ!真面目にするんだぞ!」(モブモブモブ 輝空 君がどうしてもボクの威光に惹かれて、また会いたいと願うなら、 輝空 事情聴取を受け外出許可証を得た上で、街中で遭遇する運命を願うがいいさ。 モブ導警察 「そんなに不安がらないで…お話を聞くだけなので…お話できるか分からないですけど…」(まともなモブが1人いるっぽい 謎の美少女 å©ãã¦! ãå©ãã¦ãã ããï¼! 待って! 悪いことしてない!(壁に追い詰められ、輝空に向かって何か叫んでるが 謎の美少女 ã©ãèãã¦ããã®äººéã¢ãï¼! 助けて運命の人(仮)! 輝空 魔導警察の諸君。宜しく頼むよ。得体の知れない宇宙人ではあるが、危害を加える様子はない。 輝空 一先ずは外敵としてではなく、漂流者として丁重に扱ってくれ。 モブ導警察 「イエッサーです!リーズベルト氏!」 輝空 まあ。このようなこと。わざわざボクが言うまでもなく実行するであろうがね。 謎の美少女 éå½ã®äºº――!! どうして――!! 輝空 (叫ぶ宇宙人に背を向けて 輝空 さようならだ。君がキアシスの治安を脅かさない事を祈っているよ。 輝空 もし君が、輝く空から、一際輝く空 ボク を見つけて落ちてきた星ならば、 輝空 また、会うこともあるだろうね。 輝空 (とかなんとか言いながら爪先立ち歩きで街中へと消えていく モブ導警察 「さ、行きましょう。ココは冷えますしね…!」(モブモブと美少女を囲って 謎の美少女 ãã¶ãéå½ã®äºº――! 運命の人(仮)――!(何か喚いているが捕らえられ 輝空さんが退室しました モブ導警察 (丁重に警察署へとご案内。魔導警察はなんだかんだクリーンだから!御丁重に!御丁重に! 謎の美少女 ã©ãã㦅… オョョ…… モブ導警察 (あと、仕事振りはなんだかんだ真面目だから!ご丁寧に!ご丁寧に! 謎の美少女 (連行されていく モブ導警察 (美少女。ロケット。どちらも真面目に調査致しました。―――成果の程は……… モブ導警察さんが退室しました 謎の美少女 (取り調べを受けるが どこの言語でもない そもそもこれは言葉なのか―― 謎の美少女 (魔術者であっても、読解者であっても、解けることは無かったという 謎の美少女 ãªãã§!? どうして!? 謎の美少女さんが退室しました
https://w.atwiki.jp/suchikeisan/pages/8.html
1.Newton@Ì f(x), df(x) ÌÖW𦷮ð±¯D restart; y-f(x1) = (df(x1))(x-x1); y - f(x1) = df(x1)(x - x1) y = (df(x1))(x-x1)+f(x1); y = df(x1)(x - x1) + f(x1) y=0; df(x1)(x-x1)+f(x1)=0; df(x1)(x-x1)=-f(x1); x-x1=-f(x1)/df(x1); x=-f(x1)/df(x1)+x1; y = 0 df(x1)(x - x1) + f(x1) = 0 df(x1)(x - x1) = -f(x1) f(x1) x - x1 = - ------ df(x1) f(x1) x = - ------ + x1 df(x1) ñª@Ej [g@ restart; Digits = 10; f = unapply(exp(-x)-2*exp(-2*x), x); plot(f(x), x = -1 .. 2); x - exp(-x) - 2 exp(-2 x) ñª@ x1 =-1.0 x2 =2.0 f1 =f(x1) f2 =f(x2) epsilon =0.001 delta =0.001 for i from 1 to 100 do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.0) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; printf("%10.5f, %10.5f, %10.5f, %10.5f\n", x1, x2, f1, f2); if (abs(x2-x1) delta) then break; end if; if (abs(f2-f1) epsilon) then break; end if; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 0.50000, 2.00000, -0.12923, 0.09870 0.50000, 1.25000, -0.12923, 0.12233 0.50000, 0.87500, -0.12923, 0.06931 0.68750, 0.87500, -0.00285, 0.06931 0.68750, 0.78125, -0.00285, 0.03861 0.68750, 0.73438, -0.00285, 0.01938 0.68750, 0.71094, -0.00285, 0.00866 0.68750, 0.69922, -0.00285, 0.00301 0.68750, 0.69336, -0.00285, 0.00011 0.69043, 0.69336, -0.00136, 0.00011 0.69189, 0.69336, -0.00063, 0.00011 11åç®, x = 0.69189, f(x) = -0.00063 Newton@ x4 = -1.0; df = unapply(diff(f(x), x), x); epsilon = 0.1e-2; delta = 0.1e-2; for i to 100 do x5 = x4-f(x4)/df(x4); printf("%10.5f,~%10.5f\n", x5, f(x5)); if abs(x5-x4) delta then break end if; if abs(f(x5)-f(x4)) epsilon then break end if; x4 = x5 end do; printf("%d aring; ccedil; reg;,~x~=~%.5f,~f(x)~=~%.5f\n", i, x4, f(x4)); x - -exp(-x) + 4 exp(-2 x) -0.55064, -4.28169 -0.13485, -1.47479 0.22538, -0.47607 0.49737, -0.13151 0.64833, -0.02397 0.69032, -0.00142 0.69314, -0.00001 0.69315, 0.00000 8 aring; ccedil; reg;, x = 0.69314, f(x) = -0.00001 absÍâÎl 3(2004ú) 3(a) restart; with(plots) with(plottools) Digits =10 f =unapply(exp(-x)-x^2,x); plot(f(x),x=-1..1); 2 x - exp(-x) - x x1 =-1.0 x2 =1.0 f1 =f(x1) f2 =f(x2) epsilon =0.001 delta =0.001 for i from 1 to 100 do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; printf("%10.5f, %10.5f, %10.5f, %10.5f\n", x1, x2, f1, f2); if (abs(x2-x1) delta) then break; end if; if (abs(f2-f1) epsilon) then break; end if; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 0.00000, 1.00000, 1.00000, -0.63212 0.50000, 1.00000, 0.35653, -0.63212 0.50000, 0.75000, 0.35653, -0.09013 0.62500, 0.75000, 0.14464, -0.09013 0.68750, 0.75000, 0.03018, -0.09013 0.68750, 0.71875, 0.03018, -0.02924 0.70312, 0.71875, 0.00065, -0.02924 0.70312, 0.71094, 0.00065, -0.01425 0.70312, 0.70703, 0.00065, -0.00679 0.70312, 0.70508, 0.00065, -0.00307 0.70312, 0.70410, 0.00065, -0.00121 11åç®, x = 0.70410, f(x) = -0.00121 x4 =-1.0 df =unapply(diff(f(x),x),x); epsilon =0.001 delta =0.001 for i from 1 to 100 do x5 =x4-f(x4)/df(x4); printf("%10.5f, %10.5f\n", x5, f(x5)); if (abs(x5-x4) delta) then break; end if; if (abs(f(x5)-f(x4)) epsilon) then break; end if; x4 =x5; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x4, f(x4)); x - -exp(-x) - 2 x 1.39221, -1.68973 0.83509, -0.26353 0.70983, -0.01214 0.70348, -0.00003 0.70347, -0.00000 5åç®, x = 0.70348, f(x) = -0.00003 3(b) äºåæ³ x1 =-1.0 x2 =1.0 f1 =f(x1) f2 =f(x2) x_f =0.7034674225; lst1 =[] for i from 1 to 100 do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst1 =[op(lst1), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 0.7034674225 101åç®, x = 0.70347, f(x) = 0.00000 Newtonæ³ x4 =-1.0 df =unapply(diff(f(x),x),x); x_f =0.7034674225; lst2 =[] for i from 1 to 100 do x5 =x4-f(x4)/df(x4); lst2 =[op(lst2), log(abs(x_f-x5))]; x4 =x5; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x4, f(x4)); x - -exp(-x) - 2 x 0.7034674225 101åç®, x = 0.70347, f(x) = 0.00000 l1 =listplot(lst1,linestyle=solid,legend="äºåæ³"") l2 =listplot(lst2,linestyle=dashdot,legend="Newtonæ³"") display([l1,l2]); 4(2008Nx) restart; Digits =10 f =unapply(x^4-x-0.12,x); plot(f(x),x=0..2); 4 x - x - x - 0.12 x1 =0.0 x2 =2.0 f1 =f(x1) f2 =f(x2) epsilon =1.0*10^(-8) delta =1.0*10^(-8) for i from 1 to 1000 do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; if (abs(x2-x1) delta) then break; end if; if (abs(f2-f1) epsilon) then break; end if; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 28åç®, x = 1.03717, f(x) = 0.00000 6(ü@) restart; Digits =10 func =x- x^2-4*x+1; x1 =0 x2 =2 f1 =func(x1) f2 =func(x2) plot({(z-x1)*(f1-f2)/(x1-x2)+f1,func(z)},z=0..2); x_f =fsolve(func(x)=0,x)[1]; 2 x - x - 4 x + 1 0.2679491924 with(plots) with(plottools) f =func N =40 å²ç·æ³ x1 =0 x2 =2 f1 =f(x1) f2 =f(x2) lst1 =[] for i from 1 to N do y1 =((f2-f1)/(x2-x1))*(x-x1)+f1 x3 =fsolve(y1=0,x); f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst1 =[op(lst1), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 41åç®, x = 0.26795, f(x) = 0.00000 äºåæ³ x1 =0.0 x2 =2.0 f1 =f(x1) f2 =f(x2) lst2 =[] for i from 1 to N do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst2 =[op(lst2), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 41åç®, x = 0.26795, f(x) = -0.00000 Newtonæ³ x4 =0.0 df =unapply(diff(f(x),x),x) lst3 =[] for i from 1 to N do x5 =x4-f(x4)/df(x4); lst3 =[op(lst3), log(abs(x_f-x5))]; x4 =x5; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x4, f(x4)); 41åç®, x = 0.26795, f(x) = 0.00000 l1 =listplot(lst1,linestyle=solid,legend="ü@") l2 =listplot(lst2,linestyle=dashdot,legend="ñª@") l3 =listplot(lst3,linestyle=spacedash,legend="Newton") display([l1,l2,l3]); 7(2009Nx) restart; Digits =10 with(plots) with(plottools) f =unapply(cos(x)-x^2,x); x_f =fsolve(f(x)=0,x); N =50 plot(f(x),x=0..1); 2 f = x - cos(x) - x x_f = 0.8241323123 äºåæ³ x1 =0.0 x2 =1.0 f1 =f(x1) f2 =f(x2) lst1 =[] for i from 1 to N do x3 =(x1+x2)/2; f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst1 =[op(lst1), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 51åç®, x = 0.82413, f(x) = -0.00000 å²ç·æ³ x1 =0.0 x2 =1.0 f1 =f(x1) f2 =f(x2) lst2 =[] for i from 1 to N do y1 =((f2-f1)/(x2-x1))*(x-x1)+f1 x3 =fsolve(y1=0,x); f3 =f(x3); if (evalf(f1*f3) =0.00) then f1 =f3; x1 =x3; else f2 =f3; x2 =x3; end if; lst2 =[op(lst2), log(abs(x_f-x3))]; end do printf("%dåç®, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 51åç®, x = 0.82413, f(x) = -0.00000 l1 =listplot(lst1,linestyle=solid,legend="ñª@") l2 =listplot(lst2,linestyle=dashdot,legend="ü@") display([l1,l2]); 8(2010Nx) restart; Digits =10 f =unapply(x^3-3*x+3,x); plot(f(x),x=-5..5); 3 x - x - 3 x + 3 x1 =-3; df =unapply(diff(f(x),x),x); epsilon =0.001 delta =0.001 for i from 1 to 100 do x2 =x1-f(x1)/df(x1); printf("%10.5f, %10.5f\n", x2, f(x2)); if (abs(x2-x1) delta) then break; end if; if (abs(f(x2)-f(x1)) epsilon) then break; end if; x1 =x2; end do printf("%dñÚ, x = %.5f, f(x) = %.5f\n", i, x1, f(x1)); -3 2 x - 3 x - 3 -2.37500, -3.27148 -2.14001, -0.38046 -2.10458, -0.00801 -2.10380, -0.00000 4ñÚ, x = -2.10458, f(x) = -0.00801 x3 =2; df =unapply(diff(f(x),x),x); epsilon =0.001 delta =0.001 for i from 1 to 100 do x4 =x3-f(x3)/df(x3); printf("%10.5f, %10.5f\n", x3, f(x3)); if (abs(x4-x3) delta) then break; end if; if (abs(f(x4)-f(x3)) epsilon) then break; end if; x3 =x4; end do printf("%dñÚ, x = %.5f, f(x) = %.5f\n", i, x3, f(x3)); 2 2 x - 3 x - 3 2.00000, 5.00000 1.44444, 1.68038 0.92887, 1.01482 3.39447, 31.92927 2.38301, 9.38346 1.71449, 2.89625 1.21672, 1.15109 0.41805, 1.81891 1.15275, 1.07357 0.06453, 2.80668 1.00400, 1.00005 -40.56380, -66619.86560 -27.05959, -19729.42266 úlðx=-3Ƶ½êÍ4ñÚÅðªoܵ½ªúlðx=2Ƶ½êÍvZʪoܹñŵ½B
1 2