function?[?output_args?]?=?Untitled1（?input_args?）
%UNtitleD1?Summary?of?this?function?goes?here
%??Detailed?explanation?goes?here
clc;
clear;
h=0.1;
a=0;
b=1;
ya=1;
N=（b-a）/h;
T=linspace（abN+1）;
Y=zeros（1N+1）;
Y（1）=ya;

%Euler法
for?j=1:N
????Y（j+1）=Y（j）+h*（Y（j）-T（j）*Y（j）^2）;
end
plot（TY‘g‘）;grid?on;
hold?on;

%Euler改進法
for?j=1:N
???Y（j+1）=Y（j）+（h/2）*（Y（j）-T（j）*Y（j）^2+??（Y（j）+h*（Y（j）-T（j）*Y（j）^2））-T（j+1）*（Y（j）+h*（Y（j）-T（j）*Y（j）^2））^2??）;
end
plot（TY‘r‘）;
hold?on;

%精確解
Y=1./（T-1+2*exp（-T））;
plot（TY）;
xlabel（‘x‘）;
ylabel（‘f（x）‘）;
title（‘Euler（綠色）、Euler改進法（紅色）及精確解（藍色）的比較‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????611??2010-01-03?22:41??常微分方程數值解法\常微分方程數值解法.m

?????目錄??????????0??2010-01-09?23:06??常微分方程數值解法

-----------?---------??----------?-----??----

??????????????????611????????????????????2