MATLAB: How to make the function accept vectors instead of scalars(I’m Really New to Matlab)

for loopnewbietaylor seriesvectors

I'm trying to make a taylor series expression for e^x = sum(x^n/n!) for n=0:50 and for x^n/n!>0.01 So far I have this:

function [ xs ] = myExpfuntion(x)
x=input('Enter x: ');
xl=zeros(51,1);                          %preallocation of the storage
for n= 1:50
    if ((x^n)/factorial(n))>0.01
        xl(n+1) = ((x^n)/factorial(n)); %save every iteration result
    end
end
xs=1+sum(xl) ;                            %calculate the sum of the iteration results
end

I want to delete the line "x=input('Enter x: ');" and just place my value for x in the function i.e. myExpfunction(x) and also make the function work with x being a row vector i.e x=[1 2 3]

My current function wont allow me to do so and I don't get what I'm doing wrong. Please help!

Best Answer

Hi,

I modified you code so that vectors can be used as input.

However, I am not sure, if I did it in the right way. It should help you though to get what you want:

function [ xs ] = myExpfuntion(x)
xs=zeros(1,length(x));
xl=zeros(51,length(x));                          %preallocation of the storage
for n= 1:50
    if ((x.^n)/factorial(n))>0.01
        xl(n+1,:) = ((x.^n)/factorial(n)); %save every iteration result
    end
end
for i = 1:length(x)
    xs(i)=1+sum(xl(:,i)) ;  %calculate the sum of the iteration results
end
end

This code results in:

>> xs=myExpfuntion([1,2,3,4])
xs =
      2.7083    7.0000   16.3750   34.3333

Related Solutions

MATLAB: Summation with looping

for n=0:16 is the correct syntax, you are also making another error by doing the summation in every iteration, you should instead find all values and do the summation after the loop ends

x = 1.8
xl=zeros(17,1);                 %preallocation of the storage
for n=0:16
xl(n+1) = ((x^n)/factorial(n)); %save every iteration result
end
xs=sum(xl)                      %calculate the sum of the iteration results

it can also be done without loops

n=0:16;
xs = sum(((x.^n)./factorial(n))); %notice the dots

MATLAB: Finalizing the Code for taylor series, gets error

Japoe25 - the error message is telling you that the assignment at

 xl(n+1,:) = ((x.^n)/factorial(n));

is failing due to trying to assign too many elements to x1. Note that x1 is a 51x1 matrix whereas x is a 1x4 - so you are trying to squeeze four elements into one. You need to resize x1 as

 xl=zeros(51,length(x));

The above will allow you to continue and get an answer that is almost correct. Note the condition for the if statement

 if ((x.^n)/factorial(n))>0.01

Remember that x is a 1x4 array and so this condition will be false when at least one element of the four element ((x.^n)/factorial(n)) is less than 0.01. You want to continue updating xl so long as at least one element of this matrix satisfies the condition. You can do this with the any function as

 if any(((x.^n)/factorial(n))>0.01)

Adding this will give you results that are similar to what you are looking for. When I tried this, I found the answer to be

 2.7183       7.3891       20.085       54.594

which is almost what you are looking for. You get the correct answers if you pass in each element individually as

 myExpFunction(1)
 myExpFunction(2) 
 myExpFunction(3)
 myExpFunction(4)

The difference is due to the any - as long as one element in x produces a result that is greater than 0.01 then we continue computing the expansion for all values in x. To get around this you could do

 for n= 1:50
    temp = ((x.^n)/factorial(n)); %save every iteration result
    temp(temp <0.01) = 0;
    xl(n+1,:) = temp;
 end

In the above, we create a temporary array and then set all of its values that are less than 0.01 to zero before the assignment to xl. This will give you the answer that you are looking for.

Best Answer

Related Solutions

MATLAB: Summation with looping

MATLAB: Finalizing the Code for taylor series, gets error

Related Question