I have solved an ode using ode45. The solution is correct but when I try to use the solution on other functions the program only takes the initial conditions.

My code is attached. It runs but the temperature plot is wrong. Can someone tell me how to fix this?

[Walter Roberson]

Undefined function or variable ao in

Aavg=@ (t,a)  a(1)*A1+a(2)*A2+(ao(a))*A0; %Average planetary albedo

Reminder: we do not know what the temperature plot should look like, so it would help if you posted an example of what you are hoping for, or describe what difference you observe.

Answer 1

[Ezor]

That is strange that your ao function is defined after the Aavg one which uses a0. Same thing for vector a. Is it done on purpose?