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

Polynomials, write the zeros of underroot3power2 -8x+4underroot 3

write the zeros of underroot3power2 -8x+4underroot 3

Linear algebra, solve for k such that the system 4x+ky=6 kx+y=-3

solve for k such that the system 4x+ky=6 kx+y=-3

#title.simpal harmonic motion., #questionShow that the system oscillates in...

#questionShow that the system oscillates in simple harmonic motion demonstrated by; , for which the general solution where X = (x – x0)..

Graph for the sequence - sequences and series, Graph for the Sequence F...

Graph for the Sequence First we wish to think about the term graphing a sequence. To graph the sequence {a n } we plot the points {n, a n } as n ranges over every possible valu

Calculate the limit of f (-4), Let's take a look at one more example to ens...

Let's take a look at one more example to ensure that we've got all the ideas about limits down that we've looked at in the last couple of sections. Example: Given the below gr

Numerical methods for ordinary differential equationsordinay, #k1=f(Tn, Xn)...

#k1=f(Tn, Xn), k2=f (Tn + H.Y,Xn + H.Y.k1) Xn+1=Xn + H(a.k1+ b.k2) Find a relation between Y,a and b so that the method is second order consistent.

Area related to circle, If ABCD isaa square of side 6 cm find area of shad...

If ABCD isaa square of side 6 cm find area of shaded region

Evaluate the log function, Evaluate the log function: Calculate 3log 1...

Evaluate the log function: Calculate 3log 10 2. Solution: Rule 3.             log  (A n ) = nlog b   A 3log 10  2 = log 10 (2 3 ) = log 10   8 = 0.903

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