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

Find the time required for an enlargement, 1. The polynomial G(x) = -0.006x...

1. The polynomial G(x) = -0.006x4 + 0.140x3 - 0.53x2 + 1.79x measures the concentration of a dye in the bloodstream x seconds after it is injected. Does the concentration increase

Creating input- output tables, How do I create an input-output table with t...

How do I create an input-output table with the rule -12

Iit-jee questions, can u tell me a website for iit-jee questions?

can u tell me a website for iit-jee questions?

E is irrational, If e were rational, then e = n/m for some positive integer...

If e were rational, then e = n/m for some positive integers m, n. So then 1/e = m/n. But the series expansion for 1/e is 1/e = 1 - 1/1! + 1/2! - 1/3! + ... Call the first n v

Geography, How do you find the maxima or minima on a parabolic graph?

How do you find the maxima or minima on a parabolic graph?

Continuous uniform distribution, Continuous Uniform Distribution Consid...

Continuous Uniform Distribution Consider the interest earned on a bank deposit. Let X equal the value after the decimal point. (Assume no rounding off to the nearest paise.) Fo

Geometry, in right angle triangle BAC.

in right angle triangle BAC.

Modeling , A plastic manufacturer has 1200 boxes of transparent wrap in sto...

A plastic manufacturer has 1200 boxes of transparent wrap in stock at one factory and 1000 boxes at his second factory.The manufacturer has order for this product from 3 different

Fundamental theorem of integral facts formulasproperties, Fundamental Theor...

Fundamental Theorem of Calculus, Part I If f(x) is continuous on [a,b] so, g(x) = a ∫ x f(t) dt is continuous on [a,b] and this is differentiable on (a, b) and as,

Launching of a new product, Launching a new product (Blackberry Cube) Analy...

Launching a new product (Blackberry Cube) Analysis (target market) Product features Promotions and advertisement sample design (location)

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