Question

6. [15] In this question you will write a user-defined function, which contains one sub-function. Create a vector t with exactly 20001 components starting at 0 and ending with 200. Now you will use this vector t and MATLAB's plot3 to create a curve in 3D given by the parametric equations: X(t) = t. y(t) = 1- tz(t) = A for-loop must then be used to start on the first point (1 = 0) and move along the curve till the last point = 200. At each point your program must test if the value of the x-coordinate at that point is simultaneously divisible without remainder by two positive integers n, and nz. These values of n, and n, are the inputs to the Main function. If a particular x-coordinate is in fact divisible by both n, and 12, a red star must be placed on the curve at the point with that particular x-coordinate. After a red star is plotted, use MATLAB 5 pause function to let MATLAB wait a very short while (for example a tenth of a second) before moving on to the next iteration of the loop. Doing this will make it appear as if the red stars are popping up on the curve like in an animation. With n = 3 and na = 4 the final image will look like below (box on was used). 1 0 100 200 150 -100 10 50 -200 0

Please use Matlab for coding and good commentary please.
Thank you.

Answer 1

Hi,

Please find the answer below;

%%%Driver Function

function [] = main()
%main function
%declare variables
n1=3;
n2=4;
figure
%Invoking the function with paramters
graphFunction(n1,n2);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [] = graphFunction(n1,n2)
%animation3
%t vector with points
t=linspace(0,200,2001);
%for loop to plot
for i = 1:length(t)
    %calculate x,y,z
    x=t(i);
    y=1-t(i);
    z=t(i)^2;
    %check if x is divisible by n1 and n2
    if rem(x,n1)==0 && rem(x,n2)==0
        plot3(x,y,z,'-o','Color','r','MarkerSize',15);
        hold on;
        grid on;
        pause(0.001);
    end
    plot3(x,y,z,'-o','Color','b');
    hold on;
    grid on;
    %plot labels
    xlabel('x(t)')
    ylabel('y(t)')
    zlabel('z(t)')
end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

D Х Figure 2 File Edit View Insert Tools Desktop Window Help Search Documentation Sign In N CU *104 X 4 3 z(t) 1 0 50 De w 200 -50 150 Na -100 100 - 150 50 y(t) -200 0 x(t) slanci 23 - ylabel ('y(t)') Command Window Now to MATLAR? Seo recources for Getting Started X