武汉理工大学数值分析欧拉

#include<iostream>
#include<cmath>
using namespace std;
double x[10],y[10];
double f(double a,double b)
{
  double c;
  c=b-a*2.0/b;
  return (c);   
}
int main()
{
double x0,y0,h,N,n,x1,y1,yp,yc;
cout<<"x0:";
cin>>x0;
cout<<"y0:";
cin>>y0;
cout<<"h:";
cin>>h;
cout<<"N:";
cin>>N;
n=1;
while(n!=N)
{
   x1=x0+h;
   y1=y0+h*f(x0,y0);
   cout<<"x1:"<<x1<<" "<<"y1:"<<y1<<endl;
   n=n+1;
   x0=x1;
   y0=y1;