阿当姆斯方法求解微分方程初值问题：四阶龙格库塔提供出发值（C语言）

#include<stdio.h>

float f1(float x,float y)  //第一问
{
    return y*y;//方程
}

float f2(float x,float y)//第二问
{
    return 0.1*(x*x*x+y*y);
}

//四阶龙格库塔提供出发值
void runge(float(*f)(float x,float y),float a,float b,float y0,int N,float yy[])
{
    float x=a,y=y0,K1,K2,K3,K4;
    float h=(b-a)/N;