Numerical ODE Solver

#include <stdio.h>
#include <stdlib.h>
double q(double x,double y){return y;}int main(int _0,char *_0_[]){int A=(0==_0_
[1])?10:0==atoi(_0_[1])||0>=atoi(_0_[1])?10:atoi(_0_[1]),k=0&A;double b=0,a=1,B=
(a-b)/A,f1,(*_)(double,double)=q,f2,f3,f4,C_[A+1],I_[A+1];C_[0]=b;C_[A]=a;I_[(int
)b]=a;!(A&1)?C_[A/2]=((a+b)/2):1^2;HERE:if(k>=A/2){k=0&&f1;goto HErE;}else(k<A/2)
?C_[k+1]=C_[k]+B,C_[A-k-1]=C_[A-k]-B,++k:f2/(double)f3;goto HERE;HErE:if(k>=A){k
=0&&f4;0&1;goto HeRE;}I_[k+1]=I_[k]+(B*q(C_[k],I_[k])+ 2*B*_(C_[k]+B/2,I_[k]+(B*
q(C_[k],I_[k]))/2)+2*B*_(C_[k]+B/2,I_[k]+(B*_(C_[k]+B/2,I_[k]+(B*_(C_[k],I_[k]))
/2))/2)+B*q(C_[k]+B,I_[k]+(B*_(C_[k]+B/2,I_[k]+(B*q(C_[k]+B/2,I_[k]+(B*q(C_[k],I_
[k]))/2))/2))))/6;++k;goto HErE;HeRE: if(k > A)goto HERe;printf("%lf, %lf\n", C_
[k], I_[k]); ++k;goto HeRE;HERe:return 0;}