MATLAB: Help in numerical exercise

n_p

Hi everyone,

I have a vector y2 starting from 1 to 100 with properties of y2(1)=1 and y2(i+1)-y2(i)=i+1

Knowing that i=1,…,99 what it will be y2(100)?

Thank you very much!!!

Best Regards.

Best Answer

If you look at this site in the section called Anecdotes" you will see a short cut to obtaining your answer. It doesn't require matlab at all.

http://en.wikipedia.org/wiki/Carl_Friedrich_Gauss

Related Solutions

MATLAB: Failure in initial user-supplied objective function evaluation. FSOLVE cannot continue

Hi,

this runs - not sure if it is really efficient, but it works:

x0=[1,1,1,1,1];
fsolve(@F,x0)
function m=F(x)
syms  y1 y2 y3 y4 z1 k11 k12 k13 k14 l11 
h=1/1000;
f1(y1,y2,y3,y4,z1)=(-(z1^3)/(y3^2))*(3*(y2-y1+y3^(-1)-z1/10)^2+(1/5)*(y2-y1+y3^(-1)-(z1)/10))-y4;
f2(y1,y2,y3,y4,z1)=1/10*z1-y4;
f3(y1,y2,y3,y4,z1)=(z1^3)*(3*(y2-y1+y3^(-1)-z1/10)^2+(1/5)*(y2-y1+y3^(-1)-(z1)/10));
f4(y1,y2,y3,y4,z1)=y1-y3^(-1);
g(y1,y2,y3,y4,z1)=(y1-y3^(-1))^2+y4^2-1/10*z1;
J1(y1,y2,y3,y4,z1)=jacobian([f1,f2,f3,f4],[y1,y2,y3,y4]);
J1=double(J1(2,2,1,0,10));
J2(y1,y2,y3,y4,z1)=jacobian([f1,f2,f3,f4],z1);
J2=double(J2(2,2,1,0,10));
J3(y1,y2,y3,y4,z1)=jacobian(g,[y1,y2,y3,y4]);
J3=double(J3(2,2,1,0,10));
J4(y1,y2,y3,y4,z1)=jacobian(g,z1);
J4=double(J4(2,2,1,0,10));
k11=x(1);
k12=x(2);
k13=x(3);
k14=x(4);
l11=x(5);
k1=[k11;k12;k13;k14];
m(1)=double(h*(f1(2,2,1,0,10)+0.44*(J1(1,:)*k1 +J2(1,:)*l11))-k11);
m(2)=double(h*(f2(2,2,1,0,10)+0.44*(J1(2,:)*k1 +J2(2,:)*l11))-k12);
m(3)=double(h*(f3(2,2,1,0,10)+0.44*(J1(3,:)*k1 +J2(3,:)*l11))-k13);
m(4)=double(h*(f4(2,2,1,0,10)+0.44*(J1(4,:)*k1 +J2(4,:)*l11))-k14);
m(5)=double(g(2,2,1,0,10)+0.44*(J3*k1 +J4(1,:)*l11));
end

Best regards

Stephan

MATLAB: How to merge multiple graphs/vectors in one!

To begin with, your data are apparently all row vectors, so I would concatenate them as cell arrays:

xall = {fliplr(x1); fliplr(x2); x3; x4; fliplr(x5); fliplr(x6); x7; x8};
yall = {fliplr(y1); fliplr(y2); y3; y4; fliplr(y5); fliplr(y6); y7; y8};
figure(1)
hold on
for k1 = 1:size(xall,1)
    plot(xall{k1}, yall{k1}, 'x')
end
hold off

I would not bother with unique, since the points will simply overplot each other, since you’re plotting them as discrete points and not as lines. In R2014b, the colours will cycle automatically. Otherwise, you will have to cycle through the colours yourself, but that is not difficult.

Best Answer

Related Solutions

MATLAB: Failure in initial user-supplied objective function evaluation. FSOLVE cannot continue

MATLAB: How to merge multiple graphs/vectors in one!

Related Question