1樓:匿名使用者
改進的尤拉方法
1.子函式
function[x,y]=eulerpro(fun,x0,xfinal,y0,n)
if nargin<5
n=50;
endh=(xfinal-x0)/n;%步長內
x(1)=x0;y(1)=y0;
for i=1:n
x(i+1)=x(i)+h;
y1=y(i)+h*feval(fun,x(i),y(i));
y2=y(i)+h*feval(fun,x(i+1),y1);
y(i+1)=(y1+y2)/2;
endend
2.子函式
function f=doty(x,y)
f=cos(x*y);
end3.主函
容數呼叫
[x,y]=eulerpro('doty',0,1,1,10)
用matlab求解常微分方程求初值問題,原題是xdy+(x2–y)dx=0.
2樓:我行我素
先變形為dy/dx=y/x-x,再用dsolve求通解或ode45求數值解。如:
syms y(x)
y=dsolve(diff(y)==y/x-x)結果是:
y =- x^2 + c1*x
3樓:時間變奏曲
用solve函式可解常微分方程,具體引數參考百科
用matlab求解常微分方程求初值問題原題是
先變形為dy dx y x x，再用dsolve求通解或ode45求數值解。如 syms y x y dsolve diff y y x x 結果是 y x 2 c1 x 用solve函式可解常微分方程，具體引數參考百科 大學本科學習matlab數學軟體的一般是什麼專業的學生？資訊與計算科學專業，有...
一道微分方程問題求解,求解一道常微分方程習題,提如下圖
原來題中函式baix 1處無定義,du求x趨於1 與zhix趨於1 時的極限按題意理應分別採dao 用各自的式子求專,但最後題目又屬要求補充定義後使得定義域為全體實數且連續,則要求x 1處極限存在 左右極限相等且等於該點函式值 我沒看見你昭下答案。原理上將一階微分方程的求解公式用進去就可以求出的 求...
什麼是高階常微分方程，高數 常微分方程 高階微分方程，有三道題，求大神幫忙解答！
如果在一個微分方程中出現的未知函式只含一個自變數,這個方程就叫做常微分方程,也可以簡單地叫做微分方程.高階常微分方程就是自變數的次數大於一次的常微分方程了.很高興為你解答有用請採納 高數 常微分方程 高階微分方程，有三道題，求大神幫忙解答！第一題的問題 f 1 2隱含著的條件是，f 1 2 所以,f...