#include<stdio.h> float func(float x,float y); int main() { float m1,m2,m3,m4,m,h=0.3; float x0 = 0.3, y0 = 5, xn; scanf("%f",&xn); while(x0<xn) { m1=func(x0,y0); m2=func((x0+h/2.0),(y0+m1*h/2)); m3=func((x0+h/2.0),(y0+m2*h/2)); m4=func((x0+h),(y0+m3*h)); m=((m1+2*m2+2*m3+m4)/6); y0=y0+m*h; x0=x0+h; } printf("y=%f",y0); // Final output return 0; } float func(float x,float y) { float m; m=(x*(x+1)-3*y*y*y)/10; return m; }