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

Sqares, Recently I had an insight regarding the difference between squares ...

Recently I had an insight regarding the difference between squares of sequential whole numbers and the sum of those two whole numbers. I quickly realized the following: x + (x+1)

Graph y = cos ( x ) - common graph, Graph y = cos (x) Solution: There ...

Graph y = cos (x) Solution: There actually isn't a whole lot to this one.  Given the graph for -4 ? ≤ x ≤ 4 ? . Note that we can put all values of x in cosine (that wo

Multiplication of two like terms with opposite signs, The product of -7ab a...

The product of -7ab and +3ab is (-7 x 3) a 2  b 2  = -21a 2  b 2 . In other words, a term with minus sign when multiplied with a term having a positive sign, gives a product having

Find inverse laplace transform, Question: Find Inverse Laplace Transfor...

Question: Find Inverse Laplace Transform of the following (a) F(s) = (s-1)/(2s 2 +8s+13)     (b) F(s)= e -4s /(s 2 +1) + (1/s 3 )

Limit comparison test - sequences and series, Limit Comparison Test Ass...

Limit Comparison Test Assume that we have two series ∑a n and ∑b n with a n , b n   ≥ 0 for all n. Determine, If c is positive (i.e. c > 0 ) and is finite (i.e. c

Standardizing normal variables, Standardizing Normal Variables Suppose ...

Standardizing Normal Variables Suppose we have a normal population. We can represent it by a normal variable X. Further, we can convert any value of X into a corresponding valu

Set builder notation, A={2,3,5,7,11} B={1,3,5,7,9} C={10,20,30,40,......100...

A={2,3,5,7,11} B={1,3,5,7,9} C={10,20,30,40,......100} D={8,16,24,32,40} E={W,O,R,K} F={Red,Blue,Green} G={March,May} H={Jose,John,Joshua,Javier} I={3,6,9,12,15}

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