>
> In function ex1ode, you will have to pass the source terms at position x to the calling program - so there should be no loop over the grid points in your code:
>
> function dydx=ex1ode(x,y)
> % The source terms (sco2,smea,Tg,TL) are calculated as functionS of y1-y4
> dydx=[sco2/ug
> Tg
> smea/uli
> TL];
>
> Best wishes
> Torsten.
------------------------------------------------
Removing the loop doesn't produce the desired results.
How do i pass the source terms at position x to the calling program? Please could you give me an idea what you mean?
Thank you for your help.