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

Surface area of prisms , Can you help me find out how to find the surface a...

Can you help me find out how to find the surface area of a prism

L''hospital''s rule, L'Hospital's Rule Assume that we have one of the g...

L'Hospital's Rule Assume that we have one of the given cases, where a is any real number, infinity or negative infinity.  In these cases we have, Therefore, L'H

Series - convergence or divergence, Series - Convergence/Divergence In ...

Series - Convergence/Divergence In the earlier section we spent some time getting familiar with series and we briefly explained convergence and divergence.  Previous to worryin

Area under curve, w/ You could use this sample code to test your C function...

w/ You could use this sample code to test your C functions // Please make appropriate changes to use this for C++. // Following main function contains 3 representative test cases

Harold used a 3% iodine solution and a 20% iodine solution, Harold used a 3...

Harold used a 3% iodine solution and a 20% iodine solution to make a 95- ounce solution in which was 19% iodine. How many ounces of the 3% iodine solution did he use? Let x = t

Solving a quadratic equation, In polynomials you have seen expressi...

In polynomials you have seen expressions of the form x 2 + 3x - 4. Also we know that when an expression is equated to zero or some other expression, we cal

Quadratic equation, find a quadratic equation whose roots are q+1/2 and 2p-...

find a quadratic equation whose roots are q+1/2 and 2p-1 with p+q=1

Complex numbers, Complex Numbers In the radicals section we noted that...

Complex Numbers In the radicals section we noted that we won't get a real number out of a square root of a negative number.  For example √-9 isn't a real number as there is no

Real numbers on every line, Make a file called "testtan.dat" which has 2 li...

Make a file called "testtan.dat" which has 2 lines, with 3 real numbers on every line (some negative, some positive, in the range from-1 to 3).  The file can be formed from the edi

Polynomials, find a quadratic polynomial whose zeroes are 2 and -6.verify t...

find a quadratic polynomial whose zeroes are 2 and -6.verify the relationship between the coefficients and zeroes of the polynomial

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