用4阶经典R-K法编写matlab程序

h=1/128

tf=3

t = 1:h:tf

x = zeros(1,length(t))

x(1) = 2

F_tx = @(t, x)(t.*x-x.^2)./t.^2

for i=1:(length(t)-1)

k_1 = F_tx(t(i),x(i))

k_2 = F_tx(t(i)+0.5*h,x(i)+0.5*h*k_1)

k_3 = F_tx((t(i)+0.5*h),(x(i)+0.5*h*k_2))

k_4 = F_tx((t(i)+h),(x(i)+k_3*h))

x(i+1) = x(i) + (1/6)*(k_1+2*k_2+2*k_3+k_4)*h

end

t0 = t(1)x0 = x(1)

xspan = [t0 tf]

[x_ode45, y_ode45] = ode45(F_tx,xspan,x0)

subplot(211)

plot(x_ode45,y_ode45,'--')

xlabel('t')ylabel('x')

legend('Exact')

subplot(212)

plot(t,x,'-')

xlabel('t')ylabel('x')

legend('Approximation')

和。表达式：I=I+IIR分流原理：

MATLAB 7这个版本很经典，而且也不大。

好多专业认识都用matlab7这个版本。

关键是把各项功能都熟悉了才是最重要的，版本高的是对程序进行升级和改进。

I am a t aab 2008 m ataab，2010哪个比较经典呢？现在经典呢，那就是当然就是这个22008的这是比较经典的，所以说那个2008年那个时候啊，就是它是最经典的，哎呀，我的妈呀。

---