I have a problem with ODE 45 - Programmers Heaven

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Categories

Welcome to the new platform of Programmer's Heaven! We apologize for the inconvenience caused, if you visited us from a broken link of the previous version. The main reason to move to a new platform is to provide more effective and collaborative experience to you all. Please feel free to experience the new platform and use its exciting features. Contact us for any issue that you need to get clarified. We are more than happy to help you.

I have a problem with ODE 45

jee1141jee1141 Posts: 2Member
hey Guys
i have a problem with ODE i already wrote my program but it doesnt work :(. who can help me
i need to use the TEMP Q and R in every single steps of ode45 and say to it to do that with these new values in the next period of time and also i need to force ode45 to run the function just one time during the period and jump out of function and start with my new values

Global TEMP R Q
dt=0.00001;
time =0:dt:0.05;
gt=linspace(0,0.05,(0.05/dt)+1);
g = gt;

for k=2:length(time)

t=time(k)
i = 1;
[t,y] = ode45(@(t,y) F(t,y,gt,g),[time(k-1) :dt:time(k)],yy0(i,:),options);
yy0(i,:)=y(end,:)
i = i+1;
plot(time(k),real(y(end,1)));

end
----------------------------------------------------
Function fx = F(t,y,gt,g)
Global TEMP R Q
g=interp1(gt,g,t);
fx = zeros(4,1);

if y(1)-y(2)< 20
R=0
Q=0
TEMP=100
Else
Q = f(last(R),last(y1),last(y2))
TEMP = f(last(R),Q)
R = f(last(T))
End

fx=R.y+cos(g)

Sign In or Register to comment.