降圧チョッパ回路における負荷電圧波形が見たいです

Question

前提・実現したいこと

降圧チョッパ回路の負荷電圧波形のプロットデータをC言語で出力したいです。
数値計算は4次のルンゲクッタ法を用いて、算出したいです。

発生している問題・エラーメッセージ

中々、求める波形が出力されないです。
PSIM上で同じ回路を同様の規格でシミュレーションしたものと同じにならないです。
間違えている箇所が分からないので、指摘して頂けると嬉しいです。

該当のソースコード

#include<stdio.h>
#include<math.h>

double func_ON_1(double IL,double VRout) {

	double C = 1 * pow(10, -6);
	double R = 2.27* pow(10, -3);
	//double L = 0.00001;

	double U = IL;
	double Vo = VRout;

	long double A = ((-1 * pow( R * C ,-1)) * Vo) + ( U * pow( C ,-1));
	printf("\n\n A = %lf\n", A);

	return A;

}

double func_ON_2( double VRout) {


	double L = 0.01;// *pow(10, -3);

	double Vo = VRout;

	long double B = ((200 - Vo) * 100);// pow(L, -1));
	printf("\n\n B = %lf\n", B);

	return B;

}

double func_OFF_2( double VRout) {


	double L = 0.01;// *pow(10, -3);

	double Vo = VRout;

	long double C = ((-1 * Vo) * 100);// pow(L, -1));
	printf("\n\n A = %lf\n", C);

	return C;

}

int main(void) {

	double VRout, t, step,IL, end,T;
	double k1, k2, k3, k4;
	double l1, l2, l3, l4;
	int error,n;
	FILE *fi;

	//R = 2.27*(0.1)*(0.1)*(0.1);
	//C = 0.0005;
	n = 0;
	T = 0.0;
	step = 0.00002;
	VRout = 0.0;
	IL = 0.0;
	end = 20;

	if ((error = fopen_s(&fi, "DCDC_Curcuit_Simulation_Vr.txt", "w")) != 0) {

		printf("\nThis file can't opened \n\n");

		return 0;

	}

	for (t = 0.0; t <= end; t = t + step) {

		//fprintf(fi, "%f %f\n", t, VRout);
		n = (int)(t / 0.002);
		T = (double)(t - (n * 0.002));

		if ((T > 0) && (T < 0.0015)) {

			k1 = step * func_ON_1(IL, VRout);
			printf("\n k1 = %lf\n", k1);

			l1 = step * func_ON_2( VRout);
			printf("\n l1 = %lf\n", l1);

			k2 = step * func_ON_1(IL +( step * (l1*0.5)), VRout + (step * (k1*0.5)));
			//k2 = step * func_ON_1(IL, VRout + (k1*0.5));
			printf("\n k2 = %lf\n", k2);

			l2 = step * func_ON_2( VRout + step * (l1*0.5));
			printf("\n l2 = %lf\n", l2);
			
			k3 = step * func_ON_1(IL +( step * (l2*0.5)), VRout + (step * (k2*0.5)));
			//k3 = step * func_ON_1(IL, VRout + (k2*0.5));
			printf("\n k3 = %lf\n", k3);

			l3 = step * func_ON_2( VRout + ( step * (l2*0.5)));
			printf("\n l3 = %lf\n", l3);
			
			k4 = step * func_ON_1(IL + k3, VRout + k3);
			//k4 = step * func_ON_1(IL, VRout + k3);
			printf("\n k4 = %lf\n", k4);
			
			l4 = step * func_ON_2( VRout + l3);
			printf("\n l4 = %lf\n", l4);

			VRout = VRout + (step *((k1 + (2.0*k2) + (2.0*k3) + k4) / 6.0));
			IL = IL + (step * ((l1 + (2.0*l2) + (2.0*l3) + l4) / 6.0));

			fprintf(fi, "%f %f\n", t, VRout);

			k1 = 0, l1 = 0, k2 = 0, l2 = 0, k3 = 0, l3 = 0, k4 = 0, l4 = 0;

		}
		else {

			k1 = step * func_ON_1(IL, VRout);
			printf("\n k1 = %lf\n", k1);

			l1 = step * func_OFF_2(VRout);
			printf("\n l1 = %lf\n", l1);
			
			k2 = step * func_ON_1(IL + (step * (l1*0.5)), VRout + (step * (k1*0.5)));
			//k2 = step * func_ON_1(IL, VRout + (k1*0.5));
			printf("\n k2 = %lf\n", k2);
			
			l2 = step * func_OFF_2(VRout + (step *(l1*0.5)));
			printf("\n l2 = %lf\n", l2);
			
			k3 = step * func_ON_1(IL + (step * (l2*0.5)), VRout + (step *(k2*0.5)));
			//k3 = step * func_ON_1(IL, VRout + (k2*0.5));
			
			printf("\n k3 = %lf\n", k3);
			l3 = step * func_OFF_2(VRout + (step *(l2*0.5)));
			printf("\n l3 = %lf\n", l3);
			
			k4 = step * func_ON_1(IL + k3, VRout + (step * k3));
			//k4 = step * func_ON_1(IL, VRout + k3);
			printf("\n k4 = %lf\n", k4);

			l4 = step * func_OFF_2(VRout + (step * l3));
			printf("\n l4 = %lf\n", l4);

			VRout = VRout + (step *((k1 + (2.0*k2) + (2.0*k3) + k4) / 6.0));
			IL = IL + (step * ((l1 + (2.0*l2) + (2.0*l3) + l4) / 6.0));

			fprintf(fi, "%f %f\n", t, VRout);

			k1 = 0, l1 = 0, k2 = 0, l2 = 0, k3 = 0, l3 = 0, k4 = 0, l4 = 0;


		}


	}

	fclose(fi);

	return 0;
}

試したこと

LCの値を変化させたのですが、変わりませんでした。
。