matlab代写代考

clear all Npoints=50; alpha = 5; %BC1 y(x=-1)=5 beta = 3; %BC2 y(x=1)=3 %xpoints=linspace(-1,1,Npoints)’; [xpoints,w,P]=lglnodes(Npoints); xpoints=-xpoints; tpoints=(xpoints+3)/2; y = @(t) t+4./t.^2; yexact = y(tpoints); n=length(xpoints)-1; %order of polynomial D=DerivMatrix(xpoints,n); D2 = D*D; A=D2+diag(2./(xpoints+3))*D-diag(2./(xpoints+3).^2); Amatrix=A(2:end-1,2:end-1); Cvector=-alpha*A(2:end-1,1)-beta*A(2:end-1,end); sol = Amatrix\Cvector; yapprox = [alpha;sol;beta]; hold off plot(tpoints,yapprox,’bo’) hold on plot(tpoints,yexact,’k-‘) xlabel(‘t’);ylabel(‘y’); function D=DerivMatrix(x,n) D=zeros(n+1,n+1); for i=1:n+1 num(i)=1.0; for […]

clear all close all eps=1.0; p=ones(2,1); %Guess value for p while eps>1.0e-14 [gfunc,jacobian]=ExampleFunctionWithJacobian(p); eps=max(abs(gfunc)); temp=jacobian\(-gfunc); %this is more efficient than inv(jacobian)*(-gfunc) p=temp+p; end options = optimoptions(@fsolve,’Display’,’iter’,’SpecifyObjectiveGradient’,true); [pmatlab,F,exitflag,output,JAC] = fsolve(@ExampleFunctionWithJacobian,[1;1],options); %Check to see if same as solution given by MATLAB fsolve() abs(p-pmatlab) function [gfunc,jac]=ExampleFunctionWithJacobian(p) gfunc=[2*p(1)-2*p(2)-exp(-p(1)); -p(1)+2*p(2)-exp(-p(2))]; jac=[2+exp(-p(1)) -2;-1 2+exp(-p(2))]; end

[Content_Types].xml _rels/.rels matlab/document.xml matlab/output.xml metadata/coreProperties.xml metadata/mwcoreProperties.xml metadata/mwcorePropertiesExtension.xml metadata/mwcorePropertiesReleaseInfo.xml MATLAB Root Finding Methods In this livescript, you will learn how to Plot graph of a function Use the inbuilt MATLAB \texttt{roots()} function to find roots of a polynomial function Use the inbuilt MATLAB \texttt{fzero()} function to find roots of a general function Suppose you would like

ENGR20005 Numerical Methods in Engineering Workshop 3 Part A: MATLAB Livescripts 3.1 The livescript ENGR20005 Workshop3p1.mlx runs through the use of functions in MATLAB. (a) Read through the livescript and make sure you understand what each line of code does. (b) Convert your .m files from workshop 2 to functions that output the root xr.

[Content_Types].xml _rels/.rels matlab/document.xml matlab/output.xml metadata/coreProperties.xml metadata/mwcoreProperties.xml metadata/mwcorePropertiesExtension.xml metadata/mwcorePropertiesReleaseInfo.xml The Spectral Collocation Method for Solving Boundary Value Problems In this livescript, you will learn how To solve boundary value problems using the spectral collocation method. Consider the linear two point boundary value problem \frac{d^{2}y}{dx^{2}}+2\frac{dy}{dx}-6y=0 with the boundary conditions y(-1)=1 an y'(1)=1 . y0 = 1; yN

ENGR20005 Numerical Methods in Engineering Workshop 2 Part A: MATLAB Livescripts 2.1 The livescript ENGR20005 Workshop2p1.mlx shows you how to use while loops and conditional statements by looking at the sample problem f(x) = x − 2. (2.1) Read through and experiment with the livescript and make sure you understand each of the commands in

clear all xpoints=linspace(-1,1,5); %[xpoints,w]=lgwt(25,-1,1); %[xpoints,w,P]=lglnodes(2); ypoints=1./(1+25*xpoints.^2); n=length(xpoints)-1; %order of polynomial xint=-1:0.01:1; yint=LagrangePolynomial(xpoints,ypoints,n,xint); yreal=1./(1+25*xint.^2); hold off plot(xpoints,ypoints,’ko’,’MarkerSize’,20,’MarkerFaceColor’,’r’) hold on plot(xint,yint,’b-‘,’Linewidth’,4) plot(xint,yreal,’k-‘,’Linewidth’,4) xlabel(‘x’) ylabel(‘y’) function yint=LagrangePolynomial(x,y,n,xint) yint=zeros(size(xint)); for i=1:n+1 Lix=1.0; for j=1:n+1 if j ~= i Lix=Lix.*(xint-x(j))/(x(i)-x(j)); end end yint=yint+y(i).*Lix; end end function [x,w]=lgwt(N,a,b) % lgwt.m % % This script is for computing definite integrals using Legendre-Gauss