Initial value problems, math, Mathematics

Assignment Help:
Write a Matlab function MyIVP that solves an initial-value problem (IVP) for a system of ordinary
differential equations (ODEs) of the form
x ?(t) = f (t, x(t)),
where f : R × Rn ? Rn is an arbitrary function with one one-dimensional input (for time t) and one n-dimensional input x, and n-dimensional output. The function should implement a Runge-Kutta formula (for example, the rk34 formula or the Dormand & Prince formula).
The initial value x0 is provided by the user of MyIVP. The first line of MyIVP (saved in a file MyIVP.m) should look like this
function [xend,t,xt]=MyIVP(f,x0,tspan,N) Inputs
• f: function defining the right-hand side of the ODE. f should accept two arguments: t (a number) and x (an n-dimensional vector). The function f should return an n-dimensional vector y (the time derivative). Typical calling sequence: y=f(t,x), returning the value of f at time t in position x.
• x0: initial value where integration starts from (n-dimensional vector).
• tspan: Starting time and end time for integration. Integration has to run from time t =tspan(1)
to time t =tspan(2).
• N: number of steps for integration. The integration stepsize h=(tspan(2)-tspan(1))/N should
be small.
Outputs
• xend: result of integration at t =tspan(2).
• t: vector of times at which intermediate values have been computed (this should have N + 1
entries).
• xt: intermediate values (n × (N + 1)-array). xt(:,k) should be the solution at t(k).
You can check the built-in variable nargout inside your function to see if the user wants to get three outputs or only the end value xend. If nargout==1 you don’t need to store the intermediate values.

Related Discussions:- Initial value problems, math

Calculate the price of the horseracing track, There are five horseracing tr...

There are five horseracing tracks in Kentucky. The Kentucky legislature allows only one track to be open at a time. How does this restriction affect the price the track can charge

DETERMINANT, IF 7 AND 2 ARE TWO ROOTS OF THE EQUATION |X 3 7 2 X 2 7 6 X...

IF 7 AND 2 ARE TWO ROOTS OF THE EQUATION |X 3 7 2 X 2 7 6 X |=0 THEN FIND THE THIRD ROOT IS

What is the continuously compounded forward rate, At time t an investor s...

At time t an investor shorts a $1 face value zero coupon bond that matures at time T = t and uses the entire proceeds to purchase a zero coupon bond that matures at time

Queuing Theory, A telephone exchange has two long distance operators.The te...

A telephone exchange has two long distance operators.The telephone company find that during the peak load,long distance calls arrive in a poisson fashion at an average rate of 15 p

Division of complex number, Division of complex number Now, we gave thi...

Division of complex number Now, we gave this formula a long with the comment that it will be convenient while it came to dividing complex numbers so let's look at a couple of e

Find out the radius of convergence, Example: Find out the radius of conver...

Example: Find out the radius of convergence for the following power series. Solution : Therefore, in this case we have, a n = ((-3) n )/(n7 n+1 )   a n+1 = (

Nun, how do you identify area ??

how do you identify area ??

Activity example of one to one correspondence learning, Devise one activity...

Devise one activity each to help the child understand 'as many as' and 'one-to-one correspondence'. Try them out on a child/children in your neighbourhood, and record your observat

What is the probability in which she gets heads further flip, Jennifer ?ipp...

Jennifer ?ipped a coin three times and got heads each time. What is the probability in which she gets heads on the further ?ip? The probability of heads does not modify based o

Write Your Message!

Captcha
Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd