ZKX's LAB

上古卷轴5爱欧拉 爱欧拉代码

2021-03-08知识4

谁能提供一下欧拉方法的matlab代码 %步长h=0.01时,误差会减小。function eulerclc;clear;h=0.1;x0=1;y0=0;k=0;x=x0;y=y0;xfinal=1.9;f=inline('1-y')while xk=k+1;y=y+h*feval(f,y);Y(:,k)=y;x=x+h;endt=x0:h:xfinalY=[0,Y]y=dsolve('Dy=1-y','y(1)=0');解析解,即真值y=subs(y)error_ab=Y-yerror_re=(Y-y)./y结果:t=1.0000 1.1000 1.2000 1.3000 1.4000 1.5000 1.6000 1.7000 1.8000 1.9000Y=0 0.1000 0.1900 0.2710 0.3439 0.4095 0.4686 0.5217 0.5695 0.6126y=0.0000 0.0952 0.1813 0.2592 0.3297 0.3935 0.4512 0.5034 0.5507 0.5934error_ab=0.0000 0.0048 0.0087 0.0118 0.0142 0.0160 0.0174 0.0183 0.0189 0.0191error_re=1.0000 0.0508 0.0482 0.0456 0.0431 0.0408 0.0385 0.0363 0.0343 0.0323

关于matlab function[x_star,index,it]=Newton(fun,x,ep,it_max)%求解非线性方程的Newton法,其中,%fun(x)为需要求根的函数;第一个分量是函数值,第2个分量是导数值;x为初始点;ep为精度,。

谁能告诉上古卷轴5爱欧拉的代码 谁告诉我上古卷轴5爱欧拉的代码 带着爱欧拉(吃人的那个)进山洞 自己跑路没管她 等我传出天际发现她丢了 自动存档自己顶来顶去顶没了 这丫的。

欧拉公式推导及实现代码,欧拉方法的基本思想是在小区间上用数值微分的前差公式代替方程左端的导数,右端函数fx,yx中x取中的某一点。

#爱欧拉代码

随机阅读

qrcode
访问手机版