Решить дифференциальное уравнение методом Рунге-Кутта 4 порядка (ode45)

e=15;
V1=10;
t=(1:0.1:30);
V2=1;
f=10000;
omega=2*pi*f;
k=1.36*10-23;
Lk=0.000001;
R=1600;
T=222;
B=222;
Vc=0.92837;
R=1;
Vn=sqrt(4*k*R*T*B);
function q(2)=solverDE(t,y);
f=@(t,q)([q(2); V1+V2*sin*(omega*t)-Lk*q(2)-R*q(1)-Vn-Vc*sin((2*pi*q(1))/2*e)/Lk]);

function tempq 
[T,Q] = ode45(@loreq, [0,3], [0 0 0]); 
plot(T,Q) 
legend('q(t)', 'I(t)','I2(t)' )
end 
 
function dq = loreq(t,q) 
R = 2; 
V=4; 
dq=zeros(3,1) 
dq(1)=V/R;
dq(2)=q(2);
dq(3)=q(3);
end

dq(1)=q(2)
dq(2)=V/R ;

function temp
[T, Q] = ode45(@loreq, [0, 3],[0]); 
dq = loreq(T,Q); % считаем первую производную
plot(T,Q,'b', T,dq,'.r')
legend('q(t)','dq/dt')
end
 
 
function dq=loreq(t, q) 
R=2; V=4; 
dq = V/R; 
end

function dx=loreq(t, x) 
dx=x+3;
end

function dx=loreq(t, x) 
dx = t + 3;
end

function tempq 
[ T , Y ]=ode45(@loreq , [0,3] , [0 1] ) ; 
ds=loreq(T,Y)
plot (T,Y,Y,ds) 
legend ('funkciya y(x)', 'y''(x)') 
end 
 
function dy=loreq(t , y, x) 
x=5;
dy=zeros(2,1)
dy(1)=y(2);
dy(2)=y(1)+x+1;
end

ds=loreq(T,Y)

function tempq 
[ X , Y]=ode45(@loreq , [0,3] , [0 1] ) ; 
ds=loreq(Y,X);
dq=loreq(Y,ds);
plot (Y,dq);
legend ('funkciya y(x)') 
end 
 
function dy=loreq(x , y) 
dy=[y(2); 5]
end

function tempq
[T, Y] = ode45(@loreq, [0,3], [0 1] ) ;
 
for i = 1:length(T) % считаем производные
    DY(i,:) = loreq(T(i),Y(i,:));
end
 
plot(T,Y(:,1),T,DY)
legend('y(x)','dy/dx','d^2y/dx^2')
end
 
function dy = loreq(x, y)
dy = zeros(2,1);
dy(1) = y(2);
dy(2) = y(1)+x+1;
end

function tempqwertru
[T, Q] = ode45(@zod, [0,5.7*10^-12],[1 0]);
 
for i = 1:length(T)
    DQ(i,:) = zod(T(i),Q(i,:));
end
subplot(4,1,1); plot(T,Q(:,1));
legend('q(t)')
subplot(4,1,4); plot(T, V(T));
legend('I(V)')
subplot(4,1,2); plot(T,DQ(:,1));
legend('dqdt(t)')
subplot(4,1,3); plot(T,DQ(:,2));
legend('d2qdt2(t)')
end
 
function dq = zod(t,q)
L=1*10^-7; 
R=1*10^6;
dq=zeros(2,1);
dq(1)=q(2);
dq(2)=1/L*(V(t)-R*q(1));
end
 
function y = V(t) 
V=1*10^-6;
y = V*t;
end