分别用改进的欧拉法和四阶龙格-库塔公式求解微分方程初值问题(1) Y'=Y-2X/Y,Y(0)=1,X=[0,1],H=0.1(2) Y'=X2+Y2.Y(0)=0,X=[0,1] H=0.1 注释;X2为X的平方,Y2为Y的平方,然后用C语音编程的帮我编程一下
来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/18 01:30:54
分别用改进的欧拉法和四阶龙格-库塔公式求解微分方程初值问题(1) Y'=Y-2X/Y,Y(0)=1,X=[0,1],H=0.1(2) Y'=X2+Y2.Y(0)=0,X=[0,1] H=0.1 注释;X2为X的平方,Y2为Y的平方,然后用C语音编程的帮我编程一下
分别用改进的欧拉法和四阶龙格-库塔公式求解微分方程初值问题
(1) Y'=Y-2X/Y,Y(0)=1,X=[0,1],H=0.1
(2) Y'=X2+Y2.Y(0)=0,X=[0,1] H=0.1 注释;X2为X的平方,Y2为Y的平方,
然后用C语音编程的帮我编程一下
分别用改进的欧拉法和四阶龙格-库塔公式求解微分方程初值问题(1) Y'=Y-2X/Y,Y(0)=1,X=[0,1],H=0.1(2) Y'=X2+Y2.Y(0)=0,X=[0,1] H=0.1 注释;X2为X的平方,Y2为Y的平方,然后用C语音编程的帮我编程一下
求解就自己算吧,程序如下,其中RK4与改进欧拉法分别写成一个模块,方便阅读:
#include
#include
double f1(double x,double y) {//定义方程1
\x05return (y-2*x/y);
}
double f2(double x,double y) {//定义方程2
\x05return (x*x+y*y);
}
double RK4(double xn,double yn) {
\x05double h=0.1,k1=0,k2=0,k3=0,k4=0,d=1,h2=h/2;//RK4解方程2
\x05printf(" x2 y2\n");
\x05printf("%3.1lf%10.6lf\n",xn,yn);
\x05while(xn