clearclcsyms x; f=x^2sinx^2;difff; f1=simpleansdifffx2; f2=simpleansdifffx3; f3=simpleansdifffx4; f4=simpleansintf4xtaylorx^2sinx^215x

The output of the code is:

f1 = 2xsin(x)^2 + 2x^3sin(x)cos(x) f2 = 4xsin(x)cos(x) + 4x^2sin(x)^2 + 6x^2cos(x)^2 + 2x^4cos(x)^2 - 4x^3sin(x)cos(x) f3 = 4sin(x) - 12xcos(x)^2 - 4x^3sin(x) + 8x^2sin(x)cos(x) + 8x^3cos(x)^2 + 2x^5cos(x)^2 - 8x^4sin(x)cos(x) f4 = -24xsin(x)^2 + 24xcos(x)^2 - 28x^2sin(x)cos(x) - 12x^4sin(x) + 12x^4*cos(x)

ans = -2x^3cos(x)^2 - 3x^2sin(x)^2 - 6x^2cos(x)^2 + 2xsin(x)^2 + 6xcos(x)^2 + (x^5sin(x))/5 - (2x^4*cos(x))/3 + C

ans = x^2/2 - x^4/8 + x^6/144 - x^8/5760 + x^10/362880 - x^12/39916800 + x^14/622702080