Hi Godas,
The solution to this problem can be achieved in 2 ways:
If you have Symbolic Toolbox, you can do this easily:
>> syms t;
>> f = t*sin(t);
>> result = int( int(f), [0 1])
>> vpa(result)
If you do not have symbolic toolbox, this can be done in following way:
>> format long
>> f = @(x,t) t.*sin(t)
>> y2 = @(x) x;
>> result = integral2(f2,0,1,0,y2);
Important Note: In second method, as you can see that I have put the lower limit of inside integral as 0 which is working for this case but if the result of first integral yields 
then it will become 1. Similarly, if it yields 
, then you have to use 
. Keep this in mind.
then it will become 1. Similarly, if it yields , then you have to use . Keep this in mind.Hope it helps
