hello,
Thanks for your reply,but I tried the same yesterday which did not yield any result, I will paste here the code
Fo=@ (t,i) F(i).*sin((2*pi()/T)*t);
Fo(1,2)
tspan=[0:100];
u0=zeros(202,1);
dudt0=zeros(202,1);
y0=[u0;dudt0];
[t,y]=ode45(@f,tspan,y0);
could you please give...