Conversion to double from sym is not possible.

syms t
td=5;
p0=3;
m=5;
w=4;
k=m*w^2;
t1=0:1:10;
for n = 1 : length(t1)
  if t1(n)< td;
    y(n) = p0*(t1(n)/td);
    else
        y(n) = p0;
  end
end
x=int(y.*sin(w*(t-t1))/(m*w),0,t);
plot(t,x);

??? Error using ==> plot
Conversion to double from sym is not
possible.

Error in ==> Untitled2 at 16
plot(t,x);

Answer 1

[Andrei Bobrov]

td = 5;
p0 = 3;
m = 5;
w = 4;

x = @(z)arrayfun(@(t)integral(@(t1)p0*min(t1,td)/td.*sin(w*(t-t1)),0,t)/m/w,z);

plot(0:.01:20,x(0:.01:20));