Fast på en deluppgift av ett projekt. (b). Jag får inte till RK4 metoden och dessutom ingen plot. (ljudvågor)

såhär långt har jag kommit hittills:

b0= 12.3;
c=[5000; tand(b0)]

p1=-20.2090;
p2=17.3368;
p3= 272.9057;
p4=0.7528;

Z = @(z) 4800+ p1 + p2*(z/1000)+ p3*(exp(-p4*z/1000));
Zderv = @(z) p2/1000- (p3*p4*exp(-p4*z/1000)/1000);
q0 =(Z(5000)/cos(12.3).^2);
S= @(x,c) [c(2); -q0.*(Zderv(1)/Z(c(1)).^3)]
xtol= 30*6076;
h =10;

while x < xtol
k1= S(x,c);
k2=S.*[x+0.5.*h, c+0.5*h.*k1]
k3=S.*[x+0.5.*h, c+0.5*h*k2]
k4=S.*[x+h, c+h*k3]


end

skulle verkligen uppskatta lite hjälp!