This statement is not inside any function. (It follows the END that terminates the definition of the function "trap_wfp610".)
13 views (last 30 days)
Show older comments
%% function [ w ] = trap_wfp610( n ) % Calculates the work done using numerical integration and the trapezoidal rule % Detailed explanation goes here h=(5*10^-2)/n; a=0; b=5*10^-2; f=@(a)(5*(10^5)*a+10^5)*pi*(5*10^-2)^2; g=@(b)(5*(10^5)*b+10^5)*pi*(5*10^-2)^2; w=h/2*(f(a)+g(b)); end
%% syms x k1=500000 k2=10^5 p=(k1*x+k2)*pi*.05^2 W=int(p,x,0,5*10^-2)
%% error=abs((w-W)/W)*100
0 Comments
Accepted Answer
Mischa Kim
on 11 Dec 2014
Kevin, it's hard to tell what exactly you would like to do and what is going wrong. One thing you could do is do both computations within the same function and return all results:
function [ w,Ws,error ] = trap_wfp610( n )
% Calculates the work done using numerical integration and the trapezoidal rule
% Detailed explanation goes here
h = (5*10^-2)/n;
a = 0;
b = 5*10^-2;
f = @(a)(5*(10^5)*a+10^5)*pi*(5*10^-2)^2;
g = @(b)(5*(10^5)*b+10^5)*pi*(5*10^-2)^2;
w = h/2*(f(a)+g(b));
syms x
k1 = 500000
k2 = 10^5
p = (k1*x+k2)*pi*.05^2
Ws = int(p,x,0,5*10^-2)
error = abs((w-Ws)/Ws)*100;
end
3 Comments
Mischa Kim
on 11 Dec 2014
Have you assigned a value to n? In other words try
[ w,Ws,error ] = trap_wfp610( 3 )
or
n = 3;
[ w,Ws,error ] = trap_wfp610( n )
in the MATLAB command window.
More Answers (0)
See Also
Categories
Find more on Number Theory in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!