Matlab - Metodo de Runge Kutta orden 4 para sistema de ecuación de orden superior. AYUDA

Buenas noches.
Necesito su ayuda para este mini proyecto que me encargaron, pero no entiendo 100% al parecer... Se me da una ecuación diferencial del comportamiento de un sistema y las condiciones de borde, y arreglandolo en un sistema de ecuaciones y mediante RK-4 debo hallar los resultados al comportamiento del sistema. (Anexo abajo una imagen detallada del problema)




Anexo aqui mi Pseudo-Codigo para la resolucion del problema, el cual por alguna razon me diverge el primer termino. Segun indicaciones, debo probar aleatoriamente la variable que indique como z01 para hallar el valor que cumpla las condiciones de borde, pero si sigue estando ese termino que diverge no puedo iterar varias veces para encontrarlo...

%Programa para resolucion de un problema mediante RK-4


clc

clear all

close all


%Se Definen los parametros de entrada

T0=200;

Tf=100;

h=0.01;

%Planteamos nuestro sistema de ecuaciones


dz=inline('((10^(-7))*(T+273)^4)-4*(150-T)','T');

dt=inline('z','z');


%Se hace la suposición del valor inicial de Z para aplicar RK-4


z01=150;


for x=0:h:0.5


k11=dz(T0);

k21=dz(T0+(k11*h)/2);

k31=dz(T0+(k21*h)/2);

k41=dz(T0+k31*h);



Z1=z01+(h/6)*(k11+k21+k31+k41);


k12=dt(z01);

k22=dt(z01+k12*h/2);

k32=dt(z01+k22*h/2);

k42=dt(z01+k32*h);


T1=T0+(h/6)*(k12+k22+k32+k42);


T0=T1;

z01=Z1;


end

T0
Z01

Cualquier ayuda sera bien recibida, Gracias.