Question

write a mfile with matlab which takes a and n from user and gives the n term s of taylor series of sin(x) around x=a

Answer 1

i have difficulties with math part
@experimentX

Answer 2

we must print something like this\[a_0+a_1(x-a)+a_2(x-a)^2+...\]

Answer 3

man i have nothing :(

Answer 4

Woops ... !! probably write a function do do that.

Answer 5

syms x
diff(sin(x), 'n times differentiate', 'around a');

Answer 6

@phi

Answer 7

man would u mind to join us :)

Answer 8

Woops!! lol

Answer 9

lol

Answer 10

nth derivative of sin(x)\[\sin(\frac{n\pi}{2}+x)\]

Answer 11

wrong !!!

Answer 12

Woops ... i need to differentiate it and substitute it.

Answer 13

clc
syms x
center = input('Enter the center of expansion: ');
for i=1:10
disp(subs(diff(sin(x), i), x, center));
end

Answer 14

exper can we show that like what i mentioned above?

Answer 15

for example with center 3

-0.9900-0.1411(x-3)+0.9900(x-3)^2+...

Answer 16

hold on ... I'm working.

Answer 17

I don't like Matlab's string handing !!

Answer 18

me too

Answer 19

This does not do the input. But it does give an idea
a=pi/3, n=4; str2='';
for ii=0:(n-1)

switch mod(ii,4)
case 0, ff= sin(a);
case 1, ff= cos(a);
case 2, ff= -sin(a);
case 3, ff= -cos(a);
end;
if (ff<0) sg= '-'
else sg= '+';
end;
if ii>0
str= sprintf('%2s %7.5f(x-%2.4f)^%d',sg, abs(ff)/factorial(ii),a,ii);
else
str= sprintf('%7.5f',ff);
end;
str2= [str2 str];
end;
str2

Answer 20

its output
str2 =

0.86603 + 0.50000(x-1.0472)^1 - 0.43301(x-1.0472)^2 - 0.08333(x-1.0472)^3

Answer 21

I really hate matlab for string manipulation
clc
syms x
center = input('Enter the center of expansion: ');

for i=1:3
disp( [subs(diff(sin(x), i), x, center), strcat('(x-', num2str(center), ')^', num2str(i))]);
end

Answer 22

I got this nasty code ... looks pretty close
clc
syms x
center = input('Enter the center of expansion: ');
to_out = '';
for i=1:10
to_out = strcat(to_out, '+' , strcat(char(subs(diff(sin(x), i), x, center)), strcat('(x-', num2str(center), ')^', num2str(i))));
end
disp(to_out);

Answer 23

the out put is ...for r=4 ... maybe tweak a little.
+cos(3)(x-3)^1+-sin(3)(x-3)^2+-cos(3)(x-3)^3+sin(3)(x-3)^4

Answer 24

thank u guys