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

What is the smallest possible number 3, What is the smallest possible numbe...

What is the smallest possible number in which can be created along with four decimal places using the numbers 3, 5, 6, and 8? Place the smallest number in the largest place val

How many pounds should nicole put in every basket, Nicole is forming 20 gif...

Nicole is forming 20 gift baskets. She has 15 pounds of chocolates to distribute equally between the baskets. If each basket gets the similar amount of chocolates, how many pounds

Complex number, If z=re i ? ,find the value of |e iz | Solution)   z=r(c...

If z=re i ? ,find the value of |e iz | Solution)   z=r(cos1+isin1) |e iz |=|e ir(cos1+isin1) |=|e -rsin1 |=e -rsin1

GEOMETRY, DIFFERENCE BETWEEN RIGHT ANGLE AND SCALENE

DIFFERENCE BETWEEN RIGHT ANGLE AND SCALENE

Determine the actual viewing area, Computer monitors are calculated by thei...

Computer monitors are calculated by their diagonals. If a monitor is advertised to be 19 in, Determine the actual viewing area, considerthe screen is square? (Round to the nearest

Proof of limit comparison test - sequences and series, Proof of Limit Compa...

Proof of Limit Comparison Test As 0  Now, as   we know that for large enough n the quotient a n /b n should be close to c and thus there must be a positive integer

.gradient, Draw the graph of y=x^2-4x from x=-1 to x=5.use the scale of 2cm...

Draw the graph of y=x^2-4x from x=-1 to x=5.use the scale of 2cm on the x axis and 1cm on the y axis.Estimate the gradient at point:x=4, x=2 and x=0

Diffrence between rational and irrational numbers, Q. Diffrence between Rat...

Q. Diffrence between Rational and Irrational Numbers? Ans. A number which is not rational is called irrational. The word "irrational" sounds not quite right...as though th

Laura paid $17 for jeans what was original price of jeans, Laura paid $17 f...

Laura paid $17 for a pair of jeans. The ticketed price was 20% off the original price plus the sign on the rack said, "Take an additional 15% off the ticketed price." What was the

Trigonometry, sin^2alpha *sec^2beta +tan^2 beta *cos^2alpha=sin^2alpha+tan^...

sin^2alpha *sec^2beta +tan^2 beta *cos^2alpha=sin^2alpha+tan^2 beta

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