function R=Rungkuta4(f,a,b,N,ya)h=(b-a)/N;x=zeros(1,N+1);y=zeros(1,N+1);x=a:h:b;y(1)=ya;for i=1:N k1=feval(f,x(i),y(i)); k2=feval(f,x(i)+h/2,y(i)+(h/2)*k1); k3=feval(f,x(i)+h/2,y(i)+(h/2)*k2); k4=feval(f,x(i)+h,y(i)+h*k3) y(i+1)=y(i)+(h/6)*(k1+2*k2+2*k3+k4);end这个是M文件,如果我要计算a=0,b=1,h=0.1,N=10,f=x^2-y,ya=1的情况,是需要再建立个M文件,还是在MATLAB命令窗直接输入,该怎么输?老是格式不对。。。。高手帮忙!先谢了~
|