Writing a sin function program dealing with factorials, help in editing and error?
15 次查看(过去 30 天)
显示 更早的评论
I am trying to write a specific program involving the sine function and factorials.
The problem states that:
sin(x) = x - x^3 / 3! + x^5 / 5! - x^7 / 7! + x^9 / 9!
I am trying to get the error rate very low, basically having the number of terms it takes to get enough accuracy by 10^-5, so that:
abs[sin(x) - sk(x)] / abs[sin(x)] = 10^-5
I know the sk equation should be written as:
sk = x-x^3 / 3! + .......x^(2k-1) / (2k-1)!
I am starting out my program with:
for s = 1: maxs
nsteps(s) = (ones(s)
end
I am confused on how to add on to my code to create my desired result, I know intuitively that when plotting this, the #'s should increase as x increases.
[Merged information from duplicate question]
I am trying to write a specific program involving the sine function and factorials.
The problem states that:
sin(x) = x - x^3 / 3! + x^5 / 5! - x^7 / 7! + x^9 / 9!
I am trying to get the error rate very low, basically having the number of terms it takes to get enough accuracy by 10^-5, so that:
abs[sin(x) - sk(x)] / abs[sin(x)] = 10^-5
I know the sk equation should be written as:
sk = x-x^3 / 3! + .......x^(2k-1) / (2k-1)!
I am starting out my program with:
for s = 1: maxs
nsteps(s) = (ones(s)
end
I am confused on how to add on to my code to create my desired result, I know intuitively that when plotting this, the #'s should increase as x increases.
Here is my code:
function [s] = fact(t)
s = 1:max(s);
for s = 1:max(s)
nsteps(s) = ones(s);
end
for t = (1:2:2*k-1);
(-1)^(0:k-1) .* x^t ./ factorial(t);
end
I receive this error when performing it:
??? Undefined function or variable 's'.
??? Undefined function or variable "s".
Error in ==> fact at 2
s = 1:max(s);
HELP PLEASE!
1 个评论
回答(2 个)
Walter Roberson
2012-4-23
t = (1:2:2*k-1);
(-1).^(0:k-1) .* x.^t ./ factorial(t)
Or something like that.
(Note: corrected -- ^ operators had to be .^ )
7 个评论
Walter Roberson
2012-4-23
"=" is not permitted in a "while". "while" has to be given a logical test.
Warning: You initialize s to 0, so max(s) is 0, and 1:0 is the empty vector, so your "for" loop will not do any iterations.
What is the purpose of your nsteps vector? Why are you trying to initialize it to a square matrix such as
1 1 1
1 1 1
1 1 1
(i.e., ones(3))
Richard Brown
2012-4-23
I take it this is a programming exercise, not a straight mathematics question - you can work out the number of terms pretty easily using Taylor's theorem with remainder ...
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!