I computed my results using ode45. The main function m file is given below.
function qprime=memcapfinal(t,q);
e = 1.6*10^(-19);
h = 6.63*10^(-34);
phi = (0.3)*1.6*10^(-19);
R=1;
Epsilon0 = 8.85*10^-12;
k1 = 100;
k2 = 10;
d=10*10^(-9);
s = 6*10^(-9);
m = 9.11*10^(-31);
A = 10^(-4);
b = -4*3.14*s/h;
c0=(A*Epsilon0*k1)/d;
f=10000;
v0=7.5;
T=5/f;
v=v0*sin(2*3.14*f*t);
qprime=[((v(:)/R)-(1/R)*(((d*q(1))+(s*q(2)))/(c0*d)));(sqrt(2*m*phi)/s*(e/h)^2*(((q(1)+q(2))/(2*A*Epsilon0*k2))*s)*exp(b*sqrt(2*m*phi)))*A];
AND calling m file is give as follows.
clear all;
clc;
f1=10000;
v01=7.5;
t1=0:0.0000001:0.0005;
v1=v01*sin(2*3.14*f1*t1);
q0=[10^(-20) 10^(-20)];
tspan1=[0,0.0005];
[t,q]=ode45(@memcapfinal,tspan1,q0);
subplot(2,2,2)
plot(t,q(:,1))
subplot(2,2,3)
plot(t,q(:,2))
subplot(2,2,4)
plot(t,v)
Error: ??? Undefined function or variable 'v'.
Error in ==> ode45call at 21 plot(t,v)
Comment : Although I have defined the v in memcapfinal.m file. Moreover I am getting the length(t) of t as 133 which I don't why is it so? My goal is to plot q(1) vs v and this q(1) is calculated by ode45. for that both must have the same lengths but q(1) length is 1 and v is undefined.