untitled3.m

4阶龙格--库塔法解常微分方程,比较好用,希望大家尝试.

close all
clear all
f=@f2;
a=0;
b=1;
N=10;
ya=1;
R=Rungkuta4(f,a,b,N,ya);
y=solve(f2(a:(b-a)/N:b));
m=[R,y']