04.cpp

来自「龙格-库塔方法的使用..这是计算方法的一个实验..希望对大家有帮助」· C++ 代码 · 共 42 行

#include <stdio.h>
#include <iostream.h>
#include <math.h>

float f(float x,float y)
{
	float f,f2,f3;
	//f=y-2*x/y;
    f=y;
	return f;
}

void main()
{   
	int n,N;
	float x0,y0,h,x1,y1,k1,k2,k3,k4;
    cout<<"************************欢迎使用(龙格-库塔)方法****************************"<<endl<<endl;
	cout<<"请输入x0的值：";
	cin>>x0;
    cout<<"请输入y0的值：";
	cin>>y0;
	cout<<"请输入h的值：";
	cin>>h;
	cout<<"请输入n的值：";
	cin>>N;
	n=1;
	do
	{
		x1=x0+h;
		k1=f(x0,y0);
		k2=f(x0+h/2,y0+h/2*k1);
		k3=f(x0+h/2,y0+h/2*k2);
		k4=f(x0+h,y0+h*k3);
		y1=y0+h/6*(k1+2*k2+2*k3+k4);
		cout<<"它的结果是：x1="<<x1<<" y1="<<y1<<endl;
	    n++;
		x0=x1;
		y0=y1;
	}
	while(n!=N);
    cin>>n;
}