0
|
5602 Views
|
5 Replies
|
0 Total Likes
View groups...
Share
GROUPS:

# Iterative method for a simple function

Posted 9 years ago
 Take the function V = 10 (1-Exp[-0.1 t])  I want to program a code that will give me this function as  t=t+tau V=V+D[V,t] tau  via an iterative method for say tau=0.01 and I=1 to 1000. This seems simple but I have so far had little success. I would like to be able to plot the results of this computation.
5 Replies
Sort By:
Posted 9 years ago
 I just tried your code and found it to work well. Thanks.,Ray
Posted 9 years ago
 Hi Bill, I will try your code .Thanks again. I have produced a code that does work. It goes as follows: Clear[t,V], V=10 (1-Exp[-.1 t]) ,D[V,t]=Exp[-.1 t],t=0,V=0,points={t,V} then  points=Table [t=t+.01; V=V+Exp[-.1 t] .01,{I,1,5000}], ListPlot[points] My success has probably given me some false confidence! I am happy to try your code as well, Regards, Ray
Posted 9 years ago
 See if you can put this into your Mathematica and run it without error. Clear[t]; tau = .1; V[t_] := 10 (1 - Exp[-0.1 t]); dV = D[V[t], t]; v = 0; points = Table[{time, v = v + dV /. t -> time}, {time, 0., 10, tau}]; ListPlot[points] I may not have understood exactly the function you wished to create. If this code does not do what you need then perhaps you can modify it or you can try to explain your question in another way and I will see if I can modify this for you.
Posted 9 years ago
 Thanks for reply. I put in what I thought to be your code and it did not work. Dear Bill,I put in what I thought to be your code suggestion and found that this did not work. I do not know much about programming and would have to see a detailed code that works, Then I could analyze and try to understand it. If you can provide this, I would be grateful. If not, thanks anyway,,Ray
Posted 9 years ago
 It seems like you are using several of your variable names for more than one thing. That may confuse Mathematica.Can you see if you can understand this idea and perhaps adapt it to what you want to do? Clear[t];(*Make sure t has no assigned value*) V[t_] := 10 (1 - Exp[-0.1 t]);(*Define V function of t*) dV = D[V[t], t];(*Define derivative of V at any t*) v = 0;(*Start the iteration value with v=0*) (*Now create a list of {time,v} values where t is replaced with time in dV*) points = Table[{time, v = v + dV /. t -> time}, {time, 0., 10, tau}] That tries to use the function V[t] for only that purpose, the derivative dV for only that purpose, the iteration time for only that purpose and the iteration value for only that purpose.Then you can plot the result using this ListPlot[points] If I have misunderstood and you can explain then perhaps I will be able to correct this for you