This statement is not inside any function. (It follows the END that terminates the definition of the function "trap_wfp610".)

조회 수: 14 (최근 30일)
%% 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

채택된 답변

Mischa Kim
Mischa Kim 2014년 12월 11일
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
Mischa Kim
Mischa Kim 2014년 12월 11일
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.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Calculus에 대해 자세히 알아보기

태그

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by