下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
1、Vc+6.0 實(shí)現(xiàn)微分方程求解,程序如下:#include <iostream.h>#include <math.h>double f1(double t, double x, double y, double z) return -x + 6 * y + 2 * z;double f2(double t, double x, double y, double z) return -y + 3 * z + 2 * sin(t);double f3(double t, double x, double y, double z) return -z + pow(t,2) *
2、 exp(-t) + cos(t);double RungeKutta(double x0, double y0, double z0, double a, double b, int n, double x, double y, double z, double (*f1)(double,double,double,double), double (*f2)(double,double,double,double), double (*f3)(double,double,double,double) double k1, k2, k3, k4 ; double m1, m2, m3, m4
3、; double n1, n2, n3, n4 ; double h; double t;int count;count =0;h = (b - a) / n; x0 = x0; y0 = y0; z0 = z0; for(int i = 1;i<=n;i+) t = a + (i - 1) * h; k1 = f1(t, xi - 1, yi - 1, zi - 1); m1 = f2(t, xi - 1, yi - 1, zi - 1); n1 = f3(t, xi - 1, yi - 1, zi - 1); k2 = f1(t+h/2.0, xi - 1 + k1 * h / 2.
4、0, yi - 1 + m1 * h / 2.0, zi - 1 + n1 * h / 2.0); m2 = f2(t+h/2.0, xi - 1 + k1 * h / 2.0, yi - 1 + m1 * h / 2.0, zi - 1 + n1 * h / 2.0); n2 = f3(t+h/2.0, xi - 1 + k1 * h / 2.0, yi - 1 + m1 * h / 2.0, zi - 1 + n1 * h / 2.0); k3 = f1(t+h/2.0, xi - 1 + k2 * h / 2.0, yi - 1 + m2 * h / 2.0, zi - 1 + n2 *
5、 h / 2.0); m3 = f2(t+h/2.0, xi - 1 + k2 * h / 2.0, yi - 1 + m2 * h / 2.0, zi - 1 + n2 * h / 2.0); n3 = f3(t+h/2.0, xi - 1 + k2 * h / 2.0, yi - 1 + m2 * h / 2.0, zi - 1 + n2 * h / 2.0); k4 = f1(t+h/2.0, xi - 1 + k3 * h / 2.0, yi - 1 + m3 * h / 2.0, zi - 1 + n3 * h / 2.0); m4 = f2(t+h/2.0, xi - 1 + k3
6、 * h / 2.0, yi - 1 + m3 * h / 2.0, zi - 1 + n3 * h / 2.0); n4 = f3(t+h/2.0, xi - 1 + k3 * h / 2.0, yi - 1 + m3 * h / 2.0, zi - 1 + n3 * h / 2.0); xi = xi - 1 + h * (k1 + 2 * k2 + 2 * k3 + k4) / 6.0; yi = yi - 1 + h * (m1 + 2 * m2 + 2 * m3 + m4) / 6.0; zi = zi - 1 + h * (n1 + 2 * n2 + 2 * n3 + n4) / 6.0; +count;cout<<"t= "<<t+0.05<<" "<<"x= "<<xi<<" "<<"y= "<<yi<<" "<<"z= "<<zi<<endl;cout<<"共計(jì)算的次數(shù)為:"<<count<<endl;return 0; vo
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 采購年終工作總結(jié)3篇(企業(yè)采購年終工作總結(jié))
- 人教版數(shù)學(xué)八年級(jí)上冊13.1.2.1 線段的垂直平分線的性質(zhì)和判定教案
- 高一上學(xué)期體育與健康人教版 竹竿舞 力量素質(zhì) 教案
- 【核心素養(yǎng)目標(biāo)】第9課 辛亥革命 教案
- 第四課“巨人”輩出-文藝復(fù)興美術(shù)教案
- 2024年甲氨蝶呤項(xiàng)目資金籌措計(jì)劃書代可行性研究報(bào)告
- 粵教版初中化學(xué)九年級(jí)上冊 2.3 物質(zhì)構(gòu)成的奧秘(Ⅱ)-原子和離子 教案
- 【湘美版】四下 第7課 剪紙故事 教案
- 2023年秋人教版物理九年級(jí)上學(xué)期15.4電流的測量 教案
- 4.2《憐憫是人的天性》課件+2024-2025學(xué)年統(tǒng)編版高中語文選擇性必修中冊
- 寫桂花的作文指導(dǎo)課件.ppt
- 體育與健康課程《體室內(nèi)理論課》課件PPT
- 2021年成人高考高起本理化綜合考試模擬真題及答案
- 小學(xué)數(shù)學(xué)三年級(jí)下冊 《認(rèn)識(shí)分?jǐn)?shù)》說課課件
- 水電站水輪發(fā)電機(jī)組A級(jí)檢修施工組織設(shè)計(jì)[詳細(xì)]
- 瀝青及水穩(wěn)拌合站建設(shè)標(biāo)準(zhǔn)(附CAD圖)
- GA∕T 1737-2020 公安移動(dòng)信息網(wǎng)技術(shù)要求
- 現(xiàn)代漢語詞典附錄中的我國歷代紀(jì)元表19頁
- 高職課程標(biāo)準(zhǔn)改革存在的問題及對(duì)策
- 康復(fù)醫(yī)學(xué)概論康復(fù)醫(yī)學(xué)是由理療學(xué)物理醫(yī)學(xué)逐漸展形成一門新學(xué)科課件
- 噴涂自動(dòng)生產(chǎn)線安全技術(shù)操作規(guī)程
評(píng)論
0/150
提交評(píng)論