Метод Рунге-Кутты

[peace on you]

Люди разбирающиеся, помогите пожалуйста.
Есть такое уравнение : y'=2xy^3/(1-x^2y^2).
При таких то заданных НУ.
Пытаюсь решить методом Рунге-Кутты.
Проинтегрировав данное уравнение и подставив в код получаю ошибку(выделение синей строки f1 = (-y0+k1/2)*log(1-(x0+h/2)*(x0+h/2)*(y0+k1/2)*(y0+k1/2));)

#include <stdio.h>
#include <conio.h>
#include <math.h>
void main()
{
float h=0.05,x0=2,xk=2.5,y0=1,k1,k2,k3,k4,f0,f1,f2,f3,x,y;
printf("x0=%f y0=%f\n",x0,y0);
met1:{ x=x0+h;
f0=-y0*log(1-x0*x0*y0*y0);
k1=h*f0;
f1=(-y0+k1/2)*log(1-(x0+h/2)*(x0+h/2)*(y0+k1/2)*(y0+k1/2));
k2=h*f1;
f2=(-y0+k2/2)*log(1-(x0+h/2)*(x0+h/2)*(y0+k2/2)*(y0+k2/2));
k3=h*f2;
f3=(-y0+k3)*log(1-(x0+h)*(x0+h)*(y0+k3)*(y0+k3));
k4=h*f3;
y=-y0+(k1+2*k2+2*k3+k4)/6;
printf("x=%f y=%f\n",x,y);}
 
if (x<=xk) {x0=x, y0=y;goto met1;}
getch();}

[peace on you]

чуть переделал, ошибка сместилась до строчки к4=h*f3 :

#include <stdio.h>
#include <conio.h>
#include <math.h>
void main()
{
float h=0.05,x0=2,xk=2.5,y0=1,k1,k2,k3,k4,f0,f1,f2,f3,x,y;
printf("x0=%f y0=%f\n",x0,y0);
met1:{ x=x0+h;
f0=y0*log(x0*x0*y0*y0);
k1=h*f0;
f1=(y0+k1/2)*log(x0+h/2)*(x0+h/2)*(y0+k1/2)*(y0+k1/2);
k2=h*f1;
f2=(y0+k2/2)*log(x0+h/2)*(x0+h/2)*(y0+k2/2)*(y0+k2/2);
k3=h*f2;
f3=(y0+k3)*log(x0+h)*(x0+h)*(y0+k3)*(y0+k3);
k4=h*f3;
y=y0+(k1+2*k2+2*k3+k4)/6;
printf("x=%f y=%f\n",x,y);}
if (x<=xk) {x0=x, y0=y;goto met1;}
getch();}

[peace on you]

спасибо всем кто ответил, лабусдал.