Need help solving this problem with newtons method

[James Tursa]

Newton's Method or Euler's Method?

[Rena Berman]

(Answers Dev) Restored edit

Answer 1

[James Tursa]

To match how you are calling it, you need to switch the order of the arguments in this function handle:

F = @(p, t) k*p*(L-p)-a*sin((2*pi*t)/b);

to this

F = @(t, p) k*p*(L-p)-a*sin((2*pi*t)/b);

Also, you should strongly consider decreasing your step size by quite a bit. Maybe you could experiment with N in the 100's to see what works well. E.g., start with N = 100, then N = 200, etc. until you get answers that seem to be stable.

[T]

Thanks for the help!!