I have written a script to determine total energy out of a system.
The energy equation does not work when I attempt to use an interal function for two temperatures.
I have attached the script, from lines 82 to 122 containes the issues I am having.
I am attempting to use the equation C1 + C2*((C3/T)/(sinh(C3/T))^2)+C4*((C5/T_4)/(cosh(C5/T_4))^2 to calculate the specific heat capacity.
Where C1-C5 are given.
The T variable is the integral between two limits.
This calculation of the specific heat capacity is need to be done 6 times.
I have tried setting all the C values to values with different numbers e.g C11-C15 and then the next component to C21-C25.
With every alteration i get this same result.
Any help would be much appreciated
Undefined function 'function_Cp_1' for input arguments of type 'double'.
Error in integralCalc/iterateScalarValued (line 314)
fx = FUN(t);
Error in integralCalc/vadapt (line 132)
[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);
Error in integralCalc (line 75)
[q,errbnd] = vadapt(@AtoBInvTransform,interval);
Error in integral (line 87)
Q = integralCalc(fun,a,b,opstruct);
Error in Matlab_Assignment_2 (line 90)
Cp_C4H6O2_i=integral(@function_Cp_1,Ti,T1);
T2=(Exit_temperature_degrees_Celcius+273.15);
function_Cp_1=@(T_1) (C1+(C2*((C3/T_1)/(sinh(C3/T_1)))^2)+(C4*((C5/T_1)/(cosh(C5/T_1)))^2));
Cp_C4H6O2_i=integral(@function_Cp_1,Ti,T1);
Cp_C4H6O2_f=integral(@function_Cp_1,Ti,T2);
function_Cp_2=@(T_2) (C1+(C2*((C3/T_2)/(sinh(C3/T_2)))^2)+(C4*((C5/T_2)/(cosh(C5/T_2)))^2));
Cp_O2_i=integral(@function_Cp_2,Ti,T1);
Cp_O2_f=integral(@function_Cp_2,Ti,T2);
function_Cp_3=@(T_3) (C1+(C2*((C3/T_3)/(sinh(C3/T_3))^2)+(C4*((C5/T_3)/(cosh(C5/T_3)))^2)));
Cp_CO2=integral(@function_Cp_3,Ti,T1);
function_Cp_4=@(T_4) (C1+(C2*((C3/T_4)/(sinh(C3/T_4))^2)+(C4*((C5/T_4)/(cosh(C5/T_4)))^2)));
Cp_H2O=integral(@function_Cp_4,Ti,T2);
Undefined function 'function_Cp_1' for input arguments of type 'double'.
Error in integralCalc/iterateScalarValued (line 314)
Error in integralCalc/vadapt (line 132)
[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);
Error in integralCalc (line 75)
[q,errbnd] = vadapt(@AtoBInvTransform,interval);
Error in integral (line 87)
Q = integralCalc(fun,a,b,opstruct);
Error in Matlab_Assignment_2 (line 90)
Cp_C4H6O2_i=integral(@function_Cp_1,Ti,T1);
Invalid expression. When calling a function or indexing a variable, use parentheses. Otherwise, check for mismatched delimiters.
